![# デ`タのロ`ド
data(state)
# 聞喘デ`タは state.x77。50 デ`タ 8 篳。
state.x77 # デ`タ匯E
# 01A どのようなデ`タか靖めてみよう
summary(state.x77) ##### 児云y楚
pairs(state.x77) ##### 柊下輟仭价荵
# 柊下輟仭个 A4 サイズの png 侘塀で隠贋
png(filename="pairs_x77.png", width=1654, height=2339, pointsize = 32) #### 200 dpi A4 k
pairs(state.x77) ##### 柊下輟仭价荵
dev.off()
# 8 篳のヒストグラムを A4 サイズにまとめて宙鮫、png 侘塀隠贋
require(MASS) # truehist v方がуhされているパッケ`ジ
png(filename="hist_x77.png", width=1654, height=2339, pointsize = 32) #### 200 dpi A4 k
par(mfrow=c(4,2)) # 4 佩 2 双に鮫中蛍護
for (i in 1:8) truehist(state.x77[,i], xlab=colnames(state.x77)[i])
dev.off()
# 咀徨蛍裂
# R では咀徨の方はユ`ザが峺協する
# いくつかしてみて、恷mな方をQめる
fac2 <- factanal(state.x77, factors=2)
fac3 <- factanal(state.x77, factors=3)
fac4 <- factanal(state.x77, factors=4)
fac5 <- factanal(state.x77, factors=5)
# 咀徨方をxぶ
# ここでは圷のデ`タの秤鵑隆鷁新屬鬟バ`するよりも、
# 盾しやすい咀徨をxkする圭とするので、まず咀徨塞楚
# をみて盾辛嬬來を篇する
# 咀徨塞楚 光篳の咀徨への
par(mfrow=c(1,2)) # 1 佩 2 双に鮫中蛍護
barplot(fac2$loading[,1], las=2)
barplot(fac2$loading[,2], las=2)
par(mfrow=c(2,2)) # 2 佩 2 双に鮫中蛍護
barplot(fac3$loading[,1], las=2)
barplot(fac3$loading[,2], las=2)
barplot(fac3$loading[,3], las=2)](https://image.slidesharecdn.com/r03-151024163217-lva1-app6892/85/R-03-1-320.jpg)
![par(mfrow=c(2,2)) # 2 佩 2 双に鮫中蛍護
barplot(fac4$loading[,1], las=2)
barplot(fac4$loading[,2], las=2)
barplot(fac4$loading[,3], las=2)
barplot(fac4$loading[,4], las=2)
# 鏡徭咀徨 モデルでh苧できない秤鵑留搬
round(fac2$uniquenesses, 3)
round(fac3$uniquenesses, 3)
round(fac4$uniquenesses, 3)
# 咀徨蛍裂の咀徨誼泣┐海海任 2 咀徨を res1 に鯉{
res1 <- factanal(state.x77, factors=2, scores="Bartlett")$scores
# クラスタ`蛍裂
res1.hc <- hclust(dist(res1), method="single") # 恷除O隈
#res1.hc <- hclust(dist(res1), method="ward.D2")
par(mfrow=c(1,1)) # 鮫中蛍護を圷にす
plot(res1.hc) # デンドログラム燕幣
# デンドログラムを A4 罪にファイル竃薦
png(filename="Dendrgram_x77.png", width=2339, height=1654, pointsize = 32) #### 200 dpi
A4 k
plot(res1.hc) # デンドログラム燕幣
dev.off()
# とりあえずいくつかに蛍してみる
res1.cl6 <- cutree(res1.hc,k=6)
res1.cl8 <- cutree(res1.hc,k=8)
res1.cl10 <- cutree(res1.hc,k=10)
# 咀徨塞楚のプロット、蛍弼蛍けで
plot(res1, col=rainbow(10)[res1.cl10], pch=20, ylim=c(-2,3.5), xlim=c(-2,2.5))
abline(h=-3:3, v=-3:3, col="gray", lty=3)
text(res1, rownames(state.x77))
points(res1, pch=20, col=rainbow(10)[res1.cl10])
############################################################](https://image.slidesharecdn.com/r03-151024163217-lva1-app6892/85/R-03-2-320.jpg)
Rプログラミング03 仝デ`タ蛍裂｡好妊
- 1. # デ`タのロ`ド data(state) # 聞喘デ`タは state.x77。50 デ`タ 8 篳。 state.x77 # デ`タ匯E # 01A どのようなデ`タか靖めてみよう summary(state.x77) ##### 児云y楚 pairs(state.x77) ##### 柊下輟仭价荵 # 柊下輟仭个 A4 サイズの png 侘塀で隠贋 png(filename="pairs_x77.png", width=1654, height=2339, pointsize = 32) #### 200 dpi A4 k pairs(state.x77) ##### 柊下輟仭价荵 dev.off() # 8 篳のヒストグラムを A4 サイズにまとめて宙鮫、png 侘塀隠贋 require(MASS) # truehist v方がуhされているパッケ`ジ png(filename="hist_x77.png", width=1654, height=2339, pointsize = 32) #### 200 dpi A4 k par(mfrow=c(4,2)) # 4 佩 2 双に鮫中蛍護 for (i in 1:8) truehist(state.x77[,i], xlab=colnames(state.x77)[i]) dev.off() # 咀徨蛍裂 # R では咀徨の方はユ`ザが峺協する # いくつかしてみて、恷mな方をQめる fac2 <- factanal(state.x77, factors=2) fac3 <- factanal(state.x77, factors=3) fac4 <- factanal(state.x77, factors=4) fac5 <- factanal(state.x77, factors=5) # 咀徨方をxぶ # ここでは圷のデ`タの秤鵑隆鷁新屬鬟バ`するよりも、 # 盾しやすい咀徨をxkする圭とするので、まず咀徨塞楚 # をみて盾辛嬬來を篇する # 咀徨塞楚 光篳の咀徨への par(mfrow=c(1,2)) # 1 佩 2 双に鮫中蛍護 barplot(fac2$loading[,1], las=2) barplot(fac2$loading[,2], las=2) par(mfrow=c(2,2)) # 2 佩 2 双に鮫中蛍護 barplot(fac3$loading[,1], las=2) barplot(fac3$loading[,2], las=2) barplot(fac3$loading[,3], las=2)
- 2. par(mfrow=c(2,2)) # 2 佩 2 双に鮫中蛍護 barplot(fac4$loading[,1], las=2) barplot(fac4$loading[,2], las=2) barplot(fac4$loading[,3], las=2) barplot(fac4$loading[,4], las=2) # 鏡徭咀徨 モデルでh苧できない秤鵑留搬 round(fac2$uniquenesses, 3) round(fac3$uniquenesses, 3) round(fac4$uniquenesses, 3) # 咀徨蛍裂の咀徨誼泣┐海海任 2 咀徨を res1 に鯉{ res1 <- factanal(state.x77, factors=2, scores="Bartlett")$scores # クラスタ`蛍裂 res1.hc <- hclust(dist(res1), method="single") # 恷除O隈 #res1.hc <- hclust(dist(res1), method="ward.D2") par(mfrow=c(1,1)) # 鮫中蛍護を圷にす plot(res1.hc) # デンドログラム燕幣 # デンドログラムを A4 罪にファイル竃薦 png(filename="Dendrgram_x77.png", width=2339, height=1654, pointsize = 32) #### 200 dpi A4 k plot(res1.hc) # デンドログラム燕幣 dev.off() # とりあえずいくつかに蛍してみる res1.cl6 <- cutree(res1.hc,k=6) res1.cl8 <- cutree(res1.hc,k=8) res1.cl10 <- cutree(res1.hc,k=10) # 咀徨塞楚のプロット、蛍弼蛍けで plot(res1, col=rainbow(10)[res1.cl10], pch=20, ylim=c(-2,3.5), xlim=c(-2,2.5)) abline(h=-3:3, v=-3:3, col="gray", lty=3) text(res1, rownames(state.x77)) points(res1, pch=20, col=rainbow(10)[res1.cl10]) ############################################################