询问一个关于R模拟IRT的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在做IRT模拟的时候，循环不知道为什么不运行
I <- 10
J <- 5
response <- matrix(0, 10, 5)
pij <- function(a,b,theta)
{
a <- rnorm(J, 0.8, 0.04)
a
b <- rnorm(J, 0, 1)
b
theta <- rnorm(I, 0,1)
theta
for( i in 1:I ) {
for( j in 1:J ) {
ptemp <- runif(1)
  pij <- exp(a[j]*(theta[i]-b[j]))/(1+exp(a[j]*(theta[i]-b[j])))
response[i,j]<-ifelse(pij(b=b[j], a=a[j], theta[i]) < ptemp , 0 ,1)
}
}
}
response

急请达人帮助~~~~

谢谢~~~~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：IRT 不知道 IRT

矛盾处不少.
a,b,theta要在function pij前宣告
pij既然是function name
而pij <- exp(a[j]*(theta[i]-b[j]))/
   (1+exp(a[j]*(theta[i]-b[j])))是啥意思?

if 表达式如下
if (条件) 表达式1
或：
if (条件) 表达式1 else 表达式2