Thanks a ton @dariober.
• 0 views
•
link
written
by
housebie •
How to adjust data for "continuous" covariates ?
We know packages like SVA and ComBat, which can be used to adjust data for batch effects or biological artifacts. However, in most cases, all adjusting variables are treated as factors. So I want to know whether anyone could suggest me any package where in I could adjust for the "continuous" variables/covariates. I have multiple of them to adjust for. Could anyone suggest me an easy-to-use method?
I have also read about "cbcbSEQ" with a modified ComBat function (Post: How To Adjust Data For Confounding Covariates Before Pca?). Has anyone tried it for continuous variables?
• 16,761 views
•
link
updated
by
Ram
•
written
by
housebie •
•
written
by
housebie •
2 answers
So maybe fit a linear model through your response variable as a function of the covariates you want to correct for. Then extract the residuals from the model and use those as your new response variable. E.g.:
## Example data:
n<- 100
dat<- data.frame(y= c(rnorm(n= n, mean= 1), rnorm(n= n, mean= 2)), sex= rep(c('A', 'B'), each= n), age= c(rnorm(n= n, mean= 1), rnorm(n= n, mean= 5)))
## Response y varies with age:
cor.test(dat$y, dat$age)
## Fit model:
lm1<- lm(y ~ sex + age, data= dat)
## Now residuals are not correlated with age:
cor.test(lm1$residuals, dat$age)
## Residuals vs original
plot(lm1$residuals, dat$y)
## ... Now use lm1$residuals instead of dat$y
Does it make sense?
@dariober I just have another question here. In this case, sex is the phenotype and not the covariate we want to correct for. So does this mean that, residuals that we are able to get here, are independent of the phenotype (in this case: sex)? And only depends on the covariates?
So in other words, giving another phenotype here, instead of sex, could have not altered the residuals at all??
• 0 views
•
link
written
by
housebie •
## ... Now use lm1$residuals instead of dat$y
Would you please be able to explain Why this residuals are used as new phenotype and how?
• 0 views
•
link
written
by
Bioinformatics_NewComer
https://bioinformatics.stackexchange.com/questions/9172/how-are-residuals-used-as-new-phenotype
This may be of help to you.
• 0 views
•
link
written
by
Bioinformatics_NewComer
I am not sure if your answer if correct.
Please see
https://stats.stackexchange.com/questions/286850/linear-regression-confounder
• 0 views
•
link
written
by
fansili2013 •
Log in to answer this question.
More posts like this
-
ComBat for Multiple Batch Effects in Methylation Data (Chip, Row, Plate)
written by Eliveri •
-
How to handle Single-Level Batch Variable for Batch Correction in RNA-Seq Data
written by Nelo •
-
How should I adjust covariates in eQTL analysis with an interaction term?
written by maximal_life •
-
How to choose the biological variables if you want to keep all the biological signal from your data…
written by ev97 •
-
Remove batch effects in meta-analysis
written by Morteza Hadizadeh •
-
Advice required for remove batch effect which perfectly collinear with contrast of interest
written by boymin2020 •
-
Batch effects : ComBat or removebatcheffects (limma package) ?
written by lessismore
-
JunctionSeq: using continuous covariates, is it possible?
written by dhibar •
-
Batch ID in TCGA database
written by Na Sed
-
How To Calculate Power For Genome Wide Methylation Array Data, Like Illumina 27K Or 450K?
written by housebie •
In R you can just include the values of the continuous variable in your design matrix
@russ_hyde.. thanks for your reply. But I do not quite understand how do I do that. I do not want to include these covariates in my linear model. I would like to adjust for them separately in my data and then work with the adjusted residual matrix.