multiNet: A function that generates samples for a multivariate fixed...

In netcmc: Spatio-Network Generalised Linear Mixed Models for Areal Unit and Network Data

A function that generates samples for a multivariate fixed effects and network model.

Description

This function that generates samples for a multivariate fixed effects and network model, which is given by

Y_{i_sr}|μ_{i_sr} \sim f(y_{i_sr}| μ_{i_sr}, σ_{er}^{2}) ~~~ i=1,…, N_{s},~s=1,…,S ,~r=1,…,R,

g(μ_{i_sr}) = \boldsymbol{x}^\top_{i_s} \boldsymbol{β}_{r} + ∑_{j\in \textrm{net}(i_s)}w_{i_sj}u_{jr}+ w^{*}_{i_s}u^{*}_{r},

\boldsymbol{β}_{r} \sim \textrm{N}(\boldsymbol{0}, α\boldsymbol{I})

\boldsymbol{u}_{j} = (u_{1j},…, u_{Rj}) \sim \textrm{N}(\boldsymbol{0}, \boldsymbol{Σ}_{\boldsymbol{u}}),

\boldsymbol{u}^{*} = (u_{1}^*,…, u_{R}^*) \sim \textrm{N}(\boldsymbol{0}, \boldsymbol{Σ}_{\boldsymbol{u}}),

\boldsymbol{Σ}_{\boldsymbol{u}} \sim \textrm{Inverse-Wishart}(ξ_{\boldsymbol{u}}, \boldsymbol{Ω}_{\boldsymbol{u}}),

σ_{er}^{2} \sim \textrm{Inverse-Gamma}(α_{3}, ξ_{3}).

The covariates for the ith individual in the sth spatial unit or other grouping are included in a p \times 1 vector \boldsymbol{x}_{i_s}. The corresponding p \times 1 vector of fixed effect parameters relating to the rth response are denoted by \boldsymbol{β}_{r}, which has an assumed multivariate Gaussian prior with mean \boldsymbol{0} and diagonal covariance matrix α\boldsymbol{I} that can be chosen by the user. A conjugate Inverse-Gamma prior is specified for σ_{er}^{2}, and the corresponding hyperparamaterers (α_{3}, ξ_{3}) can be chosen by the user.

The R \times 1 vector of random effects for the jth alter is denoted by \boldsymbol{u}_{j} = (u_{j1}, …, u_{jR})_{R \times 1}, while the R \times 1 vector of isolation effects for all R outcomes is denoted by \boldsymbol{u}^{*} = (u_{1}^*,…, u_{R}^*), and both are assigned multivariate Gaussian prior distributions. The unstructured covariance matrix \boldsymbol{Σ}_{\boldsymbol{u}} captures the covariance between the R outcomes at the network level, and a conjugate Inverse-Wishart prior is specified for this covariance matrix \boldsymbol{Σ}_{\boldsymbol{u}}. The corresponding hyperparamaterers (ξ_{\boldsymbol{u}}, \boldsymbol{Ω}_{\boldsymbol{u}}) can be chosen by the user.

The exact specification of each of the likelihoods (binomial, Gaussian, and Poisson) are given below:

\textrm{Binomial:} ~ Y_{i_sr} \sim \textrm{Binomial}(n_{i_s}, θ_{i_sr}) ~ \textrm{and} ~ g(μ_{i_sr}) = \textrm{ln}(θ_{i_sr} / (1 - θ_{i_sr})),

\textrm{Gaussian:} ~ Y_{i_sr} \sim \textrm{N}(μ_{i_sr}, σ_{er}^{2}) ~ \textrm{and} ~ g(μ_{i_sr}) = μ_{i_sr},

\textrm{Poisson:} ~ Y_{i_sr} \sim \textrm{Poisson}(μ_{i_sr}) ~ \textrm{and} ~ g(μ_{i_sr}) = \textrm{ln}(μ_{i_sr}).

Usage

multiNet(formula, data, trials, family, W, numberOfSamples = 10, burnin = 0, 
thin = 1, seed = 1, trueBeta = NULL, trueURandomEffects = NULL, 
trueVarianceCovarianceU = NULL, trueSigmaSquaredE = NULL, 
covarianceBetaPrior = 10^5, xi, omega, a3 = 0.001, b3 = 0.001, 
centerURandomEffects = TRUE)

Arguments

formula	A formula for the covariate part of the model using a similar syntax to that used in the lm() function.
data	An optional data.frame containing the variables in the formula.
trials	A vector the same length as the response containing the total number of trials n_{i_sr}. Only used if \texttt{family}=“binomial".
family	The data likelihood model that must be “gaussian", “poisson" or “binomial".
W	A matrix \boldsymbol{W} that encodes the social network structure and whose rows sum to 1.
numberOfSamples	The number of samples to generate pre-thin.
burnin	The number of MCMC samples to discard as the burn-in period.
thin	The value by which to thin \texttt{numberOfSamples}.
seed	A seed for the MCMC algorithm.
trueBeta	If available, the true value of \boldsymbol{β}_{1}, …, \boldsymbol{β}_{R}.
trueURandomEffects	If available, the true values of \boldsymbol{u}_{1}, …, \boldsymbol{u}_{J}, \boldsymbol{u}^{*}.
trueVarianceCovarianceU	If available, the true value of \boldsymbol{Σ}_{\boldsymbol{u}}.
trueSigmaSquaredE	If available, the true value of σ_{e1}^{2}, …, σ_{eR}^{2}. Only used if \texttt{family}=“gaussian".
covarianceBetaPrior	A scalar prior α for the covariance parameter of the beta prior, such that the covariance is α\boldsymbol{I}.
xi	The degrees of freedom parameter for the Inverse-Wishart distribution relating to the network random effects ξ_{\boldsymbol{u}}.
omega	The scale parameter for the Inverse-Wishart distribution relating to the network random effects \boldsymbol{Ω}_{\boldsymbol{u}}.
a3	The shape parameter for the Inverse-Gamma distribution relating to the error terms α_{3}. Only used if \texttt{family}=“gaussian".
b3	The scale parameter for the Inverse-Gamma distribution relating to the error terms ξ_{3}. Only used if \texttt{family}=“gaussian".
centerURandomEffects	A choice to center the network random effects after each iteration of the MCMC sampler.

Value

call	The matched call.
y	The response used.
X	The design matrix used.
standardizedX	The standardized design matrix used.
W	The network matrix used.
samples	The matrix of simulated samples from the posterior distribution of each parameter in the model (excluding random effects).
betaSamples	The matrix of simulated samples from the posterior distribution of \boldsymbol{β}_{1}, …, \boldsymbol{β}_{R} parameters in the model.
varianceCovarianceUSamples	The matrix of simulated samples from the posterior distribution of \boldsymbol{Σ}_{\boldsymbol{u}} in the model.
uRandomEffectsSamples	The matrix of simulated samples from the posterior distribution of network random effects \boldsymbol{u}_{1}, …, \boldsymbol{u}_{J}, \boldsymbol{u}^{*} in the model.
sigmaSquaredESamples	The vector of simulated samples from the posterior distribution of σ_{e1}^{2}, …, σ_{eR}^{2} in the model. Only used if \texttt{family}=“gaussian".
acceptanceRates	The acceptance rates of parameters in the model from the MCMC sampling scheme .
uRandomEffectsAcceptanceRate	The acceptance rates of network random effects in the model from the MCMC sampling scheme.
timeTaken	The time taken for the model to run.
burnin	The number of MCMC samples to discard as the burn-in period.
thin	The value by which to thin \texttt{numberOfSamples}.
DBar	DBar for the model.
posteriorDeviance	The posterior deviance for the model.
posteriorLogLikelihood	The posterior log likelihood for the model.
pd	The number of effective parameters in the model.
DIC	The DIC for the model.

Author(s)

George Gerogiannis

netcmc documentation built on Nov. 10, 2022, 5:11 p.m.