#SOSCC algorithm #George M. Garrity and Timothy G. Lilburn #Copyright Michigan State University 2004, all rights reserved "step10" = function(taxon.subset, dis.table) { all.genus <- substring(as.character(taxon.subset[, 4]), 1, regexpr("[^a-zA-Z]", as.character(taxon.subset[, 4])) - 1) taxon.level <- unique(all.genus) taxon.seq.new <- vector("list", length(taxon.level)) names(taxon.seq.new) <- taxon.level misidentified <- NULL for(i in 1:length(unique(taxon.level))) { if(sum(all.genus == taxon.level[[i]]) >= 2) { temp <- dis.table[all.genus == taxon.level[[i]], all.genus == taxon.level[[i]]] misidentified <- c(misidentified, dimnames(temp)[[1]][abs(apply(scale(temp), 1, max)) > 2.5]) } } return(misidentified) }