R에서 다음과 같은 복잡한 합계를 계산하는 방법은 무엇입니까? [닫은]

Question

 a1=c(1,2,3, 4)
 a2=c(4,5,6, 7)
 a3=c(7,8,9,10)
 x <- cbind(a1, a2, a3)
  y <- c(1,5,7,9)
   z<-c(2,3,5)

다음 합계를 계산하고 싶습니다.

Answer 1

다음을 수행 할 수 있습니다.

diag(crossprod(x, drop(x %*% z) - x * rep(z, each = NROW(x))))
#R> [1]  630 1080  756

또는 sapply/ vapply버전의 루프 ekoam이 주석에 표시됩니다.

vapply(seq_along(z), function(j) drop(x[, j] %*% (x[, -j] %*% z[-j])), 
       numeric(1))
#R> [1]  630 1080  756