Entering edit mode
I want to use BMI as a continuous outcome variable. Not to divide in groups. To do this, I use limma. But I am not sure limma is suitable to directly use continuous variable as outcome based on no divide groups.
If i cann't do this, and need consider for eaxmple linear regression model, what input data should i put in linear regression model? just use CPM data or use calcNormFactors() like limma package do, to adjust count then transform to CPM?
My code is as below:
design <- model.matrix(~ BMI+ age + sex + activities +
SBP + TG + smoking + HDL + Glu + WBC + Batch))
rownames(design) = colnames(count)
y <- voom(d, design, plot = F)
fit <- lmFit(y, design)
head(coef(fit))
tmp <- contrasts.fit(fit, coef = 2) # test "pH" coefficient
tmp <- eBayes(fit)
summary(decideTests(fit))
top.table <- topTable(tmp, sort.by = "P", n = Inf, adjust.method = "BH", coef = 2)
head(top.table, 20)
length(which(top.table$adj.P.Val < 0.05))))
rownames(design) = colnames(count)
y <- voom(d, design, plot = F)
fit <- lmFit(y, design)
head(coef(fit))
tmp <- contrasts.fit(fit, coef = 2) # test "pH" coefficient
tmp <- eBayes(fit)
summary(decideTests(fit))
top.table <- topTable(tmp, sort.by = "P", n = Inf, adjust.method = "BH", coef = 2)
head(top.table, 20)
length(which(top.table$adj.P.Val < 0.05))
sessionInfo( )
Thanks for your reply, because my outcome is a continuous variable, If i want do use linear regression model, what input data should i put in linear regression model? just use CPM data or use calcNormFactors() like limma package do, to adjust count then transform to CPM or use vst in DESeq2 to adjust that?
Sorry, I can't make sense of what you are saying.
You say that BMI is an outcome variable, but your code shows that BMI is a covariate in the linear model, so it is a predictor and not an outcome variable.
You are using voom so you must be working with sequence read counts, which means you do not have a continuous outcome variable.
What exactly is the problem?
When i try to just use linear regression model, the code is as below, BMI is my outcome and its a continious variable
That doesn't have anything to do with limma.
May I add an in-between point of view? For each gene, the outcome is the reads count in the samples of the experiment. The design matrix specifies the covariates associated to the samples. The covariates are supposed to affect the reads count of some genes. If I simplified, I suspect that what you want to achieve is :
design <- model.matrix(~BMI)in order to find genes for which the coefficient associated to BMI is statistically significant. You call BMI the outcome but it is a covariate when modelling the reads count. Hope this shortcut will help.yes ,but i don't want to use limma, try to use linear regression or negative biominal to caculate significant genes, continous variable BMI is my outcome, i am not sure which i should use, cpm or normalization count , or TPM
yes, the follow question isn't related to limma
If you aren't asking about Bioconductor packages, you are off-topic for this site. You might try over on biostars.org instead.