>DEGList<-DGEList(counts=cbind(M1$counts,M2$counts,M3$counts,M4$counts,S1$counts,S2$counts,S3$counts,S4$counts), genes=M1$annotation)
> DEGList<-sumTechReps(DEGList)  
> colnames(DEGList) <- c("M1","M2","M3","M4","S1","S2","S3","S4")
> ENSG_Annot <- readLines("Homo_sapiens.gene_info", n = 1)
> ENSG.colname <- unlist(strsplit(substring(ENSG_Annot, 10, 234), ''))
> ENSG <- read.delim("Homo_sapiens.gene_info", skip=1, header=FALSE, stringsAsFactors=FALSE)
> colnames(ENSG) <- ENSG.colname
> m <- match(DEGList$genes$GeneID, ENSG$GeneID)
> DEGList$genes$Chr <- ENSG$chromosome[m]
> DEGList$genes$Symbol <- ENSG$Symbol[m]
> DEGList$genes$Strand <- NULL
> isexpr <- rowSums(cpm(DEGList) > 1) >=3
> DEGList <- DEGList[isexpr,,keep.lib.sizes=FALSE]
> DEGList <- calcNormFactors(DEGList)
> Batch <- factor(c(1,2,3,4,1,2,3,4))
> DEGList$samples$group=c(rep("Malade", 4), rep("Sain", 4))
> Pa <- factor(DEGList$samples$group,levels=c("Malade","Sain"))
> design <- model.matrix(~ Batch + Pa)
> v <- voom(DEGList,design,plot=TRUE)
> fit <- lmFit(v, design)
> fit.de <- eBayes(fit, robust=TRUE)
> topTablefit.de, coef=ncol(fit))
> summarydecideTestsfit.de))
> ex <- diffSplice(fit[,"PaSain"], geneid = "GeneID", exonid = "Start")

Total number of exons:  16772 
Total number of genes:  16772 
Number of genes with 1 exon:  16772 
Mean number of exons in a gene:  1 
Max number of exons in a gene:  1 
Error in diffSplice(fit[, "PaSain"], geneid = "GeneID", exonid = "Start") : 
  No genes with more than one exon