We will be hosting mothur and R workshops throughout 2018. Learn more.


From mothur
Jump to: navigation, search

The corr.axes command will calculate the correlation coefficient for each column in a shared/relabund file to the axes displayed in a pcoa file.

Default settings

The corr.axes command reads a shared, relabund or metadata file as well as an axes file. If a shared file is given, the relabund is calculated. A metadata file has the same number of rows as samples, but the column(s) are floats that describe the samples (e.g. temperature, weight, etc.)

mothur > corr.axes(shared=genus.pool.shared, axes=genus.pool.thetayc.genus.lt.pcoa)


mothur > corr.axes(relabund=genus.pool.relabund, axes=genus.pool.thetayc.genus.lt.pcoa)


 mothur > corr.axes(metadata=genus.pool.metadata, axes=genus.pool.thetayc.genus.lt.pcoa)



The groups parameter allows you to specify which of the groups you would like included. The group names are separated by dashes. By default all groups in the shared or relabund file are used.


The label parameter allows you to select what distance level you would like used, if none is given the first distance is used.


The method parameter allows you to select what method you would like to use. Options are pearson, spearman and kendall. Default=pearson - http://en.wikipedia.org/wiki/Pearson_correlation_coefficient

To use the spearman method: - http://en.wikipedia.org/wiki/Spearman's_rank_correlation_coefficient

mothur > corr.axes(metadata=genus.pool.metadata, axes=genus.pool.thetayc.genus.lt.pcoa, method=spearman)

or to use the kendall method: - http://en.wikipedia.org/wiki/Kendall_tau_rank_correlation_coefficient

mothur > corr.axes(metadata=genus.pool.metadata, axes=genus.pool.thetayc.genus.lt.pcoa, method=kendall)


The numaxes parameter allows you to select the number of axes you would like to use. Default=3.

mothur > corr.axes(metadata=genus.pool.metadata, axes=genus.pool.thetayc.genus.lt.pcoa, numaxes=2)


  • 1.24.0 - added p_values.
  • 1.25.0 - fixed OTU labels.