使用MATLAB的fitlm函数进行线性回归

今天在做《数理统计》关于线性回归的作业,本来用R已经做出来了,但是由于最近使用matlab很多,所以也想看看用matlab怎么做。

matlab中有很多函数可以做各种各样的回归,也有cftool工具箱可以可视化的做回归,很方便。这里选用fitlm做回归,由于多元回归和一元回归基本思想是差不多的,操作也只是参数个数的问题,所以这里用一元线性回归做例子,记录下来以备后用。

数据选用R中的自带数据:cars数据集,是一个关于汽车速度和距离的数据,50*2的矩阵。

采用一元线性回归模型进行回归,公式这里就不说了,dist为因变量,speed为自变量。 代码来了:

clear
data = xlsread('D:\Program Files\RStudio\cars.xlsx');
x = data(:, 1);
y = data(:, 2);
tb = table(x, y, 'VariableNames', {'speed', 'dist'});
model = fitlm(tb, 'dist~speed');
plot(model);

model里含有模型的各种参数,估计值,R2值,p值等等等等。

model = 


Linear regression model:
    dist ~ 1 + speed

Estimated Coefficients:
                   Estimate      SE        tStat       pValue  
                   ________    _______    _______    __________

    (Intercept)    -17.579      6.7584    -2.6011      0.012319
    speed           3.9324     0.41551      9.464    1.4898e-12


Number of observations: 50, Error degrees of freedom: 48
Root Mean Squared Error: 15.4
R-squared: 0.651,  Adjusted R-Squared 0.644
F-statistic vs. constant model: 89.6, p-value = 1.49e-12

可以看出输入形式和输出形式与R的输出及其相似。

最后plot画出回归图。

好了,该吃饭去了。

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏程序生活

斯坦福tensorflow教程-tensorflow 实现线性回归代码结果

1182
来自专栏书山有路勤为径

Creating a Filter, Edge Detection

Below, you've been given one common type of edge detection filter: a Sobel opera...

811
来自专栏Python小屋

Python图像处理模块pillow子模块Image用法精要

Image是pillow库中一个非常重要的模块,提供了大量用于图像处理的方法。使用该模块时,首先需要导入。 >>> from PIL import Image ...

3354
来自专栏生信小驿站

R 相关性分析

1864
来自专栏数据处理

图片RGB三维矩阵的表示

5951
来自专栏数值分析与有限元编程

ANSYS里的对称与反对称约束

首先回顾一下结构力学里的概念:在平面内绕对称轴旋转180度,荷载的作用点重合,作用方向相反便是反对称荷载,如果荷载的作用点重合,作用方向相同,便是正对称荷载。通...

3234
来自专栏Petrichor的专栏

图像处理: 无损地旋转图像

2383
来自专栏Python

随机验证码

Python生成随机验证码,需要使用PIL模块.python3则是pillow 安装: ? 1 pip3 install pillow ...

2626
来自专栏杂七杂八

Matplotlib 绘3D图

Matplotlib 也可以绘制 3D 图像,与二维图像不同的是,绘制三维图像主要通过 mplot3d 模块实现。但是,使用 Matplotlib 绘制三维图像...

3497
来自专栏专知

【论文推荐】最新5篇信息抽取(IE)相关论文—开放信息抽取、不完整信息、主动学习、越南语、依存分析

【导读】专知内容组整理了最近五篇信息抽取(Information Extraction)相关文章,为大家进行介绍,欢迎查看! 1. Assertion-base...

4125

扫码关注云+社区