library(vegan)
Ityphilus_sp6 = c(0, 0, 0, 0, 0, 1, 0)
Schendylops_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Schendylops_sp3 = c(0, 0, 0, 0, 0, 0, 1)
Lamyctes_sp = c(0, 0, 1, 0, 0, 0, 0)
Scolopendra_viridicornis = c(0, 0, 0, 1, 0, 0, 0)
Sphendononema_guildingii = c(0,	0, 1, 1, 1, 0, 1)
Polydesmida_sp1 = c(0, 0, 0, 0, 0, 1, 1)
Polydesmida_sp3 = c(0, 0, 0, 1, 0, 0, 0)
Pseudonannolene_imbirenses = c(0, 0, 1, 0, 1, 0, 0)
Pseudonannolenidae_sp = c(1, 0, 1, 0, 0, 0, 0)
Pseudonannolenidae_sp2 = c(0, 0, 0, 1, 1, 0, 0)
Pseudonannolenidae_sp3 = c(0, 0, 0, 1, 0, 1, 1)
Pseudonannolenidae_sp4 = c(0, 0, 0, 1, 1, 0, 0)
Pseudonannolenidae_sp5 = c(0, 1, 0, 0, 1, 0, 1)
Pseudonannolenidae_sp6 = c(0, 0, 0, 1, 0, 0, 1)
Pseudonannolenidae_sp7 = c(0, 0, 0, 0, 0, 0, 1)
Spirostreptidae_sp12 = c(0, 0, 1, 0, 0, 1, 1)
Siphonophorida_sp1 = c(0, 0, 0, 1, 0, 0, 0)
Blattaria_sp1 = c(1, 0, 1, 0, 1, 1, 1)
Blattaria_sp2 = c(0, 1, 1, 1, 1, 0, 0)
Blattaria_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Blattaria_sp5 = c(0, 0, 1, 0, 0, 0, 0)
Blableridae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Blattellidae_sp1 = c(0, 1, 1, 1, 0, 0, 1)
Blattellidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Blattellidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Polyphagidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_indet = c(0,	0, 1, 0, 0, 0, 1)
Coleoptera_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Coleoptera_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp3 = c(0, 0, 0, 0, 1, 0, 0)
Coleoptera_sp4 = c(0, 0, 0, 0, 1, 0, 1)
Coleoptera_sp5 = c(0, 1, 1, 0, 0, 0, 0)
Coleoptera_sp6 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp7 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp8 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp9 = c(0, 0, 0, 1, 0, 0, 0)
Coleoptera_sp10 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp11 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp12 = c(0, 0, 1, 0, 0, 0, 0)
Coleoptera_sp13 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp14 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp15 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp16 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp17 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp18 = c(0, 0, 0, 0, 0, 0, 1)
Coleoptera_sp20 = c(0, 0, 0, 0, 0, 1, 0)
Cantharidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Cantharidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Cantharidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Carabidae_indet1 = c(1, 0, 0, 0, 0, 0, 0) 
Pentagonica_sp1 = c(0, 0, 1, 0, 0, 1, 1)
Platynus_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Pterostichini_sp2 = c(0, 0, 0, 0, 0, 0, 1)
Agonina_sp1 = c(0, 0, 1, 0, 0, 1, 1)
Agonina_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Loxandrus_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Paratachys_sp1 = c(0, 0, 1, 1, 0, 0, 1)
Paratachys_sp2 = c(0, 0, 1, 0, 0, 0, 1)
Pericompsus_sp1 = c(0, 0, 1, 0, 0, 1, 1)
Polyderis_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Tachys_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Carabidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Carabidae_sp2 = c(1, 0, 0, 0, 0, 0, 0)
Carabidae_sp3 = c(1, 0, 0, 0, 0, 0, 0)
Zuphiini_sp2 = c(1, 0, 0, 0, 0, 1, 0)
Chrysomelidae_sp1 = c(1, 0, 0, 0, 0, 0, 0)
Chrysomelidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Chrysomelidae_sp4 = c(0, 0, 1, 0, 0, 0, 1)
Curculionidae_sp1 = c(1, 0, 1, 0, 0, 0, 0)
Curculionidae_sp2 = c(0, 0, 0, 0, 0, 0, 1)
Dytiscidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Dryopidae_sp1 = c(1, 0, 1, 0, 0, 0, 0)
Dryopidae_sp2 = c(1, 0, 0, 0, 0, 0, 0)
Macrelmis_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Macrelmis_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Hexacylloepus_sp = c(0, 0, 1, 0, 0, 0, 0)
Lampyridae_larva_indet = c(0, 0, 1, 1, 1, 0, 0)
Ptilodactylidae_sp2 = c(0, 0, 1, 1, 0, 0, 1)
Canthon_sp = c(0, 0, 0, 0, 0, 0, 1)
Atheucus_sp = c(1, 0, 0, 0, 0, 0, 0)
Scarabaeidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Scydmaenidae_sp1 = c(0, 0, 0, 0, 0, 0, 1)
Aleocharinae_sp = c(0, 0, 1, 0, 0, 0, 1)
Holotrochus_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Oxytelinae_sp = c(0, 0, 1, 0, 0, 0, 0)
Oxytelinae_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Oxytelinae_sp2 = c(0, 0, 1, 1, 0, 0, 0)
Paederinae_sp4 = c(0, 0, 0, 0, 0, 0, 1)
Paederinae_sp5 = c(0, 0, 1, 0, 0, 0, 1)
Homeotarsus_sp9 = c(0, 0, 0, 1, 0, 0, 0)
Pselaphinae_sp1 = c(0, 0, 0, 1, 1, 0, 0)
Staphylinini_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Staphylinini_sp2 = c(0, 0, 0, 0, 0, 0, 1)
Xantholinini_sp7 = c(0, 0, 1, 0, 0, 0, 0)
Stenus_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Stenus_sp2 = c(0, 0, 0, 0, 0, 1, 0)
Tachyporinae_sp1 = c(0, 0, 0, 0, 0, 0, 1)
Tachyporinae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Tachyporinae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Staphylinidae_sp = c(1, 0, 1, 0, 0, 0, 1)
Staphylinidae_sp4 = c(0, 0, 1, 0, 0, 0, 1)
Staphylinidae_sp8 = c(0, 0, 0, 0, 1, 0, 1)
Staphylinidae_sp13 = c(0, 0, 1, 1, 1, 0, 1)
Staphylinidae_sp14 = c(0, 0, 0, 0, 1, 0, 0)
Staphylinidae_sp15 = c(0, 0, 1, 0, 0, 0, 0)
Staphylinidae_sp16 = c(0, 0, 1, 0, 0, 0, 0)
Staphylinidae_sp18 = c(0, 0, 1, 0, 0, 0, 1)
Staphylinidae_sp19 = c(0, 0, 0, 0, 0, 0, 1)
Staphylinidae_sp20 = c(0, 0, 0, 0, 0, 1, 1)
Staphylinidae_sp21 = c(0, 0, 0, 0, 0, 0, 1)
Tenebrionidae_sp1 = c(0, 1, 0, 0, 0, 1, 0)
Diptera_indet = c(0, 0, 1, 0, 0, 0, 0)
Diptera_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Diptera_sp2 = c(0, 1, 1, 0, 0, 0, 0)
Diptera_sp3 = c(0, 1, 1, 0, 0, 0, 0)
Diptera_sp4 = c(0, 0, 0, 1, 0, 0, 0)
Diptera_sp5 = c(0, 0, 0, 1, 0, 0, 0)
Diptera_sp6 = c(0, 0, 1, 1, 0, 0, 0)
Diptera_sp7 = c(0, 0, 1, 0, 0, 0, 0)
Diptera_sp9 = c(0, 0, 1, 0, 0, 0, 0)
Acalyptratae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Cecidomyiidae_sp1 = c(0, 1, 1, 0, 1, 0, 0)
Cecidomyiidae_sp2 = c(0, 1, 0, 0, 1, 0, 1)
Cecidomyiidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Ceratoponidae_sp1 = c(0, 0, 1, 0, 0, 1, 0)
Ceratoponidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Chironomidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Chironomidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Chironomidae_sp3 = c(0, 0, 1, 0, 0, 0, 1)
Chironomidae_sp4 = c(0, 0, 1, 0, 0, 0, 0)
Chironomidae_sp5 = c(0, 0, 1, 0, 0, 0, 0)
Culicidae_sp2	= c(0, 0, 0, 0, 0, 0, 1)
Drosophila_sp1 = c(1, 1, 0, 0, 0, 0, 1)
Drosophila_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Mycetophilidae_sp1 = c(0, 1, 0, 0, 0, 0, 0)
Mycetophilidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Mycetophilidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Muscidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Phoridae_sp1 = c(0, 0, 0, 0, 0, 1, 0)
Phoridae_sp2 = c(0, 0, 0, 0, 0, 0, 1)
Phoridae_sp4 = c(0, 0, 1, 0, 0, 0, 0)
Psychodidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Phlebotominae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Sciaridae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Sciaridae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Sciaridae_sp3 = c(0, 0, 1, 0, 0, 0, 1)
Sciaridae_sp4 = c(0, 0, 1, 0, 0, 0, 0)
Leptophlebidae_sp= c(0, 0, 1, 0, 0, 1, 0)
Leptophlebidae_sp1 = c(0, 1, 0, 0, 0, 0, 0)
Miroculis_sp= c(0, 0, 1, 0, 0, 0, 0)
Ephemeroptera_indet = c(1, 0, 0, 0, 0, 0, 0)
Hemiptera_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Hemiptera_sp3 = c(0, 0, 1, 1, 0, 0, 0)
Hemiptera_sp4 = c(0, 0, 1, 1, 0, 0, 0)
Cicadellidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Cicadellidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Cixiidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Zicca_sp= c(0, 0, 1, 0, 0, 0, 0)
Gelastocoridae_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Gelastocoridae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Hydrometridae_sp= c(0, 0, 0, 0, 0, 0, 1)
Lygaeidae_sp= c(0, 0, 1, 0, 0, 0, 1)
Lygaeidae_sp1 = c(0, 0, 0, 0, 0, 0, 1)
Lygaeidae_sp2 = c(0, 0, 0, 0, 0, 0, 1)
Lygaeidae_sp3 = c(0, 0, 0, 0, 0, 0, 1)
Notonectidae_sp= c(0, 0, 0, 1, 0, 0, 0)
Reduviidae_sp1 = c(0, 0, 1, 0, 0, 0, 0) 
Reduviidae_sp2 = c(0, 1, 1, 1, 0, 0, 0)
Reduviidae_sp4 = c(0, 0, 0, 0, 1, 0, 0)
Reduviidae_sp5 = c(0, 0, 1, 0, 0, 0, 0)
Zelurus_sp= c(0, 0, 1, 1, 0, 0, 1)
Emesinae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Emesinae_sp2 = c(0, 0, 1, 0, 0, 0, 1)
Emesinae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Harpactocorinae_indet = c(0, 0,	1, 0, 0, 0, 0)
Veliidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp4 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp5 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp6 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp8 = c(0, 0, 1, 0, 0, 0, 0)
Hymenoptera_sp12 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp4 = c(0, 0, 1, 0, 0, 0, 1)
Formicidae_sp5 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp6 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp9 = c(0, 0, 1, 0, 0, 0, 1)
Formicidae_sp10 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp11 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp13 = c(0, 0, 1, 0, 0, 0, 0)
Formicidae_sp14 = c(0, 0, 0, 0, 0, 0, 1)
Ichneumonoidea_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Vespidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Nasutitermitinae_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Noctuidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Tineidae_sp1 = c(0, 1, 1, 0, 0, 0, 1)
Megaloptera_larva_indet = c(0, 0, 0, 1, 0, 0, 0)
Myrmeleontidae_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Odonata_jovem_indet = c(0, 0, 0, 0, 0, 0, 1)
Orthoptera_sp1 = c(0, 0, 1, 0, 0, 1, 0)
Phalangopsidae_sp1 = c(0, 0, 1, 1, 1, 0, 0)
Phalangopsidae_sp2 = c(0, 0, 1, 0, 1, 0, 0)
Endecous_sp= c(0, 1, 0, 1, 0, 0, 1)
Endecous_sp1 = c(1, 0, 1, 0, 0, 0, 0)
Endecous_spn= c(0, 0, 1, 0, 0, 0, 0)
Eidmanacris_sp= c(0, 0, 1, 0, 0, 0, 0)
Cyrtoxipha_sp= c(0, 0, 0, 1, 0, 0, 0)
Perlidae_sp1 = c(0, 0, 0, 1, 0, 1, 0)
Perlidae_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Perlodidae_sp= c(0, 0, 1, 0, 0, 0, 0)
Psocoptera_indet = c(0, 0, 1, 0, 0, 0, 0)
Psocoptera_sp1 = c(0, 0, 1, 0, 0, 0, 0)
Psocoptera_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Psocoptera_sp3= c(0, 0, 1, 0, 0, 0, 0)
Lepismatidae_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Lepismatidae_sp2 = c(0, 0, 1, 0, 0, 0, 1)
Lepismatidae_sp3 = c(0, 0, 1, 0, 0, 0, 0)
Trichoptera_sp1 = c(0, 0, 1, 0, 0, 0, 1)
Trichoptera_sp2 = c(0, 0, 1, 0, 0, 0, 0)
Trichoptera_sp3 = c(0, 1, 1, 0, 0, 0, 1)
Hydrobiosidae_sp1_larva = c(0, 0, 1, 0, 0, 0, 0)
Hydroptilidae_sp1_larva = c(0, 0, 1, 0, 0, 0, 0)
Hydroptilidae_sp2_larva = c(0, 0, 1, 0, 0, 0, 0)
Hydropsychidae_sp1_larva = c(0, 0, 1, 0, 0, 0, 0)
Armadillidae_sp34 = c(0, 0, 1, 1, 0, 0, 0)
Dubioniscidae_sp36 = c(0, 1, 0, 0, 0, 0, 0)	
Dubioniscidae_sp42 = c(0, 0, 0, 0, 0, 1, 1)
Dubioniscidae_sp43 = c(0, 0, 1, 0, 0, 1, 1)
Platyarthridae_sp37 = c(0, 0, 1, 0, 0, 0, 1)
Styloniscidae_sp35 = c(0, 0, 1, 0, 0, 0, 0)

community=rbind(Ityphilus_sp6, Schendylops_sp1, Schendylops_sp3, Lamyctes_sp,  Scolopendra_viridicornis, Sphendononema_guildingii, Polydesmida_sp1, Polydesmida_sp3, Pseudonannolene_imbirenses, Pseudonannolenidae_sp, Pseudonannolenidae_sp2, Pseudonannolenidae_sp3, Pseudonannolenidae_sp4, Pseudonannolenidae_sp5, Pseudonannolenidae_sp6, Pseudonannolenidae_sp7, Spirostreptidae_sp12, Siphonophorida_sp1, Blattaria_sp1, Blattaria_sp2, Blattaria_sp3, Blattaria_sp5, Blableridae_sp1, Blattellidae_sp1, Blattellidae_sp2, Blattellidae_sp3, Polyphagidae_sp1, Coleoptera_indet, Coleoptera_sp1, Coleoptera_sp2, Coleoptera_sp3, Coleoptera_sp4, Coleoptera_sp5, Coleoptera_sp6, Coleoptera_sp7, Coleoptera_sp8, Coleoptera_sp9, Coleoptera_sp10, Coleoptera_sp11, Coleoptera_sp12, Coleoptera_sp13, Coleoptera_sp14, Coleoptera_sp15, Coleoptera_sp16, Coleoptera_sp17, Coleoptera_sp18, Coleoptera_sp20, Cantharidae_sp1, Cantharidae_sp2, Cantharidae_sp3, Carabidae_indet1, Pentagonica_sp1, Platynus_sp1, Pterostichini_sp2, Agonina_sp1, Agonina_sp2, Loxandrus_sp1, Paratachys_sp1, Paratachys_sp2, Pericompsus_sp1, Polyderis_sp1, Tachys_sp1, Carabidae_sp1, Carabidae_sp2, Carabidae_sp3, Zuphiini_sp2, Chrysomelidae_sp1, Chrysomelidae_sp3, Chrysomelidae_sp4, Curculionidae_sp1, Curculionidae_sp2, Dytiscidae_sp1, Dryopidae_sp1, Dryopidae_sp2, Macrelmis_sp1, Macrelmis_sp2, Hexacylloepus_sp, Lampyridae_larva_indet, Ptilodactylidae_sp2, Canthon_sp, Atheucus_sp, Scarabaeidae_sp2, Scydmaenidae_sp1, Aleocharinae_sp, Holotrochus_sp2, Oxytelinae_sp, Oxytelinae_sp1, Oxytelinae_sp2, Paederinae_sp4, Paederinae_sp5, Homeotarsus_sp9, Pselaphinae_sp1, Staphylinini_sp1, Staphylinini_sp2, Xantholinini_sp7, Stenus_sp1, Stenus_sp2, Tachyporinae_sp1, Tachyporinae_sp2, Tachyporinae_sp3, Staphylinidae_sp, Staphylinidae_sp4, Staphylinidae_sp8, Staphylinidae_sp13, Staphylinidae_sp14, Staphylinidae_sp15, Staphylinidae_sp16, Staphylinidae_sp18, Staphylinidae_sp19, Staphylinidae_sp20, Staphylinidae_sp21, Tenebrionidae_sp1, Diptera_indet, Diptera_sp1, Diptera_sp2, Diptera_sp3, Diptera_sp4, Diptera_sp5, Diptera_sp6, Diptera_sp7, Diptera_sp9, Acalyptratae_sp3, Cecidomyiidae_sp1, Cecidomyiidae_sp2, Cecidomyiidae_sp3, Ceratoponidae_sp1, Ceratoponidae_sp3, Chironomidae_sp1, Chironomidae_sp2, Chironomidae_sp3, Chironomidae_sp4, Chironomidae_sp5, Culicidae_sp2, Drosophila_sp1, Drosophila_sp2, Mycetophilidae_sp1, Mycetophilidae_sp2, Mycetophilidae_sp3, Muscidae_sp1, Phoridae_sp1, Phoridae_sp2, Phoridae_sp4, Psychodidae_sp2, Phlebotominae_sp1, Sciaridae_sp1, Sciaridae_sp2, Sciaridae_sp3, Sciaridae_sp4, Leptophlebidae_sp, Leptophlebidae_sp1, Miroculis_sp, Ephemeroptera_indet, Hemiptera_sp2, Hemiptera_sp3, Hemiptera_sp4, Cicadellidae_sp1, Cicadellidae_sp2, Cixiidae_sp1, Zicca_sp, Gelastocoridae_sp1, Gelastocoridae_sp2, Hydrometridae_sp, Lygaeidae_sp, Lygaeidae_sp1, Lygaeidae_sp2, Lygaeidae_sp3, Notonectidae_sp, Reduviidae_sp1, Reduviidae_sp2, Reduviidae_sp4, Reduviidae_sp5, Zelurus_sp, Emesinae_sp1, Emesinae_sp2, Emesinae_sp3, Harpactocorinae_indet, Veliidae_sp1, Hymenoptera_sp2, Hymenoptera_sp3, Hymenoptera_sp4, Hymenoptera_sp5, Hymenoptera_sp6, Hymenoptera_sp8, Hymenoptera_sp12, Formicidae_sp4, Formicidae_sp5, Formicidae_sp6, Formicidae_sp9, Formicidae_sp10, Formicidae_sp11, Formicidae_sp13, Formicidae_sp14, Ichneumonoidea_sp1, Vespidae_sp1, Nasutitermitinae_sp1, Noctuidae_sp1, Tineidae_sp1, Megaloptera_larva_indet, Myrmeleontidae_sp1, Odonata_jovem_indet, Orthoptera_sp1, Phalangopsidae_sp1, Phalangopsidae_sp2, Endecous_sp, Endecous_sp1, Endecous_sp, Eidmanacris_sp, Cyrtoxipha_sp, Perlidae_sp1, Perlidae_sp2, Perlodidae_sp, Psocoptera_indet, Psocoptera_sp1, Psocoptera_sp2, Psocoptera_sp3, Lepismatidae_sp1, Lepismatidae_sp2, Lepismatidae_sp3, Trichoptera_sp1, Trichoptera_sp2, Trichoptera_sp3, Hydrobiosidae_sp1_larva, Hydroptilidae_sp1_larva, Hydroptilidae_sp2_larva, Hydropsychidae_sp1_larva, Armadillidae_sp34, Dubioniscidae_sp36, Dubioniscidae_sp42, Dubioniscidae_sp43, Platyarthridae_sp37, Styloniscidae_sp35) 
colnames(community)=c("Passa_Tres", "S.Bernardo", "Angelica", "Bezerra", "S.Mateus", "T.Ronca_I", "T.Ronca_II")

beals = beals(community, reference = community, type = 0, include = FALSE) 

presences = which(community == 1) 
beals.presences = beals[presences] 
beals.presences= as.matrix(beals.presences)
length(beals.presences) # N=341

absences = which(community == 0) 
beals.absences = beals[absences]  
beals.absences= as.matrix(beals.absences)
length(beals.absences) # N=1276

beals.presences=as.numeric(beals.presences) 
pres= beals.presences
h=hist(pres) 
h$counts= cumsum(h$counts)/341 
beals.absences=as.numeric(beals.absences)

abs= beals.absences
i=hist(abs)
i$counts= cumsum(i$counts)/1276

rank.abs = rank (c(beals.absences, beals.presences))
summary(rank.abs)#  threshold= median (809.5)
threshold=beals.absences[809.5]#
threshold

isopoda= beals[c("Armadillidae_sp34", "Dubioniscidae_sp36", "Dubioniscidae_sp42", "Dubioniscidae_sp43", "Platyarthridae_sp37", "Styloniscidae_sp35"),]
comm.isopoda=community[c("Armadillidae_sp34", "Dubioniscidae_sp36", "Dubioniscidae_sp42", "Dubioniscidae_sp43", "Platyarthridae_sp37", "Styloniscidae_sp35"),]
dark.div = ifelse(isopoda>threshold & comm.isopoda==0,TRUE,FALSE)