[R-br] Comparação de modelos não lineares via teste de Vuong [nonnest2]

[Walmes Zeviani]

Recomendo proceder da seguinte forma para testar o efeito de sistema sobre
a curva ajustada, ou seja, se as curvas são iguais independente dos níveis
do fator sistema. O teste usado abaixo é o de razão de verossimilhanças. Eu
desconheço a teoria que sustenta o teste de Vuong que você mencinou, mas
não me parece correto porque você está usando dados diferentes e ajustando
o mesmo modelo. Eu creio que ele seja recomendado para avaliar a igualdade
de dois modelos aplicados ao mesmo conjunto dados. Segue CMR.

library(nlme)

n0 <- gnls(diversidade ~ a * parA^b,
           data = dados.ajuste,
           params = a + b ~ 1,
           start = list(50, -1))
summary(n0)

n1 <- gnls(diversidade ~ a * parA^b,
           data = dados.ajuste,
           params = a + b ~ sistema,
           start = list(50, 0, -1, 0))
summary(n1)

# LRT.
anova(n0, n1)

anova(n1, test = "marginal")

​À
​ disposição.

Walmes​.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20171107/e048bb89/attachment.html>