# defaut linear regression model Y=aX+b
mod <- model()
# BaRatin model for a single-control rating curve Y=a(X-b)^c
mod <- model(ID='BaRatin',nX=1,nY=1,
par=list(parameter('a',10,prior.dist='LogNormal',prior.par=c(log(10),0.1)),
parameter('b',-1,prior.dist='Gaussian',prior.par=c(-1,1)),
parameter('c',5/3,prior.dist='Gaussian',prior.par=c(5/3,0.05))),
xtra=xtraModelInfo(object=matrix(1,nrow=1,ncol=1)))
Run the code above in your browser using DataLab