Dear all,
I have a question concerning the RNA-seq experiment.
My file targets is as follows:
Filename Group Time
1 A2780_1 Control 0
2 A2780_2 Control 0
3 4PTX_1 Treat 1
4 4PTX_2 Treat 1
5 8PTX_1 Treat 2
6 8PTX_2 Treat 2
7 16PTX_1 Treat 3
8 16PTX_2 Treat 3
9 32PTX_1 Treat 4
10 32PTX_2 Treat 4
11 64PTX_1 Treat 5
12 64PTX_2 Treat 5
13 128PTX_1 Treat 6
14 128PTX_2 Treat 6
and I'm trying to find genes that are significant across these time points and groups. Then I'm using code from limma guide for time course experiment with many time points and at the end I'm getting some warnings. Could someone help me what is wrong in my code...
> X <- ns(targets$Time, df=5)
> Group <- factor(targets$Group)
> design <- model.matrix(~Group*X)
> x <- DGEList(counts=counts$counts[,c(9,41,13:17,45:49,65,66)],
+ genes=counts$annotation[,c("GeneID","Length")])
> isexpr <- rowSums(cpm(x) > 5) >= 2
> x <- x[isexpr,]
> y <- voom(x, design)
Coefficients not estimable: GroupTreat:X1 GroupTreat:X2 GroupTreat:X3 GroupTreat:X4 GroupTreat:X5
Partial NA coefficients for 11610 probe(s)
> fit <- lmFit(y, design)
Coefficients not estimable: GroupTreat:X1 GroupTreat:X2 GroupTreat:X3 GroupTreat:X4 GroupTreat:X5
Partial NA coefficients for 11610 probe(s)
fit <- eBayes(fit)
Estimation of var.prior failed - set to default value
> topTableF(fit, adjust="BH")
F-statistics not found in fit
I will be grateful for any help.
Best regards, Joanna