# NOT RUN {
set.seed(2342)
bmat <- sim_bmat(n_items = 10, k = 2)$bmat
theta1 <- rnorm(100)
theta2 <- rnorm(100, mean = -1)
dat1 <- sim_data(bmat = bmat, theta = theta1)
dat2 <- sim_data(bmat = bmat, theta = theta2)
# estimate each model with fixed-effects and k = 0
fmp0_1 <- fmp(dat = dat1, k = 0, em = FALSE)
fmp0_2 <- fmp(dat = dat2, k = 0, em = FALSE)
# Stocking-Lord linking
# }
# NOT RUN {
sl_res <- sl_link(bmat1 = fmp0_1$bmat[1:5, ],
bmat2 = fmp0_2$bmat[1:5, ],
k_theta = 0)
hb_res <- hb_link(bmat1 = fmp0_1$bmat[1:5, ],
bmat2 = fmp0_2$bmat[1:5, ],
k_theta = 0)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab