## simulate DINA and DINO data according to a tetrachoric correlation
# define Q-matrix for 4 items and 2 attributes
q.matrix <- matrix(c(1,0,0,1,1,1,1,1), ncol = 2, nrow = 4)
# Slipping parameters
slip <- c(0.2,0.3,0.4,0.3)
# Guessing parameters
guess <- c(0,0.1,0.05,0.2)
set.seed(1567) # fix random numbers
dat1 <- sim.din(N = 200, q.matrix, slip, guess,
# Possession of the attributes with high probability
mean = c(0.5,0.2),
# Possession of the attributes is weakly correlated
Sigma = matrix(c(1,0.2,0.2,1), ncol=2), rule = "DINA")$dat
head(dat1)
set.seed(15367) # fix random numbers
res <- sim.din(N = 200, q.matrix, slip, guess, mean =
c(0.5,0.2), Sigma = matrix(c(1,0.2,0.2,1), ncol=2),
rule = "DINO")
# extract simulated data
dat2 <- res$dat
# extract attribute patterns
head( res$alpha )
## [,1] [,2]
## [1,] 1 1
## [2,] 1 1
## [3,] 1 1
## [4,] 1 1
## [5,] 1 1
## [6,] 1 0
# simulate data based on given attributes
# -> 5 persons with 2 attributes -> see the Q matrix above
alpha <- matrix( c(1,0,1,0,1,1,0,1,1,1) ,
nrow=5,ncol=2, byrow=TRUE )
sim.din( q.matrix = q.matrix , alpha = alpha )
Run the code above in your browser using DataLab