Open ghost opened 8 years ago
From stephane.
require("coda") fem_joint_A=read.csv("HB_female_joint_WTA_A.csv") geweke.diag(fem_joint_A$beta_sal_own_WTA, frac1=0.1, frac2=0.5) geweke.diag(fem_joint_A$beta_sal_par_WTA, frac1=0.1, frac2=0.5)
Does this work? A is a mean of some parameter samples and not samples themselves
From stephane.
require("coda") fem_joint_A=read.csv("HB_female_joint_WTA_A.csv") geweke.diag(fem_joint_A$beta_sal_own_WTA, frac1=0.1, frac2=0.5) geweke.diag(fem_joint_A$beta_sal_par_WTA, frac1=0.1, frac2=0.5)