Question: RAxML for each gene family
0
gravatar for spaul8505
2.5 years ago by
spaul850510
United States
spaul850510 wrote:

I would like to obtain the RAXML ML tree for each gene family. I have around 800 .phy trimmed alignments. I tried looping the .phy files to produce one tree each but somehow only one tree is produced in the end. Should I create sub-directories to store them or is there a better way to do this? I am specifically wanting to make single trees and hence the question.

So each gene family is named as

Ortho1.phy Ortho2.Phy ....

I tried using

for f in /home/Single_trees/trimmed_alignment/*.phy; do g=$(echo $f | sed 's/.phy$//g'); raxmlHPC -f a -x 100 -m PROTGAMMAAUTO -p 100 -s $f -N 100 -n $g.tree; done;

I wanted to name each tree based on the gene family name now but htis always gives me an error saying $ symbol is not allowed.

raxml • 658 views
ADD COMMENTlink modified 2.5 years ago by Caesar10 • written 2.5 years ago by spaul850510
1
gravatar for Caesar
2.5 years ago by
Caesar10
Sweden/Uppsala
Caesar10 wrote:

Try this:

for f in /home/Single_trees/trimmed_alignment/*.phy
do 
pathName=${f%.*}
raxmlHPC -f a -x 100 -m PROTGAMMAAUTO -p 100 -s $f -N 100 -n ${pathName}.tree
done
ADD COMMENTlink written 2.5 years ago by Caesar10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1077 users visited in the last hour