====== 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)