Displaying 2 results from an estimated 2 matches for "quasi2".
Did you mean:
quasi
2005 Jun 16
1
mu^2(1-mu)^2 variance function for GLM
...on(mu) all(mu > 0) && all(mu < 1)
dev.resids <- function(y, mu, wt) 2 * wt * ((2 * y - 1) *
(log(ifelse(y == 0, 1, y/mu)) - log(ifelse(y == 1, 1,
(1 - y)/(1 - mu)))) - 2 + y/mu + (1 - y)/(1 - mu))
I've modified the `quasi' function accordingly (into `quasi2' given
below) and my results are very much in line with the ones cited by
McCullagh and Nelder on p.330-331:
> data(leafblotch, package = "faraway")
> summary(fit <- glm(blotch ~ site + variety,
+ family = quasi2(link = "logit", variance = "mu^2(1-mu...
2019 Apr 26
1
Error in glm(..., family=quasi(..., variance=list(...)))
...arse() breaks lines at width.cutoff = 60L characters. In stepping through quasi() during debug, setting vtemp <- (vtemp, collapse=" ") on line 576 avoids the error.
A workaround from https://tolstoy.newcastle.edu.au/R/help/05/06/6795.html appears to be to define one's own complete quasi2() function with the desired variance function pre-stored.
Is this known/expected, or should I file a bug?
Many thanks and best regards
Daniel
> sessionInfo()
R version 3.6.0 (2019-04-26)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 14393)
Matrix products: d...