[BioC] beadarray_1.7.6. problem/error with readBeadSummaryData()
md392 at cam.ac.uk
Tue Mar 18 17:23:28 CET 2008
This script was a bit outdated I'm afraid. The package was recently
changed so that it would use the ProbeID column in the input files to
give a unique row name in the resulting ExpressionSet. Unfortunately,
the file provided with the package did not have this column, hence R
gives an error when we try and assign a vector of length 0 to the
names of a matrix.
If you add the argument ProbeID= "TargetID" to readBeadSummaryData it
should work. Thanks to Vince Carey for also pointing this error out
and the solution.
There is also a problem in the code where we try and use BiomaRT to
get information about a set of probes. The vector illuids passed to
getBM is never defined. However, if you set to something appropriate
(maybe the rownames of the topTable output) then this command should
work fine also.
I will update the script shortly with these changes.
Hope this helps,
On 18 Mar 2008, at 15:34, Christian Kohler wrote:
> Dear Bioconductor-list,
> I recently started to go in for analysing Illumina (BeadChip) datasets
> with the help of the 'beadarray' package.
> As an Example R-script (within BeadSummaryExample) is provided, the
> introduction to this topic is quite straightforward.
> Until today, I used 'beadarray' in version 1.7.0 and the code snippets
> of the example script worked out fine.
> With 'beadarray' version 1.7.6. I experienced the following problem:
>> BSData <- readBeadSummaryData(dataFile, qcFile=qcFile,
> sampleSheet=sampleSheet,skip=7, columns=list(exprs="AVG_Signal",
> se.exprs="BEAD_STDEV", >NoBeads="Avg_NBEADS"),
> qc.columns=list(exprs="AVG.Signal", se.exprs="SeqVAR"),qc.sep=","
> ,sep=",", qc.skip=7, annoPkg="illuminaHumanv1")
>> Error in `row.names<-.data.frame`(`*tmp*`, value = "ProbeID") :
>> invalid 'row.names' length
> Now, I am just wondering if other 'beadarray'-users came across
> this error .
> R version 2.6.2 (2008-02-08)
> attached base packages:
>  tools stats graphics grDevices utils datasets
>  base
> other attached packages:
>  beadarray_1.7.6 affy_1.16.0 preprocessCore_1.0.0
>  affyio_1.6.0 geneplotter_1.16.0 lattice_0.17-4
>  annotate_1.16.0 xtable_1.5-1 AnnotationDbi_1.0.4
>  RSQLite_0.6-3 DBI_0.2-3 Biobase_1.16.1
>  limma_2.12.0
> loaded via a namespace (and not attached):
>  KernSmooth_2.22-22 RColorBrewer_1.0-1 grid_2.6.2
>  tcltk_2.6.2
> Thanks a lot for any help/suggestions.
> Christian Kohler
> Institute of Functional Genomics
> Computational Diagnostics
> University of Regensburg (BioPark I)
> D-93147 Regensburg (Germany)
> Tel. +49 941 943 5055
> Fax +49 941 943 5020
> christian.kohler at klinik.uni-regensburg.de
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives: http://news.gmane.org/
More information about the Bioconductor