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...

Abraço,

---
Fernando de Pol Mayer
Laboratório de Estatística Ambiental - LEA
Instituto de Matemática, Estatística e Física - IMEF
Universidade Federal do Rio Grande - FURG
URL: http://fernandomayer.github.io
e-mail: fernando.mayer [@] {gmail.com, furg.br}

2014-04-15 12:24 GMT-03:00 Wagner Bonat <wbonat@gmail.com>:

Ok ! Mas vocês não acharam estranho não precisar fazer a mesma coisa para o link ?

Vc teria alguma explicação para isso ?

Obrigado !!

Em 15 de abril de 2014 14:20, Éder Comunello <comunello.eder@gmail.com> escreveu:

Wagner, bom dia!

Supondo que o vetor de valores possíveis é da classe 'expression':

### <code r>

link = "log"

variance = "mu"

quasi(link = link, variance = as.expression(variance))

# 

# Family: quasi 

# Link function: log 

variance = "mu(1-mu)"

quasi(link = link, variance = as.expression(variance))

#

#Family: quasi 

#Link function: log 

### </code>

Éder Comunello <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

--
Wagner Hugo Bonat
LEG - Laboratório de Estatística e Geoinformação
UFPR - Universidade Federal do Paraná

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.