Cost function:
J(p) = 1/2*[(Yo-Y(p))'Co ^(-1) (Yo-Y(p))+(p-pb) ' Pb^(-1) (p-pb)]
where Yo represents the data vector, Y(p) is the model output vector, Co and Pb the error covariance matrix on data and model parameters, respectively. p is the parameter vector. pb is a priori values of p.
pb = [0.06, 1600, 0.05, 1];
My question is how to estimate the parameter (a, b, c, d) using WinBUGS?