多元统计分析及R语言建模（第五版）——第3章

表3-2是2004年广东省各市高新技术产品情况。试对该资料按本章介绍的多元图示方法做直观分析。

地区	工业总产值	工业增加值	产品销售收入	产品出口销售收入
广州	1486.05	478.55	1454.67	31.75
韶关	28.55	7.56	27.75	0.77
深圳	3266.52	940.86	3120.18	350.6
珠海	465.07	74.6	447.08	37.74
汕头	147.73	43.8	140.01	5.39
佛山	930.49	191.99	891.08	45.63
江门	251.23	52.86	235.79	6.9
湛江	75.66	21.38	72.7	1.34
茂名	25.37	10.37	23.26	0.1
肇庆	70.06	16.9	65.64	3.17
惠州	525.67	109.43	495.99	40.52
梅州	18.09	4.72	16.97	0.26
汕尾	26.21	7.66	25.79	1.44
河源	2.41	0.68	2.17	0.01
阳江	29.52	8.78	26.52	0.51
清远	18.93	3.08	18.34	1.02
东莞	679.87	128.34	673.82	51.2
中山	414.64	77.25	390.16	25.87
潮州	40.98	8.21	40.19	1.94
揭阳	26.34	6.88	25.05	0.41
云浮	18.77	5.85	16.08	0.64

library(openxlsx)
T2=read.xlsx("E:/多元统计分析及R语言建模/mvexer5.xlsx",sheet="E3.2",rowNames=TRUE)
head(T2)
msa.X=function(df){
  X=df[,-1]
  rownames(X)=df[,1]
  X=return(X)
}
t2=msa.X(T2)
summary(t2)
attach(t2)
rm=apply(t2,1,mean);rm
barplot(rm,las=3)#按行做均值图
cm=apply(t2,2,mean);cm
barplot(cm,horiz = T,las=1)#按列做均值图
barplot(apply(t2,1,median))#按行做中位数条图
barplot(apply(t2,2,median),col=1:4)#按列取色
boxplot(t2)#按列做箱尾图
boxplot(t2,col=c("red","blue","green","pink"),notch=T)
stars(t2,key.loc=c(8,1))#全角星相图
stars(t2,full=F,key.loc=c(8,1))#半角星相图
stars(t2,full=T,draw.segments = T,key.loc = c(10,1))#扇形星相图
stars(t2,full=F,draw.segments = T,key.loc = c(10,1))#半圆星相图
library(aplpack)
faces(t2,ncol.plot=6)#彩色脸谱图
faces(t2[2:4],ncol.plot=6)#去掉第一个变量的彩色脸谱图
library("TeachingDemos")
faces2(t2,ncol.plot=6)#黑白脸谱图
library(andrews)
andrews(t2)#调和曲线