如何从R中的glmnet中提取Cox模型的损失函数值?
我通过glmnet R软件包和我的 例如:如何从R中的glmnet中提取Cox模型的损失函数值?,r,glmnet,cox-regression,R,Glmnet,Cox Regression,我通过glmnet 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 i、 例如,目标函数=损失函数+惩罚函数。
library(fastcox);data(FHT);attach(FHT) #
library(glmnet)
library(survival)
fit = glmnet(x,Surv(y,status),family="cox",alpha=1)
从帮助文档中,我们知道glmnet适合以下受惩罚的模型
-loglik/nobs + λ*penalty
i、 例如,目标函数=损失函数+惩罚函数
。
我想获取-loglik/nobs
(损失函数值,
拟合模型的负偏对数似然
还是两届
从fit
对象开始的对数可能性的泰勒级数展开。
有什么想法吗?Tks
顺便说一句,我们也试过了
fit0 = glmnet(x,Surv(y,status),family="cox",alpha=1,lambda=0)
根据-loglik/nobs+λ*惩罚
,但它显示错误