标签: r function glmnet
我需要创建一个函数,将glm泊松模型作为输入,并返回定量偏差/自由度。我用glmnet包尝试了这个:
functionGLM <- function(x,y,z) { x<- glm(y~.,data=z,family="poisson") return(deviance(x)/x$df.residual) }
然而,它产生的比率不正确,如模型(x),它也采用y变量。你知道我怎么能解决这个问题,或者有没有更好的方法来创建这个功能?