ex1.dag.data

10000 observations simulated from a DAG with 10 variables from Poisson, Bernoulli and Gaussian distributions.

datasets

Bayesian network analysis is a form of probabilistic graphical models which derives from empirical data a directed acyclic graph, DAG, describing the dependency structure between random variables.
An additive Bayesian network model consists of a form of a DAG where each node comprises a generalized linear model, GLM. Additive Bayesian network models are equivalent to Bayesian multivariate regression using graphical modelling, they generalises the usual multivariable regression, GLM, to multiple dependent variables.
'abn' provides routines to help determine optimal Bayesian network models for a given data set, where these models are used to identify statistical dependencies in messy, complex data. The additive formulation of these models is equivalent to multivariate generalised linear modelling (including mixed models with iid random effects).
The usual term to describe this model selection process is structure discovery.
The core functionality is concerned with model selection - determining the most robust empirical model of data from interdependent variables. Laplace approximations are used to estimate goodness of fit metrics and model parameters, and wrappers are also included to the INLA package which can be obtained from <https://www.r-inla.org>.
The computing library JAGS <https://mcmc-jags.sourceforge.io> is used to simulate 'abn'-like data.
A comprehensive set of documented case studies, numerical accuracy/quality assurance exercises, and additional documentation are available from the 'abn' website <http://r-bayesian-networks.org>.

Matteo Delucchi

Modelling Multivariate Data with Additive Bayesian Networks

Reinhard Furrer

Gilles Kratzer

Fraser Iain Lewis

Marta Pittavino

Kalina Cherneva

ex1.dag.data function

A data frame, binary variables are factors.
The relevant formulas are given below (note these do not give parameter
estimates just the form of the relationships, like in glm(),
e.g. logit()=1+p1 means a logit link function and comprises of an
intercept term and a term involving p1).<ul>
<li>b1binary, logit()=1</li>
<li>p1poisson, log()=1</li>
<li>g1gaussian, identity()=1</li>
<li>b2binary, logit()=1</li>
<li>p2poisson, log()=1+b1+p1</li>
<li>b3binary, logit()=1+b1+g1+b2</li>
<li>g2gaussian, identify()=1+p1+g1+b2</li>
<li>b4binary, logit()=1+g1+p2</li>
<li>b5binary, logit()=1+g1+g2</li>
<li>g3gaussian, identity()=1+g1+b2</li>
</ul>

Format

Synthetic validation data set for use with abn library examples — ex1.dag.data

Synthetic validation data set for use with abn library examples

Last chance! 50% off unlimited learning

ex1.dag.data: Synthetic validation data set for use with abn library examples

Description

Usage

Arguments

Format

Examples