[BioC] arrayQualityMetrics error with MAList
Wolfgang Huber
whuber at embl.de
Thu Mar 1 22:28:40 CET 2012
Dear Karen
in fact, there is no coercion method from MAList to NChannelSet, yet
there is one from MAList to ExpressionSet (and all that actually makes
sense). Thus, can you please try
library(convert)
MA.avg_TEST1_eSet <- as(MA.avg_TEST1, 'ExpressionSet')
arrayQualityMetrics(MA.avg_TEST1_eSet, outdir="Test1_report")
Let me know if that works. Sometimes it happens that this fails with an
error, since the metadata in limma MALists tends to be less strictly
organised than what is expected by Biobase. In that case, send me the
offending object, and I will try to tweak the coercion method to deal
with it.
Best wishes
Wolfgang
Mar/1/12 9:48 PM, Wolfgang Huber scripsit::
> Dear Karen
>
> what happens when you run
>
> library("convert")
>
> before
>
> MA.avg_TEST1_nChannel<-as(MA.avg_TEST1, 'NChannelSet')
>
> If that does not help, can you please send me (offline) the offending
> object 'MA.avg_TEST1'?
>
> I suspect I need to fix the NAMESPACE of arrayQualiyMetrics, sorry for
> that.
>
> best wishes
> Wolfgang
>
> Mar/1/12 3:45 PM, karen power scripsit::
>> Hi,
>>
>> I'm trying to use arrayQualityMetrics to investigate 2 colour Agilent
>> microarrays. However I get the following error;
>>
>>> arrayQualityMetrics(MA.avg_TEST1, outdir="Test1_report")
>> The report will be written into directory 'Test1_report'.
>> Error in as(expressionset, "NChannelSet") :
>> no method or default for coercing "MAList" to "NChannelSet"
>> Error in prepdata(expressionset, intgroup = intgroup, do.logtransform =
>> do.logtransform) :
>> Argument 'expressionset' is of class 'MAList', and its automatic
>> conversion
>> into 'NChannelSet' failed. Please try to convert it manually.
>>
>>
>> I have no idea how to succesfully turn an MAList into NChannelSet as
>> "MA.avg_TEST1_nChannel<-as(MA.avg_TEST1, 'NChannelSet')"
>> doesn't work.
>>
>> Is it something I'm doing wrong with the arrayQualityMetrics command?
>> Or else is
>> there a way around it? Session is as follows;
>>
>>> sessionInfo()
>> R version 2.14.2 (2012-02-29)
>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
>>
>> locale:
>> [1] en_IE.UTF-8/en_IE.UTF-8/en_IE.UTF-8/C/en_IE.UTF-8/en_IE.UTF-8
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods base
>>
>> other attached packages:
>> [1] arrayQualityMetrics_3.10.0 BiocInstaller_1.2.1 limma_3.10.2
>>
>> loaded via a namespace (and not attached):
>> [1] affy_1.32.1 affyio_1.22.0 affyPLM_1.30.0
>> annotate_1.32.1 AnnotationDbi_1.16.18 beadarray_2.4.1 Biobase_2.14.0
>> Biostrings_2.22.0
>> [9] Cairo_1.5-1 cluster_1.14.2 DBI_0.2-5
>> genefilter_1.36.0 grid_2.14.2 Hmisc_3.9-2 hwriter_1.3
>> IRanges_1.12.6
>> [17] lattice_0.20-0 latticeExtra_0.6-19 preprocessCore_1.16.0
>> RColorBrewer_1.0-5 RSQLite_0.11.1 setRNG_2009.11-1 splines_2.14.2
>> survival_2.36-12
>> [25] SVGAnnotation_0.9-0 tools_2.14.2 vsn_3.22.0 XML_3.9-4
>> xtable_1.7-0 zlibbioc_1.0.0
>>
>>
>>
>> Regards,
>>
>> Karen
>>
>> Karen Power BSc PhD
>> Room 039,
>> UCD Centre for Food Safety,
>> School of Public Health, Physiotherapy and Population Science
>> Veterinary Sciences Centre,
>> University College Dublin,
>> Belfield,
>> Dublin 4,
>> Ireland.
>>
>>
>>
>>
>> [[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
>
>
--
Best wishes
Wolfgang
Wolfgang Huber
EMBL
http://www.embl.de/research/units/genome_biology/huber
More information about the Bioconductor
mailing list