Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 1 addition & 2 deletions scripts/Analysis/DAS/EDGER.R
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
4 changes: 1 addition & 3 deletions scripts/Analysis/DE/EDGER.R
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
5 changes: 2 additions & 3 deletions scripts/Analysis/DEU/EDGER.R
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
Loading