#単回帰分析
dat1=read.csv("dat1.csv")
dat1
res1=lm(HbA1c~Glucose,data=dat1)
summary(res1)
cor(dat1$HbA1c,dat1$Glucose)
plot(dat1)
abline(res1)

#重回帰分析
dat2=read.csv("dat2.csv")
dat2
res2=lm(Cmax~RBC+Ht+Alb+gammaGTP+LDH, data=dat2)
summary(res2)

#変数選択
install.packages("psych")
library(psych)
cor(dat2)
plot(dat2)
res2b=lm(Cmax~RBC+LDH, data=dat2)
summary(res2b)
res2b$fitted.values
plot(dat2$Cmax,res2b$fitted.values)
abline(0,1)