如何从R中的glmnet中提取Cox模型的损失函数的值?

时间:2017-07-29 07:48:42

标签: r glmnet cox-regression

我使用Cox模型通过glmnet R包和我的 小R的例子是:

    library(fastcox);data(FHT);attach(FHT)  #   
    library(glmnet)
    library(survival)
    fit = glmnet(x,Surv(y,status),family="cox",alpha=1)

从帮助文档中,我们知道glmnet适合像

这样的惩罚模型
    -loglik/nobs + λ*penalty

即,objective function = loss function + penalty function。 我想获取-loglik/nobs(损失函数值, 拟合模型的负部分对数似然 或两个学期 来自fit对象的泰勒级数对数似然扩展。 任何的想法? TKS

顺便说一下,我们也试过

fit0 = glmnet(x,Surv(y,status),family="cox",alpha=1,lambda=0)

根据-loglik/nobs + λ*penalty,但它显示错误。

0 个答案:

没有答案