#ho17_복습1 y=c(43,63,71,61,81,43,58,71,72,67,64,67,69,68,77,81,74,65,65,50, 50,64,53,40,63,66,78,48,85,82) x1=c(51,64,70,63,78,55,67,75,82,61,53,60,62,83,77,90,85,60,70,58, 40,61,66,37,54,77,75,57,85,82) x2=c(30,51,68,45,56,49,42,50,72,45,53,47,57,83,54,50,64,65,46,68, 33,52,52,42,42,66,58,44,71,39) x3=c(39,54,69,47,66,44,56,55,67,47,58,39,42,45,72,72,69,75,57,54, 54,62,50,58,48,63,74,45,71,59) x4=c(61,63,76,54,71,54,66,70,71,62,58,59,55,59,79,60,79,55,75,64, 45,66,63,50,66,88,80,51,77,64) x5=c(92,73,86,84,83,49,68,66,83,80,67,74,63,77,77,54,79,80,85,78, 64,80,80,57,75,76,78,83,74,78) x6=c(45,47,48,35,47,34,35,41,31,41,34,41,25,35,46,36,63,60,46,52, 33,41,37,49,33,72,49,38,55,39) #p459 x1=c(10,8,12,7,12,8,12,5,8,5,11,12,6,10,10) x2=c(8,6,9,16,17,15,8,10,4,16,7,12,6,4,4) x3=c(5.5,2.5,8,3,2.9,3,8,9,4,6.5,5.5,5,6,5,3.5) y=c(79,200,163,200,178,146,31,292,160,339,160,86,237,107,155) x=cbind(x1,x2,x3) ####설명변수 개수에 따라 달라짐 mydata=data.frame(y,x1,x2,x3) ####설명변수 개수에 따라 달라짐 pairs(x) cor(x) vif=diag(solve(cor(x))) #분산팽창인자 vif eigv=eigen(cor(x))$values ci=sqrt(eigv[1]/eigv) #상태지수 ci ls=regsubsets(y~x1+x2+x3, data=mydata) #모든 모형 고려--(나) ####설명변수 개수에 따라 달라짐 rs=summary(ls) rs #각 p에 따라 최적모형 제시 cp=rs$cp #선택된 모형의 cp--(다) cp adjr2=rs$adjr2 #선택된 모형의 adjr2--(다) adjr2 fit=lm(y~x1+x2+x3,data=mydata) ####설명변수 개수에 따라 달라짐 step=stepAIC(fit, direction="both") #(라), AIC 구하기 #17_1 y=c(78.5,74.3,104.3,87.6,95.9,109.2,102.7,72.5,93.1,115.9,83.8,113.3,109.4) x1=c(7,1,11,11,7,11,3,1,2,21,1,11,10) x2=c(26,29,56,31,52,55,71,31,54,47,40,66,68) x3=c(6,15,8,8,6,9,17,22,18,4,23,9,8) x4=c(60,52,20,47,33,22,6,44,22,26,34,12,12) #17_2peru #y=blpr x1=age x2=year x3=wt x4=abdo y=c(170,120,125,148,140,106,120,124,134,116,114,130,118,138,134, 120,120,114,124,114,136,126,124,128,128,112,118,134,128, 140,138,118,110,142,134,116,132,152) x1=c(21,22,24,24,25,27,28,31,32,33,33,34,35,35,36,36,37,37, 38,38,38,38,39,39,39,41,41,41,42,43,43,43,44,44,45,47,50,54) x2=c(1,6,5,1,1,19,5,6,13,13,10,15,18,2,12,15,16,17,10,18,11,11, 21,24,14,25,32,5,12,25,26,10,19,18,10,1,43,40) x3=c(71,56.5,56,61,65,62,53,65,57,66.5,59.1,64,69.5,64,56.5,57, 55,57,58,59.5,61,57,57.5,74,72,62.5,68,63.4,68,69,73,64,65, 71,60.2,55,70,87) x4=c(12.7,8,4.3,4.3,20.7,5.7,8,10,6,8.3,10.3,7,7,6.7,11.7,6,7,11.7, 13,7.7,4,3,5,15.7,13.3,8,11.3,13.7,10.7,6,5.7,7,7.7,4.3,3.3,4, 11.7,11.3) #17_3river #y=평균 질소 농도 x1=농지 비율 x2=산림지 비율 #x3=주거지 비율 x4=상업/산업 용지 비율 y=c(1.10,1.01,1.90,1.00,1.99,1.42,2.04,1.65,1.01,1.21,1.33,0.75, 0.73,0.80,0.76,0.87,0.80,0.87,0.66,1.25) x1=c(26,29,54,2,3,19,16,40,28,26,26,15,6,3,2,6,22,4,21,40) x2=c(63,57,26,84,27,61,60,43,62,60,53,75,84,81,89,82,70,75,56,49) x3=c(1.2,0.7,1.8,1.9,29.4,3.4,5.6,1.3,1.1,0.9,0.9,0.7,0.5,0.8,0.7,0.5, 0.9,0.4,0.5,1.1) x4=c(0.29,0.09,0.58,1.98,3.11,0.56,1.11,0.24,0.15,0.23,0.18,0.16, 0.12,0.35,0.35,0.15,0.22,0.18,0.13,0.13)