[R-br] Ajuda quasi likelihood

[Rodrigo Sant'Ana]

Wagner,

  não sei se é isso que você quer...

mas a solução é só atribuir um valor para mu. Pois a variance tem que ser
numerica.
E com isso tua função pode ser genérica. Olha o exemplo que fiz, talvez
seja isso, não?

link <- "log"
mu <- 2
quasi(link = link, variance = mu)

Assim, teu mu pode assumir quaisquer valor conforme a necessidade da
aplicação da função.
Espero ter te ajudado...

Abs,
RS

___________________________________________________
*Rodrigo Sant'Ana*
--
Mestre em Ciência e Tecnologia Ambiental - MCTA/UNIVALI
Graduado em Oceanografia - CTTMar/UNIVALI
Universidade do Vale do Itajaí - UNIVALI
--



Em 15 de abril de 2014 05:08, Wagner Bonat <wbonat em gmail.com> escreveu:

> Caros,
>
> Alguem tem ideia do que acontece nesta situacao
>
> quasi(link = "log", variance = "mu")
>
> Family: quasi
> Link function: log
>
> Resultado perfeito, agora preciso deixar generico para usar a funcao quasi
> dentro de outra funcao para obter valores iniciais, entao pensei algo assim
>
> link = "log"
> variance = "mu"
> quasi(link = link, variance = variance)
> Error in quasi(link = link, variance = variance) (from quasi.R em 3553CaW#69)
> :
>   'variance' "NA" is invalid: possible values are "mu(1-mu)", "mu",
> "mu^2", "mu^3" and "constant"
>
> Alguma ideia ??
>
> Obrigado
>
>
>
> --
> Wagner Hugo Bonat
> LEG - Laboratório de Estatística e Geoinformação
> UFPR - Universidade Federal do Paraná
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140415/a50e66a9/attachment.html>