LMest (version 3.1.2)

drawLMbasic: Draw samples from the basic LM model

Description

Function that draws samples from the basic LM model.

Usage

drawLMbasic(piv, Pi, Psi, n, est = NULL, format = c("long","matrices"), seed = NULL)

Value

Y

matrix of response configurations unit by unit

S

matrix of distinct response configurations

yv

corresponding vector of frequencies

piv

vector of initial probabilities of the latent Markov chain

Pi

set of transition probabilities matrices (k x k x TT)

Psi

array of conditional response probabitlies (mb x k x r)

n

sample size

est

object of class LMbasic (LMbasic-class)

Arguments

piv

vector of initial probabilities of the latent Markov chain

Pi

set of transition probabilities matrices (k x k x TT)

Psi

array of conditional response probabitlies (mb x k x r)

n

sample size

est

object of class LMbasic (LMbasic-class)

format

character string indicating the format of final responses matrix

seed

an integer value with the random number generator state

Author

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

Examples

if (FALSE) {

# draw a sample for 1000 units and only one response variable
n <- 1000
TT <- 6
k <- 2
r <- 1 #number of response variables
mb <- 3 #maximum number of response categories

piv <- c(0.7,0.3)
Pi <- matrix(c(0.9,0.1,0.1,0.9), k, k)
Pi <- array(Pi, c(k, k, TT))
Pi[,,1] <- 0
Psi <- matrix(c(0.7,0.2,0.1,0.5,0.4,0.1), mb, k)
Psi <- array(Psi, c(mb, k, r))

out <- drawLMbasic(piv, Pi, Psi, n = 1000)


data("data_SRHS_long")
SRHS <- data_SRHS_long[1:2400,]

SRHS$srhs <- 5 - SRHS$srhs

est <- lmest(responsesFormula = srhs ~ NULL,
             index = c("id","t"),
             data = SRHS,
             k = 3)

out1 <- drawLMbasic(est = est, format = "matrices", seed = 4321, n = 100)
}