majorityModel: Combine Model-based Recursive Partitioning with a Majority...

In schiffner/locClass: Collection of Local Classification Methods

Description

Combine Model-Based Recursive Partitioning with a Majority Classifier.

Usage

majorityModel

\method{reweight}{majorityModel}(object, weights, ...)

\method{deviance}{majority}(object, ...)

\method{estfun}{majority}(x, ...)

predict.majorityModel(object, out = c("class", "posterior"), ...)

Arguments

object	An object of class "majorityModel" and "majority", respectively.
weights	A vector of observation weights.
x	An object of class "majority".
out	Should class labels or posterior probabilities be returned?
...	Further arguments.

Format

An object of class StatModel of length 1.

Details

This page lists all ingredients to combine a Majority Classifier with Model-Based Recursive Partitioning (mob from package party). See the example for how to do that.

majorityModel is an object of class StatModel-class implemented in package modeltools that provides an infra-structure for an unfitted majority model.

Moreover, methods for majority and majorityModel objects for the generic functions reweight, deviance, estfun, and predict are provided.

Value

reweight: The re-weighted fitted "majority Model" object.
deviance: The value of the deviance for the Majority Classifier extracted from object, i.e. the log-likelihood.
estfun: The empirical estimating (or score) function for the Majority Classifier, i.e. the derivatives of the log-likelihood with respect to the parameters, evaluated at the training data.
predict: Either a vector of predicted class labels or a matrix of class posterior probabilities.

References

Zeileis, A., Hothorn, T. and Kornik, K. (2008), Model-based recursive partitioning. Journal of Computational and Graphical Statistics, 17(2) 492–514.

See Also

reweight, deviance, estfun, predict.

Examples

library(benchData)

data <- vData(500)
x <- seq(0,1,0.05)
grid <- expand.grid(x.1 = x, x.2 = x)

fit <- mob(y ~ x.1 + x.2 | x.1 + x.2, data = data, model = majorityModel,
control = mob_control(objfun = deviance, minsplit = 20))

## predict posterior probabilities
pred <- predict(fit, newdata = grid, out = "posterior")
post <- matrix(0, length(pred), 2)
colnames(post) = 1:2
for (i in seq_along(pred))
    post[i, colnames(pred[[i]])] = pred[[i]]

image(x, x, matrix(as.numeric(post[,1]), length(x)), xlab = "x.1", ylab = "x.2")
contour(x, x, matrix(as.numeric(post[,1]), length(x)), levels = 0.5, add = TRUE)
points(data$x, pch = as.character(data$y))

## predict node membership
splits <- predict(fit, newdata = grid, type = "node")
contour(x, x, matrix(splits, length(x)), levels = min(splits):max(splits), add = TRUE, lty = 2)

## training error
mean(predict(fit) != as.numeric(data$y))

schiffner/locClass documentation built on May 29, 2019, 3:39 p.m.