新しい変数を 1 つ追加するたびに、いくつかの Cox モデルを実行しています。
私が実行した主なモデルは
model1 <- coxph(Surv(time1sec, time2sec, event) ~ gain + Buy + Lev + TP + gain:TP + SL + gain:SL + InstrumentID, data)
weekこれは問題なく機能しましたが、取引が開始された時間 (または週) を考慮に入れる固定効果を追加したいと考えました。だから私は以下を実行しました：
model2 <- coxph(Surv(time1sec, time2sec, event) ~ gain + Buy + Lev + TP + gain:TP + SL + gain:SL + InstrumentID + week, data)
このモデルをラップトップで実行すると、R がクラッシュします。クラスターで実行すると、次のエラーが発生します。
*** caught segfault ***
address 0xfffffffc00000000, cause 'memory not mapped'
Traceback:
1: docount(y, x, weight)
2: survConcordance.fit(y, lp, strata, weights)
3: fitter(X, Y, strats, offset, init, control, weights = weights, method = method, row.names(mf))
4: coxph(Surv(time1sec, time2sec, event) ~ gain + Buy + Lev + TP + gain:TP + SL + gain:SL + InstrumentID + week, data)
aborting ...
/var/spool/gridscheduler/execd/node074/job_scripts/27524: line 28: 37078 Segmentation fault (core dumped)
model1 では 8 GB の RAM を使用し、正常に動作しました。model2 では、クラスターで 128 GB に増やしてみましたが、動作しませんでした。
メモリとは関係なく、代わりにweek変数を追加する必要があると思います。このエラーは、週変数が取引の開始日と終了日であるtime1secとに何らかの形で関連付けられているために発生していますか?time2sec
更新グラフに示すように、weekに対して
プロットしました。time1secこれら 2 つの変数が完全に相関していることはわかっていますが、持続時間も time2sec に基づいているため、固定効果として週を含めることはできませんか? それとも、これは time1sec 変数を介してモデルによって考慮されますか?
