Torus para leg x nleg

torus.entorno.09set.r

Na linha “dados.s1$area.basal=as.numeric(area.basal)” dá erro

• Error in `$←.data.frame`(`*tmp*`, “area.basal”, value = c(8171.28, 4071.5, : replacement has 39967 rows, data has 39966
  

Fiz o match encontrei NA no dbh_09
Continua dando erro!!! só que dados.s1 tem 2 elementos a menos agora
str(area.basal)
num [1:39967(1d)] 8171 4072 139205 16837 7698 …
str(dados.s1)
'data.frame': 39965 obs. of 11 variables:
- attr(*, “dimnames”)=List of 1 ..$ : chr [1:39967] “2” “3” “4” “8” …
?(1d) vc sabe o que significa este parenteses com 1d qdo peço str(area.basal)?
Fiz diversos match e não consigo corrigir…
Imagino que no conjunto dados.s1 (no qual selecionamos apenas os fustes 1) tenha menos elementos porque alguma (duas no caso) das árvores no conjunto de dados não possuem fuste 1 mas possuem fuste 2 por exemplo, por algum erro… Mas não consigo consertar, é através do match, porque quero saber quais são os dois elementos que estão a mais no area.basal… Mas não consigo!

02out10.rdata Dentro do Rdata o objeto peic10 é a compilação de todos os dados e o peic10.uso é um subsetting no qual tirei colunas que não me interessavam, só pra vc saber…

• Função conferida e funcional… problemas com grafico ainda só…
• Modifiquei a função conta.res: utiliza densidade relativa… considero, assim, que se a densidade relativa observada (pra habitat leg) for mais extrema que 97.5% das densidades simuladas a espécie é associada
• Positiva = proporção em que simulada foi menor que observada for maior ou igual a 0.95
• Negativa = proporção em que simulada foi maior que observada for maior ou igual a 0.95
• Assim, resolvemos o problema do igual.obs…
• modificações acima baseadas no Harms 2001

Script com resultados pras 3 legs....

set.rdata

BALIZIA - Positivamente associadas = Ourapa (0.963) (mas, 13 indivs apenas) - Negats = Poutbe (0.95 - 282 indivs.) e ternbr (0.98 - 652 indivs.)
ORMOSIA - Posits. = Ternbr (0.958) - isso pode ser indício de um modelo nulo fraco? pode ser também que esteja negativamente com balipe porque está positivamente com ormosia, certo? - Negats. = nenhuma!
ANDIRA - Positd. = Blepsa (0.969), Guatau (0.994) e Ocotve (0.983) - Negats. = Peragl (0.966)

torus.julia.r torus.entorno09ago2010ii.r