2008-08-23から1日間の記事一覧

オブジェクトからの値の取り出し

R glm

resultを取り出すには・・・ summary summary$coefficients, もしくは、summary[12]に以下の内容がリストで格納されているので Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -0.5898 0.2627 -2.246 0.02473 * G 0.5073 0.1880 2.699 0.0…

漸近近似ベイズ因子

こちらにGWASでのSNP関連検定に漸近近似ベイズ因子を用いる論文について書いた。 以下、20080821と重複になる部分もあるが、漸近近似に中心を移して書く。 事前確率分布のおき方によっては、解析的にベイズ因子を算出するのは難しい。 適当な事前確率分布を…