norm_PLN: Runs PLN function from the PLNmodels package and normalized...

In Rmomal/nestor: Network inference from Species counTs with missing actORs

Description

Runs PLN function from the PLNmodels package and normalized the outputs

Usage

norm_PLN(Y, X = NULL, O = NULL)

Arguments

Y	Count dataset (n x p).
X	Matrix of covariates (n x d).
O	Matrix of offsets (n x p).

Value

• MO: normalized means (n x p).

• SO: normalized marginal variances (n x p).

• sigma_O: correlation matrix (p x p) computed from the variance-covariance matrix estimated by PLNmodels, and corresponding to observed variables.

• theta: matrix of the covariates regression coefficients (p x (d+1) including the intercept).

Examples

 n=100 # 100 samples
p=10  # 10 species
r=1 # 1 missing actor
data=generate_missing_data(n=n,p=p,r=r,type="scale-free", plot=TRUE)
X=data.frame(X1=rnorm(n), X2=runif(n))
normPLNfit<-norm_PLN(data$Y,X)
str(normPLNfit)

Rmomal/nestor documentation built on Dec. 16, 2020, 7:54 p.m.