Displaying 1 result from an estimated 1 matches for "gaprior".
2011 Sep 19
2
Poisson-Gamma computation (parameters and likelihood)
...r rates which has a conjugate Gamma
prior.
1) The first question is precisely how I work out the parameters.
#Suppose I assign values to theta with *seq()*
*theta<-seq(0,1,len=500)*
#Then I try out the parameters that seem to fit with a certain prior idea on
theta (see next)
*a=182*
*b=3530*
*gaprior<-dgamma(theta,a,b)*
*
*
It should work by trial-and-error (according to "Bayesian Computation with
R") , but how can I check the parameters turned out well: should I just look
at the plot, or evaluate it through the *1 - pgamma(x,a,b) *function, having
knowledge of the 5th percentile (...