Re: [R] lme question

From: Dieter Menne <dieter.menne_at_menne-biomed.de>
Date: Thu 09 Feb 2006 - 00:33:48 EST

Mahdi Osman <m_osm <at> gmx.net> writes:

>
> Hi list,
>
> I am fitting microarray data (intensity) model using the lme package in R
> environment. I have 5 fixed variables in the model. One of the fixed
> variables is genes. I am trying to get p-values for different genes. But I
> am getting only one p-value for all genes together. I can get a list of
> p-value when I run lm. Why can't this work in lme?
>
> I was wondering if some one can help me solve this problem. That is getting
> a list of p-value for each gene in the model using the lme.
>

The bioconductor list might be a better place for this question, and an example would haven been welcome. So my only guess is: could it be that you forgot to make genes a factor variable? When genes is a number, you get a linear regression which is probably not what you want.

Dieter



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Feb 09 01:08:41 2006

This archive was generated by hypermail 2.1.8 : Fri 10 Feb 2006 - 21:35:16 EST