ExpectedValue_EstParmFunc: Obtaining the value of the Dirichlet parameters, the expected...

In CoDaLoMic: Compositional Models to Longitudinal Microbiome Data

Obtaining the value of the Dirichlet parameters, the expected value and the variance.

Description

This function calculates the value of the Dirichlet parameters, the expected value and the variance for the Dirich-gLV model.

Usage

ExpectedValue_EstParmFunc(Param.Estimates, especie)

Arguments

Param.Estimates	Vector with the estimated parameters. This value is the output of the Estimate_Param_EstParmFunc function.
especie	Matrix that contains at row i the bacterial taxa of bacteria i at all time points. The bacteria placed in the last row of this matrix is the one used as reference in the alr transformation that the model applies.

Value

Returns a list with:

• Dirichlet.Param: Matrix. Matrix that contains at row i the Dirichlet parameter of the bacteria i at all time points.

• Expected.Value: Matrix. Matrix that contains at row i the expected value of the bacteria i at all time points. The bacterias are placed at the same orden than in especies.

• Variance.Value: Matrix. Matrix that contains at row i the variance of the bacteria i at all time points. The bacterias are placed at the same orden than in especies.

References

Creus-Martí, I. and Moya, A. and Santonja, F. J. (2018). A Statistical Model with a Lotka-Volterra Structure for Microbiota Data. Lucas Jodar, Juan Carlos Cortes and Luis Acedo, Modelling or engineering and human behavior 2018, Instituto Universitario de Matematica Multidisciplinar. ISBN: 978-84-09-07541-6

Examples

especie1=cbind(c(0.5,0.3,0.2), c(0.1,0.3,0.6))
tau1=0.4
parms1= cbind(c(0.1,0.2),c(-0.2,0.1),c(0.3,0.2))
parms11=c(tau1,as.vector( parms1))

ExpectedValue_EstParmFunc(parms11,especie1)

CoDaLoMic documentation built on April 12, 2025, 2:18 a.m.