現在,R言語でパーティクルフィルタを用いて,ある問題の重み付け平均を算出しようとしています.
その問題というのが「りんごとみかんを1個ずつ買ったときの価格が100円.これまでの経験からりんごは60円程度,みかんは30円程度である.りんごとみかんはいくらか」というものです.
以下スクリプトの現状なのですが,"lik"で尤度の計算を行なっております.
ここからどのように,1000個のサンプルに重みを付ければ良いのでしょうか?
n <- 1000 x1 <- rnorm(n,60,20);x2 <- rnorm(n,30,10) rmat <- diag(n)*10^2 t <- (100-x1-x2)/10 lik <- (1/((2*pi)^(n/2)*(rmat^0.5)))*exp(-0.5*t^2)
あなたの回答
tips
プレビュー