https://github.com/cran/CollocInfer
Tip revision: f072e42c0357a24ae8f48e3596c5f02c0ca4dd2d authored by Giles Hooker on 25 April 2010, 00:00:00 UTC
version 0.1.1
version 0.1.1
Tip revision: f072e42
trymcmc.r
source('Devel/Colloc.MCMC.R')
walk.var = diag(c(0.005,0.005,0.05))
cscale = 1
prior = function(pars){ return( sum(dnorm(pars/c(1,1,10),log=TRUE))) }
nstep = 50
hmm = Colloc.MCMC(FhNtimes,FhNdata,FhNpars,Ires2$coefs,srklik,srkproc,prior,walk.var,cscale,nstep,in.meth='SplineEst',control.in=control.in)