### 第4章 数据可视化 ###
## 4.1 用R的基础绘图系统作图
## 4.2 用ggplot2包作图
getwd()
setwd("C:/我的/广西医科大学/专业核心课/R语言课程/5.27")
## 4.1 用R的基础绘图系统作图 ##
# 函数plot() 对于数值型数据 → 散点图
# 分类数据 → 箱线图
# 统计模型 → 相应图形,如:生存分析→生存曲线
#示例数据,表示某病病人对2种药物,5个剂量(dose)水平上的响应情况
dose <- c(20,30,40,45,60)
drugA <- c(16,20,27,40,60)
drugB <- c(15,18,25,31,40)
plot(dose, drugA)
plot(dose, drugA, type = "b") # type默认为"p",得到散点图。type="b",得到点线图
#用低水平作图函数
plot(dose, drugA, type = "b",
xlab = "Dosage", ylab = "Response",
lty = 1, pch = 15)
lines(dose, drugB, type = "b", lty = 2, pch = 17) #添加线
legend("topleft", title = "Drug Type",
legend = c("A","B"),
lty = c(1,2),
pch = c(15,17)) #添加图例。注意lyt和pch要与上面命令相对应
# 直方图和密度曲线图
#用于展示连续型变量分布,用途:探索分布
#数据集anorexia,关于年轻女性厌食症患者体重变化的研究
data(anorexia, package = "MASS") #调用数据集
str(anorexia) #显示对象内部结构
attach(anorexia) #将数据框添加到搜索路径
hist(Prewt) #绘制直方图
plot(density(Prewt)) #绘制密度曲线
hist(Prewt, freq = FALSE, col = "red",
xlab = "体重",
main = "治疗前体重分布直方图",
las = 1)
lines(density(Prewt), col = "blue", lwd = 2) #颜色、线条宽度
rug(Prewt) #添加轴须图
detach(anorexia) #将对象移出搜索