토목공학과

#1번
data=read.csv("mid_bucheon.csv", fileEncoding = "euc-kr")
head(data)
install.packages("psych")
library(psych)
pairs.panels(data)
data2=data[ , 2:5]
head(data2)
data2$불량률=with(data2, max(불량률)-불량률)
head(data2)
h_pca=prcomp(data2, scale=TRUE)
h_pca
summary(h_pca)
plot(h_pca, type='l')
cor(data$최종평가,h_pca$x[,1])
biplot(h_pca, choices=c(1,2))

#2번
data("attitude")
attitude
write.csv(attitude, file="atti.csv")
data=read.csv("atti.csv", fileEncoding = "euc-kr")
head(data)
data2=data[  , 2:7]
head(data2)
fa1=factanal(data2,3)
print(fa1,digits=2,sort=T)
install.packages("GPArotation")
library(GPArotation)
fa3=factanal(data2, 3, rotation="none")
print(fa3, digits=2, sort=T)

#3번
mid=read.csv("mid_bucheon.csv", fileEncoding = "euc-kr")
head(mid)
mid2=mid [ , 2:5]
head(mid2)
hc2=hclust(dist(mid2,method="euclidean"),method="complete")
plot(hc2)

hc3=hclust(dist(mid2,method="euclidean"),method="single")
plot(hc3)

install.packages("mclust")
library(mclust)

USA_scale=scale(mid2)
mc=Mclust(USA_scale)
summary(mc)