Penggandaan vektor dalam BUGS dan JAGS

9

Dalam R, c (3,1,0) * c (2,0,1) == c (6,0,0). Ini bukan produk titik dan bukan produk lintas. Pertama, apa nama untuk produk ini, dan kedua, apakah itu berfungsi di WinBUGS, OpenBUGS, dan / atau JAGS?

bugs jags

— Jack Tanner
sumber

2

ini hanya penggandaan elemen-bijaksana. Saya tidak yakin apakah itu berfungsi di WinBUGS, OpenBUGS, atau JAGS.

— Diasumsikan normal

3

Tidak seperti JAGS, WinBUGS dan OpenBUGS tidak melakukan bentuk vektorisasi ini; Anda harus menulis satu lingkaran, dan menghitung setiap elemen 'dengan tangan', seperti dijelaskan di atas.

— tamu
sumber

4

Martyn Plummer menunjukkan bahwa ini diimplementasikan dalam JAGS, yang saya lewatkan saat membaca manual. Dari Ch 5:

Fungsi skalar yang mengambil argumen skalar secara otomatis di-vektor-kan. Mereka juga bisa dipanggil ketika argumen array dengan dimensi yang sesuai, atau skalar. Jadi, misalnya, skalar dapat ditambahkan ke matriks menggunakan $c$ $A$
B <- A + c
bukannya bentuk yang lebih verbose
D <- dim(A)
for (i in 1:D[1])
  for (j in 1:D[2]) {
    B[i,j] <- A[i,j] + c
  }
}

— Jack Tanner
sumber

2

Untuk melakukan perkalian elemen-bijaksana, Anda bisa membuat for for loop dalam bahasa tersebut dan hanya itu! Saya telah menggunakan loop di WinBUGS tanpa masalah.

— Ingin tahu
sumber

Apa pertanyaan dari alamat balasan ini? Tampaknya tidak relevan di sini.

— whuber

@ Whubber, mengapa? Ini sangat relevan. Ok, saya mengubah posting sedikit menjadi lebih jelas.

— Penasaran

Yup, for for adalah apa yang telah saya lakukan sejauh ini; Saya hanya ingin tahu apakah versi vektor itu mungkin.

— Jack Tanner

Saya telah mengirimkan permintaan fitur ke JAGS: sourceforge.net/tracker/…

— Jack Tanner

Terima kasih, Tomas. Sekarang saya melihat hubungannya: Anda tidak menjawab pertanyaan seperti yang dinyatakan, tetapi Anda menawarkan solusi.

— whuber

2

Kebetulan, perkalian elemen-bijaksana dari dua vektor dengan panjang yang sama disebut produk Hadamard (alias produk Schur).

— Anders Gorm
sumber