plot.tucker3 {rrcov3way}

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

Author(s)

Valentin Todorov [e-mail] and Maria Anna Di Palma [e-mail] and Michele Gallo [e-mail]

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