我想在RMARK中运行模型,以便使用因子和协变量进行已知的命运分析。
我使用的代码是:
do_analysis=function()
{
S.dot=list(formula=~1)
S.time=list(formula=~time)
S.Sex=list(formula=~Sex)
S.S_ESH=list(formula=~S_ESHerb)
model.list=create.model.list("Known")
Surv.results=mark.wrapper(model.list,data=dp,ddl=ddl,
invisible=FALSE,threads=1)
return(Surv.results)
}
Surv.results=do_analysis()
Surv.results$model.table=model.table(Surv.results,use.lnl=TRUE,model.name=FALSE)
Surv.results
我收到了错误:
S.Sex
Error in make.mark.model(data.proc, title = title, parameters = model.parameters, :
The following individual covariates are not allowed because they are factor variables: Sex
Error in mark(model.parameters = model.parameters, initial = initial, :
Misspecification of model or internal error in code
如何使用因子运行模型?预先感谢您的帮助。
答案 0 :(得分:0)
我认为您需要在模型中将sex
定义为组变量。
示例:survival.processed <- process.data(data, model="Known",groups=c("sex", "Age"))