昨日は混合比の推定も含めた、混合分布推定をBUGSでやった(こちら)ある一次元データが得られたとき、その個々の観測を中心としたカーネル関数の混合分布として背景の分布を推定することがあるカーネル関数として正規分布を取ることも多く、そのときに問題となるのは、正規分布の分散をいくつにしたらよいかを調整する必要があるということこれをBUGSでやってみようまずはモデル観測値がyというベクトルで観測数がNそのyのそれぞれの値を中心とし、共通のtauによって定められる正規分布(N個の正規分布)から、値を採取する(この中心の値にはmmという変数を使って書いてある。yと同じベクトル)N個の分布のどれにするかは、等確率で取る(dcat(piで指定。このpiに、等確率ベクトルを与える)。混合分布だが、N個の混合分布。g=1:NとしてN個のグループを与える。GWAS/GWASメタアナリシスとその周辺解析用のアプリケーションへのリンク複合遺伝性疾患とは----------News----------20万アクセス2008/07/2040万アクセス2010/09/12-------------------------はてなブログをはじめよう!ryamada22さんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。