[BioC] class label

Martin Morgan mtmorgan at fhcrc.org
Sat Dec 18 17:02:23 CET 2010


On 12/18/2010 07:56 AM, Peevi Ijkl wrote:
> Hello 
> ok my phenodata(target.txt) looks something like this
> 
> 1.CEL  A
> 1.CEL  A
> 3.CEL  A    
> 4.CEL  H
> 5.CEL  H
> 6.CEL  H
> 
>  i loaded the above into my R using 
> library(affy)
> library(simpleaffy)
> eset<-read.AnnotatedDataFrame("target.txt", header=T, row.names=1, sep="\t" )
> 
> to view my data:
> pData(eset)
> 
> then i did normalization on it
> 
> Rmaeset<-rma(eset)
> 
> now i want to perform multtest..but first i have to create class label..can 
> anyone tell me how do i do tht??

Hi peevi -- You need to be more helpful to us, and stick with your
original question. Likely the problem is that your pd$Diabetic variable
is not what you think it is, e.g., a factor when you're expecting it to
be a character, or of the wrong length, or ???. But how are we to know?
You have to show it to us. And you can't really say 'multtest doesn't
work', because we know that

  example(multtest)

*does* work. And since even simple emails like this can take five
minutes to write, you have to make it easy for us to see what your
problem might be, so no sense giving us a dozen lines of code -- reduce
your real problem to something that shows what the issue is, in just a
few lines, and nothing more.


Martin

> 
> peevi
> 
> 
> 
>       
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor


-- 
Computational Biology
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109

Location: M1-B861
Telephone: 206 667-2793



More information about the Bioconductor mailing list