通常の頻度統計で行う。 library(AER) # CPS1985 データを含むパッケージ data("CPS1985", package = "AER") fit.freq <- lm(wage ~ age, data = CPS1985) print(fit.freq) 結果。 Call: lm(formula = wage ~ age, data = CPS1985) Coefficients: (Intercept) age 6.16747 0.07755 ベイズ統計での計算 単回帰のstanファイルを作る data { int<lower=0> N; // データポイントの数 ve…