/* Semi-conjugate prior with normal likelihood (Section 9.4)
 */

#include "random.h"
#include <iomanip.h>
#include <iomanip.h> */
#include <fstream.h> */
using namespace std;

const int iter=10;   /* Number of iterations of the EM algorithm  */

const int m=500;     /* Number of replications                    */

const int t=10;      /* Number of iterations                      */

long idum=65537;     /* Seed for random number generator          */

int main (void)
{
  ofstream ans("semiconj.ans");
  int j,s,nu,nu0,nu1,n,n0,n1;
  double xbar,sxx,s0,s1;
  double theta,theta0,theta1,thetabar,thetass,thetavar;
  double phi,phi0,phi1,phibar,phiss,phivar;
  n=12;
  xbar=119;
  sxx=13045;
  s0=2700;
  nu0=11;
  n0=15;
  theta0=110;
  phi0=s0/(n0*(nu0-2));
  thetabar=0;
  phibar=0;
  thetass=0;
  phiss=0;
  ans << endl;
  ans << "Data quoted in P M Lee, `Bayesian Statistics: "
      << "An Introduction'," << endl;
  ans << "Arnold 1989, Section 2.13.  "
      << "Taking n=12, xbar=139, S=13,045 and"< */