Modifiquei a função e rodei análise de torus para as espécies que apresentam pelo menos 50 indivíduos na parcela. A tabela de resultados é composta pela contagem dos indivíduos para cada hábitat analisado (Neo, Trans e Esp Aren) e as probabilidades associadas ao número de vezes que o valor observado é maior, menor ou igual (Gr, Ls, Eq) à simulação torus.

Tabela de resultados brutos:

tab_torus.txt

** Densidade observada maior para o habitat Neossolo **

torus.peic[(torus.peic[,2]/512)>0.95,2]/512
   Astrac    Calobr    Euteed    Gomisc    Guatau    Heissi    Jacapu    Marlra    Myrcra    Ormoar 
0.9921875 0.9746094 0.9980469 0.9511719 0.9570312 0.9960938 0.9531250 0.9980469 0.9824219 0.9531250 

## NEO associação negativa      ("Ls.Hab.neo")

torus.peic[(torus.peic[,3]/512)>0.95,3]/512
 
## Blepsa    Calyco    Cluscr    Maytro    Myrcgl    Pimecf    Poutbe    Psidca    Ternbr 
##0.9550781 0.9980469 0.9843750 0.9648438 0.9804688 0.9882812 0.9511719 0.9960938 0.9921875 


##  Trans Positivo ("Gr.Hab.trans")


rownames(torus.peic)[(torus.peic[,6]/512)>0.95]
#[1] "Psidca"

torus.peic[(torus.peic[,6]/512)>0.95,6]/512
#[1] 0.9628906



## Trans NEgativo  ("Ls.Hab.trans")

torus.peic[(torus.peic[,7]/512)>0.95,7]/512

##   Euteed    Guatau 
##0.9960938 0.9882812 

## Esp.aren Positivo ("Gr.Hab.esp.aren")

torus.peic[(torus.peic[,10]/512)>0.95,10]/512

## Blepsa    Micocu    Myrcgl    Peragl    Psidca    Ternbr 
##0.9531250 0.9960938 0.9863281 0.9960938 0.9667969 0.9804688 



## Esp.aren Negativo ("Ls.Hab.esp.aren")

torus.peic[(torus.peic[,11]/512)>0.95,11]/512

## Alchtr    Anibvi    Astrac    Euteed    Garcga    Gomisc    Heissi   Marlra    Matagu    Myrcgr    Myrcra 
##0.9921875 0.9902344 0.9843750 0.9902344 0.9941406 0.9824219 0.9785156 0.9921875 0.9902344 0.9980469 0.9980469


## Positivo Esp Hist ("Gr.Hab.esp.hist")

torus.peic[(torus.peic[,14]/512)>0.95,14]/512

##  Alchtr    Anibvi    Byrsli    Garcga    Matagu    Ocotpr 
## 0.9843750 0.9863281 0.9980469 0.9882812 0.9648438 0.9960938 

## Negativo Esp Hist ("Ls.Hab.esp.hist")

torus.peic[(torus.peic[,15]/512)>0.95,15]/512

##Andian    Balipe    Ilexth    Myrcsp    Schean 
##0.9707031 0.9804688 0.9667969 0.9765625 0.9667969 

1. torus07set09.rdata
2. script.torus.ale07set09.r