問題タブ [rstan]

私は openCPU を使用して API を作成し、rstan でかなり基本的なモデルを推定しています。

openCPU から呼び出されるのではなく、通常の R コンソール セッションから呼び出されたときに、関数が自分の環境で機能することを確認しました。

しかし、openCPU 経由で呼び出すと、応答は次のようになります。

rstan のソースにチェックポイントを挿入することで、これを R/cxxfunplus.R の 1 つの呼び出しに追跡しました ( github の行へのリンクはこちら):

この呼び出しで、失敗する可能性のある接続を読み取ろうとしているのは何ですか?

私のデータ (m) の尤度関数は次のようになります。

STAN コード内に記述したユーザー定義の対数尤度関数は次のようになります。

. 私はそれを正しく行っているかどうか知りたいですか？ユーザー定義関数内で「exponential_log」や「normal_cdf_log」などの STAN 関数を使用できますか?

github.com/stan-dev/example-models のバージョンを使用して、古いバグの例「leuk」のスタン バージョンを動作させようとしています。そこにある leuk.stan のバージョンは正常に動作しますが、「<-」、「increment_log_prob」、および「poisson_log」はすべて非推奨であるという警告が表示されます。「<-」をすべて「=」に編集し、「increment_log_prob」を「target +=」に編集しましたが、「poisson_log」の正しい置き換えが何であるかわかりません。Stan Reference 2.14.0 を確認したところ、490 ページで、「poisson_log」に相当する新しいものは「poisson_log_lpmf」であるという提案が見つかりました (「,」を「|」に変更)。これは構文的に正しく実行されますが、間違った答えが返されます。「poisson_log」の正しい置き換えを教えてもらえますか? よろしくお願いします。

私はスタンと確率的プログラミングが初めてです。非線形成長モデルを構築しようとしています。私はモデルを構築することができましたNLS

私が使用したNLS式は次のとおりです。Trump_Pct ~ alpha - beta * lambda^Population

私のNLSの要約は次のとおりです。

言い換えれば、基本的な指数減衰曲線です。私はスタンと複製しようとしています。

私のデータは次のとおりです。

データセットにN観測があります。予測子は郡の人口 (「人口」) であり、予測された Y はトランプへの投票率「Trump_Pct」です。

このモデルを構築する 2 つの方法を試しました。

1. 1 つは、データの各コンポーネントをモデルにベクトルとして渡します。

2. もう 1 つの方法では、各データ コンポーネントをリストとして残し、各データ ポイントを使用しようとします。

どちらの場合も、モデルを正常に実行することはできません。

ここに私のモデルがあります：

ケース 1:

これは、このモデルの適応です。

ここでは、Trump_Pct 列と Population 列のベクトル化されたバージョンを作成しました。

このモデルは、指数の行で次の理由で失敗します。

`構文エラー、パーサーからのメッセージ:

^ への引数はプリミティブ (実数または整数) でなければなりません。block=local` では実数をベクトルで累乗できません

使用してみましpow()たが、先に進む方法が見つかりません。任意のヒント？

ケース 2:

ケース 2 では、パラメーターの推定値を範囲内に保つことができません。

"Informational Message: The current Metropolis proposal is about to be rejected because of the following issue:"
[2] "Exception thrown at line 21: normal_log: Location parameter[2873] is -inf, but must be finite!"

私の式の単純な非線形モデルについてアドバイスをくれる人はいますか?