myplotDD: Distance-Distance Plot

In robustfa: Object Oriented Solution for Robust Factor Analysis

Distance-Distance Plot

Description

"myplotDD" is a revised version of ".myddplot" in "plot-utils.R" in the package "rrcov". In "myplotDD", id.n and ind are printed out.

Usage

myplotDD(x, cutoff, id.n)

Arguments

x	An S4 object of class "CovRobust".
cutoff	The cutoff value used. If missing, cutoff <- sqrt(qchisq(0.975, p)) by default.
id.n	Number of observations to identify by a label. If not supplied, the number of observations with robust distance larger than cutoff is used.

Details

Distance-Distance Plot: Plot the vector y=rd (robust distances) against x=md (mahalanobis distances). Identify by a label the id.n observations with largest rd. If id.n is not supplied, calculate it as the number of observations larger than cutoff. Use cutoff to draw a horisontal and a vertical line. Draw also a dotted line with a slope 1.

"myplotDD(x)" is equivalent to "plot(x, which="dd")". which: indicate what kind of plot. If which = "dd", then a distance-distance Plot.

Value

A distance-distance plot is shown. Return a list with components:

cutoff	The cutoff value used. If missing, cutoff <- sqrt(qchisq(0.975, p)) by default.
id.n	Number of observations to identify by a label. If not supplied, the number of observations with robust distance larger than cutoff is used.
sort.y	A list containing the sorted values of y (the robust distance)
ind	The indices of the largest id.n observations whose robust distances are larger than cutoff.

Author(s)

Ying-Ying Zhang (Robert) [e-mail]

References

Zhang, Y. Y. (2013), An Object Oriented Solution for Robust Factor Analysis.

See Also

plot, qchisq, CovClassic, getDistance

Examples

data(stock611)
covMcd=CovRobust(x=scale(stock611[,3:12]), control="mcd"); covMcd

## "myplotDD" shows id.n and ind.
## Note: id.n and ind change each time due to covMcd changes each time!
## However, the ind of largest robust distances do not change.
result = myplotDD(x=covMcd); result

## "myplotDD" is equivalent to "plot(x=covMcd, which="dd")".
plot(x=covMcd, which="dd")

robustfa documentation built on April 16, 2023, 5:18 p.m.