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.