请选择 进入手机版 | 继续访问电脑版
楼主: Sohnia
2586 13

[问答] 关于项目反应理论程序的问题 [推广有奖]

  • 1关注
  • 6粉丝

硕士生

13%

还不是VIP/贵宾

-

威望
0
论坛币
1893 个
通用积分
3.9948
学术水平
15 点
热心指数
15 点
信用等级
15 点
经验
3318 点
帖子
46
精华
1
在线时间
90 小时
注册时间
2017-1-18
最后登录
2021-4-17

Sohnia 发表于 2018-3-1 21:14:23 |显示全部楼层 |坛友微信交流群
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
在做项目反应理论模拟的时候,我写程序模拟的是a和b基于正态分布,公式是 1/(1 + exp(-1.702 * z))(其中z=a*theta - a*b),所以a是有负值的。但是我用anote和winbugs软件估计的时候,却发现a都是正值。大家遇到过这种情况吗?
下面附上我的进行模拟作答的python程序,有懂的朋友可以看一下。
[code]import numpy as np
import os
import pickle
import warnings
'''
creat by zhu
date:2018/2/23
'''
class bi_simulater(object):
    def __init__(self,n,issave= False,ispickle=False,savepath=None,randseed=100):
        #检查文件夹是否存在,不存在新建一个
        if not os.path.exists(savepath):
            os.makedirs(savepath)[code]#设置
二维码

扫码加我 拉你入群

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

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

关键词:项目反应理论 Winbugs软件 warnings simulate warning

simulater.zip

3.28 KB

python模拟二分作答

Sohnia 发表于 2018-3-1 21:15:12 |显示全部楼层 |坛友微信交流群
代码没法给全这个论坛,真的希望能有个好点的专门为心理测量提供交流的论坛。

使用道具

Sohnia 发表于 2018-3-1 21:20:41 |显示全部楼层 |坛友微信交流群
[code]#设置

使用道具

刘彦楼 发表于 2018-3-6 08:53:38 |显示全部楼层 |坛友微信交流群
软件在估计参数值的时候加约束了?另外你的a模拟用正态不合适

使用道具

Sohnia 发表于 2018-3-6 15:28:06 |显示全部楼层 |坛友微信交流群
刘彦楼 发表于 2018-3-6 08:53
软件在估计参数值的时候加约束了?另外你的a模拟用正态不合适
刘老师,您好。所有的参数估计的脚本都是用的手册里的,a基于正态分布是这篇文章《基于神经网络集成的IRT参数估计_余嘉元》中提到的,但他的文章里面用于参数估计的软件是BILOG-MG。

使用道具

作答矩阵那块是刻意的加入了混淆数据?

一般来说直接np.random.binomial(1, p)就可以了

使用道具

Sohnia 发表于 2018-3-7 22:48:07 |显示全部楼层 |坛友微信交流群
chinachrist 发表于 2018-3-7 16:18
作答矩阵那块是刻意的加入了混淆数据?

一般来说直接np.random.binomial(1, p)就可以了
对,你那个语句正解,我是怕日后再用到它改成其他分布容易。作答矩阵的原理是生成n*n维的随机[0,1]矩阵,然后作答概率减去这个矩阵,大于等于0答对,小于0答错。

使用道具

Sohnia 发表于 2018-3-7 22:48
对,你那个语句正解,我是怕日后再用到它改成其他分布容易。作答矩阵的原理是生成n*n维的随机[0,1]矩阵, ...
改分布?改二项分布?

使用道具

Sohnia 发表于 2018-3-9 12:28:32 |显示全部楼层 |坛友微信交流群
刘彦楼 发表于 2018-3-6 08:53
软件在估计参数值的时候加约束了?另外你的a模拟用正态不合适
已经解决,他的另一篇文章中a分布是均匀分布[0,2],这个得到了比较好的估计结果。
已有 1 人评分经验 论坛币 收起 理由
刘彦楼 + 100 + 100 鼓励积极发帖讨论

总评分: 经验 + 100  论坛币 + 100   查看全部评分

使用道具

Sohnia 发表于 2018-3-9 12:30:27 |显示全部楼层 |坛友微信交流群
chinachrist 发表于 2018-3-8 13:05
改分布?改二项分布?
已经解决,他的另一篇文章中a分布是均匀分布[0,2],这个得到了比较好的估计结果。得把分布改成均匀分布[0,2]用你的那个参数估计的python程序得到了比较好的估计结果。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 21:55