[BioC] Affy normalized data using cdf env and rma
Junshi Yazaki
jyazaki at salk.edu
Tue May 24 08:18:21 CEST 2005
Hi all,
I have constructed cdf environment on my Mac. My environment and
normalization by rma looks work under below working flow (******).
But I can not get whole results. I got normalized data from each
control gene spots but all spots data from genes has merged all in
one line like "data ex.". copy9_.... means control spots. non_at and
random_at means genes. I do not know why I can not get normalized
data from all each spots. Could you please help me someone?
"data ex."
.....many control genes...........
copy9_AFFX-r2-At-Ubq-5_s_at 10.8117692745457
10.9959999081142 10.1008637507680
10.0693470115082
copy9_AFFX-r2-At-Ubq-5_x_at 11.9637427616359
12.0014521293535 11.8231131697170
11.9823112800146
non_at 2.77632266000964 3.06545224338101
2.77787004643169 3.07748076274474
random_at 4.2258095000546 3.72907697958284
4.35455022393967 3.71650013156617
******
> make.cdf.package("1-1.CDF", packagename="49chip14micronmercurycdf")
Reading CDF file.
Creating CDF environment
Wait for about 6 dots.......
Creating package in
/Users/junshiyazaki/Desktop/R_test_data/49chip14micronmercurycdf
[1] "49chip14micronmercurycdf"
> cel.files=list.files(pattern=".CEL$")
> data=ReadAffy(filenames=cel.files)
> chip14micronmercurycdf=env
> data at cdfName<- "chip14micronmercurycdf"
> temp=rma(data)
Background correcting
Normalizing
Calculating Expression
> write.exprs(temp, file="mydata8.txt")
******
--
***********************************************************
***********************************************************
Junshi Yazaki, Ph D
The Salk Institute for Biological Studies
Plant Biology Laboratories
10010 North Torrey Pines Road
La Jolla, CA 92037
Phone (858) 453-4100 x1533
FAX (858) 558-6379
Email: jyazaki at salk.edu
Web addresses:
http://signal.salk.edu
http://qtlpc.salk.edu/pbio/Web/ecker.html
More information about the Bioconductor
mailing list