我使用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
,但它显示错误。