14.03.2014 Views

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 6 Fitting Generalized Linear Models 181<br />

Examples of Generalized Linear Models<br />

Table 6.1 Examples of Generalized Linear Models<br />

Model Response Variable Distribution Canonical Link<br />

Function<br />

Traditional Linear<br />

Model<br />

continuous Normal identity, g(μ) = μ<br />

Logistic Regression<br />

a count or a binary<br />

r<strong>and</strong>om variable<br />

Binomial<br />

logit,<br />

g( μ)<br />

=<br />

log -----------<br />

μ<br />

<br />

1 – μ<br />

<br />

Poisson Regression in<br />

Log Linear Model<br />

a count Poisson log, g(μ) = log(μ)<br />

Exponential<br />

Regression<br />

positive continuous<br />

Exponential<br />

--<br />

1<br />

μ<br />

JMP fits a generalized linear model to the data by maximum likelihood estimation of the parameter vector.<br />

There is, in general, no closed form solution for the maximum likelihood estimates of the parameters. JMP<br />

estimates the parameters of the model numerically through an iterative fitting process. The dispersion<br />

parameter φ is also estimated by dividing the Pearson goodness-of-fit statistic by its degrees of freedom.<br />

Covariances, st<strong>and</strong>ard errors, <strong>and</strong> confidence limits are computed for the estimated parameters based on the<br />

asymptotic normality of maximum likelihood estimators.<br />

A number of link functions <strong>and</strong> probability distributions are available in JMP. The built-in link functions<br />

are<br />

identity: g(μ) = μ<br />

logit:<br />

g( μ)<br />

=<br />

log -----------<br />

μ<br />

<br />

1 – μ<br />

<br />

probit: g(μ) = Φ -1 (μ), where Φ is the st<strong>and</strong>ard normal cumulative distribution function<br />

log: g(μ) = log(μ)<br />

1<br />

reciprocal: g(μ) = --<br />

μ<br />

<br />

power: g( μ)<br />

μ λ =<br />

if ( λ ≠ 0)<br />

<br />

log( μ)<br />

if λ = 0<br />

complementary log-log: g(m) = log(–log(1 – μ))<br />

When you select the Power link function, a number box appears enabling you to enter the desired power.<br />

The available distributions <strong>and</strong> associated variance functions are<br />

normal: V(μ) = 1<br />

binomial (proportion): V(μ) = μ(1 – μ)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!