https://github.com/cran/CollocInfer
Tip revision: 41bf713b4d66e5c1d58114455c66a35ff8f8f7c3 authored by Giles Hooker on 29 July 2011, 00:00:00 UTC
version 0.1.7
version 0.1.7
Tip revision: 41bf713
trymcmc2.r
source('Devel/Colloc.MCMC.R')
walk.var = diag(c(0.001,0.001,0.01))
cscale = 0.000001
prior = function(pars){ return( sum(dnorm(pars/c(1,1,10),log=TRUE))) }
nstep = 100
hmm = Colloc.MCMC(times,data,pars+c(-0.1,0.1,-1),coefs,lik,proc,prior,walk.var,cscale,nstep,in.meth='SplineEst',control.in=control.in)