Rで散布図（仮）

#csvファイルの読み込み、１列目をrow.namesに設定

d <- read.csv("f:/Program Files/R/R-2.9.1/data1.csv",header=T,row.names=1)
d
library(lattice)
xyplot(d[,1] ~ d[,3])
plot(d)

plot(d[,1]*100,d[,3]*100,type="n",xlab="住宅シェア（％）",ylab="総預金シェア（％）")
text(d[,1]*100,d[,3]*100,rownames(d))

plot(d[,3]*100,d[,1]*100,type="n",ylab="住宅シェア（％）",xlab="総預金シェア（％）")
text(d[,3]*100,d[,1]*100,rownames(d))

plot(d[,3]*100,d[,4],type="n",ylab="預金顧客数（人）",xlab="預金シェア（％）")
text(d[,3]*100,d[,4],rownames(d))

plot(d[,4],d[,2],type="n",xlab="預金顧客数（人）",ylab="住宅顧客数（人）")
text(d[,4],d[,2],rownames(d))

#df1　へ預金・貸金の契約者のみ移行する
df2 <- data.frame(cbind(d[,2],d[,4]))
#colnames(df2) <- c(colnames(d[2]),colnames(d[4]))
df2
cor(df2)
plot(df2)

df2.lm <-lm(X2~X1,data=df2)
summary(df2.lm)
plot(df2)
abline(df2.lm,lwd=2)
plot(df2.lm)


cor(d[,4],d[,2])
d.lm <-lm(d[,4]~d[,2],data=d)
summary(d.lm)
plot(d)

plot(d[,4],d[,2])
abline(df2.lm,lwd=2)


#df1　へ預金・貸金のシェアのみ移行する
df1 <- data.frame(cbind(rownames(d),d[,1],d[,3]),row.names=1)
colnames(df1) <- c(colnames(d[1]),colnames(d[3]))
df1

df1.d <-dist(df1)
round(df1.d)
df1.d
df1.hc<- hclust(df1.d)
summary(df1.hc)
plot(df1.hc)
df1.hc<- hclust(df1.d,method="ward")
summary(df1.hc)
plot(df1.hc)

#df1　へ預金契約者、預金・貸金のシェアのみ移行する
df1 <- data.frame(cbind(rownames(d),d[,1],d[,3]),row.names=1)
colnames(df1) <- c(colnames(d[1]),colnames(d[3]))
df1

df1.d <-dist(df1)
round(df1.d)
df1.d
df1.hc<- hclust(df1.d)
summary(df1.hc)
plot(df1.hc)
df1.hc<- hclust(df1.d,method="ward")
summary(df1.hc)
plot(df1.hc)