PartitionFinder

Descrição

De acordo com a página do PartitionFinder, PartitionFinder é um software livre e open-source para selecionar esquemas de partição e modelos de evolução molecular para análises filogenéticas.

Versões Disponíveis

  • partitionfinder/2.1.1 (default)

Carregando o Módulo

# Carregar PartitionFinder
module load partitionfinder/2.1.1

# Variáveis disponíveis
echo $PartitionFinder           # PartitionFinder.py
/opt/gridunesp/dist/partitionfinder/2.1.1/bin/PartitionFinder.py

echo $PartitionFinderProtein    # PartitionFinderProtein.py
/opt/gridunesp/dist/partitionfinder/2.1.1/bin/PartitionFinderProtein.py

echo $PartitionFinderMorphology # PartitionFinderMorphology.py
/opt/gridunesp/dist/partitionfinder/2.1.1/bin/PartitionFinderMorphology.py

Submissão de Jobs Seriais

submit_partitionfinder.sh
#!/bin/bash
#SBATCH -J partitionfinder
#SBATCH -N 1
#SBATCH -n 1
#SBATCH -t 23:00:00
#SBATCH --mem=8G

export INPUT="examples/nucleotide"
export OUTPUT="pf_results/"

module load partitionfinder/2.1.1

job-nanny python $PartitionFinder examples/nucleotide

Exemplo de Arquivo de Configuração

partition_finder.cfg
## ALIGNMENT FILE ##
alignment = my_alignment.nex;

## BRANCHLENGTHS ##
branchlengths = linked;

## MODELS OF EVOLUTION ##
models = all;

## MODEL SELECCTION ##
model_selection = aicc;

## DATA BLOCKS ##
[data_blocks]
gene1 = 1-1000;
gene2 = 1001-2000;

## SCHEMES ##
[schemes]
scheme1 = gene1, gene2;

## USER SCHEMES ##
[user schemes]

Job Array para Múltiplos Alinhamentos

submit_pf_array.sh
#!/bin/bash
#SBATCH -J pf_array
#SBATCH --array=1-10
#SBATCH -N 1
#SBATCH -n 1
#SBATCH -t 24:00:00
#SBATCH --mem=8G

ALIGNMENTS=(
    "gene1.nex"
    "gene2.nex"
    "gene3.nex"
    "gene4.nex"
    "gene5.nex"
    "gene6.nex"
    "gene7.nex"
    "gene8.nex"
    "gene9.nex"
    "gene10.nex"
)

ALIGNMENT=${ALIGNMENTS[$SLURM_ARRAY_TASK_ID-1]}
export INPUT="$ALIGNMENT partition_finder.cfg"
export OUTPUT="pf_${ALIGNMENT%.nex}/"

module load partitionfinder/2.1.1

mkdir -p pf_${ALIGNMENT%.nex}
cd pf_${ALIGNMENT%.nex}
cp ../$ALIGNMENT .
cp ../partition_finder.cfg .

job-nanny python $PartitionFinder . > pf_output.log

Análise de Resultados

analyze_pf.sh
#!/bin/bash
#SBATCH -J analyze_pf
#SBATCH -N 1
#SBATCH -n 1
#SBATCH -t 00:30:00
#SBATCH --mem=2G

# Extrair melhor esquema
grep "Best scheme" */*.log > best_schemes.txt

# Extrair scores AICc
grep "AICc score" */*.log | awk '{print $4}' > aicc_scores.txt

# Extrair modelos selecionados
grep "Model for partition" */*.log > models.txt

# Gerar resumo
echo "Resumo das análises:" > summary.txt
echo "=====================" >> summary.txt

for dir in pf_*/; do
    echo "Diretório: $dir" >> summary.txt
    grep "Best scheme" ${dir}partition_finder.log >> summary.txt
    grep "AICc score" ${dir}partition_finder.log >> summary.txt
    echo "" >> summary.txt
done

Referências

Ver também