dad (version 4.1.6)

interpret.dstatis: Scores of the dstatis function vs. moments of the densities

Description

Applies to an object of class "dstatis", 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

# S3 method for dstatis
interpret(x, nscore = 1, moment=c("mean", "sd", "var", "cov", "cor",
    "skewness", "kurtosis"), ...)

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.

Arguments

x

object of class "dstatis" (returned by the dstatis.inter function).

nscore

numeric. Selects the column of the data frame x$scores consisting of a score vector.

Note that since dad-4, nscore can only be a single value (in earlier versions, it could be a vector of length > 1).

Warning: nscore cannot be greater than the nb.factors argument in the call of the dstatis.inter 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)

...

Arguments to be passed to methods.

Author

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

Details

A graphics device can contain up to 9 graphs. If there are too many (more than 36) graphs for each score, one can display the graphs in a multipage PDF file.

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

References

Lavit, C., Escoufier, Y., Sabatier, R., Traissac, P. (1994). The ACT (STATIS method). Computational Statistics & Data Analysis, 18 (1994), 97-119.

See Also

dstatis.inter; plot.dstatis.

Examples

data(roses)
rosesf <- as.folder(roses[,c("Sha","Den","Sym","rose")])

# Dual STATIS on the covariance matrices
if (FALSE) {
result <- dstatis.inter(rosesf, group.name = "rose")
interpret(result)
interpret(result, moment = "var")
interpret(result, moment = "cor")
interpret(result, nscore = 2)
}