ArrayVision and limma
1
0
Entering edit mode
Isaac Mehl ▴ 60
@isaac-mehl-417
Last seen 10.2 years ago
I just received some wonderful help from Jean on reading ArrayVision output files into marrayInput. Now i would like to try and analyze with limma. if i understand correctly there is no way to do this unless source = "ArrayVision" which doesn't exist. Read: https://www.stat.math.ethz.ch/pipermail/bioconductor/2003-August/00218 2.html which leads me to believe that the only way to accomplish this is for Gordon to write a separate function. Is this the only way to read in data? Won't this be a pain with many, many different software packages making spot calls?? Isn't it possible to simply specify which column contains the needed information (RF, Gf, Rb, etc) in a matrix?? I am attaching an example of a .csv output from ArrayVision if anyone would like a crack at it. -- -isaac mehl gene expression lab (gele) salk institute 10010 n. torrey pines rd. la jolla ca. 92037 http://genex.salk.edu -------------- next part -------------- A non-text attachment was scrubbed... Name: 12371987B.csv Type: text/x-comma-separated-values Size: 404677 bytes Desc: not available Url : https://www.stat.math.ethz.ch/pipermail/bioconductor/attachments /20030826/3e07d249/12371987B-0001.bin
• 1.1k views
ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 50 minutes ago
WEHI, Melbourne, Australia
At 12:12 PM 26/08/2003, Isaac Mehl wrote: >I just received some wonderful help from Jean on reading ArrayVision >output files into marrayInput. Now i would like to try and analyze with >limma. if i understand correctly there is no way to do this unless >source = "ArrayVision" which doesn't exist. Read: >https://www.stat.math.ethz.ch/pipermail/bioconductor/2003-August/0021 82.html >which leads me to believe that the only way to accomplish this is for >Gordon to write a separate function. Is this the only way to read in >data? Won't this be a pain with many, many different software packages >making spot calls?? Isn't it possible to simply specify which column >contains the needed information (RF, Gf, Rb, etc) in a matrix?? I don't see the problem. If you've read in the data using marrayInput and then normalized using marrayNorm, you can use any of the linear model functions in limma to analyse your data, as explained in the section called "Using limma with the marray Packages" in the LIMMA User's Guide. The easiest function to use would be 'lmFit', in the development version of limma, because it accepts 'marrayNorm' objects directly. If for some reason you want to read in your data using limma instead of marrayNorm, you could do so by using the 'columns' argument of 'read.maimages' to specify which columns of your data files correspond to Rf, Gf, Rb and Gb. This option has always been available in limma, although the development version of limma will do a better job of handling column names containing spaces, like the names in your file. As the help document for 'read.maimages' says, "At present, the image analysis programs SPOT, GenePix, QuantArray and Imagene are supported. Other image analysis programs can potentially be accommodated though a user-specified 'columns' argument." The reason that Imagene data required special intervention from me was (i) because Imagene is almost unique in that it writes the red and green channel information to separate files and (ii) because neither limma nor marrayNorm supported imagene explicitly. Gordon >I am attaching an example of a .csv output from ArrayVision if anyone >would like a crack at it. >-- >-isaac mehl >gene expression lab (gele) >salk institute >10010 n. torrey pines rd. >la jolla ca. 92037 >http://genex.salk.edu
ADD COMMENT
0
Entering edit mode
Sorry, I have to withdraw my statement that 'read.maimages' can read ArrayVision files. Looking at the attachment from Isaac, it seems that ArrayVision writes column names over two rows of the file, and in addition these column names include non-ascii characters. 'read.maimages' can't handle a non-standard feature like this without special programming. However, since Jean has apparently just written a function to read ArrayVision data into marray, the obvious solution is to use her function. Gordon At 03:50 PM 26/08/2003, Gordon Smyth: >At 12:12 PM 26/08/2003, Isaac Mehl wrote: >>I just received some wonderful help from Jean on reading ArrayVision >>output files into marrayInput. Now i would like to try and analyze with >>limma. if i understand correctly there is no way to do this unless >>source = "ArrayVision" which doesn't exist. Read: >>https://www.stat.math.ethz.ch/pipermail/bioconductor/2003-August/002 182.html >>which leads me to believe that the only way to accomplish this is for >>Gordon to write a separate function. Is this the only way to read in >>data? Won't this be a pain with many, many different software packages >>making spot calls?? Isn't it possible to simply specify which column >>contains the needed information (RF, Gf, Rb, etc) in a matrix?? > >I don't see the problem. If you've read in the data using marrayInput and >then normalized using marrayNorm, you can use any of the linear model >functions in limma to analyse your data, as explained in the section >called "Using limma with the marray Packages" in the LIMMA User's Guide. >The easiest function to use would be 'lmFit', in the development version >of limma, because it accepts 'marrayNorm' objects directly. > >If for some reason you want to read in your data using limma instead of >marrayNorm, you could do so by using the 'columns' argument of >'read.maimages' to specify which columns of your data files correspond to >Rf, Gf, Rb and Gb. This option has always been available in limma, >although the development version of limma will do a better job of handling >column names containing spaces, like the names in your file. As the help >document for 'read.maimages' says, > > "At present, the image analysis programs SPOT, GenePix, QuantArray > and Imagene are supported. Other image analysis programs can > potentially be accommodated though a user-specified 'columns' argument." > >The reason that Imagene data required special intervention from me was (i) >because Imagene is almost unique in that it writes the red and green >channel information to separate files and (ii) because neither limma nor >marrayNorm supported imagene explicitly. > >Gordon > >>I am attaching an example of a .csv output from ArrayVision if anyone >>would like a crack at it. >>-- >>-isaac mehl >>gene expression lab (gele) >>salk institute >>10010 n. torrey pines rd. >>la jolla ca. 92037 >>http://genex.salk.edu
ADD REPLY

Login before adding your answer.

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