 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
计算机自适应化考试中预查项目校准的方法与比较
下载积分：
内容提示：计算机自适应化考试中预查项目校准的方法与比较,对比,校准,项目,校准方法,方法,计算机化,自适应考试,计算机化自适应,考试中的,自适应,计算机,考试的,方法的比较,恋爱考试中字,考试中,考试中的作文,在体育考试中,在某次考试中,在一次考试中
文档格式：PDF|
浏览次数：1|
上传日期： 13:07:17|
文档星级：
该用户还上传了这些文档
官方公共微信
下载文档:计算机自适应化考试中预查项目校准的方法与比较.PDF 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
计算机自适应化考试中预查项目校准的方法与比较
下载积分：
内容提示：计算机自适应化考试中预查项目校准的方法与比较,某学生在n次考试中,在某次考试中,在一次考试中,一次考试中,在体育考试中对部分,驾驶证考试中的路考,张含韵 考试中,银行考试中宏观经济
文档格式：PDF|
浏览次数：1|
上传日期： 04:41:43|
文档星级：
该用户还上传了这些文档
官方公共微信
下载文档:计算机自适应化考试中预查项目校准的方法与比较.PDF[精品文档]基于COS方法的参数估计研究
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
[精品文档]基于COS方法的参数估计研究
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
（概率论与数理统计专业论文）计算机自适应化考试中预查项目校准的方法与比较
下载积分：
内容提示：（概率论与数理统计专业论文）计算机自适应化考试中预查项目校准的方法与比较
文档格式：PDF|
浏览次数：1|
上传日期： 13:45:46|
文档星级：
该用户还上传了这些文档
官方公共微信
下载文档:（概率论与数理统计专业论文）计算机自适应化考试中预查项目校准的方法与比较.PDF查看: 233|回复: 14
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
本帖最后由 超级大菜鸟 于
17:17 编辑
我编了一个单参数IRT的联合极大释然估计函数JMLE.1P（），但是报错：Error in JMLE.1P(data1$U_Matrix) : object 'theta.temp' not found&&。其中data1$U_Matrix是作答矩阵。请问我是哪里出错了呢？谢谢！求大家指点！其中CMLE.theta.1P（）和CMLE.beta.1P()这两个函数是我已经编好的能力参数和难度参数的条件极大似然估计函数。JMLE.1P()会调用上面两个函数。#1PLM/Rasch moedel JMLE
#作答矩阵需事先排除满分和零分的被试数据，以及全对和全错的项目数据
if (!exists(&CMLE.theta.1P&)) source(&CMLE.theta.1P.R&)
if (!exists(&CMLE.beta.1P&))&&source(&CMLE.beta.1P.R&)
JMLE.1P&- function(Resp.Matrix,n2=20,tol2=0.005){
&&N&- nrow(Resp.Matrix)
&&m&- ncol(Resp.Matrix)
&&T&- rowSums(Resp.Matrix)
&&mean_T&- mean(T)
&&sd_T&- sd(T)
&&p&- numeric()
&&for (j in 1:m) p[j]&- sum(Resp.Matrix[,j])/N
&&q&- 1-p
&&
&&#取在题目上得分的被试的卷面总分平均分，一共有m道题，mean_Tp就有m个元素
&&mean_Tp&- numeric()
&&for (j in 1:m) mean_Tp[j]&- sum(Resp.Matrix[,j]*T)/sum(Resp.Matrix[,j])
&&
&&#在题目上失分的被试的卷面总分平均分
&&mean_Tq=numeric()
&&for (j in 1:m) mean_Tq[j]&- sum((1-Resp.Matrix[,j])*T)/sum(1-Resp.Matrix[,j])
&&
&&#设定beta初始值
&&#有题目通过率极端计算对应的概率密度，用于将后面的点二列相关转换为二列相关，进而设定项目参数估计初值
&&y&- dnorm(qnorm(p,lower.tail=F))
&&r&- (mean_Tp-mean_Tq)*(p*q)/(sd_T*y) #二列相关
&&beta.start&- qnorm(p,lower.tail=F)/r
&&beta.temp&- beta.start
&&theta.start&- log(T/(m-T))
&&P.start&- matrix(,N,m)
&&for (i in 1:N) {
& & for(j in 1:m) {
& && &P.start[i,j]&- 1/(1 + exp(-1.7 * (theta.start[i] - beta.start[j])))
& & }
&&}
&&logLik.start&-&&sum(Resp.Matrix*log(P.start)+(1-Resp.Matrix)*log(1-P.start))
&&
&&logLikDiff&-
&&iterN&- 1
&&while (logLikDiff & tol2 & iterN & n2){
& & cmle.theta.temp&- apply(Resp.Matrix,1,CMLE.theta.1P,beta=beta.temp)
& & theta.temp&- numeric()
& & for (i in 1:N) theta.temp[i]&- cmle.theta.temp[[i]]& JMLE.1P(data1$U_Matrix)
Error in JMLE.1P(data1$U_Matrix) : object 'theta.temp' not found
point estimator'
& &
& & cmle.beta.temp&- apply(Resp.Matrix,2,CMLE.beta.1P,theta=theta.temp,k=20)
& & beta.temp&- numeric()
& & for(j in 1:m) beta.temp[j]&- cmle.beta.temp[[j]]& JMLE.1P(data1$U_Matrix)
Error in JMLE.1P(data1$U_Matrix) : object 'theta.temp' not found
point estimator'
& &
& & P.temp&- matrix(,N,m)
& & for (i in 1:N) {
& && &for(j in 1:m) {
& && &&&P.temp[i,j]&- 1/(1 + exp(-1.7 * (theta.temp[i] - beta.temp[j])))
& && &}
& & }
& & logLik.temp&- sum(Resp.Matrix*log(P.temp)+(1-Resp.Matrix)*log(1-P.temp))
& & logLikDiff&- abs(logLik.start-logLik.temp)
& & iterN&- iterN+1
&&}
&&Results&- list(theta.estimator=theta.temp, beta.estimator=beta.temp)
&&return(Results)
17:17:27 上传
支持楼主：、
购买后，论坛将把您花费的资金全部奖励给楼主，以表示您对TA发好贴的支持
载入中......
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
超级大菜鸟 发表于
我编了一个单参数IRT的联合极大释然估计函数JMLE.1P（），但是报错：Error in JMLE.1P(data1$U_Matrix) : o ...自己顶，各位大神帮帮忙啊！！！！
阅读权限26威望0 级论坛币1333 个学术水平2 点热心指数3 点信用等级1 点经验2107 点帖子292精华0在线时间420 小时注册时间最后登录
积分 718, 距离下一级还需 82 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 205 天连续签到: 1 天[LV.7]常住居民III
Resp.Matrix 和 data1$U_Matrix 是同样的吗？
CMLE.theta.1P（）和CMLE.beta.1P(） 这两个函数不在，所以很难运行你的Code 帮你查问题出在哪里
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
calsunny 发表于
Resp.Matrix 和 data1$U_Matrix 是同样的吗？
CMLE.theta.1P（）和CMLE.beta.1P(） 这两个函数不在，所 ...CMLE.beta.1P&- function(resp.vect,D=1.7,theta,beta.start=0,k=20,n1=50,tol1=0.0005,...){
&&beta&- beta.start
&&d1.logLikDiff&-
&&it.log&- matrix(ncol=3,nrow=0)
&&colnames(it.log)&- c(&beta_estimator&,&d1.logLik&,&d2.logLik&)
&&
&&#将某题的作答向量与theta合并成一个数据框
&&A&- data.frame(resp.vect,theta)
&&A&- A[order(A$theta),] #以A中的theta排序
&&A$theta_space&- cut(A$theta,breaks=k)
&&
&&#将theta分成k组后，取每组的中位数代表该组的能力水平
&&med_theta_space&- aggregate(A$theta,list(A$theta_space),median)
&&names(med_theta_space)&- c(&theta_space&,&theta_level&)
&&
&&#B是每个能力水平答对r和答错(f-r)组成的矩阵
&&B&- table(A$resp.vect,A$theta_space)
&&f&- colSums(B)
&&r&- B[2,]
&&
&&iterN&- 1 #iteration number
&&
&&while(d1.logLikDiff & tol1 & iterN& n1) {
& & probs.1&- 1/(1+exp(-D*(med_theta_space$theta_level - beta)))
& & probs.0&- 1-probs.1
& &
& & #First derivarive
& & d1.logLik&- sum(D*(f*probs.1-r))
& &
& & #Second derivative
& & d2.logLik&- -sum(D*f*probs.1*probs.0)
& & it.log&- rbind(it.log, c(beta, d1.logLik, d2.logLik))
& &
& & #Newton Raphson
& & beta&- beta - (d1.logLik/d2.logLik)
& &
& & #Convergence&&critieria& &
& & if(nrow(it.log) & 1) {
& && &d1.logLikDiff&- abs(d1.logLik-it.log[(nrow(it.log) - 1 ), 2])
& & }
& & iterN&- iterN+1
&&}
&&
&&
&&#Output of the CMLE.beta funtion
&&Results&- list(beta, iterN, it.log)
&&names(Results) &- c(&point estimator&, &iteration Num&,&&&iteration log&)
&&return(Results)
}复制代码#单个theta的条件极大似然估计；作答矩阵需事先排除满分和零分的被试数据，以及全对和全错的项目数据
CMLE.theta.1P&-function(resp.vect,D=1.7,beta,theta.start=0,n1=50,tol1=0.0005,...){
&&theta&- theta.start
&&d1.logLikDiff&-
&&it.log&- matrix(ncol=3,nrow=0)
&&colnames(it.log)&- c(&theta_estimator&,&d1.logLik&,&d2.logLik&)
&&iterN=1 #iteration number
&&
&&while(d1.logLikDiff&tol1 & iterN&n1){
& & probs.1&- 1/(1+exp(-D*(theta-beta)))
& & probs.0&- 1-probs.1
& &
& & #First derivative
& & d1.logLik&- sum(D * (resp.vect - probs.1))
& &
& & #Second derivative
& & d2.logLik&- -D*sum(probs.1 * probs.0)
& & it.log&- rbind(it.log,c(theta,d1.logLik,d2.logLik))
& &
& & #Newton-Raphson
& & theta&- theta-(d1.logLik / d2.logLik)
& &
& & #Convergence critieria
& & if(nrow(it.log) & 1){
& && &d1.logLikDiff&- abs(d1.logLik - it.log[(nrow(it.log) - 1), 2])
& & }
& & iterN&- iterN+1
&&}
&&
&&#Standard error
&&SE.theta&- 1/sqrt(-1*d2.logLik)
&&
&&#Output of the CMLE.theta function
&&Results &- list(theta,SE.theta,iterN,it.log)
&&names(Results) &- c(&point estimator&,&standard error&, &iteration Num&, &iteration log&)
&&return(Results)
}复制代码
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
calsunny 发表于
Resp.Matrix 和 data1$U_Matrix 是同样的吗？
CMLE.theta.1P（）和CMLE.beta.1P(） 这两个函数不在，所 ...Resp.Matrix是参数，data1$U_Matrix是具体的作答矩阵。麻烦你帮我看看，谢谢了！
阅读权限26威望0 级论坛币1333 个学术水平2 点热心指数3 点信用等级1 点经验2107 点帖子292精华0在线时间420 小时注册时间最后登录
积分 718, 距离下一级还需 82 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 205 天连续签到: 1 天[LV.7]常住居民III
本帖最后由 calsunny 于
03:43 编辑 超级大菜鸟 发表于
Resp.Matrix是参数，data1$U_Matrix是具体的作答矩阵。麻烦你帮我看看，谢谢了！OK. I will. But it takes time.
Do you have any particular reasons for not using Winstep or Bilog but writing the functions by yourself?
阅读权限26威望0 级论坛币1333 个学术水平2 点热心指数3 点信用等级1 点经验2107 点帖子292精华0在线时间420 小时注册时间最后登录
积分 718, 距离下一级还需 82 积分
权限: 自定义头衔, 签名中使用图片
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡, 匿名卡, 抢沙发下一级可获得
权限: 隐身
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 205 天连续签到: 1 天[LV.7]常住居民III
Is Resp.Matrix 是具体的作答矩阵 while data1$U_Matrix is 参数, right?
Could you double check where data1$U_Matrix is from?
I believe the problem is here:
theta.temp[i]&- cmle.theta.temp[[i]]& JMLE.1P(data1$U_Matrix)
I did not get where data1$U_Matrix is from.
This might be the reason you got the error message:
Error in JMLE.1P(Resp.Matrix = resp.vect, n2 = 20, tol2 = 0.005) :
&&object 'theta.temp' not found
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
calsunny 发表于
Is Resp.Matrix 是具体的作答矩阵 while data1$U_Matrix is 参数, right?
Could you double check where ...data1$U_Matrix是我事先就模拟出来的作答矩阵。
之前忘了发模拟这个矩阵的函数。Resp.Matrix.1P&- function(N,m,D=1.7) {
&&#True value of Difficulty
&&b_T&- rnorm(m)
&&for (j in 1:m) {
& & while (b_T[j] & (-3) | b_T [j] & 3) b_T[j]&- rnorm(1)
&&}
&&
&&#True value of Theta
&&theta_T&- rnorm(N)
&&for (i in 1:N) {
& & while (theta_T[i] & (-3) | theta_T[i] & 3) theta_T[i]&- rnorm(1)
&&}
&&
&&#Rand_Matrix: uniform distribution u~(0,1)
&&#P_Matrix: probability matrix
&&#U_Mattix: generated response matrix
&&Rand_Matrix&- matrix(runif(N*m,0,1),nrow=N,ncol=m)
&&P_Matrix&- matrix(0,N,m)
&&U_Matrix&- matrix(0,N,m)
&&for (i in 1:N) {
& & for(j in 1:m) {
& && &P_Matrix[i,j]&- 1/(1 + exp(-D * (theta_T[i] - b_T[j])))
& && &if (Rand_Matrix[i,j] &= P_Matrix[i,j]) U_Matrix[i,j]&- 1
& && &else U_Matrix[i,j]&- 0&&
& & }
&&}&&
&&Results&- list(N&-N,m&-m,b_T&-b_T,theta_T&-theta_T,P_Matrix&-P_Matrix,U_Matrix&-U_Matrix)
&&names(Results)&- c(&N&,&m&,&b_T&,&theta_T&,&P_Matrix&,&U_Matrix&)
&&return(Results)
}复制代码data1就是这个函数的结果data1&- Resp.Matrix.1P()复制代码
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
calsunny 发表于
OK. I will. But it takes time.
Do you have any particular reasons for not using Winstep or Bi ...我是想试着练习一下，方便以后进一步学习IRT。
阅读权限20威望0 级论坛币59 个学术水平0 点热心指数0 点信用等级0 点经验1082 点帖子31精华0在线时间49 小时注册时间最后登录
积分 105, 距离下一级还需 40 积分
权限: 自定义头衔
道具: 彩虹炫, 雷达卡, 热点灯, 雷鸣之声, 涂鸦板, 金钱卡, 显身卡下一级可获得
道具: 匿名卡
购买后可立即获得
权限: 隐身
道具: 金钱卡, 雷鸣之声, 彩虹炫, 雷达卡, 涂鸦板, 热点灯
开心签到天数: 13 天连续签到: 1 天[LV.3]偶尔看看II
calsunny 发表于
OK. I will. But it takes time.
Do you have any particular reasons for not using Winstep or Bi ...我是想试着练习一下，方便以后进一步学习IRT。
论坛好贴推荐}