search for: nulldf

...xnames[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, fitted.values = mu, 192 effects = fit$effects, R = Rmat, rank = fit$rank, qr = fit[c("...

A non-text attachment was scrubbed... Name: not available Type: text Size: 1502 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980106/e79288f9/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1502 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980106/e79288f9/attachment.pl

...ag[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 = G$sp,df=G$df,nsdf=G$nsdf,Vp=G$Vp,mgcv.conv=G$conv,gcv.ubre=G$gcv.ubre, hat=G$hat) } diff -ur c:/usr/temp/wo...