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

View source: R/RidgeParamEst.R

Maximum likelihood estimation of the ridge parameter by cross-validation

1 2 3 |

`dat` |
the data matrix. |

`X` |
the design matrix. |

`weights` |
weights on the cases of the design matrix. |

`refs` |
a vector specifying validation group membership. Default is to
construct |

`tol` |
the sensitivity in calculations near zero. |

`only.ridge` |
logical, whether only the ridge Parameters should be passed back or additionally the Cross Validation penalised likelihood. |

`doPlot` |
logical, whether a plot of -2logL vs a candidate for the ridge parameter should be drawn. |

`col` |
color of Plot symbols. |

`type` |
type of Plot symbols. |

`...` |
further plot arguments. |

This function estimates the ridge parameter when applying ridge regularization to a sample correlation matrix of residuals. The ridge parameter is estimated to maximize the normal likelihood as estimated via cross validation (Warton 2008).

A list with the following component:

`ridgeParameter` |
the estimated ridge parameter |

If `only.ridge=FALSE`

the returned list additionally contains the element:

`minLL` |
the minimum of the negative log-likelihood |

.

David Warton <[email protected]> and Ulrike Naumann.

Warton D.I. (2008). Penalized normal likelihood and ridge regularization of
correlation and covariance matrices. *Journal of the American
Statistical Association* 103, 340-349.

1 2 3 4 5 | ```
data(spider)
spiddat <- mvabund(spider$abund)
X <- spider$x
ridgeParamEst(dat = spiddat, X = model.matrix(spiddat~X))
``` |

mvabund documentation built on May 29, 2017, 9:50 a.m.

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.