setwd('C:\\Users\\LJL\\Desktop\\统计软件\\data')
cq<-read.csv(file.choose(),encoding = "UTF-8")
D<-cq[1:58,3:15]
D<-as.matrix(D) #转化为矩
D<-scale(D) #标准化数据
mean_D<-colMeans(D)#计算每列均值
cov_D<-cov(D)#协方差
cor_D<-cor(D)#相关系数
mean_D
cov_D
cor_D
install.packages('psych')
library('psych')
describe(D)
corr.test(D,use='complete',adjust = 'none')#皮尔逊相关性检验
print(corr.test(D,adjust = 'none'),short = F)
#轮廓图
outline<-function(x, txt=TRUE){if (is.data.frame(x)==TRUE)#轮廓图
x<-as.matrix(x);m<-nrow(x); n<-ncol(x)
plot(c(1,n), c(min(x),max(x)), type="n", main="数据的轮廓图",xlab="变量",ylab="取值")
for(i in 1:m){lines(x[i,], col=i)
if (txt==TRUE){k<-dimnames(x)[[1]][i]
text(1+(i-1)%%n, x[i,1+(i-1)%%n], k)}}}
outline(D)
windows()
#星图(每个样本单列)
x<-as.character(cq[,2]) #股票名称
stars(D,draw.segments = T,col.segments = rainbow(7),main = '重庆板块股票星图',frame.plot=T,labels=x,flip.labels=F,family="Courier")
#半幅星图
y<-as.character(cq[1:20,2]) #股票名称
stars(D[1:20,],full=FALSE, draw.segments = TRUE,key.loc = c(5,0.5), mar = c(2,0,0,0),labels=y)#半幅星图
#雷达图或蛛网图
stars(D, locations = c(0, 0), scale = TRUE, radius = FALSE, key.loc = c(0, 0), main = "重庆股票雷达图")
#脸谱图
library(aplpack)
library(tcltk)
faces(D,face.type=2,main='重庆板块股票脸谱图',labels=x)
#调和曲线图
unison<-function(x){ if (is.data.frame(x)==TRUE)
x<-as.matrix(x); t<-seq(-pi, pi, pi/30); m<-nrow(x); n<-ncol(x); f<-array(0, c(m,length(t)))
for(i in 1:m){f[i,]<-x[i,1]/sqrt(2)
for( j in 2:n){ if (j%%2==0) f[i,]<-f[i,]+x[i,j]*sin(j/2*t)
else f[i,]<-f[i,]+x[i,j]*cos(j%/%2*t) } }
plot(c(-pi,pi), c(min(f),max(f)), type="n", main="三角多项式图", xlab="t", ylab="f(t)")
for(i in 1:m) lines(t, f[i,] , col=i)}
unison(D)
# 气泡图
D1<-cq[1:58,3:17]
within(D1,{
r<-sqrt(D1$换手./pi) #气泡半径
symbols(D1$总量,D1$总金额,circle=r,inches=0.5,fg="white",bg="lightslateblue",main="气泡图(点的大小代表换手率)",xlab="总量",ylab="总金额")
text(D1$总量,D1$总金额,D1[,15],cex=0.6,pos=1)
}
)
##相关系数矩阵图
install.packages("PerformanceAnalytics")
library(PerformanceAnalytics)#加载包
chart.Correlation(D, histogram=TRUE, pch=36)
#词云图(总金额)
library(wordcloud2)
c<-cq[,c(2,17)]
wordcloud2(c,size=2,fontFamily = '微软雅黑',color = 'random-light',minRotation = -pi/6,maxRotation = -pi/6)
星图绘图结果如上,中文字体出现乱码,请问该如何解决?