-
Notifications
You must be signed in to change notification settings - Fork 0
/
full_analysis_JC.Rev
36 lines (32 loc) · 1.13 KB
/
full_analysis_JC.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
################################################################################
#
# 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 = "pps_example"
data <- readDiscreteCharacterData(inFile)
model_name = "JC"
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 = "pps_example"
data <- readDiscreteCharacterData(inFile)
model_name = "JC"
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 = "pps_example"
data <- readDiscreteCharacterData(inFile)
model_name = "JC"
model_file_name = "scripts/"+model_name+"_Model.Rev"
source("scripts/PosteriorPredictive_MCMC.Rev")
clear()
q()