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

This function plots Cattell's scree-test or the BIC selection, using parameters coming from `hdda` or `hddc` functions.

1 2 |

`x` |
A 'hdc' class object obtained using |

`method` |
The method used to select the intrinsic dimension. It can be "BIC" or "Cattell". By default it takes the method used when obtaining the parameters using |

`threshold` |
The threshold used in Cattell's Scree-Test. By default it takes the threshold in the argument x, if none, the default value of the threshold is 0.2. |

`noise.ctrl` |
This parameter avoids to have a too low value of the 'noise' parameter b. It garantees that the dimension selection process do not select too many dimensions (which leads to a potential too low value of the noise parameter b). When selecting the intrinsic dimensions using Cattell's scree-test or BIC, the function doesn't use the eigenvalues inferior to noise.ctrl, so that the intrinsic dimensions selected can't be higher or equal to the order of these eigenvalues. |

`...` |
Arguments based from or to other methods. |

- If method = "Cattell"
The plot of the eigen values and of the sequential differences of the eigen values. The dimension to retain is the one before the last fall of the eigenvalues' differences below the threshold.

- If method = "BIC"
The BIC related to the dimension for each class. It stops after the first fall of the BIC.

Laurent Berge, Charles Bouveyron and Stephane Girard

Bouveyron, C. Girard, S. and Schmid, C. (2007) “High Dimensional Discriminant Analysis”, *Communications in Statistics: Theory and Methods*, vol. **36** (14), pp. 2607–2623

Bouveyron, C. Girard, S. and Schmid, C. (2007) “High-Dimensional Data Clustering”, *Computational Statistics and Data Analysis*, vol. **52** (1), pp. 502–519

Berge, L. Bouveyron, C. and Girard, S. (2012) “HDclassif: An R Package for Model-Based Clustering and Discriminant Analysis of High-Dimensional Data”, *Journal of Statistical Software*, **46**(6), 1–29, url: http://www.jstatsoft.org/v46/i06/

`hdda`

, `hddc`

, `predict.hdc`

.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ```
# Example 1 :
data(wine)
a <- wine[,-1]
z <- wine[,1]
prms1 <- hdda(a, z, model="AkBkQkDk", scaling=TRUE, d_select="bic")
#the plot related to the selection that has been done: BIC
plot(prms1)
#it shows the plot of Cattell's scree-test, with a threshold of .3
plot(prms1,"Cattell",0.3)
prms2 <- hdda(a, z, model="AkBkQkD", scaling=TRUE, d_select="cattell")
#the plot related to the selection that has been done: Cattell's scree-test
plot(prms2)
#the plot of the BIC
plot(prms2,"b")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.