rrcov3way (version 0.5-0)

plot.tucker3: Plot a parafac or a tucker3 object

Description

Different plots for the results of Parafac or Tucker3 analysis, stored in a Parafac or a tucker3 object, see Details.

Usage

# S3 method for tucker3
plot(x, which = c("dd", "comp", "allcomp", "jbplot", 
    "tjplot", "all"), ask = (which == "all" && dev.interactive(TRUE)), id.n, ...)
    # S3 method for parafac
plot(x, which = c("dd", "comp", "percomp", "allcomp", 
    "all"), ask = (which == "all" && dev.interactive(TRUE)), id.n, ...)

Arguments

x

A tucker3 object

which

Which plot to select (see Details). Default is dd, distance-distance plot.

ask

Generates all plots in interactive mode

id.n

Number of items to highlight

...

Other parameters to be passed to the lower level functions.

Author

Valentin Todorov and Maria Anna Di Palma and Michele Gallo

Details

Different plots for a tucker3 or parafac object will be produced. Use the parameter which to select which plot to produce:

dd

Distance-distance plot

comp

Paired components plot

percomp

Per-component plot - only for Parafac

allcomp

All components plot

jbplot

Joint biplot - only for Tucker3

tjplot

Trajectory plot - only for Tucker3

References

Kiers, H.A. (2000).Some procedures for displaying results from three-way methods. Journal of Chemometrics. 14(3): 151-170.

Kroonenberg, P.M. (1983).Three-mode principal component analysis: Theory and applications (Vol. 2), DSWO press.

Examples

#############
##
## Example with the UNIDO Manufacturing value added data

data(va3way)
dim(va3way)

## Treat quickly and dirty the zeros in the data set (if any)
va3way[va3way==0] <- 0.001

##
res <- Tucker3(va3way)
res
print(res$fit)
print(res$A)

## Print the core matrix
print(res$GA)

## Distance-distance plot
plot(res, which="dd", main="Distance-distance plot")

## Paired component plot, mode A
plot(res, which="comp", main="Paired component plot (mode A)")

## Paired component plot, mode B
plot(res, which="comp", mode="B", main="Paired component plot (mode B)")

## Joint biplot
plot(res, which="jbplot", main="Joint biplot")

## Trajectory
plot(res, which="tjplot", choices=c(1:4), arrows=FALSE, main="Trajectory biplot")