[BioC] localization of mm values in affybatch exprs matrix
Karin Lagesen
karin.lagesen at medisin.uio.no
Thu Jan 25 16:26:47 CET 2007
I have a custom affy chip that I read into R using ReadAffy():
> newdata = ReadAffy()
> newdata
AffyBatch object
size of arrays=754x754 features (17777 kb)
cdf=E_colia530222N (11378 affyids)
number of samples=4
number of genes=11378
annotation=ecolia530222n
>
I now want to look at different values in this object.
For instance, some pm values:
> pm(newdata)[1:5,]
chip1.CEL chip2.CEL chip3.CEL chip4.CEL
[1,] 1855.0 2180.8 1444.0 2932.0
[2,] 2812.0 3451.0 2276.5 3406.0
[3,] 4162.3 4301.0 2996.0 5088.0
[4,] 1608.5 1758.0 1123.0 1987.0
[5,] 2290.0 3189.0 2474.5 2838.3
>
I now also look at the values in the affybatch exprs matrix:
> newdata at exprs[1:5,]
chip1.CEL chip2.CEL chip3.CEL chip4.CEL
[1,] 942.0 776.0 281 1475
[2,] 24422.0 26071.0 8914 21826
[3,] 1024.5 908.8 227 1594
[4,] 26267.0 27674.0 16199 22104
[5,] 130.0 193.0 168 145
>
I also notice that the dimension of the exprs matrix is such that
there is one column for each chip, and as many rows as there are pm
plus mm values.
Are the first half of rows the pm values, with the mm values
following, or are the pm values every other row with the corresponding
mm value below, or is this set up in some other way? Is there any way
for me to look at a value in the exprs matrix and find out which entry
in the pm/mm value list it is?
TIA,
Karin
--
Karin Lagesen
More information about the Bioconductor
mailing list