R에서이 합계를 계산하는 방법은 무엇입니까?
Nov 26 2020
데이터 세트가있는 경우 $x_{ij}$크기 4 x 3 행렬, i = 1,2,3,4 및 j = 1,2,3. y_i는 벡터이고 z_ {ij}도 4x3 행렬입니다.
a1=c(1,0,1, 2)
a2=c(2,2,1, 2)
a3=c(1,3,2,3)
a4=c(-1,0,-1,-1)
x <- cbind(a1, a2, a3)
y <- c(1,2,3,2)
z<-cbind(a2, a3,a4)
답변
1 BenjaminChristoffersen Nov 26 2020 at 21:26
나는 이것이 할 것이라고 생각한다.
z1 <- z[1, ]
vapply(
seq_along(z1), function(j)
drop((y - x[, -j] %*% z1[-j])) %*% x[, j],
numeric(1))
#R> [1] 10 15 -15
for 루프 버전은 다음과 같습니다.
out <- numeric(3)
for(j in 1:3)
for(i in 1:4){
f1 <- y[i]
for(k in setdiff(1:3, j))
f1 <- f1 - x[i, k] * z[1, k]
out[j] <- out[j] + f1 * x[i, j]
}
out
#R> [1] 10 15 -15