当前位置: 首页> 健康> 知识 > 绘制拟合联合密度分布

绘制拟合联合密度分布

时间:2025/7/11 9:31:30来源:https://blog.csdn.net/qq_42830971/article/details/140690756 浏览次数:0次

绘制拟合联合密度分布


import numpy as np
import matplotlib.pyplot as plt
from sklearn.neighbors import KernelDensity# 生成模拟数据
np.random.seed(42)
mean = [0, 0]
cov = [[1, 0.5], [0.5, 1]]
data = np.random.multivariate_normal(mean, cov, size=200)
a = data[:, 0]
y = data[:, 1]# 创建一个网格用于绘制密度函数
x_grid = np.linspace(a.min(), a.max(), 100)
y_grid = np.linspace(y.min(), y.max(), 100)
X, Y = np.meshgrid(x_grid, y_grid)
positions = np.vstack([X.ravel(), Y.ravel()]).T# 创建KDE模型
kde = KernelDensity(kernel='gaussian', bandwidth=0.2).fit(data)# 计算每个点上的log-density
log_dens = kde.score_samples(positions)
dens = np.reshape(np.exp(log_dens), X.shape)# 绘制等高线图
plt.figure(figsize=(8, 6))
plt.contourf(X, Y, dens, cmap='viridis')
plt.colorbar(label='Density')
plt.scatter(a, y, color='red', label='Data Points')
plt.xlabel('a')
plt.ylabel('y')
plt.title('Joint Density of a and y')
plt.legend()
plt.show()

在这里插入图片描述

关键字:绘制拟合联合密度分布

版权声明:

本网仅为发布的内容提供存储空间,不对发表、转载的内容提供任何形式的保证。凡本网注明“来源:XXX网络”的作品,均转载自其它媒体,著作权归作者所有,商业转载请联系作者获得授权,非商业转载请注明出处。

我们尊重并感谢每一位作者,均已注明文章来源和作者。如因作品内容、版权或其它问题,请及时与我们联系,联系邮箱:809451989@qq.com,投稿邮箱:809451989@qq.com

责任编辑: