(3)5.3.2項の続きで、アルバイトが好きかどうかごとにYを集計せよ元データはこんなもの。data-attendance-3.txt PersonID,A,Score,Weather,Y 1,0,69,B,1 1,0,69,A,1 1,0,69,C,1 1,0,69,A,1 1,0,69,B,1 1,0,69,B,1 1,0,69,C,0 1,0,69,B,1 1,0,69,A,1コードは…

(1)model5-3.stan(P.57)をRから実行して得られるMCMCサンプルを使って、それぞれのnにおけるε[n]=Y[n]-μ[n]のMCMCサンプルを算出するRコードを書けインプットするデータは以下のようなもの（冒頭の１０行ほど）data-attendance-1.txt A,Score,Y 0,69,0.286 1…

t検定に相当することをStanで行う。以下のRコードで生成されたデータを用いて、グループ１（Y1）の平均パラメーターμ1とグループ２（Y2）の平均パラメーターμ2に差があるか、Prob[μ1 set.seed(123) N1 <- 30 N2 <- 20 Y1 <- rnorm(n=N1, mean=0, sd=5) Y2 <-…