在做项目反应理论模拟的时候,我写程序模拟的是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]#设置