Agilent's one color array
2
0
Entering edit mode
Weiyin Zhou ▴ 220
@weiyin-zhou-1970
Last seen 9.6 years ago
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20070614/ 2b1040c4/attachment.pl
• 644 views
ADD COMMENT
0
Entering edit mode
Weiwei Shi ★ 1.2k
@weiwei-shi-1407
Last seen 9.6 years ago
can you just use newData = as.matrix(data[, 2:5]) rownames(newData) = data[, 1] So your rownames contain ID you need. HTH, Weiwei On 6/14/07, Weiyin Zhou <weiyin.zhou at="" exonhit-usa.com=""> wrote: > Dear listers, > > > > I am analyzing Agilent's one-color array. I use limma package to read > data then do quantile normalization between arrays on RG$G and log2 the > data. Then I try to convert the data to exprSet object so I can > analyzing it the way similar to Affmetrix array. > > > > > data[1:5,] > > ID File1 File2 File3 File4 File5 File6 > > 1 A_15_P100001 9.92 9.91 9.86 10.07 10.03 9.86 > > 2 A_15_P100002 9.40 9.20 8.81 8.82 8.89 8.94 > > 3 A_15_P100003 4.08 3.67 3.32 2.47 2.89 2.91 > > 4 A_15_P100004 2.30 2.41 2.51 2.81 2.68 2.57 > > 5 A_15_P100005 11.29 11.38 11.25 10.94 11.06 11.29 > > > > > pd <- read.phenoData("pData.txt", header=TRUE, row.names=1) > > > show(pData(pd)) > > Target > > File1 group1 > > File2 group1 > > File3 group1 > > File4 group2 > > File5 group2 > > File6 group2 > > > > > eset <- new("exprSet", exprs = as.matrix(data), phenoData = pd) > > Error in validObject(.Object) : invalid class "exprSet" object: 1: > number of exprs columns different from number of pData rows > > invalid class "exprSet" object: 2: sampleNames different from names of > phenoData rows > > > > I know the problem is because the ID column. Is there any way I can > only read data starting from second column(File1) but also associate the > result with ID later on? > > > > > sessionInfo() > > R version 2.4.0 (2006-10-03) > > i386-pc-mingw32 > > > > locale: > > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > > > > attached base packages: > > [1] "splines" "tools" "methods" "stats" "graphics" > "grDevices" > > [7] "utils" "datasets" "base" > > > > other attached packages: > > biomaRt RCurl XML GOstats Category genefilter > survival > > "1.8.1" "0.8-0" "1.4-0" "2.0.4" "2.0.3" "1.12.0" > "2.30" > > KEGG RBGL annotate GO graph affy > affyio > > "1.12.0" "1.10.0" "1.12.0" "1.14.1" "1.12.0" "1.12.2" > "1.2.0" > > Biobase limma > > "1.12.2" "2.9.1" > > > > > > Thanks in advance, > > > > Weiyin > > > > > > > > Weiyin Zhou > > Senior Research Associate > > ExonHit Therapeutics, Inc. > > 217 Perry Parkway, Building # 5 > > Gaithersburg, MD 20877 > > > > email: Weiyin.zhou at exonhit-usa.com > > phone: 240.404.0184 > > fax: 240.683.7060 > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. "Did you always know?" "No, I did not. But I believed..." ---Matrix III
ADD COMMENT
0
Entering edit mode
Weiyin Zhou ▴ 220
@weiyin-zhou-1970
Last seen 9.6 years ago
It works. Thanks, Weiwei. Weiyin -----Original Message----- From: Weiwei Shi [mailto:helprhelp@gmail.com] Sent: Thursday, June 14, 2007 4:36 PM To: Weiyin Zhou Cc: bioconductor at stat.math.ethz.ch Subject: Re: [BioC] Agilent's one color array can you just use newData = as.matrix(data[, 2:5]) rownames(newData) = data[, 1] So your rownames contain ID you need. HTH, Weiwei On 6/14/07, Weiyin Zhou <weiyin.zhou at="" exonhit-usa.com=""> wrote: > Dear listers, > > > > I am analyzing Agilent's one-color array. I use limma package to read > data then do quantile normalization between arrays on RG$G and log2 the > data. Then I try to convert the data to exprSet object so I can > analyzing it the way similar to Affmetrix array. > > > > > data[1:5,] > > ID File1 File2 File3 File4 File5 File6 > > 1 A_15_P100001 9.92 9.91 9.86 10.07 10.03 9.86 > > 2 A_15_P100002 9.40 9.20 8.81 8.82 8.89 8.94 > > 3 A_15_P100003 4.08 3.67 3.32 2.47 2.89 2.91 > > 4 A_15_P100004 2.30 2.41 2.51 2.81 2.68 2.57 > > 5 A_15_P100005 11.29 11.38 11.25 10.94 11.06 11.29 > > > > > pd <- read.phenoData("pData.txt", header=TRUE, row.names=1) > > > show(pData(pd)) > > Target > > File1 group1 > > File2 group1 > > File3 group1 > > File4 group2 > > File5 group2 > > File6 group2 > > > > > eset <- new("exprSet", exprs = as.matrix(data), phenoData = pd) > > Error in validObject(.Object) : invalid class "exprSet" object: 1: > number of exprs columns different from number of pData rows > > invalid class "exprSet" object: 2: sampleNames different from names of > phenoData rows > > > > I know the problem is because the ID column. Is there any way I can > only read data starting from second column(File1) but also associate the > result with ID later on? > > > > > sessionInfo() > > R version 2.4.0 (2006-10-03) > > i386-pc-mingw32 > > > > locale: > > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 > > > > attached base packages: > > [1] "splines" "tools" "methods" "stats" "graphics" > "grDevices" > > [7] "utils" "datasets" "base" > > > > other attached packages: > > biomaRt RCurl XML GOstats Category genefilter > survival > > "1.8.1" "0.8-0" "1.4-0" "2.0.4" "2.0.3" "1.12.0" > "2.30" > > KEGG RBGL annotate GO graph affy > affyio > > "1.12.0" "1.10.0" "1.12.0" "1.14.1" "1.12.0" "1.12.2" > "1.2.0" > > Biobase limma > > "1.12.2" "2.9.1" > > > > > > Thanks in advance, > > > > Weiyin > > > > > > > > Weiyin Zhou > > Senior Research Associate > > ExonHit Therapeutics, Inc. > > 217 Perry Parkway, Building # 5 > > Gaithersburg, MD 20877 > > > > email: Weiyin.zhou at exonhit-usa.com > > phone: 240.404.0184 > > fax: 240.683.7060 > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > -- Weiwei Shi, Ph.D Research Scientist GeneGO, Inc. "Did you always know?" "No, I did not. But I believed..." ---Matrix III
ADD COMMENT

Login before adding your answer.

Traffic: 869 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6