plot-methods function

Plot a parafac or a tucker3 object

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

## S3 method for class 'tucker3'
plot(x, which = c("dd", "comp", "allcomp", "jbplot", 
    "tjplot", "all"), ask = (which == "all" && dev.interactive(TRUE)), id.n, ...)
    ## S3 method for class '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.

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.

Author(s)

Valentin Todorov valentin.todorov@chello.at and Maria Anna Di Palma madipalma@unior.it and Michele Gallo mgallo@unior.it

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")