この言語って何の言語か分かりますか?
my.cancor <- function(x, gr1, gr2)
{
geneig2 <- function(a, b, k, sd)
{
a <- as.matrix(a)
b <- as.matrix(b)
if (nrow(a) == 1) {
res <- list(values=a/b, vectors=as.matrix(1))
}
else {
res <- eigen(b)
g <- diag(1/sqrt(res$values))
v <- res$vectors
res <- eigen(g %*% t(v) %*% a %*% v %*% g)
res$vectors <-v %*% g %*% res$vectors
}
unstd.vectors <- (std.vectors <- res$vectors[,1:k])/sd
list(values=res$values[1:k], std.vectors=std.vectors, unstd.vectors=unstd.vectors)
}