Call: nls(formula = diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == ...
Call: nls(formula = diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == ...
dados.ajuste <- expand.grid(parA = c(0.0,0.99,2,3,4,5,6)+1, sistema = factor(c("MAFA01","MAFA04"))); perfilMAFA01; perfilMAFA04
dados.ajuste$diversidade <- c(55.00000,19.572140, 10.230871, 7.224990, 6.004212, 5.402218, 5.056630,
62.00000, 31.69523, 21.41784, 17.61259, 15.76542, 14.70719, 14.03339); dados.ajuste
expMAFA01 <- nls(diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == "MAFA01",], start = list(a = 50, b = -1))
summary(expMAFA01)
lines(fitted.values(expMAFA01))
expMAFA04 <- nls(diversidade ~ a * parA^b, data = dados.ajuste[dados.ajuste$sistema == "MAFA04",], start = list(a = 50, b = -1))
summary(expMAFA04)
lines(fitted.values(expMAFA04))
ggplot(dados.ajuste,aes(parA,diversidade))+
stat_function(fun=function(parA){coef(expMAFA01)[[1]] * parA^coef(expMAFA01)[[2]]},
aes(colour="MAFA01", linetype="MAFA01")) +
stat_function(fun=function(parA){coef(expMAFA04)[[1]] * parA^coef(expMAFA04)[[2]]},
aes(colour="MAFA04", linetype="MAFA04"))
AIC(expMAFA01,expMAFA04)
vuongtest(expMAFA01,expMAFA04)
--
Engenheiro Florestal - UFRA
Mestre em Ciência Florestal - UFVJM
Doutorando em Ciência Florestal - UFVJM