我一直在寻找使用R生成代码的方法，但是找不到适合我的学习的代码。 我迫切需要帮助！
我必须生成包含2个维度（2个因数）的20个项目的数据，例如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))