diff --git a/scripts/Analysis/DAS/EDGER.R b/scripts/Analysis/DAS/EDGER.R index ff19817e..f38a7c21 100755 --- a/scripts/Analysis/DAS/EDGER.R +++ b/scripts/Analysis/DAS/EDGER.R @@ -245,8 +245,7 @@ for (compare in comparisons[[1]]) { # } # contrast <- as.numeric(contrast[,1]) - # create files topSpliced by gene, simes and exon method - # sp <- diffSpliceDGE(fit, contrast=contrast, geneid="genes", exonid="exons", verbose=FALSE) + # create files topSpliced by gene, simes and exon method sp <- diffSpliceDGE(fit, geneid = "genes", exonid = "exons", verbose = FALSE) # add comp object to list for image diff --git a/scripts/Analysis/DE/EDGER.R b/scripts/Analysis/DE/EDGER.R index bef0a4c4..8606547e 100755 --- a/scripts/Analysis/DE/EDGER.R +++ b/scripts/Analysis/DE/EDGER.R @@ -318,9 +318,7 @@ for (contrast in comparison[[1]]) { # contrast[which(contrast[,2]==i)]<- plus # } # contrast <- as.numeric(contrast[,1]) - - ## Testing - # qlf <- glmQLFTest(fit, contrast=contrast) ## glm quasi-likelihood-F-Test + AvsB <- makeContrasts(TreatvsUntreat = paste("condition", A, sep = ""), levels = design) qlf <- glmQLFTest(fit, contrast = AvsB, prior.count = 5) ## glm quasi-likelihood-F-Test # add comp object to list for image diff --git a/scripts/Analysis/DEU/EDGER.R b/scripts/Analysis/DEU/EDGER.R index 451fc10c..763a4afa 100755 --- a/scripts/Analysis/DEU/EDGER.R +++ b/scripts/Analysis/DEU/EDGER.R @@ -239,9 +239,8 @@ for (compare in comparisons[[1]]) { # } # contrast <- as.numeric(contrast[,1]) - # quasi-likelihood F-Test - # qlf <- glmQLFTest(fit, contrast = contrast) - qlf <- glmQLFTest(fit, prior.count = 5) + # quasi-likelihood F-Test + qlf <- glmQLFTest(fit) is.de <- decideTests(qlf, p.value = 0.05) # add comp object to list for image