本期介绍了利用R语言ggplot2包绘制核密度估计图。

核密度估计图（Kernel Density Estimation, KDE），是在概率论中用来估计未知的密度函数，属于非参数检验方法之一，由Rosenblatt （1955）和Emanuel Parzen（1962）提出，又名Parzen窗（Parzen window）。分析核密度函数时主要观察其面积，而不是取值。核密度图中纵轴与横轴所围成的面积为1。

1 数据准备

数据输入格式（csv格式）：

2 R包加载及数据导入

#下载包#install.package("ggplot2")#加载包#library(ggplot2)library(reshape2)#数据载入#data  = read.table(file = 'C:/Rdata/jc/density1.csv', sep = ',', header = T) #header=T表示数据中的第一行是列名，如果没有列名就用 header=F# 转换数据成长数据格式（ggplot常用）data = melt(data)            # melt为reshape2的函数head(data)

如基因表达量数据或其他需要log函数处理可用下列代码

#data=log10(data[,2]+0.0001) #加0.0001主要为了计算有0不好计算,2列进行计算；如有2到12列进行计算则用[,2:12]#round(data,2)#data=data.frame(data)#head(data)#write.csv(data,'C:/Rdata/jc/density1.csv') #数据导出#write.table(data,'C:/Rdata/jc/density1.txt')

3 密度图绘制

#基础密度图绘制#p = ggplot(data, aes(x =value))#x轴表示基因表达值，y轴表示频率就不需要指定p + geom_density(color = "black", #线和点的颜色fill = "gray") #填充颜色

图1 基础核密度图

#美化-线条#p + geom_density(aes(color = variable))#按照不同组改变线条颜色

图2 未填充核密度图

#美化-填充#p + geom_density(aes(fill = variable), #按照不同组改变填充颜色alpha=0.5, #调整透明度linetype = 1, # 线条类型1是实线，2是虚线size=0.5    # 线条粗细)  # fill   指填充颜色# color  指线和点的颜色# colour 指图形边界颜色

图3美化+填充核密度图

好了本次分享就到这里。

关注公众号，发送“核估计密度图”获得完整代码

【R语言】-核密度估计图绘制相关推荐

R语言绘图—热图绘制
R语言绘图-热图绘制原创 R语言与医学生 2022-06-24 20:59 发表于广东热图使用颜色来表示二维图中第三个变量的变化和大小.热图能清楚直观的看到数据的变化.今天就简单介绍下R语言中热 ...
【R语言】——火山图绘制
本期介绍利用R语言筛选差异表达基因及绘制火山图. 一.什么是火山图? 火山图(volcano plot)是散点图的一种,它将统计测试中的统计显著性量度(如p-value.FDR)和变化幅度相结合,从而 ...
R语言之时间序列图绘制
目录 01 读取数据 02 赋值变量 03 时间序列值绘制 04 读取包 05 单位根检验,又称平稳性检验 06 自相关图和偏相关图 07 白噪声检验或纯随机性检验 08 模型拟合 09 显著性检验 ...
R语言使用table1包绘制（生成）三线表、使用单变量分列构建三线表、编写自定义三线表结构（将因子变量细粒度化重新构建三线图）、为不同的变量显示不同的统计信息
R语言使用table1包绘制(生成)三线表.使用单变量分列构建三线表.编写自定义三线表结构(将因子变量细粒度化重新构建三线图).为不同的变量显示不同的统计信息目录
三线表是什么？R语言使用table1包绘制（生成）三线表、使用单变量分列构建三线表、编写自定义三线表结构（将因子变量细粒度化重新构建三线图）、编写自定义函数在三线表中添加p值
三线表是什么?R语言使用table1包绘制(生成)三线表.使用单变量分列构建三线表.编写自定义三线表结构(将因子变量细粒度化重新构建三线图).编写自定义函数在三线表中添加p值目录
R语言ggplot2可视化：绘制堆叠的密度图（Stacked Area Chart）
R语言ggplot2可视化:绘制堆叠的密度图(Stacked Area Chart) 目录 R语言ggplot2可视化:绘制堆叠的密度图(Stacked Area Chart) #仿真数据 #绘制堆叠 ...
R语言使用persp函数绘制三维图像实战（3D）：自定义3D图、图像旋转、添加轴标签
R语言使用persp函数绘制三维图像实战(3D):自定义3D图.图像旋转.添加轴标签目录
R语言使用beanplot包绘制豆荚图实战：与箱图、小提琴图对比、绘制不对称的豆荚图
R语言使用beanplot包绘制豆荚图实战:与箱图.小提琴图对比.绘制不对称的豆荚图目录 R语言使用beanplot包绘制豆荚图实战 #包的安装和导入 #豆荚图与箱图
R语言可视化包ggplot2绘制甘特图（gantt chart）实战
R语言可视化包ggplot2绘制甘特图(gantt chart)实战目录 R语言可视化包ggplot2绘制甘特图(gantt chart)实战 #仿真数据

【R语言】-核密度估计图绘制