[BioC] Error by using gpr data with marray

Wolfgang Huber whuber at embl.de
Fri Mar 9 14:52:10 CET 2012


Dear Sandra

in principle, this should work:

  library("ArrayExpress")
  a = ArrayExpress("E-MTAB-109")

I just tried this. An object with data was constructed (a list of 
NChannelSets). However, errors and warnings were thrown. The resulting 
datastructure was more messy than it should be, and the metadata was 
dropped. I'll chase this up with the package maintainer and with the EBI 
people.

You can try to read to the files provided by EBI 'manually' using 
limma's read.maimages etc., but for this you will need to twiddle with 
function call parameters and file formats, and for this you might be 
best off recruiting a local R expert.

	Best wishes
	Wolfgang


Sandra Heise [guest] scripsit 03/09/2012 11:19 AM:
>
> hi,
> I am just a beginner in using R and bioconductor. I have a microarray experiment from the arrayexpress database.
>
> http://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-109
>
> There are given the target txt and row data (gpr). Now i want used
> bioconductor to generate a nomalized intensity matrix for all knockout experiments. I thought i have to use limma package and tried to read the data with marray function read.GenePix but there is still an error. Then I tried to read the info txt but there is still an error too. Maybe i just chose the wrong package for my dataset. Unfortunately I don't have a clue how i can generate the normalized datamatrix but the autors of the dataset wrote thez generate the log_fold_change matrix with Imfit limma function. So i thought there is a possibility for me.
>
> please help
>
>   -- output of sessionInfo():
>
> library(limma)
> library(marray)
>
> read.GenePix("3088.gpr")
>
>>   Error in if (skip>  0L) readLines(file, skip) :
>>   missing value where TRUE/FALSE needed
>
> a.target<- read.marrayInfo("E-MTAB-109.sdrf.txt")
>
>>   Error in names(x)<- value :
>>   'names' attribute [30] must be the same length as the vector [7]
>
>
>
> --
> Sent via the guest posting facility at bioconductor.org.
>
> _______________________________________________
> 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


-- 
Best wishes
	Wolfgang

Wolfgang Huber
EMBL
http://www.embl.de/research/units/genome_biology/huber



More information about the Bioconductor mailing list