Analyze Single Color Agilent Microarray Data
2
0
Entering edit mode
@srikanth-manda-srinivas-4590
Last seen 9.6 years ago
Hi all, I have data from single color agilent array exposed and unexposed to drug. The exposed data is from 1x44 platform and unexposed from 4x44k. When i try limma and load the files using read.maimages i get the error: Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of items to replace is not a multiple of replacement length Please help me with this. How can load and perform analysis on this data. Regards, Srikanth -- Srinivas M. Srikanth Ph.D. Student Institute of Bioinformatics Discoverer, 7th Floor, International Technology Park, Bangalore, India Mob:+917259692031 [[alternative HTML version deleted]]
• 1.2k views
ADD COMMENT
0
Entering edit mode
@sean-davis-490
Last seen 3 months ago
United States
Hi, Srikanth. Please always provide the code you are using (cut-and-paste), the error message, and sessionInfo() when asking for help. That will help the list members to more easily figure out what the problem might be. Sean On Wed, Aug 3, 2011 at 8:49 AM, Srinivas M. Srikanth <srikanth at="" ibioinformatics.org=""> wrote: > Hi all, > > I have data from single color agilent array exposed and unexposed to drug. > The exposed data is from 1x44 platform and unexposed from 4x44k. When i try > limma and load the files using read.maimages i get the error: > > Error in RG[[a]][, i] <- obj[, columns[[a]]] : > ?number of items to replace is not a multiple of replacement length > > Please help me with this. How can load and perform analysis on this data. > > Regards, > Srikanth > > -- > Srinivas M. Srikanth > Ph.D. Student > Institute of Bioinformatics > Discoverer, 7th Floor, > International Technology Park, > Bangalore, India > Mob:+917259692031 > > ? ? ? ?[[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >
ADD COMMENT
0
Entering edit mode
Hi Sean, I used following commands to load the files in R. I created a target.txt file with content as SerialNumber FileName Condition 1 Normal1.txt Normal >library(limma) >target<-readTargets() >RG<-read.maimages(target,(targets,columns = list(G = "gMedianSignal", Gb = "gBGMedianSignal", R = "gProcessedSignal",Rb = "gIsPosAndSignif"), annotation = c("Row", "Col","FeatureNum", "ControlType","ProbeName"))) I used some values of R and Rb After it reads first 18 files it throws the message: Read Normal18.txt Read Wles1.txt Read Wles2.txt Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of items to replace is not a multiple of replacement length The sessionInfo() gives: > sessionInfo() R version 2.11.1 (2010-05-31) i386-pc-mingw32 locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] limma_3.4.5 loaded via a namespace (and not attached): [1] Biobase_2.8.0 Please help. Regards, Srikanth On Wed, Aug 3, 2011 at 6:22 PM, Sean Davis <sdavis2@mail.nih.gov> wrote: > Hi, Srikanth. > > Please always provide the code you are using (cut-and-paste), the > error message, and sessionInfo() when asking for help. That will help > the list members to more easily figure out what the problem might be. > > Sean > > > On Wed, Aug 3, 2011 at 8:49 AM, Srinivas M. Srikanth > <srikanth@ibioinformatics.org> wrote: > > Hi all, > > > > I have data from single color agilent array exposed and unexposed to > drug. > > The exposed data is from 1x44 platform and unexposed from 4x44k. When i > try > > limma and load the files using read.maimages i get the error: > > > > Error in RG[[a]][, i] <- obj[, columns[[a]]] : > > number of items to replace is not a multiple of replacement length > > > > Please help me with this. How can load and perform analysis on this data. > > > > Regards, > > Srikanth > > > > -- > > Srinivas M. Srikanth > > Ph.D. Student > > Institute of Bioinformatics > > Discoverer, 7th Floor, > > International Technology Park, > > Bangalore, India > > Mob:+917259692031 > > > > [[alternative HTML version deleted]] > > > > _______________________________________________ > > Bioconductor mailing list > > Bioconductor@r-project.org > > https://stat.ethz.ch/mailman/listinfo/bioconductor > > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > > > -- Srinivas M. Srikanth Ph.D. Student Institute of Bioinformatics Discoverer, 7th Floor, International Technology Park, Bangalore, India Mob:+917259692031 [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
Thanks, Srikanth, for the details. I think that the input that you have gotten already probably sums things up--your best bet is to load things separately (the two different platforms) and then combine them in R. I do not have a script or know of a simple command to do so. Whether or not the data are going to be comparable across platforms is another question that may or may not be important, depending on your experimental design. As an aside, you are running a fairly old version of R, so you might consider updating to a newer version. In the end, I do not have a direct answer for you--sorry. Sean On Wed, Aug 3, 2011 at 11:18 AM, Srinivas M. Srikanth <srikanth at="" ibioinformatics.org=""> wrote: > Hi Sean, > > I used following commands to load the files in R. I created a target.txt > file with content as > SerialNumber FileName Condition > 1 Normal1.txt Normal > >>library(limma) >>target<-readTargets() >>RG<-read.maimages(target,(targets,columns = list(G = "gMedianSignal", Gb = > "gBGMedianSignal", R = "gProcessedSignal",Rb = "gIsPosAndSignif"), > annotation = c("Row", "Col","FeatureNum", "ControlType","ProbeName"))) > I used some values of R and Rb > > After it reads first 18 files it throws the message: > > Read Normal18.txt > Read Wles1.txt > Read Wles2.txt > Error in RG[[a]][, i] <- obj[, columns[[a]]] : > ?number of items to replace is not a multiple of replacement length > > > > The sessionInfo() gives: >> sessionInfo() > R version 2.11.1 (2010-05-31) > i386-pc-mingw32 > > locale: > [1] LC_COLLATE=English_United States.1252 > [2] LC_CTYPE=English_United States.1252 > [3] LC_MONETARY=English_United States.1252 > [4] LC_NUMERIC=C > [5] LC_TIME=English_United States.1252 > > attached base packages: > [1] stats ? ? graphics ?grDevices utils ? ? datasets ?methods ? base > > other attached packages: > [1] limma_3.4.5 > > loaded via a namespace (and not attached): > [1] Biobase_2.8.0 > > > Please help. > > Regards, > Srikanth > > On Wed, Aug 3, 2011 at 6:22 PM, Sean Davis <sdavis2 at="" mail.nih.gov=""> wrote: > >> Hi, Srikanth. >> >> Please always provide the code you are using (cut-and-paste), the >> error message, and sessionInfo() when asking for help. ?That will help >> the list members to more easily figure out what the problem might be. >> >> Sean >> >> >> On Wed, Aug 3, 2011 at 8:49 AM, Srinivas M. Srikanth >> <srikanth at="" ibioinformatics.org=""> wrote: >> > Hi all, >> > >> > I have data from single color agilent array exposed and unexposed to >> drug. >> > The exposed data is from 1x44 platform and unexposed from 4x44k. When i >> try >> > limma and load the files using read.maimages i get the error: >> > >> > Error in RG[[a]][, i] <- obj[, columns[[a]]] : >> > ?number of items to replace is not a multiple of replacement length >> > >> > Please help me with this. How can load and perform analysis on this data. >> > >> > Regards, >> > Srikanth >> > >> > -- >> > Srinivas M. Srikanth >> > Ph.D. Student >> > Institute of Bioinformatics >> > Discoverer, 7th Floor, >> > International Technology Park, >> > Bangalore, India >> > Mob:+917259692031 >> > >> > ? ? ? ?[[alternative HTML version deleted]] >> > >> > _______________________________________________ >> > Bioconductor mailing list >> > Bioconductor at r-project.org >> > https://stat.ethz.ch/mailman/listinfo/bioconductor >> > Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > >> > > > > -- > Srinivas M. Srikanth > Ph.D. Student > Institute of Bioinformatics > Discoverer, 7th Floor, > International Technology Park, > Bangalore, India > Mob:+917259692031 > > ? ? ? ?[[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >
ADD REPLY
0
Entering edit mode
Axel Klenk ★ 1.0k
@axel-klenk-3224
Last seen 16 hours ago
UPF, Barcelona, Spain
Hi Srikanth, you can of course try to read them into separate EList objects, filter those for the common oligo subset (will depend on exactly what array type you have), combine the subsets and continue with limma. This will get you around the purely technical problem of having different array layouts BUT obviously not around the more serious problem of the unfortunate study design that is confounding platform difference and treatment effect... Best, - axel Axel Klenk Research Informatician Actelion Pharmaceuticals Ltd / Gewerbestrasse 16 / CH-4123 Allschwil / Switzerland From: "Srinivas M. Srikanth" <srikanth at="" ibioinformatics.org=""> To: <bioconductor at="" stat.math.ethz.ch=""> Date: 03.08.2011 14:49 Subject: [BioC] Analyze Single Color Agilent Microarray Data Sent by: bioconductor-bounces at r-project.org Hi all, I have data from single color agilent array exposed and unexposed to drug. The exposed data is from 1x44 platform and unexposed from 4x44k. When i try limma and load the files using read.maimages i get the error: Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of items to replace is not a multiple of replacement length Please help me with this. How can load and perform analysis on this data. Regards, Srikanth -- Srinivas M. Srikanth Ph.D. Student Institute of Bioinformatics Discoverer, 7th Floor, International Technology Park, Bangalore, India Mob:+917259692031 [[alternative HTML version deleted]] _______________________________________________ Bioconductor mailing list Bioconductor at r-project.org https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor The information of this email and in any file transmitted with it is strictly confidential and may be legally privileged. It is intended solely for the addressee. If you are not the intended recipient, any copying, distribution or any other use of this email is prohibited and may be unlawful. In such case, you should please notify the sender immediately and destroy this email. The content of this email is not legally binding unless confirmed by letter. Any views expressed in this message are those of the individual sender, except where the message states otherwise and the sender is authorised to state them to be the views of the sender's company. For further information about Actelion please see our website at http://www.actelion.com
ADD COMMENT

Login before adding your answer.

Traffic: 738 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