[BioC] problem with nsFilter and hthgu133a format eset
Fraser Sim
fjsim at buffalo.edu
Wed Sep 30 19:07:43 CEST 2009
Hi,
I am having a problem running the following code which utilizes an eset
derived from hthgu133a CEL files.
>eset
ExpressionSet (storageMode: lockedEnvironment)
assayData: 13834 features, 9 samples
element names: exprs, se.exprs
phenoData
sampleNames: <<masked>> (9 total)
varLabels and varMetadata description:
sample: arbitrary numbering
CEL: NA
phenotype: NA
ShortName: NA
featureData
featureNames: 1007_s_at, 1053_at, ..., AFFX-r2-Ec-bioD-5_at (13834 total)
fvarLabels and fvarMetadata description: none
experimentData: use 'experimentData(object)'
Annotation: hthgu133a
>library(genefilter)
>nsF = nsFilter(eset, var.filter = FALSE)$eset
Error in .checkKeys(value, Lkeys(x), x at ifnotfound) :
value for "AFFX-BioB-3_at" not found'
> sessionInfo()
R version 2.9.1 (2009-06-26)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] grid stats graphics grDevices datasets utils methods
[8] base
other attached packages:
[1] RODBC_1.2-5 bioDist_1.16.0 KernSmooth_2.23-2
[4] hthgu133a.db_2.2.11 RSQLite_0.7-1 DBI_0.2-4
[7] AnnotationDbi_1.6.1 genefilter_1.24.2 RColorBrewer_1.0-2
[10] gplots_2.7.1 caTools_1.9 bitops_1.0-4.1
[13] gdata_2.4.2 gtools_2.6.1 Biobase_2.4.1
[16] rcom_2.2-1 rscproxy_1.3-1
loaded via a namespace (and not attached):
[1] annotate_1.22.0 splines_2.9.1 survival_2.35-4 xtable_1.5-5
The code works as expected using an eset derived from hgu133plus2 CEL files.
Cheers,
Fraser
More information about the Bioconductor
mailing list