我正在尝试使用 R 中的以下脚本对 9 x 484995115 矩阵运行条件逻辑回归
library(splines)
library(survival)
All <- read.table(pipe("cut -d';' -f1,2,4-9 /scratch/rvonscha/all.txt"), header=TRUE,
sep=";", na.strings="NA", dec=",", strip.white=TRUE)
model <- clogit(Alliance ~ OVB + BVB + CVC + EarlyStage + AvgVCSize + NumberVC +
strata(Strata), data=All)
summary(model)
但是在具有 256GB 内存的节点上计算 2 小时后,程序停止并显示以下错误消息:
*** caught segfault ***
address 0x7f687db38700, cause 'memory not mapped'
Traceback:
1: model.matrix.default(Terms2, m)
2: model.matrix(Terms2, m)
3: coxph(formula = Surv(rep(1, 484995115L), Alliance) ~ OVB + BVB + CVC +
EarlyStage + AvgVCSize + NumberVC + strata(Strata), method = "exact")
4: eval(expr, envir, enclos)
5: eval(coxcall, sys.frame(sys.parent()))
6: clogit(Alliance ~ OVB + BVB + CVC + EarlyStage + AvgVCSize + NumberVC +
strata(Strata))
aborting ...
这个错误的原因是什么?在计算期间从未超过内存容量。我可以在代码中进行哪些更改以避免此错误?
我真的很感激任何帮助来解决这个问题。
非常感谢您提前。