https://github.com/cran/CollocInfer
Tip revision: 0b3ff16925952344f8d6573cc0a2afafa3a28b59 authored by Giles Hooker on 08 May 2014, 00:00:00 UTC
version 1.0.1
version 1.0.1
Tip revision: 0b3ff16
Colloc.MCMC.rd
\name{Colloc.MCMC}
\alias{Colloc.MCMC}
\title{Bayesian MCMC for Collocation Inference}
\description{ An MCMC algorithm for Collocation Inference }
\usage{
Colloc.MCMC(times,data,pars,coefs,lik,proc,prior,walk.var,
nstep,in.meth='house',control.in=NULL)
}
\arguments{
\item{times}{ Vector observation times for the data.}
\item{data}{ Matrix of observed data values. }
\item{pars}{ Initial values of parameters to be estimated processes. }
\item{coefs}{ Vector giving the current estimate of the coefficients in the spline. }
\item{lik} { \code{lik} object defining the observation process. }
\item{proc}{ \code{proc} object defining the state process. }
\item{niter}{ Number of MCMC steps to take. }
\item{walk.var}{ Random walk variance for parameters. }
\item{in.meth}{ Inner optimization function to be used, currently one of 'nlminb', 'MaxNR', 'optim' or 'house'.
The last calls \code{SplineEst.NewtRaph}. This is fast but has poor convergence. }
\item{control.in}{ Control object for inner optimization functions. }
}
\value{
\item{parhist}{The Markov chain history of the parameters}
\item{coefhist}{The Markov chain history of the coefficients. }
}
}
\details{
The MCMC algorithm proceeds by taking a mean-zero Gaussian random walk step in the
parameter space, minimizes the inner optimization in coefficients and then takes a random
walk in coefficient space with variance given by the second derivative matrix of the
inner optimization.
}
\seealso{sse.setup, multinorm.setup, Profile.sse, Profile.multinorm, inneropt}
\examples{}