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