检查远程错误时出错 (val):5 个节点产生错误:未找到对象

2024-05-25

我尝试使用并行处理(parLapply)进行 10 倍交叉验证并估计联合模型的模型性能。我试图找出为什么我收到错误消息: “checkForRemoteErrors(val) 中的错误:五个节点产生错误:未找到对象‘Week’”

代码如下:

# Validation using 10-fold CV
    library("parallel")
    set.seed(123)
    V <- 10
    n <- nrow(dfC)
    splits <- split(seq_len(n), sample(rep(seq_len(V), length.out = n)))
    CrossValJM <- function (i) {
        library("JM")
        library("nlme")
        trainingData <- dfL[!dfL$ID %in% i, ]
        trainingData_ID <- trainingData[!duplicated(trainingData$ID), ]
        testingData <- dfL[dfL$ID %in% i, ]

        lmeFit <- lme(DA ~ ns(Week, 2), data = trainingData,
                           random = ~ ns(Week, 2) | ID)
        coxFit <- coxph(Surv(TT_event, Event) ~ Gender * Age, data = 
                           trainingData_ID, 
                             x = TRUE)

        jointFit <- jointModel(lmeFit, coxFit, timeVar = "Week")

        pe <- prederrJM(jointFit, newdata = testingData, Tstart = 10, 
                                                  Thoriz = 20)
        auc <- aucJM(jointFit, newdata = testingData, Tstart = 10, 
                                                  Thoriz = 20)
        list(pe = pe, auc = auc)
    }

    cl <- makeCluster(5)
    res <- parLapply(cl, splits, CrossValJM)
    stopCluster(cl)

函数本身被接受,但是当运行 Cluster 命令时,我遇到了这个错误,其中提到它无法识别函数中给出的对象..它们应该在函数本身中定义吗?或者我没有正确使用 parLapply 函数?

P.S.:数据如下(dfL 是长度 ~ 1000 的数据帧,dfC ~ 200):

dfL <- data.frame(ID = c(1, 1, 1, 2, 2, 3), DA = c(0.4, 1.8, 1.2, 3.2, 3.6, 2.8), Week = c(0, 4, 16, 4, 20, 8), Event = c(1, 1, 1, 0, 0, 1), TT_Event = c(16, 20, 8), Gender = c(0, 0, 0, 1, 1, 0), Age = c(24, 24, 24, 56, 56, 76))

dfC <- data.frame(ID = c(1, 2, 3, 4, 5, 6), DA = c(1.2, 3.6, 2.8, 2.4, 1.9, 3.4), Week = c(16, 20, 8, 36, 24, 32), Event = c(1, 0, 1, 1, 1, 0), TT_Event = c(16, 20, 8, 36, 24, 32), Gender = c(0, 1, 0, 0, 1, 1), Age = c(24, 56, 76, 38, 44, 50))

Thnx :)


Stack Overflow 上已经回答了非常相关的问题。 基本上,您有三种解决方案:

  • use clusterExport()将需要的变量导出到集群中(最常用的方法)
  • 将所有变量作为函数的参数传递CrossValJM()以便它们自动导出到集群(我更喜欢的解决方案,在编程上最正确的解决方案)
  • 使用 R 包 {future} ,它应该自动检测要导出的变量(懒惰的解决方案,但似乎也工作得很好)

参见示例this https://stackoverflow.com/q/47197566/6103040.

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

检查远程错误时出错 (val):5 个节点产生错误:未找到对象 的相关文章

随机推荐