interpret.fpcad: Scores of the 'fpcad' function vs. moments of the densities

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/interpret.fpcad.R

Description

Applies to an object of class fpcad, plots the principal scores vs. the moments of the densities (means, standard deviations, variances, correlations, skewness and kurtosis coefficients), and computes the correlations between these scores and moments.

Usage

1
2
## S3 method for class 'fpcad'
interpret(x, nscore = 1:3, moment = "mean")

Arguments

x

object of class fpcad (returned by the fpcad function).

nscore

numeric vector. Selects the columns of the data frame x$scores to be interpreted.

Warning: Its components cannot be greater than the nb.factors argument in the call of the fpcad function.

moment

characters string. Selects the moments to cross with scores:

  • "mean" (means)

  • "sd" (standard deviations)

  • "cov" (covariances)

  • "cor" (correlation coefficients)

  • "skewness" (skewness coefficients)

  • "kurtosis" (kurtosis coefficients)

  • "all" (for univariate densities only. It simultaneously considers means, standard deviations, variances and skewness and kurtosis coefficients)

Details

A new graphics device is opened for each score column. A device can contain up to 9 graphs. If there are too many (more than 36) graphs for each score, a multipage PDF file is created in the current working directory, and the graphs are displayed in it.

The number of principal scores to be interpreted cannot be greater than nb.factors of the data frame x$scores returned by the function fpcad.

Value

Returns a list including:

pearson

matrix of Pearson correlations between selected scores and moments.

spearman

matrix of Spearman correlations between selected scores and moments.

Author(s)

Rachid Boumaza, Pierre Santagostini, Smail Yousfi, Gilles Hunault, Sabine Demotes-Mainard

References

Boumaza, R., Yousfi, S., Demotes-Mainard, S. (2015). Interpreting the principal component analysis of multivariate density functions. Communications in Statistics - Theory and Methods, 44 (16), 3321-3339.

See Also

fpcad; plot.fpcad.

Examples

1
2
3
4
5
6
7
data(roses)
rosefold <- as.folder(roses[,c("Sha","Den","Sym","rose")])
result1 <- fpcad(rosefold)
interpret(result1)
interpret(result1, moment = "var")
interpret(result1, moment = "cor")
interpret(result1, nscore = 1:2)

dad documentation built on Sept. 2, 2017, 1:04 a.m.