是否可以在elrm包中设置种子?

是否可以在elrm包中设置种子?,r,logistic-regression,seeding,R,Logistic Regression,Seeding,我对r很陌生,正在使用elrm包进行精确逻辑回归(ELR)。我希望能够复制我的结果,所以从我的研究来看,似乎我需要为模拟设定一个种子 我的数据集与犯罪再次犯罪有关,与之前的记录特征、标准化评估的风险水平以及我们的计划如何影响再次犯罪率有关。因此,治疗是我主要的自变量,而累犯率(累犯总数/案件总数)是我的因变量 我试过几种不同的方法 首先,我跑了 set.seed(66672) 运行此命令后,我为ELR运行以下代码: m.treatment <- elrm(formula = RECID

我对r很陌生,正在使用elrm包进行精确逻辑回归(ELR)。我希望能够复制我的结果,所以从我的研究来看,似乎我需要为模拟设定一个种子

我的数据集与犯罪再次犯罪有关,与之前的记录特征、标准化评估的风险水平以及我们的计划如何影响再次犯罪率有关。因此,治疗是我主要的自变量,而累犯率(累犯总数/案件总数)是我的因变量

我试过几种不同的方法

首先,我跑了

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));