[BioC] data reduction in "hist" function?

[Wolfgang Huber]

Hi Hinnerk,

what's 

 length(which(is.na(exprs(Mydata.rma))))

respectively

 length(which(!is.na(exprs(Mydata.rma))))

???
Best regards

Wolfgang Huber
w.huber@dkfz.de
http://www.dkfz.de/abt0840/whuber
Tel +49-6221-424709
Fax +49-6221-42524709
DKFZ, Division of Molecular Genome Analysis 
69120 Heidelberg, Germany

-----Original Message-----
From: bioconductor-admin@stat.math.ethz.ch
[mailto:bioconductor-admin@stat.math.ethz.ch]On Behalf Of Hinnerk Boriss
Sent: Thursday, November 21, 2002 11:21 AM
To: 'Bioconductor'
Subject: [BioC] data reduction in "hist" function?


Hi all!

Can someone explain to my why the hist function applied to the
expression values returns less values, i.e. less genes, than there are
on the hu68 chip?  

Here my code fragments:

> Mydata <- ReadAffy(CDFfile="hu6800.cdf",...)
> Mydata.rma <- express(Mydata)
> dim(exprs(Mydata.rma))
[1] 7129   18
# ok, here we have 7129 expression values

> Mydata.rma.hist <- hist(exprs(Mydata.rma),,plot=F)
> sum(Mydata.rma.hist$counts)
[1] 42774
# here we are down to 42774 genes

Cheers,
Hinnerk

_______________________________________________
Bioconductor mailing list
Bioconductor@stat.math.ethz.ch
http://www.stat.math.ethz.ch/mailman/listinfo/bioconductor