Меня искали для генерации кодов с использованием R, но я не могу найти подходящий код для моего исследования.Я остро нуждаюсь в помощи!
Мне нужно сгенерировать данные, которые содержат 20 элементов с 2 измерениями (2 фактора), например, элемент 1-10 = коэффициент 1, элемент 11-20 = коэффициент 2. Элементы должны иметь категориальные значения.ответы, и я сделал каждый пункт имеет 4 категориальные переменные.
Однако, 4 категории не работают четко.Есть ли способ сбалансировать каждую категорию на сюжете?
#two dimensional categorical data simulation
Theta <- rmvnorm(10000, sigma = matrix(c(1, .5, .5, 1), 2)) #correlation of .5
set.seed(12345)
#slope matrix of 20 rows=items, a1=10 factor 1, a2=10 factor 2
a <- matrix(c(rlnorm(10,.2,.3),rep(0,20),rlnorm(10,.2,.3)),20)
#runif 4 categorical answers
d<-matrix(runif((20*4),min=-3,max=3), 20)
d <- t(apply(d, 1, sort, decreasing=TRUE))
#sort since intercepts are ordered
polytomous <- simdata(a, d, 10000, Theta=Theta, itemtype = 'gpcm')
summary(polytomous)
#mirt and plot
model1=mirt(data=polytomous, model=1, method="EM",itemtype="gpcm")
plot(model1,type="trace", theta_lim=c(-4,4))