
Plot the results of cross-validation for Bhat obtained from the rrda.cv function.
rrda.plot.RdThis function visualizes the results of cross-validation for the estimated Bhat matrix obtained from the rrda.cv function. It creates a plot of the Mean Squared Error (MSE) for each combination of rank and lambda regularization parameter, allowing for the selection of specific ranks and lambda ranges to be plotted. Error bars representing the standard error of the MSE can be displayed for the best rank.
Usage
rrda.plot(
cv_result,
nrank = NULL,
min_l = NULL,
max_l = NULL,
show_error_bar = FALSE,
title = NULL
)Arguments
- cv_result
A result list from the function
rrda.cv, containing a matrix of MSE values for each rank and lambda, and a vector of lambda values.- nrank
A numeric vector specifying the ranks of Bhat to be plotted. Default is
NULL, which plots all ranks.- min_l
Minimum lambda value to be plotted. Default is
NULL, which uses the minimum lambda value incv_result.- max_l
Maximum lambda value to be plotted. Default is
NULL, which uses the maximum lambda value incv_result.- show_error_bar
Logical value indicating if the error bar is shown on the line that gives the best MSE value.
- title
Title of the figure
Examples
if (FALSE) { # \dontrun{
set.seed(10)
simdata<-rdasim1(n = 100,p = 200,q = 200,k = 3) # data generation
X <- simdata$X
Y <- simdata$Y
cv_result<- rrda.cv(Y = Y, X = X, maxrank = 5, nfold = 5) # cv
rrda.summary(cv_result = cv_result)
rrda.plot(cv_result = cv_result)
} # }