Laplace approximation of posterior to normal

This function generates mode and variance-covariance for a normal proposal distribution for the bayesian logit.

lapl_aprx(y, x, glmobj = FALSE)

Arguments

y	the binary dependent variable y
x	the matrix of independent variables.
glmobj	logical for returning the logit glm object

Value

val A list of mode variance-covariance matrix, and scale factor for proposal draws from the multivariate normal distribution.

Examples

 y = indicat(faithful$waiting,mean(faithful$waiting))
 x = scale(cbind(faithful$eruptions,faithful$eruptions^2))
 gg<- lapl_aprx(y,x)