0
点赞
收藏
分享

微信扫一扫

(R语言)离差阵 和 一些图的画法

戴老师成长记录仪 2022-02-27 阅读 92

一、求均值向量,离差阵,协方差阵,相关阵,标准化后的协方差矩阵

原始数据长这样:

城市X1X2X3X4X5
北京4295638276
天津8396848245
石家庄124110158205
太原17459554224
呼和浩特13337140294
沈阳18357442287
长春10325942305
哈尔滨17326447303
上海6374132319
南京7365631304
杭州6385530334
合肥7395636311
福州5213821364
南昌9295833335
济南13368850223
郑州9398451230
武汉8365837309
长沙7274841309
广州7364323331
南宁8244626357
海口4112914361
重庆8395333331
成都6376441280
贵阳10184123362
昆明9264224366
拉萨7192912366
西安8418851250
兰州15477634312
西宁15366135337
银川14357236301
乌鲁木齐9367547279

X= read.csv("C:\\Users\\HuXiang\\Documents\\Tencent Files\\1085244341\\FileRecv\\Ch01_kongqizhiliang2020.csv",header=T,sep=",")
X<-data.frame(X[,2:6]) #转换为数据表
Mean_X=colMeans(X,na.rm=T) #按照列求均值向量
Dim_X=dim(X)  #求数据表的维数
Scatter_X = (Dim_X[1]-1)*cov(X) #求离差阵
Covariance_X=cov(X) #求协方差阵
Correlation_X=cor(X) #求相关系数阵
Standard_X<-scale(X,center = T,scale = T)  #样本数据的0-1标准化
Covariance_SX = cov(Standard_X)  #标准化后的样本数据的协方差矩阵

二、绘制条形图、矩阵散点图,箱型图,雷达图,星图,轮廓图,脸谱图和调和曲线图

还是使用刚才的数据

X= read.csv("C:\\Users\\HuXiang\\Documents\\Tencent Files\\1085244341\\FileRecv\\Ch01_kongqizhiliang2020.csv")
barplot(apply(X[,2:6],1,mean),names.arg=X[,1],las=3)#按行绘制均值条形图
barplot(apply(X[,2:6],1,mean),col=1:6)#按列绘制彩色均值条形图
pairs(X[,2:6]) #绘制矩阵散点图
boxplot(X[,2:6]) #按列绘制箱型图,垂直放置
boxplot(X[,2:6],horizontal = T) #按列绘制箱型图,水平放置
boxplot(X[,2:6],horizontal = T,las=1)
#绘制脸谱图
install.packages("aplpack",depend=TRUE) #安装aplpack包
library("aplpack")                      #加载aolpack包

a = faces(data.matrix(X),fill=T,which =2:6,labels = X[,1],cex=0.7 )
plot(a,face.type=1,cex=0.7)
#绘制雷达图
stars(X[,2:6],radius=T,labels = X[,1],cex=0.7)
#绘制星图
stars(X[,2:6],radius=F,labels = X[,1],cex=0.7)
#绘制轮廓图
matplot(0:30,X[,2:6],xlim=c(0.5,30),type="l",lty = 5:1,col=1:6,las=2,
        xlab="",cex.axis=0.7,cex.lab=0.7,ylab = "平均浓度或天数")
axis(1,0:30,t(X[,1]),las=2,cex.axis=0.7)
for(i in 1:5)lines(25:28,rep(230-10*(i-1),4),lty=8-i,col=i)
text(29,230-10*(i-1),colnames(X)[i+1],cex=0.7)
#绘制调和曲线图
install.packages("andrews")
library(andrews)
andrews(X[,2:6],clr=3,ymax = 3)




举报

相关推荐

0 条评论