Comparing expression of DEG gene in samples
0
0
Entering edit mode
2.5 years ago
Maryam • 0

Hi, I got the DEG table with the limma package in R, I searched for the SNAI1 gene in the DEG table and now I want to compare my gene expression in the samples. I want to know in which samples SNAI1 is up-regulated and in which samples it is reduced. How can I do that? With what code?

 library(edgeR)
library(limma)
dge <- DGEList(rawdata)
keep <- filterByExpr(dge, design = design)
keep2 <- data.frame(keep)
filt <- dge[keep, ,keep.lib.size = F]
norm <- calcNormFactors(filt, method = "TMM")
v <- voom(norm, design = design, plot = T)
E <- v$E
ex <- 2^E
Ex <- log2(ex + 1)


##DEG using edgeR and limma##

fit <- lmFit(Ex, design = design)
contrast <- makeContrasts(cancer-normal, levels = design)
fit2 <- contrasts.fit(fit = fit,contrasts = contrast)
fit2 <- eBayes(fit = fit2)  
diff <- topTable(fit = fit2, number = Inf, sort.by = "P" ) 
DEG limma R gene • 764 views
ADD COMMENT
0
Entering edit mode

have you read the manual of limma/edgeR/DEseq2 etc?

ADD REPLY
0
Entering edit mode

Well I took a quick look at limma and adgeR in doing 'Differentially Expressed Genes' but I did not read carefully. I should read carefully .Thanks

ADD REPLY

Login before adding your answer.

Traffic: 1646 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6