Hi,
I think I did not understand exactly what you want to do with your simulations, but the way DEXSeq stores the dispersion parameters is the following (demonstrated in the pasilla data package):
Code:
data("pasillaExons", package="pasilla") pasillaExons <- fitDispersionFunction(estimateDispersions( estimateSizeFactors( pasillaExons ) )) fData(pasillaExons)$dispBeforeSharing pasillaExons@dispFitCoefs fData(pasillaExons)$dispFitted fData(pasillaExons)$dispersion
Leave a comment: