如何使用RMARK中的因子进行已知的命运分析?

时间:2015-06-11 19:35:03

标签: r

我想在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

如何使用因子运行模型?预先感谢您的帮助。

1 个答案:

答案 0 :(得分:0)

我认为您需要在模型中将sex定义为组变量。

示例:survival.processed <- process.data(data, model="Known",groups=c("sex", "Age"))