results - bpp/bpp-tutorial-geneflow GitHub Wiki
mkdir -p ~/DAY-4/bf-results/model-{1..6}
cd ~/DAY-4/bf-results
# download betaweights.csv files for each model
for i in {1..6}; do cd model-$i; wget https://raw.githubusercontent.com/bpp/bpp-tutorial-geneflow/main/fourth-day/output/bf-betaweights/model-$i/betaweights.csv; cd ../; done
# calculate marginal likelihoods
for i in {1..6}; do cut -f 2,3 -d',' model-$i/betaweights.csv | tail -n +2 | awk -F ',' '{printf "%.6f\n", $1 * $2}' | paste -sd+ | bc -l | awk '{printf "%.6f\n", $1/2}'; done
| Model |
Marginal log-likelihood |
| 1 |
-452696.201629 |
| 2 |
-452697.409983 |
| 3 |
-452690.896586 |
| 4 |
-452694.108056 |
| 5 |
-452696.876828 |
| 6 |
-452695.399386 |