внешняя функция R

Question

y <- matrix(c(1,2,3,4,5,6,7,8,9),3,3) 
p <- 3 
oldgamma <- as.matrix(rep(1,3)) 
func1 <- function(i, j, y, oldgamma) { 
    return(y[i,j]/(oldgamma[i] + oldgamma[j])) 
} 
func1 <- Vectorize(func1) 
outer(1:p, 1:p, func1, y, oldgamma) 

возвращает сообщение об ошибке: Ошибка в у [I, J]: неправильное число измерений

Может кто-нибудь помочь мне решить эту проблему?

Answer 1

Может быть, это помогает

y/(oldgamma[row(y)]+oldgamma[col(y)]) 
#  [,1] [,2] [,3] 
#[1,] 0.5 2.0 3.5 
#[2,] 1.0 2.5 4.0 
#[3,] 1.5 3.0 4.5 

Или

y/with(expand.grid(rep(list(1:p),2)), 
       oldgamma[Var1]+oldgamma[Var2]) 

Или используя код

func1 <- function(i, j, y, oldgamma) { 
    return(y[cbind(i,j)]/(oldgamma[i] + oldgamma[j])) 
} 
outer(1:p, 1:p, func1, y, oldgamma) 
# [,1] [,2] [,3] 
#[1,] 0.5 2.0 3.5 
#[2,] 1.0 2.5 4.0 
#[3,] 1.5 3.0 4.5