[BioC] Problems using text to subset array information from an expression set
Benilton Carvalho
bcarvalh at jhsph.edu
Tue Apr 4 20:31:57 CEST 2006
isn't
exprs(alldata)["1007_s_at",]
exprs(alldata)[, "AA100.CEL"]
what you want?
b
On Tue, 4 Apr 2006, Jeff Lande wrote:
> I have an odd problem that I cannot seem to figure out.
>
> I have a set of CEL files in a directory, which I read using the ReadAffy()
> command. Then I run the rma command to preprocess.
>
>> Data <- ReadAffy()
>> alldata <- rma(Data)
>
> I've done this many times before without problems. However, when I try to
> use text instead of numbers for subscripting, I get an error.
>
> For example, I am able to access data from the first row and column using
> numeric subscripts
>
>> alldata[1,1]
> Expression Set (exprSet) with
> 1 genes
> 1 samples
> phenoData object with 1 variables and 1 cases
> varLabels
> sample: arbitrary numbering
>
> but using text for either subscript, I get an error.
>
>> alldata["1007_s_at",]
> Error in alldata["1007_s_at", ] : no 'dimnames' attribute for array
>> alldata[,"AA100.CEL"]
> Error in alldata[, "AA100.CEL"] : no 'dimnames' attribute for array
>
> I actually went through what I think was the same process last week (and
> many times previously) and had no problems, so I'm stumped.
>
> Here is my session information:
>
>> sessionInfo()
> R version 2.1.0, 2005-04-18, ia64-unknown-linux-gnu
>
> attached base packages:
> [1] "tools" "methods" "stats" "graphics" "grDevices" "utils"
> [7] "datasets" "base"
>
> other attached packages:
> hgu133acdf affy reposTools Biobase
> "1.4.3" "1.6.7" "1.5.19" "1.5.12"
>
> I must be missing something obvious, but I just can't figure out what is
> going wrong. Does anyone have insight into this problem?
>
> Jeff Lande
> Post-Doctoral Associate
> University of Minnesota
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>
More information about the Bioconductor
mailing list