getFitsMean {LineagePulse} | R Documentation |
Return mean model fits per gene and cell as matrix for chosen model.
getFitsMean(lsMuModel, vecGeneIDs = NULL)
lsMuModel |
(list) Object containing description of gene-wise mean parameter models. |
vecGeneIDs |
(vector of strings) Gene IDs for which mean model fits are to be extracted. |
(numeric matrix genes x cells) Mean 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 mean parameter fits on alternative model: # Use H1 model fits. vecMeanFits <- getFitsMean( lsMuModel = lsMuModelH1(objLP), vecGeneIDs = rownames(lsSimulatedData$counts)[1]) #plot(lsSimulatedData$annot$continuous, vecMeanFits)