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
trymcmc.r
source('Devel/Colloc.MCMC.R')
walk.var = diag(c(0.005,0.005,0.05))
cscale = .3
prior = function(pars){ return( sum(dnorm(pars/c(1,1,10),log=TRUE))) }
nstep = 10000
hmm = Colloc.MCMC(FhNtimes,FhNdata,FhNpars,Ires2$coefs,srklik,srkproc,prior,walk.var,cscale,nstep,in.meth='SplineEst',control.in=control.in)
hmm = Colloc.MCMC(seq(0,20,len=401),data,FhNpars,Ires2$coefs,srklik,srkproc,prior,walk.var,cscale,nstep,in.meth='SplineEst',control.in=control.in)