MOTU Cleaning and Analyses Script

This script removes cleans the raw MOTU data and runs all of the analyses in the paper including within-sample metazoan MOTU accumulation curves, sample richness estimates, site metazoan MOTU accumulation curves by sample number or read count, dissimilarity NMDS ordination, dissimilarity by distance analysis for Adelaide sites, and metazoan MOTU richness estimates standardize to 20 samples across all sites. This script was written for R version 3.4.1, requires input files "Unfiltered_18S.csv", "Unfiltered_COI.csv", "SampleIDs_18S.csv", "SampleIDs_COI.csv", and "geo dists adelaide.csv", and requires packages "vegan", "reshape2", and 'iNEXT".

本脚本用于清理原始MOTU（分子操作分类单元，Molecular Operational Taxonomic Unit）数据，并完成论文涵盖的全部分析内容，具体包括：样本内后生动物MOTU累积曲线、样本丰富度估计、按样本数量或读长计数划分的位点后生动物MOTU累积曲线、非相似性NMDS（非度量多维标度，Non-metric Multidimensional Scaling）排序、阿德莱德研究位点的距离衰减相似性分析，以及将所有位点的后生动物MOTU丰富度统一标准化至20个样本的分析。本脚本适配R 3.4.1版本，所需输入文件包括"Unfiltered_18S.csv"、"Unfiltered_COI.csv"、"SampleIDs_18S.csv"、"SampleIDs_COI.csv"及"geo dists adelaide.csv"，需依赖"vegan"、"reshape2"与"iNEXT"三个R包。