Leg + Não-leg juntas

Comparei 5 modelos: 1) O mais simples : copa.vero1←function(m=2.542,s=1.96)

{
   -sum(stats::dnorm(alo$Raio.copa.m, m, s, log=TRUE))
}

copa1←mle(copa.vero1) 2) Com a função original: copa.vero2←function(a=0,b=1,s=1.96)

{
m = a+b*((0.1*alo$Dap.mm)^0.5/2)
-sum(stats::dnorm(x=alo$Raio.copa.m,mean=m,sd=s,log=T))
}

copa2←mle(copa.vero2) 3) Com outra função, sem elevar a 0.5: copa.vero3←function(a=0,b=1,s=1.96)

{
m = a+b*((0.1*alo$Dap.mm)/2)
-sum(stats::dnorm(x=alo$Raio.copa.m,mean=m,sd=s,log=T))
}

copa3←mle(copa.vero3) 4) Elevando ao quadrado sem dividir: copa.vero4←function(a=0,b=1,s=1.96)

{
m = a+b*(0.1*alo$Dap.mm)^2
-sum(stats::dnorm(x=alo$Raio.copa.m,mean=m,sd=s,log=T))
}

copa4←mle(copa.vero4)

5) Com média e desvio variando pela mesma função: copa.vero5←function(a=1.43,b=0.002,s1=1,s2=2)

{
m = a+b*((0.1*alo$Dap.mm)/2)
s = s1+s2*((0.1*alo$Dap.mm)/2)
-sum(stats::dnorm(x=alo$Raio.copa.m,mean=m,sd=s,log=T))
}

copa5←mle(copa.vero5) PS: Quando estava passando os dados para o wiki agora, percebi que no modelo 5 tinha elevado a média ao quadrado ao invés de dividir por 2, por isso, ontem comparando AICs tinha ficado com o modelo 3, pois o modelo 5 tinha ficado com AIC semelhante ao 3, mas o 3 tem menos parâmetros # AIC #copa1 111.70791 #copa2 71.15744 #copa3 68.34604 #copa4 77.40857 #copa5 63.01979

* Como a curiosidade é grande rodei alguns nulos alterando a função r.copa para o modelo 3 e o resultado foi:

• Para Abarbr, 259 médias ficaram inferiores ao valor observado. Com a função anterior esse valor era 0.
• Balipe continuo com 0, como possui bastante dbh grande acredito que a função possa estar errando aí…
• Ormoar teve 8 valores inferiores, com a função anterior era 1 (sem muita diferença então)