1.题目

使用教材P84页表3-2，进行：
1、多元数据简单R分析：定量变量的分析（直方图、散点图）、定性变量的分析（单因素分析、多维列联表）。参考教材P45-P52页。
2、绘制均值条图、箱尾图、星相图、调和曲线图。参考教材P65-P73页，P77- P78页。

2.题目详解

2.1.多元数据简单R分析：定量变量的分析（直方图、散点图）、定性变量的分析（单因素分析、多维列联表）。参考教材P45-P52页。

读数据：

library(xlsx) #导包
data = read.xlsx("C:\\Users\\stu\\Downloads\\广东省各市高新技术产品情况.xlsx",'Sheet1',row.names=T) #读数据
data

可查看到读取的数据：

进行定量变量分析（该步骤较简单，先上代码后上图）

hist(data$工业总产值) #直方图
plot(data$工业总产值,data$工业增加值) #散点图

运行结果如下：


定性变量分析（单因素）：
步骤：
1、对工业总产值进行分组（我这里是按300/组，大家可以自行分组如150/组，200/组）
2、查看分组是否符合类别
3、绘制条形图

data$group <- cut(data$工业总产值, breaks = c(-Inf, 300, 600, 900, 1200, Inf), labels = c("0-300","300-600","600-900","900-1200","1200+"), right=FALSE) #对工业总产值进行分组
table(data$group)

绘制条形图：

barplot(table(data$group),col=1:5) #绘制条形图

绘制饼图：

pie(table(data$group)) #绘制饼图

两因素分析：
注意：该分组仅根据个人看法分组，大家可以自行选择其余特征进行分组，如产品销售收入过五百亿跟未过五百亿等，我这里按照工业增加值是否过百亿进行分组，分出两组即可

data$group1 <- cut(data$工业增加值, breaks = c(-Inf, 100, Inf), labels = c("未过百亿","过百亿"), right=FALSE) #对工业总产值进行分组
table(data$group,data$group1)

以工业总产值（group）及工业增加值分别分组绘制条形图：

barplot(table(data$group,data$group1),beside=T,col=1:5) #工业总产
barplot(table(data$group1,data$group),beside=T,col=1:2) #工业增加

2.2.绘制均值条图、箱尾图、星相图、调和曲线图。参考教材P65-P73页，P77- P78页。

均值条图：

data = read.xlsx("C:\\Users\\6\\Desktop\\R语言\\实验4\\广东省各市高新技术产品情况.xlsx",'Sheet1',row.names=T) #读数据
barplot(apply(data,1,mean))
barplot(apply(data,1,mean),las=3)

箱尾图：
按列做箱尾图：

boxplot(data) #按列作
boxplot(data,horizontal = T) #按水平做

按列作

按水平做：

星相图：（每一行对应每一个图）

stars(data,full = T,key.loc = c(10,2))
stars(data,full = F,key.loc = c(10,2))
stars(data,full = T,draw.segments = T,key.loc = c(10,2))

星象图：

半角星象图：

扇形星象图：

一般调和曲线：

install.packages("andrews")
library(andrews)
andrews(data,clr = 4,ymax=4)

改进调和曲线（msaR.R包是该书作者自己撰写的包，如果有人需要的话可以在底下留言，我把百度网盘发出来）：

source("C:\\Users\\6\\Desktop\\R语言\\实验2\\msaR.R")
msa.andrews(data)