# 7.1節　売上データの分析
candy0 <- scale(candy, center=T, scale=F) # 平均偏差行列を作る
y <- candy0[,"sales"]    # 基準変数のベクトル
X <- candy0[,2:3]        # 説明変数の行列を用意する
# 偏回帰係数の推定 
b <-solve(t(X)%*%X)%*%t(X)%*%y
# 外部ファイルに保存
attach(candy) #  データフレーム candy の変数名を使えるようにする
yhat <- mean(sales)+X%*%b
error <- sales - yhat
candyresult <- data.frame(candy,yhat,error)
write.csv(candyresult, file = "output.csv" ,row.names = F)

# Rで相関係数を求める
r23 <- cor(temp,time)  # 説明変数間の相関係数
r <- cor(y,yhat)  # 重相関係数
determ <- r^2   # 決定係数 coefficient of determination
