flip: flip

In cbmodels: Collection of "combined" models: PCA-LDA, PLS-LDA, etc.

Description

Flipping of models

lda models are invariant to flipping, i.e. the predicted posterior probabilities and classes do not change if the LD scores lda$x are multiplied by -1.

PCA (prcomp) models are invariant to flipping of both a loading and its score.

Usage

  flip(object, dims, ...)

  ## S3 method for class 'lda'
 flip(object, dims = FALSE, ...)

  ## S3 method for class 'pcalda'
 flip(object, dims = FALSE, ...)

  ## S3 method for class 'plslda'
 flip(object, dims = FALSE, ...)

  ## S3 method for class 'prcomp'
 flip(object, dims = FALSE, ...)

Arguments

object	a model to flip
dims	dimensions which should be mirrored
...	further parameters

Details

flip is a post-processing method for models. It flips (mirrors) dimensions of the model space.

Value

the flipped/mirrored model

Author(s)

Claudia Beleites

See Also

lda

pcalda

plslda

prcomp

Examples

## flip 1st dimension
model <- lda (Species ~ ., data = iris)
pred <- predict (model)
plot (pred$x, col = iris$Species, pch = 19, cex = 0.3, asp = 1)

model.flip <- flip (model, 1)
pred.flip <- predict (model.flip)
points (pred.flip$x, col = iris$Species, cex = 0.5)

## check difference between original model's predictions and rotated model's predictions
diff <- pred$posterior - pred.flip$posterior
summary (diff)
boxplot (diff)
## mirror x axis
model <- pcalda (X = iris [,1:4], grouping = iris$Species, comps=1:3)
pred <- predict (model)
plot (pred$x, col = iris$Species, pch = 19, cex = 0.3, asp = 1)

model.flip <- flip (model, 1)
pred.flip <- predict (model.flip)
points (pred.flip$x, col = iris$Species, cex = 0.5)

## check difference between original model's predictions and flipped model's predictions
diff <- pred$posterior - pred.flip$posterior
summary (diff)
boxplot (diff)
## mirror x axis of plslda scores
model <- plslda (X = iris [,1:4], grouping = iris[,5], ncomp=2)
pred <- predict (model)
plot (pred$x, col = iris$Species, pch = 19, cex = 0.3, asp = 1)

model.flip <- flip (model, 1)
pred.flip <- predict (model.flip)
points (pred.flip$x, col = iris$Species, cex = 0.5)

## check difference between original model's predictions and flipped model's predictions
diff <- pred$posterior - pred.flip$posterior
summary (diff)
boxplot (diff)
## flip 1st dimension (loading & score)
model <- prcomp (iris [-5])
model.flip <- flip (model, 1)

pairs (model$x, col = iris$Species, pch = 19, cex = 0.3, asp = 1)
pairs (model.flip$x, col = iris$Species, pch = 19, cex = 0.3, asp = 1)

matplot (model$rotation, type = "l")
abline (h = 0)
matpoints (model.flip$rotation, type = "b")

## check difference between original model's predictions and rotated model's predictions
diff <- tcrossprod (model$x, model$rotation) - tcrossprod (model.flip$x, model.flip$rotation)
summary (diff)
boxplot (diff)

cbmodels documentation built on May 31, 2017, 2:11 a.m.