#!/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-J
metadir=/data/mramseylab/16S/metadata/plaque_meta_J.tsv
visdir=/data/mramseylab/16S/visualizations/AL_2023


#!/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-A
metadir=/data/mramseylab/16S/metadata/plaque_meta_A.tsv
visdir=/data/mramseylab/16S/visualizations/AL_2023


declare -a files=("faith_pd_vector") # Replace with your filenames minus the .qza

for file in "${files[@]}"; do
  qiime evident univariate-power-analysis \
    --m-sample-metadata-file $procdir/${file}.qza \
    --m-sample-metadata-file $metadir \
    --p-data-column faith_pd \
    --p-group-column media \
    --p-alpha 0.01 0.05 0.1 \
    --p-total-observations $(seq 2 1 24) \
    --o-power-analysis-results $visdir/${file}_alpha_results-J.qza

  qiime evident plot-power-curve \
    --i-power-analysis-results $visdir/${file}_alpha_results-J.qza \
    --p-target-power 0.8 \
    --p-style alpha \
    --o-visualization $visdir/${file}_power_curve_alpha-J.qzv
done

declare -a files=("shannon_vector") # Replace with your filenames minus the .qza

for file in "${files[@]}"; do
  qiime evident univariate-power-analysis \
    --m-sample-metadata-file $procdir/${file}.qza \
    --m-sample-metadata-file $metadir \
    --p-data-column shannon \
    --p-group-column media \
    --p-alpha 0.01 0.05 0.1 \
    --p-total-observations $(seq 2 1 24) \
    --o-power-analysis-results $visdir/${file}_alpha_results-J.qza

  qiime evident plot-power-curve \
    --i-power-analysis-results $visdir/${file}_alpha_results-J.qza \
    --p-target-power 0.8 \
    --p-style alpha \
    --o-visualization $visdir/${file}_power_curve_alpha-J.qzv
done

declare -a files=("evenness_vector") # Replace with your filenames minus the .qza

for file in "${files[@]}"; do
  qiime evident univariate-power-analysis \
    --m-sample-metadata-file $procdir/${file}.qza \
    --m-sample-metadata-file $metadir \
    --p-data-column pielou_e \
    --p-group-column media \
    --p-alpha 0.01 0.05 0.1 \
    --p-total-observations $(seq 2 1 24) \
    --o-power-analysis-results $visdir/${file}_alpha_results-J.qza

  qiime evident plot-power-curve \
    --i-power-analysis-results $visdir/${file}_alpha_results-J.qza \
    --p-target-power 0.8 \
    --p-style alpha \
    --o-visualization $visdir/${file}_power_curve_alpha-J.qzv
done