getting residuals from dream

I would like to get the residuals for the variables I used in linear mixed effect model I used in dream. The variable I am interested in is continuous if that matters. I'm unsure how to do this or if it is even possible to do with dream? Here is where I am so far:

cpm.tmm.logno <- cpm(tmm.normfact.out, log=FALSE, normalized.lib.sizes = TRUE) 

form.dream <- ~ Disease + BMI  + (1|Batch) + RIN + (1|Sex) 

vobj <- voomWithDreamWeights(cpm.tmm.logno,   form.dream,   demos)

fitmm <- dream(vobj, form.dream,  demos)

toptable.out <- topTable(fitmm, coef=2, number=nrow(y), adjust.method = 'BH',p.value = 1.0)

If I wanted to then see the residuals for BMI, is there a way I could do that from this point?

dream limma residuals

• 880 views • link

written 5.3 years ago by RNAseqer

0 answers

No answers yet.

More posts like this

LIMMA: where to find the ANOVA?
written 3.3 years ago by Sebastian Hesse

I repeatedly read that LIMMA is performing an ANOVA analysis, apparently always in the background, and also that the ANOVA results are obtainable. Unfortunately, though, …
limma - voom with TPM data
written 4.5 years ago by AlexStar •

Hi, I'm doing a differential expression analysis to RNA-seq data with limma - voom. Unfortunately, I do not have acceess to the raw counts, just …
heatmaps for limma-voom DE genes
written 4.5 years ago by AlexStar •

Hi, I concducted a differential expression analysis with limma - voom for RNA-seq data. Starting with a raw RNA-seq expression data, this is my code: …
DMRcate: extractRanges error
written 5.4 years ago by macsuejacques •

Hello, I am trying to conduct a DMR analyses using the DMRcate package, and all seems to work well until the point where I attempt …
Absence of significant values in the analyzes with GEO2R
written 5.7 years ago by RafaelMP

Hi! I intend to use different array experiments, stored in GEO, to conduct a reanalysis. As a starting point, I used the GEO2R code, but …
No significant DEG: A request to double check my commands for limma.
written 7.3 years ago by RNAseqer

Hello everyone, I am currently using limma to look for Differentially Expressed Genes (DEG) in a small dataset (78 samples) Log2 Quantile Normalized gene expression …
Differential Methylation Analysis using continuous outcome
written 7.6 years ago by daniellemo •

Hi there, I am trying to do a differential methylation analysis of probes in R using "days to death" as my clinical factor. I have …
lesser fold change value using RMA background correction
written 7.9 years ago by Paul

I was trying to do background correction using RMA for my GEO dataset, since the dataset was produced using Affymetrix array (GPL570). And preprocessing using …
Not able to get DEG from FireBrowse data using R
written 8.4 years ago by lawarde.ankita1 •

Hello everyone, I'm working on LUAD (Lung adenocarcinoma) data from FireBrowse. I want to find Differentially expressed genes according to the stage of the tumor. …
Single-color Agilent array analyzing in R
written 8.7 years ago by Leite

Hello everyone, I'm working with data that was performed using Agilent-014850 Whole Human Genome Microarray 4x44K G4112F arrays and the One Color Quick Amp Labeling …