サンプル重要度リサンプリングアルゴリズムの次のコードに問題があります。未使用の引数未使用引数エラーRコード
シータの宣言のために:それは本当
ログ=返します。
T = 1e5 theta <- runif(T,0,1, log=TRUE) log.p <- function(x) dbeta(x, 3000+711, 17000+2201-711) log.s <- function(x) dunif(x, 0, 1,) w <- function(t) log.p(t)/log.s(t) HA <- sum(w(theta)%*%theta)/T
ログはrunifの引数だったという考えをどこから得ましたか? –
Rコンソールに '?runif'と入力して文書を確認してください。 'log'引数はありません。値が渡される 'n'、' min'、 'max'の3つの引数だけがあります。 –
'theta = log(runif(T、0,1))'を探していますか? –