Yeah, they just factorized the colData or whatever the equivalent is in the older version that you're using (I'll have to download that version and have a look).
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
So the first thing to try is to just reset ecs@designColumns to contain factors (if they don't already). This won't solve the original issue, since I suspect that the conversion to a dataframe is occurring after this step. The quickest way to go about this is to just look at the code for estimatelog2FoldChanges (I'll paste it below) and just execute every line manually. Then you'll know exactly what it's processing and why the duplicate row name error is happening (this will also give you the chance to avert that error then).
Code:estimatelog2FoldChanges <- function(ecs, fitExpToVar="condition", denominator="", getOnlyEffects=FALSE, averageOutExpression=TRUE, nCores=1, quiet=FALSE, file="") { stopifnot(is(ecs, "ExonCountSet")) if(any(is.na(sizeFactors(ecs)))){ stop("Please estimate sizeFactors first\n")} if(!fitExpToVar %in% ecs@designColumns){ stop("fitExpToVar parameter is not in the design columns, double check ecs@designColumns")} if(sum(is.na(featureData(ecs)$dispersion))==nrow(counts(ecs))){ stop("No dispersion parameters found, first call function estimateDispersions...\n")} frm <- as.formula(paste("count ~", fitExpToVar, "* exon")) testablegenes <- as.character(unique(fData(ecs)[which(fData(ecs)$testable),]$geneID)) geteffects <- function(geneID){ coefficients <- fitAndArrangeCoefs(ecs, geneID=geneID, frm, balanceExons=TRUE) if( is.null( coefficients ) ){ return(coefficients) } ret <- t(getEffectsForPlotting(coefficients, averageOutExpression=averageOutExpression, groupingVar=fitExpToVar)) rownames(ret) <- paste(geneID, rownames(ret), sep=":") return(ret) } if( nCores > 1 ){ if(!is.loaded("mc_fork", PACKAGE="parallel")){ stop("Please load first parallel package or set parameter nCores to 1...")} alleffects <- parallel:::mclapply( testablegenes, function(x){ geteffects(x) }, mc.cores=nCores ) }else{ alleffects <- lapply( testablegenes, function(x){geteffects(x)}) } names(alleffects) <- testablegenes alleffects <- do.call(rbind, alleffects) alleffects <- vst(exp( alleffects ), ecs) toadd <- matrix(NA, nrow=nrow(ecs), ncol=ncol(alleffects)) rownames(toadd) <- featureNames(ecs) if( getOnlyEffects ){ colnames(toadd) <- colnames(alleffects) toadd[rownames(alleffects), colnames(alleffects)] <- alleffects }else{ if( denominator == "" ){ denominator <- as.character(design(ecs, drop=FALSE)[[fitExpToVar]][1]) } stopifnot( any( colnames(alleffects) %in% denominator ) ) denoCol <- which(colnames(alleffects) == denominator) alleffects <- log2(alleffects / alleffects[,denoCol]) colnames(alleffects) <- sprintf("log2fold(%s/%s)", colnames(alleffects), denominator) colnames(toadd) <- colnames(alleffects) alleffects <- alleffects[,-denoCol, drop=FALSE] toadd <- toadd[,-denoCol, drop=FALSE] toadd[rownames(alleffects), colnames(alleffects)] <- alleffects } fData(ecs) <- cbind(fData(ecs), toadd) ecs }
Comment
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Comment