ether I should just code up a work-around for predict.gam?
Simon
--
> Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK
> +44 1225 386603 www.maths.bath.ac.uk/~sw283
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
> > get_all_vars(~z+X,dat)
>
> z X.1 X.2 X.3
> 1 26 1 6 11
> 2 27 2 7 12
> 3 28 3 8 13
> 4 29 4 9 14
> 5 30 5 10 15
>
> > dat <- data.frame(list(X=I(matrix(1:15,5,3)),z=26:30))
> > get_all_vars(~z+X,dat)
>
>
e first fit iteration only, of course. Is it worth fixing? I
found this after someone reported a problem with mgcv::gamm (which
calls MASS:glmmPQL which calls glm), which is why I'm not just supplying
my own mustart and getting on with it...
best,
Simon
>- Simon Wood, Mat
;>
> >>> A possible replacement would be ...
> >>>
> >>> initialize = expression({
> >>>n <- rep.int(1, nobs)
> >>>if (is.null(etastart) && is.null(start) && is.null(mustart) &&
"inverse" && all(y == 0)) ||
> > (family$link == "log" && all(y <= 0 stop(
> > "cannot find valid starting values: please specify some")
> >mustart <- y
> > if (family$link==&
(y <= 0 stop(
"cannot find valid starting values: please specify some")
mustart <- y
if (family$link=="log") {
iy <- y<=0
if (sum(iy)) mustart[iy] <- min(y[!iy])*.5
} else if (family$link=="inverse"
