<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">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.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:monospace,monospace">library(nlme)<br><br>n0 <- gnls(diversidade ~ a * parA^b,<br>           data = dados.ajuste,<br>           params = a + b ~ 1,<br>           start = list(50, -1))<br>summary(n0)<br><br>n1 <- gnls(diversidade ~ a * parA^b,<br>           data = dados.ajuste,<br>           params = a + b ~ sistema,<br>           start = list(50, 0, -1, 0))<br>summary(n1)<br><br># LRT.<br>anova(n0, n1)<br><br>anova(n1, test = "marginal")<br></span><br></div>​À<div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default">​ disposição.</div><div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default"><br></div><div style="font-family:trebuchet ms,sans-serif;display:inline" class="gmail_default">Walmes​.</div></div>