search for: nulldev

...age source code? .Fortran("lognet", parm=alpha,nobs,nvars,nc,as.double(x),y,offset,jd,vp,ne,nx,nlam,flmin,ulam,thresh,isd,maxit,kopt, lmu=integer(1), a0=double(nlam*nc), ca=double(nx*nlam*nc), ia=integer(nx), nin=integer(nlam), nulldev=double(1), dev=double(nlam), alm=double(nlam), nlp=integer(1), jerr=integer(1),PACKAGE="glmnet" ) Thanks for your help. Regards, Axel. [[alternative HTML version deleted]]

...mapped' Traceback: 1: .Fortran("lognet", parm = alpha, nobs, nvars, nc, as.double(x), y, offset, jd, vp, ne, nx, nlam, flmin, ulam, thresh, isd, maxit, kopt, lmu = integer(1), a0 = double(nlam * nc), ca = double(nx * nlam * nc), ia = integer(nx), nin = integer(nlam), nulldev = double(1), dev = double(nlam), alm = double(nlam), nlp = integer(1), jerr = integer(1), PACKAGE = "glmnet") 2: lognet(x, is.sparse, ix, jx, y, weights, offset, type, alpha, nobs, nvars, jd, vp, ne, nx, nlam, flmin, ulam, thresh, isd, vnames, maxit, HessianExact, family)...

...9 names(weights) <- ynames 180 names(y) <- ynames 181 names(fit$effects) <- c(xxnames[seq(fit$rank)], rep("", sum(good) - 182 fit$rank)) 183 wtdmu <- if (intercept) 184 sum(weights * y)/sum(weights) 185 else linkinv(offset) 186 nulldev <- sum(dev.resids(y, wtdmu, weights)) 187 n.ok <- nobs - sum(weights == 0) 188 nulldf <- n.ok - as.integer(intercept) 189 resdf <- n.ok - fit$rank 190 aic.model <- aic(y, n, mu, weights, dev) + 2 * fit$rank 191 list(coefficients = coef, residuals = residuals,...

...+1:direct.mesh],g =ddiag[1:m],h=ddiag[(m+1):(2*m)], e=ddiag[(2*m+1):(3*m)],iter=idiag[1],init.ok=as.logical(idiag[2]),step.fail= as.logical(idiag[3])) # unpack results back to correct place in output (which includes fixed d.f. and free d.f. terms) @@ -1247,7 +1250,7 @@ null.deviance = nulldev, iter = iter, weights = wt, prior.weights = weights, #df.residual = resdf, df.null = nulldf, y = y, converged = conv,sig2=G$sig2,edf=G$edf, - boundary = boundary,sp = G$sp,df=G$df,nsdf=G$nsdf,Vp=G$Vp,mgcv.conv=G$conv,gcv.ubre=G$gcv.ubre) + boundary = boundary,sp =...

Dear R-team As I didn't get any answer to my bug-report last week I have taken the effort and extracted a minimal data set from my data (see below) where the following bug occurs: > glm(SKR.ein.aus ~ ., family = binomial, data = bugdata, na.action = na.omit) Error in names<-.default(*tmp*, value = ynames) : names attribute must be the same length as the vector In addition: Warning