Matlab结合ChatGPT—如何计算置信区间?

​前面分享了带置信区间的折线图和带置信区间的折线散点图的绘图教程:

很多人表示,昆哥,图是很好看啦,但咱不会求置信区间啊,咋办嘞?

说实话,这种事情属于数据处理,一般都是在画图前解决的,所以我在教程中是不会体现的。

另一方面,数据处理工作涉及面太广,很多我自己用不到,也就不会深入研究,更谈不上出干货教程(但不影响画图)。

不过,最近这不学会ChatGPT了嘛,搞不定找他呗。

Matlab&ChatGPT计算置信区间

首先,为了防止他随便编个代码骗我,我从《Nature》论文中找到一幅带置信区间的折线散点图

利用之前分享的图片数据提取工具,提取其中的离散点,作为原始数据,而图片本身作为验证信息,检验ChatGPT给出的结果是否正确

然后,向ChatGPT中提问:

这是他给出的代码:

复制代码
% 数据准备data = [0.046053 0.096491 0.14912 0.19518 0.24342 0.30044 0.3443 0.39474 0.44737 0.5 0.54825 0.59649 0.70175 0.74781 0.80044 0.64912 0.85088 1.0044 0.95175 0.90132;         -1.7629 -1.6466 -1.7949 -1.8486 -1.7438 -1.754 -1.7412 -1.662 -1.731 -1.7195 -1.7412 -1.699 -1.6837 -1.7195 -1.7272 -1.4575 -1.4166 -1.3553 -1.4639 -1.4895];x = data(1, :);y = data(2, :);​% 线性拟合[p, S] = polyfit(x, y, 1);​% 计算拟合直线的y值y_fit = polyval(p, x);​% 计算标准误差n = length(y); % 样本数量y_residual = y - y_fit; % 残差SSE = sum(y_residual.^2); % 残差平方和sy = sqrt(SSE / (n - 2)); % 残差标准差​% 计算置信区间alpha = 0.05; % 置信水平 95%t_value = tinv(1 - alpha/2, n - 2); % 学生 t 分布临界值x_mean = mean(x); % x 的均值Sxx = sum((x - x_mean).^2); % x 的方差和delta = t_value * sy * sqrt(1/n + (x - x_mean).^2 / Sxx); % 置信区间​% 可视化figure;plot(x, y, 'bo'); % 原始数据点hold on;plot(x, y_fit, 'r-', 'LineWidth', 2); % 拟合直线plot(x, y_fit + delta, 'r--', 'LineWidth', 1.5); % 上置信区间plot(x, y_fit - delta, 'r--', 'LineWidth', 1.5); % 下置信区间legend('原始数据', '拟合直线', '置信区间', 'Location', 'Best');xlabel('x');ylabel('y');title('带噪声离散点的线性拟合及其置信区间');hold off;

将代码复制粘贴到Matlab,得到结果:

可以看到,除了图片比例之类的细节问题,结果图的整体趋势同验证图没啥区别,说明ChatGPT给的代码没问题。

简单调一下绘图细节:

保险起见,我又问了一下置信区间的计算原理,这是他给出的答案:

找了一下网上相关内容,基本也是大差不差,说明他还是很靠谱的。

以上。

相关推荐
Zhen (Evan) Wang几秒前
(豆包)xgb.XGBRegressor 如何进行参数调优
开发语言·python
我爱一条柴ya5 分钟前
【AI大模型】线性回归:经典算法的深度解析与实战指南
人工智能·python·算法·ai·ai编程
Qiuner11 分钟前
【源力觉醒 创作者计划】开源、易用、强中文:文心一言4.5或是 普通人/非AI程序员 的第一款中文AI?
人工智能·百度·开源·文心一言·gitcode
虾球xz15 分钟前
CppCon 2018 学习:THE MOST VALUABLE VALUES
开发语言·c++·学习
未来之窗软件服务22 分钟前
chrome webdrive异常处理-session not created falled opening key——仙盟创梦IDE
前端·人工智能·chrome·仙盟创梦ide·东方仙盟·数据调式
AI街潜水的八角39 分钟前
深度学习图像分类数据集—蘑菇识别分类
人工智能·深度学习·分类
阿蒙Amon1 小时前
C#扩展方法全解析:给现有类型插上翅膀的魔法
开发语言·c#
飞睿科技1 小时前
乐鑫代理商飞睿科技,2025年AI智能语音助手市场发展趋势与乐鑫芯片解决方案分析
人工智能
许泽宇的技术分享1 小时前
从新闻到知识图谱:用大模型和知识工程“八步成诗”打造科技并购大脑
人工智能·科技·知识图谱
坤坤爱学习2.01 小时前
求医十年,病因不明,ChatGPT:你看起来有基因突变
人工智能·ai·chatgpt·程序员·大模型·ai编程·大模型学