ProbitUpdate: EM Algorithm Update for Probit Regression

Description Usage Arguments Value Author(s) See Also Examples

View source: R/ProbitUpdate.R

Description

Function performs an EM update (both the E and M steps) of the parameters for a Probit regression model.

Usage

1
ProbitUpdate(beta.hat, X, y)

Arguments

beta.hat

A vector of length p. The current estimates of the regression parameters.

X

The n x p design matrix for the Probit regression model.

y

Vector of length n containing binary outcomes (either 0 or 1).

Value

A vector of length p - the updated parameter values.

Author(s)

Nicholas Henderson

See Also

ProbitSimulate, ProbitLogLik

Examples

1
2
3
4
5
6
7
8
n <- 200
npars <- 5
true.beta <- .5*rt(npars, df=2) + 2
XX <- matrix(rnorm(n*npars), nrow=n, ncol=npars)
yy <- ProbitSimulate(true.beta, XX)

initial.beta <- rep(0.0, npars)
new.beta <- ProbitUpdate(initial.beta, XX, yy)

daarem documentation built on Feb. 10, 2019, 1:04 a.m.