partykit (version 1.3-0)

glmtree: Generalized Linear Model Trees

Description

Model-based recursive partitioning based on generalized linear models.

Usage

glmtree(formula, data, subset, na.action, weights, offset, cluster,
  family = gaussian, epsilon = 1e-8, maxit = 25, method = "glm.fit", ...)

Value

An object of class glmtree inheriting from modelparty. The info element of the overall party and the individual

nodes contain various informations about the models.

The methodology is described in partykit::Zeileis+Hothorn+Hornik:2008.

Arguments

formula

symbolic description of the model (of type y ~ z1 + ... + zl or y ~ x1 + ... + xk | z1 + ... + zl; for details see below).

data, subset, na.action

arguments controlling formula processing via model.frame.

weights

optional numeric vector of weights. By default these are treated as case weights but the default can be changed in mob_control.

offset

optional numeric vector with an a priori known component to be included in the model y ~ x1 + ... + xk (i.e., only when x variables are specified).

cluster

optional vector (typically numeric or factor) with a cluster ID to be employed for clustered covariances in the parameter stability tests.

family, method

specification of a family and fitting method for glm.

epsilon, maxit

control parameters passed to glm.control.

...

optional control parameters passed to mob_control.

Details

Convenience interface for fitting MOBs (model-based recursive partitions) via the mob function. glmtree internally sets up a model fit function for mob, using glm.fit. Then mob is called using the negative log-likelihood as the objective function.

Compared to calling mob by hand, the implementation tries to avoid unnecessary computations while growing the tree. Also, it provides a more elaborate plotting function.

References

*

See Also

mob, mob_control, lmtree

Examples

if(require("mlbench") && require("vcd")) {

## synthetic diabetes data
data("SynthDiabetes", package = "mlbench")

## recursive partitioning of a logistic regression model
sd_tree2 <- glmtree(diabetes ~ glucose | pressure + triceps +
  insulin + mass + age, data = SynthDiabetes, family = binomial)

## printing whole tree or individual nodes
print(sd_tree2)
print(sd_tree2, node = 1)

## visualization
plot(sd_tree2)
plot(sd_tree2, tp_args = list(cdplot = TRUE))
plot(sd_tree2, terminal_panel = NULL)

## estimated parameters
coef(sd_tree2)
coef(sd_tree2, node = 5)
summary(sd_tree2, node = 5)

## deviance, log-likelihood and information criteria
deviance(sd_tree2)
logLik(sd_tree2)
AIC(sd_tree2)
BIC(sd_tree2)

## different types of predictions
sd <- head(SynthDiabetes)
predict(sd_tree2, newdata = sd, type = "node")
predict(sd_tree2, newdata = sd, type = "response")
predict(sd_tree2, newdata = sd, type = "link")

}