Arguments
- object
- verbose
logical; If
TRUE
, include more details of the component definitions. IfFALSE
, only show basic component definition information. Default:FALSE
- ...
arguments passed on to component summary functions, see
summary.component()
.- x
An object to be printed
Examples
# \donttest{
if (bru_safe_inla()) {
# Simulate some covariates x and observations y
input.df <- data.frame(x = cos(1:10))
input.df <- within(input.df, {
y <- 5 + 2 * x + rnorm(10, mean = 0, sd = 0.1)
})
# Fit a Gaussian likelihood model
fit <- bru(y ~ x + Intercept(1), family = "gaussian", data = input.df)
# Obtain summary
fit$summary.fixed
}
#> mean sd 0.025quant 0.5quant 0.975quant mode
#> x 1.945352 0.03544012 1.874535 1.945352 2.016167 1.945352
#> Intercept 5.006657 0.02505637 4.956588 5.006658 5.056724 5.006658
#> kld
#> x 5.807738e-06
#> Intercept 5.807833e-06
if (bru_safe_inla()) {
# Alternatively, we can use the bru_obs() function to construct the likelihood:
lik <- bru_obs(family = "gaussian",
formula = y ~ x + Intercept,
data = input.df)
fit <- bru(~ x + Intercept(1), lik)
fit$summary.fixed
}
#> mean sd 0.025quant 0.5quant 0.975quant mode
#> x 1.945352 0.03544012 1.874535 1.945352 2.016167 1.945352
#> Intercept 5.006657 0.02505637 4.956588 5.006658 5.056724 5.006658
#> kld
#> x 5.807738e-06
#> Intercept 5.807833e-06
# An important addition to the INLA methodology is bru's ability to use
# non-linear predictors. Such a predictor can be formulated via bru_obs()'s
# \code{formula} parameter. The z(1) notation is needed to ensure that
# the z component should be interpreted as single latent variable and not
# a covariate:
if (bru_safe_inla()) {
z <- 2
input.df <- within(input.df, {
y <- 5 + exp(z) * x + rnorm(10, mean = 0, sd = 0.1)
})
lik <- bru_obs(
family = "gaussian", data = input.df,
formula = y ~ exp(z) * x + Intercept
)
fit <- bru(~ z(1) + Intercept(1), lik)
# Check the result (z posterior should be around 2)
fit$summary.fixed
}
#> mean sd 0.025quant 0.5quant 0.975quant mode
#> z 2.007911 0.005651808 1.996618 2.007911 2.019205 2.007911
#> Intercept 5.042930 0.029760284 4.983461 5.042930 5.102396 5.042930
#> kld
#> z 5.735138e-06
#> Intercept 5.735015e-06
# }