久し振りにlinear regressionをやってみた。備忘録として。
Stata、RとWinBUGSで。ベイズ統計についてはいずれ書きたいと思います。まぁ、分析の方が主だから哲学的な所にはあまりこだわってはいないけど。ただ、ベイズ統計の方がpost-estimationがより楽しい点はあるかな。
データはStataのmroz。モデルはここを参照。
*****Stata
set more off
set memory 500m
cd "E:\"
use mroz, clear
reg lwage exper expersq educ age kidslt6 kidsge6 city
*****R
setwd("E:/")
library(foreign)
mroz <-read.dta("mroz.dta")
N <- nrow(mroz)
lwage <- mroz$lwage
exper <- mroz$experexpersq <- mroz$expersq
educ <- mroz$educ
age <- mroz$age
kidslt6 <- mroz$kidslt6
kidsge6 <- mroz$kidsge6
city <- mroz$city
model <- lm(lwage ~ exper+expersq+educ+age+kidslt6+kidsge6+city)
model
*****BUGSlibrary(R2WinBUGS)
data <-list("N", "lwage", "exper", "expersq", "educ", "age", "kidslt6", "kidsge6", "city")
params <- c("b")
inits <- function () {list(b=runif(8))}
BUGSModel <- bugs (data, inits, params, "Model.txt", n.chains=5, n.iter=10000, debug=T)