fit_poilog: Composite maximum likelihood estimation for Poisson...

In ChristofSeiler/cytoeffect: Regression Models with Multivariate Outcomes for Mass Cytometry Experiments

Composite maximum likelihood estimation for Poisson log-normal model

Description

fit_poilog uses composite maximum likelihood estimation to fit the mean and the covariance of a Poisson log-normal mixed model.

Usage

fit_poilog(Y, ncores = 1)

Arguments

Y	Matrix with proteins counts with dimensions (number of cells) x (number of markers)
ncores	Number of CPU cores to use

Value

A list containing

beta	estimated mean parameter
Sigma	estimated covariance parameter

Examples

set.seed(1)
df = simulate_data(n_cells = 10)
str(df)
Y = dplyr::select(df, names(df)[3:ncol(df)])
fit_poilog(Y)

ChristofSeiler/cytoeffect documentation built on Jan. 11, 2023, 1:05 p.m.