Phép nhân vectơ trong BUGS và JAGS

Trong R, c (3,1,0) * c (2,0,1) == c (6,0,0). Đây không phải là sản phẩm chấm và nó không phải là sản phẩm chéo. Đầu tiên, tên của sản phẩm này là gì và thứ hai, nó có hoạt động trong WinBUGS, OpenBUGS và / hoặc JAGS không?

Không giống như JAGS, WinBUGS và OpenBUGS không thực hiện hình thức vector hóa này; bạn phải viết một vòng lặp và tính toán từng phần tử 'bằng tay', như được mô tả ở trên.

Martyn Plummer chỉ ra rằng điều này được thực hiện trong JAGS, điều mà tôi đã bỏ lỡ khi đọc hướng dẫn. Từ Ch 5:

Các hàm vô hướng lấy các đối số vô hướng được tự động véc tơ. Chúng cũng có thể được gọi khi các đối số là các mảng có kích thước phù hợp hoặc vô hướng. Vì vậy, ví dụ, vô hướng có thể được thêm vào ma trận bằng cách sử dụng$c$$A$

B <- A + c

thay vì hình thức dài dòng hơn

D <- dim(A)
for (i in 1:D[1])
  for (j in 1:D[2]) {
    B[i,j] <- A[i,j] + c
  }
}

Để thực hiện phép nhân phần tử, bạn chỉ cần tạo một vòng lặp for trong các ngôn ngữ đó và đó là nó! Tôi đã sử dụng các vòng lặp trong WinBUGS mà không gặp vấn đề gì.

Ngẫu nhiên, phép nhân phần tử của hai vectơ có độ dài bằng nhau được gọi là sản phẩm Hadamard (còn gọi là sản phẩm Schur).