getFitsDispersion {LineagePulse} | R Documentation |
Return dispersion model fits per gene and cell as matrix for chosen model.
getFitsDispersion(lsDispModel, vecGeneIDs = NULL)
lsDispModel |
(list) Object containing description of gene-wise dispersion parameter models. |
vecGeneIDs |
(vector of strings) Gene IDs for which dispersion model fits are to be extracted. |
(numeric matrix genes x cells) Dispersion parameter fits.
David Sebastian Fischer
lsSimulatedData <- simulateContinuousDataSet( scaNCells = 20, scaNConst = 2, scaNLin = 2, scaNImp = 2, scaMumax = 100, scaSDMuAmplitude = 3, vecNormConstExternal=NULL, vecDispExternal=rep(20, 6), vecGeneWiseDropoutRates = rep(0.1, 6)) objLP <- runLineagePulse( counts = lsSimulatedData$counts, dfAnnotation = lsSimulatedData$annot, strMuModel = "impulse") # Get dispersion parameter fits on alternative model: # Use H1 model fits. vecDispersionFits <- getFitsDispersion( lsDispModel = lsDispModelH1(objLP), vecGeneIDs = rownames(lsSimulatedData$counts)[1])