/cons_test/update_plots.R

 1#!/usr/bin/env Rscript
 2
 3source("common.R");
 4
 5
 6# F84 model parameters:
 7f84.true.kappa<-1;
 8f84.true.base.freqs<-c(1.5/6,1.25/6,1.5/6,1.25/6);
 9f84.true.gamma.shape<-0.5;
10
11# model parameters:
12wag.true.gamma.shape<-0.1;
13
14# GY94 model parameters:
15gy94.true.omega<-2;
16gy94.true.kappa<-4;
17
18
19source("nucleotide_funcs.R");
20res<-read.table(file="results/constest_nucleotide.tab");
21pdf("results/constest_nucleotide.pdf");
22plot_dframe(res,tree);
23dev.off();
24
25source("aa_funcs.R");
26res<-read.table(file="results/constest_aa.tab");
27pdf("results/constest_aa.pdf");
28plot_dframe(res,tree);
29dev.off();
30
31source("codon_funcs.R");
32res<-read.table(file="results/constest_codon.tab");
33pdf("results/constest_codon.pdf");
34plot_dframe(res,tree);
35dev.off();
36
37
38
39