#!/bin/bash
#SBATCH -t 12:00:00
#SBATCH --nodes=2 --ntasks-per-node=2
#SBATCH --export=NONE
#SBATCH --mem=96G
#SBATCH --mail-user=alex_labossiere@uri.edu
#SBATCH --mail-type=BEGIN,END,FAIL

module load QIIME2/2023.7

procdir=/data/mramseylab/16S/proc_reads/2023_AL_Master/pseudo-trip-data/core-metric-results-An
metadir=/data/mramseylab/16S/metadata/plaque_meta_An.tsv
visdir=/data/mramseylab/16S/visualizations/AL_2023

qiime evident multivariate-power-analysis \
  --i-data $procdir/bray_curtis_distance_matrix.qza \
  --m-sample-metadata-file $metadir \
  --p-group-column media \
  --p-alpha 0.01 0.05 0.1 \
  --p-total-observations $(seq 5 1 25) \
  --o-power-analysis-results $visdir/results_evident_bray.qza 

qiime evident plot-power-curve \
    --i-power-analysis-results /$visdir/results_evident_bray.qza \
    --p-target-power 0.8 \
    --p-style alpha \
    --o-visualization $visdir/power_curve_AN2.qzv