Description Usage Arguments Details Value Author(s) See Also

View source: R/predict.laplace.R

This may be used to predict either new, unobserved instances of
*y* (called *y[new]*) or
replicates of *y* (called
*y[rep]*), and then perform posterior
predictive checks. Either *y[new]* or
*y[rep]* is predicted given an object of
class `laplace`

, the model specification, and data. This function
requires that posterior samples were produced with
`LaplaceApproximation`

.

1 2 |

`object` |
An object of class |

`Model` |
The model specification function is required. |

`Data` |
A data set in a list is required. The dependent
variable is required to be named either |

`CPUs` |
This argument accepts an integer that specifies the number
of central processing units (CPUs) of the multicore computer or
computer cluster. This argument defaults to |

`Type` |
This argument specifies the type of parallel processing to
perform, accepting either |

`...` |
Additional arguments are unused. |

Since Laplace Approximation characterizes marginal posterior
distributions with modes and variances, and posterior predictive
checks involve samples, the `predict.laplace`

function requires
the use of independent samples of the marginal posterior
distributions, provided by `LaplaceApproximation`

when
`sir=TRUE`

.

The samples of the marginal posterior distributions of the target
distributions (the parameters) are passed along with the data to the
`Model`

specification and used to draw samples from the deviance
and monitored variables. At the same time, the fourth component in the
returned list, which is labeled `yhat`

, is a vector of
expectations of *y*, given the samples, model
specification, and data. To predict *y[rep]*,
simply supply the data set used to estimate the model. To predict
*y[new]*, supply a new data set instead (though
for some model specifications, this cannot be done, and
*y[new]* must be specified in the `Model`

function). If the new data set does not have *y*, then
create `y`

in the list and set it equal to something sensible,
such as `mean(y)`

from the original data set.

The variable `y`

must be a vector. If instead it is matrix
`Y`

, then it will be converted to vector `y`

. The vectorized
length of `y`

or `Y`

must be equal to the vectorized length
of `yhat`

, the fourth component of the returned list of the
`Model`

function.

Parallel processing may be performed when the user specifies
`CPUs`

to be greater than one, implying that the specified number
of CPUs exists and is available. Parallelization may be performed on a
multicore computer or a computer cluster. Either a Simple Network of
Workstations (SNOW) or Message Passing Interface is used (MPI). With
small data sets and few samples, parallel processing may be slower,
due to computer network communication. With larger data sets and more
samples, the user should experience a faster run-time.

For more information on posterior predictive checks, see https://web.archive.org/web/20150215050702/http://www.bayesian-inference.com/posteriorpredictivechecks.

This function returns an object of class `laplace.ppc`

(where
“ppc” stands for posterior predictive checks). The returned object
is a list with the following components:

`y` |
This stores |

`yhat` |
This is a |

`Deviance` |
This is a vector of length |

`monitor` |
This is a |

Statisticat, LLC.

`LaplaceApproximation`

and
`SIR`

.

LaplacesDemon documentation built on May 19, 2017, 10:22 p.m.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.