#!/bin/sh

DATA="sample.dat"
BUGFILE="inf_mu_sigma_pred"

# number of data points (needed in .bug)
NNNN=`wc -l $DATA | awk '{print $1}'`
echo $NNNN
sed s/NNNN/$NNNN/ $BUGFILE.txt > $BUGFILE.bug

# set .in file 
echo 'source("set_init.R")' | R --slave

# set command file and run BUGS
CMDS="inf_mu_sigma_pred_cmd.txt"
BUGS="/home/giulio/w/bugs/bugs06/bugs"
NBURNIN=10000
NSAMPLE=100000

sed s/NBURNIN/$NBURNIN/ $CMDS | sed s/NSAMPLE/$NSAMPLE/  | $BUGS