用于系统发育树的循环

Question

我正在使用一个称为raxml的系统发育软件，而我wnat为每个Phylip文件构建单个树。对于具有三个Phylip文件的目录，我进行了以下

 ##files in directory
 Ortho1.phy Ortho6.Phy Ortho6.Phy
for f in /home/Single_trees/trimmed_alignment/*.phy; do raxmlHPC -f a -x 100 -m PROTGAMMAAUTO -p 100 -s $f -N 100 -n $f.tree; done;

但这总是给我一个错误，说不允许$符号。

raxmlHPC: axml.c:5236: analyzeRunId: Assertion `0' failed.

错误字符/在运行ID中不允许

有更好的方法吗？我尝试使用此链接使用此链接在此处使用作业数组https://rc.fas.harvard.edu/resources/documentation/submittit-large-numbers-numbers-obs-jobs-jobs-to-odyssey/它。

这是我尝试的数组作业提交的方法：

 #!/bin/bash -l
 #
 # raxml.sbatch
 #
 #SBATCH -J consensus       # A single job name for the array
 #SBATCH -p high # best partition for single core small jobs
 #SBATCH -n 12              # one core
 #SBATCH -N 1              # on one node
 #SBATCH -t 100:00:00         # Running time of 2 hours
 #SBATCH --mem 18000        # Memory request of 4 GB
 #SBATCH -o raxml_%A_%a.out # Standard output
 #SBATCH -e raxml_%A_%a.err # Standard error
 module load raxml
  for FILES in /home/aligned_fasta/.phy; do
   echo ${FILES}
  done;
  # grab out filename from the array exported from our 'parent' shell
  FILENAME=${FILES[$SLURM_ARRAY_TASK_ID]}
  # make & move into new directory, and run!
 mkdir ${FILENAME}_out
 cd ${FILENAME}_out
 raxmlHPC -f a -x 100 -m PROTGAMMAAUTO -p 100  -s $FILENAME -N 100 -n $FILENAME.tree 
 #Now, we grab all the appropriate files and submit them en-batch with an  array:
 # grab the files, and export it so the 'child' sbatch jobs can access it

export files =（$（LS -1 .phy））

 # get size of array
 NUMPHY=${#FILES[@]}
 # now subtract 1 as we have to use zero-based indexing (first cell is 0)
 ZBNUMPHY=$(($NUMPHY - 1))
 # now submit to SLURM
 if [ $ZBNUMPHY -ge 0 ]; then
 sbatch --array=0-$ZBNUMPHY raxml.sbatch 
 fi

我使用sbatch -array = 0-10 raxml.sh提交，但它不起作用。

Answer 1

刚刚发现了一些东西。基本上，我将重命名文件，以便它们是连续的，并且可以使用slurm。

 ls *.phy | cat -n | while read num file; do mv $file ${file/./.$num.}; done

，文件将为

 Ortho1.1.phy Ortho6.2.Phy Ortho6.3.Phy

然后，您可以按照以下方式进行操作：

#!/bin/bash -l
# SBATCH -J tree
###### Standard out and Standard Error output files with the job number in the name.
#SBATCH -o tre_%A.%a.out
#SBATCH -e tre_%A.%a.err
###### number of nodes
###SBATCH --nodes=6
###SBATCH --nodes=6
###### number of processors
#SBATCH -n 16
###SBATCH --cpus-per-task=4
###### Spread the tasks evenly among the nodes
####BATCH --ntasks-per-node=8
###### coupled with array
####SBATCH --ntasks=1-179
#SBATCH --time=300:00:00
#SBATCH -p high
#SBATCH --mem 24000 
###### Want the node exclusively
### SBATCH --exclusive
#SBATCH --array=1-3
module load raxml
for i in $SLURM_ARRAY_TASK_ID.phy
do
echo $i
done
### tree
 raxmlHPC-PTHREADS -f a -x 100 -m PROTGAMMAAUTO -p 100 -T 16 -s $i -N 100 -n $i.tree 

这将为您提供单独的树，您可以用来建造共识树。

Answer 2

我认为您可能正在引导并获得共识系统发育，在这种情况下，RAXML中有一个特殊情况可以做到这一点。当我有分钟时，我会在这里发布。