Greetings,
My questions are related to producing confidence intervals via
mcmcsamp and HPDinterval. As you can see in the output below, my
model has several fixed effect parameters. Of interest to me are each
of the "Treatment" vs. "Control" comparisons for each "site"(in each
fixed effect parameter name, these are specified by the text
immediately following "sitett"). As I understand it, if I want to
produce a CI for each "Treatment" parameter from the anova output, I
need to simply subtract the Control estimate from the Treatment
estimate for each site. The SE's are a little more complicated, but
that's not important here. This is easy enough. However, if I want
to use the mcmcsamp+HPDinterval combo to produce my CI's, I unsure how
to proceed. The HPDinterval output below lists upper and lower bounds
for each fixed effect parameter, but I have no idea how to get the
"Treatment" vs. "Control" comparison for each site when working with
just the upper and lower bounds of the CI. Any advice on how to
proceed here would be much appreciated.
Thanks,
Gus