-
Notifications
You must be signed in to change notification settings - Fork 0
/
full_analysis_GTR.Rev
37 lines (33 loc) · 1.13 KB
/
full_analysis_GTR.Rev
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
################################################################################
#
# RevBayes Example: Bayesian model testing using posterior predictive simulation
#
# authors: Lyndon M. Coghill, Sebastian Hoehna and Jeremy M. Brown
#
################################################################################
## EMPIRICAL MCMC
inFile = "data/primates_and_galeopterus_cytb.nex"
analysis_name = "primates"
data <- readDiscreteCharacterData(inFile)
model_name = "GTR"
model_file_name = "scripts/"+model_name+"_Model.Rev"
source("scripts/MCMC_Simulation.Rev")
clear()
## POSTERIOR PREDICTIVE DATA SIMULATION
inFile = "data/primates_and_galeopterus_cytb.nex"
analysis_name = "primates"
data <- readDiscreteCharacterData(inFile)
model_name = "GTR"
model_file_name = "scripts/"+model_name+"_Model.Rev"
source("scripts/PosteriorPredictive_Simulation.Rev")
clear()
## POSTERIOR PREDICTIVE MCMC
inFile = "data/primates_and_galeopterus_cytb.nex"
analysis_name = "primates"
data <- readDiscreteCharacterData(inFile)
model_name = "GTR"
model_file_name = "scripts/"+model_name+"_Model.Rev"
source("scripts/PosteriorPredictive_MCMC.Rev")
clear()
# END IT ALL
q()