R保护堆栈溢出偶数设置（表达式=5e5）：无法拆分数据帧

我正在运行R中的基因表达分析，并处理大量数据。我想对这个数据应用Cox模型（249R乘以23000C），但当我尝试这样做时，我收到一条错误消息（Protect（）：Protection Stack Overflow）。我意识到很多其他人都问过这个问题，我也尝试过使用前面这些问题的答案，比如选项（表达式=5e5）和glmnet，但后者给了我一个错误（错误在

[.data.frame

（数据，c（lhsVars，rhsVars））： 我不知道在进行Cox回归时如何使用它，而且前一个选项不起作用。另外，以我执行项目的方式，拆分数据帧对我来说不是一个好选项

我已经附上了我的代码，非常感谢任何帮助！如果我的代码凌乱，我很抱歉——我是一名学生，还有很多东西要学

totranspose <- read.delim("/Users/menon/OneDrive/Desktop/csrsef files/FilesToTranspose.txt", stringsAsFactors = F, header = T)
rownames(totranspose) <- totranspose$X
totranspose$X <- NULL
finaldatacsrsef <- t(totranspose)
colnames(finaldatacsrsef) <- as.character(unlist(finaldatacsrsef[1,]))
finaldatacsrsef <- finaldatacsrsef[-1,]

library(survival)
initData <- data.frame(finaldatacsrsef)
eventfreetime <- initData[ , c("EFST")] 
vitalstatus <- initData[ , c("Status")]
temporaryObj <- Surv(as.numeric(eventfreetime), vitalstatus == 2)
finalDataa <- data.frame(SurvObj = temporaryObj)
finalDataa <- cbind(finalDataa, initData[, 3:ncol(initData)])
options(expressions = 5e5)
finalcox <- coxph(SurvObj ~ ., data =  finalDataa)
summary(finalcox)$coefficients

---

嗨，安基塔，在调试了一些部分后，我尝试了你的代码：

tobind

未定义，

coxph（SurvObj~）

，并使用249行23000列的伪数据，但我没有得到任何错误。因此，您可能需要在一些数据上显示

dput

的输出，以便我们了解您正在建模的变量类型。只需前20行和10列就可以了。

dput（finalDataa[1:20，1:10）

@Edward--再次感谢您的帮助！我无意中包含了tobind--它在我的实际代码中被注释掉了，因为我不需要它。下面是dput的输出（整个东西对于stack来说太长了，但这里有一点）：structure（list）（SurvObj=structure）（c（166、193、189、233、45、176、141、63、75、194、182、94、191、178、27、67、228、200、136、78、3、0、0、3、3、6、3、3、3、3、0、0、0、2、3、3、0、0），.Dim=c（20L、2L），.Dimnames=list（NULL、c（“时间”、“状态”）），AADACL3=structure（c）（TARGET.30.PAAPFA=219L，TARGET.30.PACLJN=103L，TARGET.30.PACRYY=83L…你能把它放在问题中吗？@Edward是的，我已经把它放在了问题中，谢谢。我真的很抱歉，但是你能为

initData

而不是

finalDataa

再做一次吗？Surv对象的结构看起来不正确。但更重要的是，为什么您的一些状态值是3和6吗？它们应该是0或1（或1或2）。

> dput(initData[1:10, 1:10])
structure(list(EFST = structure(c(TARGET.30.PAAPFA = 166L, TARGET.30.PACLJN = 193L, 
TARGET.30.PACRYY = 189L, TARGET.30.PACRZM = 233L, TARGET.30.PACSNL = 45L, 
TARGET.30.PACSSR = 176L, TARGET.30.PACUGP = 141L, TARGET.30.PACVNB = 63L, 
TARGET.30.PACYGY = 75L, TARGET.30.PACZPX = 194L), .Label = c(" 2.00000", 
" 9.00000", "1002.00000", "1013.00000", "1022.00000", "1030.00000", 
...
"165.00000", "1659.00000", "168.00000", "1696.00000", "1707.00000", 
"1710.00000", "172.00000", "1749.00000", "1795.00000", "1836.00000", 
"1843.00000", "1880.00000", "190.00000", "1933.00000", "1964.00000", 
"2007.00000", "2025.00000", "203.00000", "205.00000", "206.00000", 
"2147.00000", "2155.00000", "217.00000", "2183.00000", "219.00000", 
...
"887.00000", "922.00000", "932.00000", "956.00000", "964.00000", 
"968.00000"), class = "factor"), Status = structure(c(TARGET.30.PAAPFA = 4L, 
TARGET.30.PACLJN = 1L, TARGET.30.PACRYY = 1L, TARGET.30.PACRZM = 4L, 
TARGET.30.PACSNL = 4L, TARGET.30.PACSSR = 7L, TARGET.30.PACUGP = 4L, 
TARGET.30.PACVNB = 4L, TARGET.30.PACYGY = 4L, TARGET.30.PACZPX = 1L
), .Label = c("   1.00000", "   2.00000", "  1.00000", "  2.00000", 
"  2.000000", " 1.00000", " 2.00000"), class = "factor"), AADACL3 = structure(c(TARGET.30.PAAPFA = 219L, 
TARGET.30.PACLJN = 103L, TARGET.30.PACRYY = 83L, TARGET.30.PACRZM = 234L, 
TARGET.30.PACSNL = 150L, TARGET.30.PACSSR = 247L, TARGET.30.PACUGP = 140L, 
TARGET.30.PACVNB = 118L, TARGET.30.PACYGY = 227L, TARGET.30.PACZPX = 99L
), .Label = c("   4.70471", "   4.71884", "   4.75278", "   4.77252", 
"   4.79924", "   4.82726", "   4.83177", "   4.85297", "   4.85348", 
...
"  5.74958", "  5.76020", "  5.77283", "  5.77490", "  5.79089", 
"  5.91416", " 5.02643", " 5.16070", " 5.20377", " 5.33849", 
" 5.43057", " 5.70252", " 5.88764"), class = "factor"), AADACL4 = structure(c(TARGET.30.PAAPFA = 201L, 
TARGET.30.PACLJN = 104L, TARGET.30.PACRYY = 75L, TARGET.30.PACRZM = 213L, 
TARGET.30.PACSNL = 112L, TARGET.30.PACSSR = 247L, TARGET.30.PACUGP = 170L, 
TARGET.30.PACVNB = 199L, TARGET.30.PACYGY = 140L, TARGET.30.PACZPX = 30L
), .Label = c("   4.37049", "   4.41956", "   4.56399", "   4.60298", 
"   4.66085", "   4.66949", "   4.67545", "   4.69108", "   4.69205", 
...
"  5.77346", " 4.84060", " 4.93969", " 5.01135", " 5.03256", 
" 5.22472", " 5.29179", " 5.70923"), class = "factor"), ACADM = structure(c(TARGET.30.PAAPFA = 182L, 
TARGET.30.PACLJN = 37L, TARGET.30.PACRYY = 86L, TARGET.30.PACRZM = 218L, 
TARGET.30.PACSNL = 238L, TARGET.30.PACSSR = 241L, TARGET.30.PACUGP = 215L, 
TARGET.30.PACVNB = 201L, TARGET.30.PACYGY = 131L, TARGET.30.PACZPX = 70L
), .Label = c("   6.25757", "   6.33278", "   6.39292", "   6.61567", 
"   6.70676", "   6.77349", "   6.82724", "   6.84948", "   6.88238", 
"   7.03772", "   7.08124", "   7.11364", "   7.16697", "   7.21509", 
"   7.21749", "   7.22232", "   7.22508", "   7.24340", "   7.26963", 
"   7.27326", "   7.28288", "   7.29225", "   7.30015", "   7.31206", 
….
" 7.94490", " 8.10245", " 8.15717"), class = "factor"), ACOT11 = structure(c(TARGET.30.PAAPFA = 211L, 
TARGET.30.PACLJN = 25L, TARGET.30.PACRYY = 35L, TARGET.30.PACRZM = 176L, 
TARGET.30.PACSNL = 196L, TARGET.30.PACSSR = 247L, TARGET.30.PACUGP = 156L, 
TARGET.30.PACVNB = 186L, TARGET.30.PACYGY = 199L, TARGET.30.PACZPX = 63L
), .Label = c("   5.93635", "   5.94244", "   5.95103", "   5.98901", 
"   6.03946", "   6.11429", "   6.11859", "   6.13903", "   6.15869", 
 ...
" 6.54330", " 6.58306", " 6.78981"), class = "factor"), ACOT7 = structure(c(TARGET.30.PAAPFA = 177L, 
TARGET.30.PACLJN = 88L, TARGET.30.PACRYY = 15L, TARGET.30.PACRZM = 227L, 
TARGET.30.PACSNL = 176L, TARGET.30.PACSSR = 243L, TARGET.30.PACUGP = 112L, 
TARGET.30.PACVNB = 234L, TARGET.30.PACYGY = 190L, TARGET.30.PACZPX = 2L
), .Label = c("   7.27372", "   7.44455", "   7.57713", "   7.63049", 
"   7.67607", "   7.67940", "   7.67984", "   7.69863", "   7.74928", 
"   7.79042", "   7.79045", "   7.79654", "   7.81007", "   7.83489", 
"   7.83490", "   7.83660", "   7.85289", "   7.86943", "   7.87172", 
"   7.87769", "   7.88242", "   7.89099", "   7.91394", "   7.91521", 
"   7.91740", "   7.93080", "   7.93191", "   7.94568", "   7.96315", 
"   7.97490", "   7.98972", "   7.99591", "   8.00916", "   8.01633", 
"   8.02621", "   8.03156", "   8.03209", "   8.04701", "   8.05113", 
"   8.05589", "   8.05917", "   8.06400", "   8.06733", "   8.08233", 
"   8.10082", "   8.10544", "   8.11094", "   8.12160", "   8.13253", 
"   8.13579", "   8.14512", "   8.14566", "   8.15169", "   8.15522", 
"   8.15733", "   8.16079", "   8.16116", "   8.16126", "   8.16969", 
"   8.17116", "   8.18370", "   8.20090", "   8.20675", "   8.20713", 
"   8.21006", "   8.21736", "   8.21749", "   8.21878", "   8.24505", 
"   8.24802", "   8.25781", "   8.27362", "   8.27436", "   8.27646", 
"   8.28431", "   8.28967", "   8.29238", "   8.29331", "   8.29626", 
"   8.31274", "   8.31330", "   8.31784", "   8.31881", "   8.32197", 
"   8.34618", "   8.35837", "   8.37840", "   8.38011", "   8.38226", 
"   8.38984", "   8.40466", "   8.44161", "   8.45682", "   8.46116", 
"   8.49594", "   8.50602", "   8.52247", "   8.53077", "   8.55986", 
"   8.57609", "   8.70593", "   8.78929", "   8.79241", "   8.83021", 
"   8.88021", "   8.89001", "   8.90448", "   9.11322", "   9.43475", 
"  7.52861", "  7.66466", "  7.68660", "  7.70969", "  7.71580", 
"  7.72014", "  7.73900", "  7.73927", "  7.74026", "  7.75380", 
"  7.75922", "  7.76969", "  7.77150", "  7.78116", "  7.78417", 
"  7.79180", "  7.79952", "  7.80261", "  7.80837", "  7.81161", 
"  7.81396", "  7.82174", "  7.82495", "  7.82584", "  7.82670", 
"  7.83582", "  7.83650", "  7.85396", "  7.85827", "  7.87575", 
"  7.88243", "  7.88348", "  7.88520", "  7.88999", "  7.90379", 
"  7.91414", "  7.91715", "  7.93460", "  7.93812", "  7.96672", 
"  7.97637", "  7.98229", "  7.98352", "  8.00525", "  8.01628", 
"  8.01696", "  8.01744", "  8.01888", "  8.01977", "  8.02221", 
"  8.02997", "  8.03639", "  8.04479", "  8.05963", "  8.07101", 
"  8.07438", "  8.07717", "  8.08261", "  8.08476", "  8.09201", 
"  8.09559", "  8.10024", "  8.10100", "  8.10383", "  8.10417", 
"  8.11497", "  8.11502", "  8.11807", "  8.12097", "  8.12552", 
"  8.13045", "  8.13261", "  8.13270", "  8.13377", "  8.13397", 
"  8.13897", "  8.14309", "  8.14993", "  8.15376", "  8.15771", 
"  8.17004", "  8.17409", "  8.17599", "  8.19405", "  8.20733", 
"  8.21219", "  8.22061", "  8.23543", "  8.23934", "  8.24026", 
"  8.24202", "  8.24563", "  8.26212", "  8.26387", "  8.26394", 
"  8.27844", "  8.27997", "  8.29101", "  8.30884", "  8.31898", 
"  8.32005", "  8.32071", "  8.32172", "  8.32829", "  8.32889", 
"  8.33430", "  8.33625", "  8.34192", "  8.34423", "  8.36700", 
"  8.37298", "  8.37877", "  8.39039", "  8.41815", "  8.42202", 
"  8.42291", "  8.43199", "  8.43645", "  8.43953", "  8.48244", 
"  8.49337", "  8.49969", "  8.52678", "  8.54624", "  8.55327", 
"  8.57020", "  8.57290", "  8.66956", "  8.67772", "  8.78580", 
"  9.096890", " 7.80377", " 7.84975", " 8.01090", " 8.18525", 
" 8.22819", " 8.25531", " 8.47783"), class = "factor"), ACTB = structure(c(TARGET.30.PAAPFA = 139L, 
TARGET.30.PACLJN = 93L, TARGET.30.PACRYY = 90L, TARGET.30.PACRZM = 218L, 
TARGET.30.PACSNL = 164L, TARGET.30.PACSSR = 237L, TARGET.30.PACUGP = 211L, 
TARGET.30.PACVNB = 230L, TARGET.30.PACYGY = 202L, TARGET.30.PACZPX = 37L
), .Label = c("  10.29300", "  10.33920", "  10.37040", "  10.51230", 
"  10.52510", "  10.55200", "  10.57640", "  10.60880", "  10.63620", 

"10.84560", "10.85970", "11.00500", "11.11320"), class = "factor"), 
    ACTL8 = structure(c(TARGET.30.PAAPFA = 212L, TARGET.30.PACLJN = 73L, 
    TARGET.30.PACRYY = 24L, TARGET.30.PACRZM = 196L, TARGET.30.PACSNL = 151L, 
    TARGET.30.PACSSR = 247L, TARGET.30.PACUGP = 148L, TARGET.30.PACVNB = 124L, 
    TARGET.30.PACYGY = 236L, TARGET.30.PACZPX = 61L), .Label = c("   5.24716", 
    "   5.37788", "   5.46886", "   5.52772", "   5.56753", "   5.57768", 
    ...
    " 6.69476"), class = "factor")), row.names = c("TARGET.30.PAAPFA", 
"TARGET.30.PACLJN", "TARGET.30.PACRYY", "TARGET.30.PACRZM", "TARGET.30.PACSNL", 
"TARGET.30.PACSSR", "TARGET.30.PACUGP", "TARGET.30.PACVNB", "TARGET.30.PACYGY", 
"TARGET.30.PACZPX"), class = "data.frame")