Description Usage Arguments Details Value Author(s) References Examples

View source: R/targetselection.R

Implements the rule of thumb proposed by Touloumis (2015) for target matrix selection. If the estimated optimal shrinkage intensities of the three target matrices are of similar magnitude, then the average and the range of the sample variances should be inspected in order to adopt the most plausible target matrix.

1 | ```
targetselection(data, centered = FALSE)
``` |

`data` |
a numeric matrix containing the data. |

`centered` |
a logical indicating if the mean vector is the zero vector. |

The rows of the data matrix `data`

correspond to variables and the
columns to subjects.

Prints the estimated optimal shrinkage intensities, the range and average of the sample variances and returns an object of the class 'targetsel' that has components:

`optimal_sphericity` |
The estimated optimal intensity for a target matrix with equal variances. |

`optimal_identity` |
The estimated optimal shrinkage intensity for the identity target matrix. |

`optimal_diagonal` |
The estimated optimal intensity for a target matrix with unequal variances. |

`range` |
The range of the sample variances. |

`average` |
The average of the sample variances. |

Anestis Touloumis

Touloumis, A. (2015) Nonparametric Stein-type Shrinkage
Covariance Matrix Estimators in High-Dimensional Settings.
*Computational Statistics & Data Analysis* **83**, 251–261.

1 2 3 4 5 6 7 8 9 10 11 12 | ```
data(colon)
normal_group <- colon[, 1:40]
targetselection(normal_group)
## Similar intensities, the range of the sample variances is small and the
## average is not close to one. The scaled identity matrix seems to be the
## most suitable target matrix for the normal group.
tumor_group <- colon[, 41:62]
targetselection(tumor_group)
## Similar intensities, the range of the sample variances is small and the
## average is not close to one. The scaled identity matrix seems to be the
## most suitable target matrix for the colon group.
``` |

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.