drawLMbasic {LMest}

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)

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

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)

Author(s)

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

Examples

## Not run: 

# 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)

## End(Not run)