getting differential expressed genes between tissues, species and time series
1
0
Entering edit mode
5.2 years ago
A ★ 3.9k

hi,

I have a microarray experiment (GSE63706) with 6 time points, tissues, three species and two replication. i am going retrieve genes with fold change more than two. 

may someone tell me a design for limma because i got confused as already..

I have a design for treatment and normal in fout time points but i can adapt for my case

the design is like below

 

X <- ns(targets$time.h, df= 2)

Group <- factor(targets$treatment)
design <- model.matrix(~Group*X)
colnames(design) <- c("Intercept", "Groupheat","Control_X1", "Control_X2","Heated_X1", "Heated_X2")
design
fit <- lmFit(data, design)
fit <- eBayes(fit)



R gene • 1.2k views
ADD COMMENT
1
Entering edit mode
5.2 years ago

You'll need to make a contrast matrix to design the actual comparisons you make. Check out the limma users guide on making contrast matrices. I'd also very strongly suggest that you read up on interaction models, if you're not from a statistical background then it can get very confusing, very quickly. Check out section 9.5.4 of the limma users guide, which gives an example of classic interaction models.

ADD COMMENT
0
Entering edit mode

Thank you, i was reading that is why i got totally confused

ADD REPLY

Login before adding your answer.

Traffic: 2025 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