全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
569 0
2023-12-13
##归因危险度
library(survival)
set.seed(123)
total <- nutrDII0[,c("DII","age","sex","Townsend_index","race","BMI","MET","pack_smoking",
                     "familiy_cancer","assess_centre","education","cancer_all","survival_all")]
total$DIIpar <- ifelse(total$DII<"-1",0,ifelse(total$DII>="-1"&total$DII<=1,1,2))
#创建生存对象
#total$survival_all <- as.numeric(total$survival_all)
#total$cancer_all <- as.factor(total$cancer_all)
surv_obj <- with(total,Surv(survival_all,cancer_all))
#拟合生存模型
fit <- survfit(surv_obj~as.factor(DIIpar)+as.numeric(age) +
                 as.factor(sex) + as.numeric(Townsend_index) +
                 as.factor(race)  + as.numeric(BMI) + as.numeric(MET) +
                 as.factor(familiy_cancer) +as.factor(education) +
                 as.factor(assess_centre)+as.factor(pack_smoking),data=total)
par_estimation <- par(fit)
print(par_estimation)
运行上述代码后出现“[size=13.3333px]Error: cannot allocate vector of size 5569.7 Gb”报错情况,怎么解决啊啊啊啊啊


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群