setwd("d:/r")

product <- read.csv("product.csv", header=T)
product

Y <- product[,1] ; LY=log(Y) # データ Y は１列目にある
K <- product[,2] ; LK=log(K) # データ K は２列目にある
L <- product[,3] ; LL=log(L) # データ L は３列目にある

library(lmtest)
library(tseries)

# -------- log(Y) = a + b*log(K) + c*log(L) + u --------------
result <- lm(LY~LK+LL) # LY=a+b*LK+c*LL+u
summary(result)
deviance(result)

# 誤差項の仮定をテストする統計量
fv <- result$fitted.values
bptest(result,~fv**2)
dwtest(result)
e <- residuals(result)
jarque.bera.test(e)
resettest(result,power=2:3)
logLik(result)
AIC(result)

plot(result)

# -------- log(Y/L) = a + b*log(K/L) + u --------------
LYL=log(Y/L)
LKL=log(K/L)
result <- lm(LYL~LKL)
summary(result)

# 誤差項の仮定をテストする統計量
fv <- result$fitted.values
bptest(result,~fv**2)
dwtest(result)
e <- residuals(result)
jarque.bera.test(e)
reset(result)
logLik(result)
AIC(result)