S1 = matrix(c(1, 0.9, 0.8, 0.5,
0.9, 1.1, 0.7, 0.4,
0.8, 0.7, 1.5, 0.2,
0.5, 0.4, 0.2, 1), ncol = 4)
S2 = t(S1)%*% S1
S2 = S2/2
explained_var(COVS = list(S1, S2),
PC = c(1,0,0,0,sqrt(2),0,0,-sqrt(2)),
k = 1,
cor = FALSE,
gamma = 0.5)
explained_var(COVS = list(cov2cor(S1), cov2cor(S2)),
PC = c(1,0,0,0,sqrt(2),0,0,-sqrt(2)),
k = 1,
cor = TRUE,
gamma = 0.5)
Run the code above in your browser using DataLab