drawLMlatent: Draw samples from LM model with covariates in the latent...

In LMest: Generalized Latent Markov Models

Draw samples from LM model with covariates in the latent model

Description

Function that draws samples from the LM model with individual covariates with specific parameters.

Usage

drawLMlatent(Psi, Be, Ga, latentFormula, data, index,
             paramLatent = c("multilogit","difflogit"), est = NULL,
             format = c("long","matrices"), fort = TRUE, seed = NULL)

Arguments

Psi	array of conditional response probabilities (mb x k x r)
Be	parameters affecting the logit for the initial probabilities
Ga	parametes affecting the logit for the transition probabilities
latentFormula	a symbolic description of the model to be fitted. Detailed description is given in lmest
data	a data frame in long format, with rows corresponding to observations and columns corresponding to variables, a column corresponding to time occasions and a column containing the unit identifier
index	a character vector with two elements indicating the name of the "id" column as first element and the "time" column as second element
paramLatent	type of parametrization for the transition probabilities ("multilogit" = standard multinomial logit for every row of the transition matrix, "difflogit" = multinomial logit based on the difference between two sets of parameters)
est	object of class LMlatent (LMlatent-class)
format	character string indicating the format of final responses matrix
fort	to use fortran routine when possible (FALSE for not use fortran)
seed	an integer value with the random number generator state

Value

Y	matrix of response configurations
U	matrix containing the sequence of latent states (n x TT)
Psi	array of conditional response probabilities (mb x k x r)
Be	parameters affecting the logit for the initial probabilities
Ga	parametes affecting the logit for the transition probabilities
latentFormula	a symbolic description of the model to be fitted. Detailed description is given in lmest
data	a data frame in long format, with rows corresponding to observations and columns corresponding to variables, a column corresponding to time occasions and a column containing the unit identifier
est	object of class LMlatent (LMlatent-class)

Author(s)

Francesco Bartolucci, Silvia Pandolfi, Fulvia Pennoni, Alessio Farcomeni, Alessio Serafini

Examples

## Not run: 
data(data_SRHS_long)
dataSRHS <- data_SRHS_long
data_SRHS_long$srhs <- 5 - data_SRHS_long$srhs
est  <- lmest(responsesFormula = srhs ~ NULL,
              latentFormula =  ~
              I(gender - 1) +
              I( 0 + (race == 2) + (race == 3)) +
              I(0 + (education == 4)) +
              I(0 + (education == 5)) +
              I(age - 50) + I((age-50)^2/100),
              index = c("id","t"),
              data = data_SRHS_long,
              k = 2,
              paramLatent = "multilogit",
              start = 0)

out <- drawLMlatent(est = est,format = "matrices",seed = 4321)

out1 <- drawLMlatent(Psi = est$Psi,Be = est$Be,Ga = est$Ga,
                     data = data_SRHS_long, index = c("id","t"),
                     latentFormula =  ~
                     I(gender - 1) +
                     I( 0 + (race == 2) + (race == 3)) +
                     I(0 + (education == 4)) +
                     I(0 + (education == 5)) +
                     I(age - 50) + I((age-50)^2/100),
                     paramLatent = "multilogit", format = "matrice",
                     seed = 4321)


## End(Not run)

LMest documentation built on Aug. 27, 2023, 5:06 p.m.