![](https://cdn-ak-scissors.b.st-hatena.com/image/square/56298f6af967f519ddf2f7692c9f0a12dc357427/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-9f5428127621718a910c8b63951390ad.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTkxNiZoPTMzNiZ0eHQ9U3RhbiUyMCVFMyU4MSVBNyVFNiU5NiVCMCVFMyU4MSU5NyVFMyU4MSU4NCVFNSU4MCVBNCVFMyU4MSVBQiVFNSVBRiVCRSVFMyU4MSU5OSVFMyU4MiU4QiVFNCVCQSU4OCVFNiVCOCVBQyVFNSU4QyVCQSVFOSU5NiU5MyVFMyU4MiU5MiVFNSU4NyVCQSVFMyU4MSU5OSUyMCUyM3JzdGF0c2omdHh0LWNvbG9yPSUyMzIxMjEyMSZ0eHQtZm9udD1IaXJhZ2lubyUyMFNhbnMlMjBXNiZ0eHQtc2l6ZT01NiZ0eHQtY2xpcD1lbGxpcHNpcyZ0eHQtYWxpZ249bGVmdCUyQ3RvcCZzPTMzMGUzNzhhMTEzOTNiNWM5NjY0ZDdiZTYxMzc5NjQ4%26mark-x%3D142%26mark-y%3D112%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTYxNiZ0eHQ9JTQwaG94b19tJnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9MzYmdHh0LWFsaWduPWxlZnQlMkN0b3Amcz01M2IzYTE1MzRjMTk5MTU1OWVkYTFkZDU2ZTNlOGQwNQ%26blend-x%3D142%26blend-y%3D491%26blend-mode%3Dnormal%26s%3D2166cd5f4a0ca2b721360a9920d29f20)
エントリーの編集
![loading...](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/common/loading@2x.gif)
エントリーの編集は全ユーザーに共通の機能です。
必ずガイドラインを一読の上ご利用ください。
Stan で新しい値に対する予測区間を出す #rstatsj - Qiita
記事へのコメント0件
- 注目コメント
- 新着コメント
このエントリーにコメントしてみましょう。
注目コメント算出アルゴリズムの一部にLINEヤフー株式会社の「建設的コメント順位付けモデルAPI」を使用しています
![アプリのスクリーンショット](https://b.st-hatena.com/bdefb8944296a0957e54cebcfefc25c4dcff9f5f/images/v4/public/entry/app-screenshot.png)
- バナー広告なし
- ミュート機能あり
- ダークモード搭載
関連記事
Stan で新しい値に対する予測区間を出す #rstatsj - Qiita
偉大なるごみ箱師に Stan で新しい値に対する予測区間を出す方法を教わったので、やってみる(ずいぶん放... 偉大なるごみ箱師に Stan で新しい値に対する予測区間を出す方法を教わったので、やってみる(ずいぶん放置していた)。 やり方は一言で言えば欠測値推定と同じやり方と言えばわかる人にはわかるだろう。 ######################################## パッケージの読み込み library(rstan) library(ggplot2) library(dplyr) ######################################## モデル作成用のデータを生成 set.seed(314) N <- 1000 x <- rnorm(N, mean = 50, sd = 10) y <- 10 + 0.8 * x + rnorm(N, mean =0, sd = 7) ggplot(data.frame(x=x, y=y), aes(x=x, y=y))