On Thu, May 19, 2011 at 10:58 PM, karena wrote:
> Thank you so much for this reply, Peter. It helps.
>
> I know this is one way to adjust for covariates. However, if what I want is
> to get the 'remaining values' after adjustment. For example, say, 'gene
> expression' value is denoted as 'ge', and
Thanks a lot for the great help, Timothy!
K
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-do-covariate-adjustment-in-R-tp3537463p3538962.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mail
Dear Karena,
x = 1:100
y = rnorm(100)
fit = lm(x~y)
# what properties does a fit have?
names(fit)
# [1] "coefficients" "residuals" "effects" "rank"
"fitted.values" "assign""qr"
# [8] "df.residual" "xlevels" "call" "terms" "model"
Thank you so much for this reply, Peter. It helps.
I know this is one way to adjust for covariates. However, if what I want is
to get the 'remaining values' after adjustment. For example, say, 'gene
expression' value is denoted as 'ge', and for each gene,
ge=a*age+b*sex+c*per_se
My question is:
On Thu, May 19, 2011 at 7:21 PM, karena wrote:
> Hi, I have a question about how to do covariate adjustment.
>
> I have two sets of 'gene expression' data. They are from two different
> tissue types, 'liver' and 'brain', respectively.
> The purpose of my analysis is to compare the pattern of the w
Hi, I have a question about how to do covariate adjustment.
I have two sets of 'gene expression' data. They are from two different
tissue types, 'liver' and 'brain', respectively.
The purpose of my analysis is to compare the pattern of the whole genome
'gene expression' between the two tissue typ
6 matches
Mail list logo