Wagner,
se olhar dentro da função quasi, no início tem isso para o link
linktemp <- substitute(link)
if (!is.character(linktemp))
linktemp <- deparse(linktemp)
if (linktemp %in% c("logit", "probit", "cloglog", "identity",
"inverse", "log", "1/mu^2", "sqrt"))
stats <- make.link(linktemp)
else if (is.character(link)) {
stats <- make.link(link)
linktemp <- link
Ou seja, passando um caracter como
link <- "log"
vai fazer com que a função use
make.link(link)
e funciona como esperado. Já para a variance, a parte relevante eh essa
vtemp <- substitute(variance)
if (!is.character(vtemp))
vtemp <- deparse(vtemp)
variance_nm <- vtemp
switch(vtemp, constant = {
...
Portanto, se fizermos
variance <- "mu"
e seguir a lógica
(vtemp <- substitute(variance))
variance
class(vtemp)
[1] "name"
vtemp ja não eh mais um character, portanto ele vai executar
vtemp <- deparse(vtemp)
vtemp
[1] "variance"
variance_nm <- vtemp
variance_nm
[1] "variance"
Como dentro das opções do switch(...) ele não vai achar nenhuma opção, ele vai para a alternativa no final que é
switch(vtemp, constant = {
...
}, variance_nm <- NA)
e variance_nm sendo NA vai cair nisso aqui
if (
is.na(variance_nm)) {
if (is.character(variance))
stop(gettextf("'variance' \"%s\" is invalid: possible values are \"mu(1-mu)\", \"mu\", \"mu^2\", \"mu^3\" and \"constant\"",
variance_nm), domain = NA)
que é o erro. Agora, se isso faz sentido ou se é uma "falha" de programação na função já não sei dizer (acredito que isso seja intencional, mas o motivo não da pra saber).
De qualquer maneira vc pode modificar a função quasi() pra fazer do jeito que vc esperaria, começando por tirar o substitute() ali no começo...