Educational dataset example

res <- ordinal_regression_simplex(
  xdata,
  ydata,
  weights = weights_orig,
  lambda2_grid = 0,
  method = "gcv",
  do_bootstrap = TRUE,
  B = 200,
  compute_opi = TRUE,
  compute_R2 = TRUE,
  compute_OCC = TRUE
)

cat("Estimated A_hat\n")

## Estimated A_hat

print(round(res$A_hat, 4))

##        [,1]   [,2]   [,3]
## [1,] 0.8814 0.0561 0.0387
## [2,] 0.1186 0.9439 0.0282
## [3,] 0.0000 0.0000 0.9331

cat("Wasserstein R2\n")

## Wasserstein R2

print(res$R2)

## [1] 0.5296394

cat("OPI\n")

## OPI

print(res$OPI)

## [1] 0.8494624

cat("OCC\n")

## OCC

print(res$OCC)

## [1] 0.8516129

cat("Frechet mean matrix\n")

## Frechet mean matrix

print(round(res$bootstrap$Frechet, 4))

##        [,1]  [,2]   [,3]
## [1,] 0.7713 0.062 0.0387
## [2,] 0.2287 0.938 0.0311
## [3,] 0.0000 0.000 0.9302