com.numericalmethod.suanshu.stats.regression.linear.glm
Interface Fitting

All Known Implementing Classes:

IWLS

public interface Fitting

This interface represents a fitting method for estimating β in a Generalized Linear Model (GLM).

John Nelder and Robert Wedderburn proposed an iteratively re-weighted least squares method for maximum likelihood estimation of the model parameters, β. Maximum-likelihood estimation remains popular and is the default method on many statistical computing packages. Other approaches, including Bayesian approaches and least squares fits to variance stabilized responses, have been developed.

See Also:

• Wikipedia: Generalized linear model
• {@link IWLS}

Method Summary

ImmutableVector	betaHat() Get the estimates of β, β^, as in E(Y) = μ = g-1(Xβ)
void	fit(GlmProblem problem, Vector beta0Initial) Fit a Generalized Linear Model.
double	logLikelihood()
ImmutableVector	mu() Get μ as in E(Y) = μ = g-1(Xβ)
ImmutableVector	weights() Get the weights to the observations.

Method Detail

fit

void fit(GlmProblem problem,
         Vector beta0Initial)

Fit a Generalized Linear Model.

This method must be called before the three get methods.

Parameters:

problem - the generalized linear regression problem to be solved

beta0Initial - initial guess for betaHat

mu

ImmutableVector mu()

Get μ as in

 E(Y) = μ = g-1(Xβ)
 

Returns:

μ

betaHat

ImmutableVector betaHat()

Get the estimates of β, β^, as in

 E(Y) = μ = g-1(Xβ)
 

Returns:

β^

weights

ImmutableVector weights()

Get the weights to the observations.

Returns:

the weights

logLikelihood

double logLikelihood()