是否可以在elrm包中设置种子?
我对r很陌生,正在使用elrm包进行精确逻辑回归(ELR)。我希望能够复制我的结果,所以从我的研究来看,似乎我需要为模拟设定一个种子 我的数据集与犯罪再次犯罪有关,与之前的记录特征、标准化评估的风险水平以及我们的计划如何影响再次犯罪率有关。因此,治疗是我主要的自变量,而累犯率(累犯总数/案件总数)是我的因变量 我试过几种不同的方法 首先,我跑了是否可以在elrm包中设置种子?,r,logistic-regression,seeding,R,Logistic Regression,Seeding,我对r很陌生,正在使用elrm包进行精确逻辑回归(ELR)。我希望能够复制我的结果,所以从我的研究来看,似乎我需要为模拟设定一个种子 我的数据集与犯罪再次犯罪有关,与之前的记录特征、标准化评估的风险水平以及我们的计划如何影响再次犯罪率有关。因此,治疗是我主要的自变量,而累犯率(累犯总数/案件总数)是我的因变量 我试过几种不同的方法 首先,我跑了 set.seed(66672) 运行此命令后,我为ELR运行以下代码: m.treatment <- elrm(formula = RECID
set.seed(66672)
运行此命令后,我为ELR运行以下代码:
m.treatment <- elrm(formula = RECID_sum / case_sum ~ TREATMENT + FDRUGVSP +
FPROPVSP + FELEVENTVSP + RISK + street_time + TREATMENT, interest =~ TREATMENT,
iter=100000, dataset=final, burnIn=100)
m.treatment不幸的是,我不这么认为。深入研究显示底层代码正在使用系统:
<>这本身有点麻烦,因为R的内置随机数生成器(可以从C++代码中访问)是.<
更糟糕的是,RNG的初始化似乎是硬编码的:
## binomMCMC.cpp, line 251
srand((unsigned int)time(NULL));
我认为您在这里的唯一选择是联系维护人员(维护人员(“elrm”)
)…使用正确的格式非常感谢您的帮助!当情感被欣赏时,StackOverflow却反对;如果这个答案有用,你可以投票给它(如果你有足够的声誉),在任何情况下,如果它满意地回答了你的问题,你被鼓励点击复选标记接受它。
## binomMCMC.cpp, line 31:
double runif = (float) rand()/RAND_MAX;
## binomMCMC.cpp, line 251
srand((unsigned int)time(NULL));