search for: tgenus

..., I am working with taxonomic data, represented as a list of classes, orders, families, genera and finally species. > class(mydata) [1] "data.frame" > mode(mydata) [1] "list" > names(mydata) [1] "tclass" "torder" "tfamily" "tgenus" "tspecies" > length(mydata$tclass) [1] 161590 The first 10 rows look like the following: > mydata[1:10,] tclass torder tfamily tgenus 1 Chlorophyta Chlorophyceae Dunaliellaceae Collodictyon 2 Chlorophyta Chlorophyceae Dunaliella...

..." ... $ tclass : Factor w/ 8 levels "Annelida","Arachnida",..: 1 5 5 5 5 ... $ torder : Factor w/ 18 levels "Achtinedida",..: 13 5 5 10 13 5 7 ... $ tfamily : Factor w/ 81 levels "","Ameletidae",..: 79 46 46 14 42 ... $ tgenus : Factor w/ 206 levels "","Acentrella",..: 1 10 10 1 140 ... $ tspecies : Factor w/ 60 levels "","aequalis",..: 1 1 1 1 1 1 1 ... $ quant : int 22 527 22 22 11 97 333 11 108 11 ... $ stream : Factor w/ 7 levels "BCrk",&...

...stions? Thanks, Olga ---------------- result<-list() for(i in 1:1000)try({ hclass<-mydata2$tclass[mydata2$tspecies%in%com[[i]][[5]][[1]]] horder<-mydata2$torder[mydata2$tspecies%in%com[[i]][[5]][[1]]] hfamily<-mydata2$tfamily[mydata2$tspecies%in%com[[i]][[5]][[1]]] hgenus<-mydata2$tgenus[mydata2$tspecies%in%com[[i]][[5]][[1]]] hspecies<-mydata2$tspecies[mydata2$tspecies%in%com[[i]][[5]][[1]]] hier<-cbind(hspecies,hgenus,hfamily,horder,hclass) s<-nrow(hier) s<-ifelse(s>1000,s/10,s) taxa<-taxa2dist(hier[sample(s,replace=F),]) abd<-t(cbind(c(1:s),rep(1,s))) taxd...