data923 <- read.csv("Table923.csv")

data923

​

attach(data923)

cor(data923)#相関行列

​

model1 <- lm(Y~X6)#公開授業の数
summary(model1)

AIC(model1)

​

model2 <- lm(Y~X7)#企画数
summary(model2)

AIC(model２)

​

model3 <- lm(Y~X6+X7)#企画数、公開授業の数
summary(model3)

AIC(model3)

​

​

step(model3)

​

​

# 標準偏回帰係数（β）を求める
z <- scale(data923)     # 得点を標準化
z <- data.frame(z)  # データフレーム形式に戻す 
summary(lm(Y~X6+X7, z))

​

 

#多重共線性のチェック
zx <- z[7:8]
r <- cor(zx)
VIF <- diag(solve(r))
tolerance = 1/VIF

data.frame(tolerance, VIF)

​