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.