Description Usage Arguments Value Examples

For models fit using MCMC only, the `log_lik`

method returns the
*S* by *N* pointwise log-likelihood matrix, where *S* is the size
of the posterior sample and *N* is the number of data points, or in the
case of the `stanmvreg`

method (when called on `stan_jm`

model objects) an *S* by *Npat* matrix where *Npat* is the number
of individuals.

1 2 3 4 5 6 7 8 9 |

`object` |
A fitted model object returned by one of the
rstanarm modeling functions. See |

`newdata` |
An optional data frame of new data (e.g. holdout data) to use
when evaluating the log-likelihood. See the description of |

`offset` |
A vector of offsets. Only required if |

`...` |
Currently ignored. |

`m` |
Integer specifying the number or name of the submodel |

`newdataLong, newdataEvent` |
Optional data frames containing new data
(e.g. holdout data) to use when evaluating the log-likelihood for a
model estimated using |

For the `stanreg`

and `stanmvreg`

methods an *S* by
*N* matrix, where *S* is the size of the posterior sample and
*N* is the number of data points. For the `stanjm`

method
an *S* by *Npat* matrix where *Npat* is the number of individuals.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | ```
roaches$roach100 <- roaches$roach1 / 100
fit <- stan_glm(
y ~ roach100 + treatment + senior,
offset = log(exposure2),
data = roaches,
family = poisson(link = "log"),
prior = normal(0, 2.5),
prior_intercept = normal(0, 10),
iter = 500 # to speed up example
)
ll <- log_lik(fit)
dim(ll)
all.equal(ncol(ll), nobs(fit))
# using newdata argument
nd <- roaches[1:2, ]
nd$treatment[1:2] <- c(0, 1)
ll2 <- log_lik(fit, newdata = nd, offset = c(0, 0))
head(ll2)
dim(ll2)
all.equal(ncol(ll2), nrow(nd))
``` |

rstanarm documentation built on Oct. 4, 2019, 1:04 a.m.

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.