I am trying to see which OTUs are common among all 6 samples using heatmap.bin. Is there a way to add a proper label on the “Y-axis” instead of them just outputing OTU#####? Thank you very much.
Below is a copy of my scripts:
*cluster.split(fasta=sb_stability.trim.contigs.good.unique.good.filter.unique.precluster.pick.pick.pick.fasta, count=sb_stability.trim.contigs.good.unique.good.filter.unique.precluster.pick.pick.count_table, taxonomy=sb_stability.trim.contigs.good.unique.good.filter.unique.precluster.pick.gg.wang.pick.pick.taxonomy, splitmethod=classify, taxlevel=6, cutoff=0.01, processors=64)*
*make.shared(list=current, count=current, label=0.01)*
*classify.otu(list=current, count=current, taxonomy=current, label=0.01)*
*get.oturep(fasta=current, count=current, list=current, method=abundance)*
* *
*count.groups(shared=current)*
* *
*sub.sample(shared=current, size=10000)*
* *
*heatmap.bin(shared=sb_stability.trim.contigs.good.unique.good.filter.unique.precluster.pick.pick.pick.opti_mcc.0.01.subsample.shared, numotu=100)*