interpret.dstatis function

Scores of the dstatis function vs. moments of the densities

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.

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

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.

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 .

Returns

Returns a list including: - pearson: matrix of Pearson correlations between selected scores and moments.

• spearman: matrix of Spearman correlations between selected scores and moments.

References

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

Author(s)

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

See Also

dstatis.inter ; plot.dstatis .

Examples

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

# Dual STATIS on the covariance matrices
## Not run:

result <- dstatis.inter(rosesf, group.name = "rose")
interpret(result)
interpret(result, moment = "var")
interpret(result, moment = "cor")
interpret(result, nscore = 2)
## End(Not run)