私はrがとても新しく、正確なロジスティック回帰(ELR)のためにelrmパッケージを使用しています。私は結果を再現することができるようにしたいので、私の研究から、シミュレーションのシードを設定する必要があるようです。elrmパッケージに種を設定することはできますか?
私のデータセットは、過去の記録特性、標準化された評価のリスクレベル、およびプログラムが再犯の割合にどのように影響するかに関する犯罪的再犯と関係があります。だから、治療は私の主要な独立変数であり、再犯率(RECID_sum/case_sum)は私の従属変数です。
私はこれをいくつかの方法で試しました。
まず、私はこれを実行した後
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 <- elrm(formula = RECID_sum/case_sum ~ TREATMENT + FDRUGVSP +
FPROPVSP + FELEVENTVSP + RISK + street_time + TREATMENT, interest =~ TREATMENT,
iter = 100000, dataset = final, burnIn = 100); set.seed(66672)
が、異なる結果に私はそれを実行するたびに受けました。私はこのコードを実行してもエラーは発生しませんでした。
私は一般的にモンテカルロシミュレーションのシードを設定するための研究を行ってきましたが、そのアプローチはここでは当てはまりそうにありません。 elrmパッケージにELRの種を設定することについての情報はありません。私はrでの私の相対的な未熟さのためにそれを見逃す可能性があり、私は私が探しているものを正確に知ることができないかもしれません。通常は私はSPSSを使用しています。これは非常に扱いにくく、残念なことにELRを意味のある方法で実行することはできません。
これは簡単な回答であれば、あなたが提供できるヘルプは大変ありがとうございます。
使用適切な形式のhttp:さらに悪い何
はRNGの初期化がハードコードされているようだということです://stackoverflow.com/help/how-to-ask –