楼主: xsx小虾米
7727 1

[经济] matlab 回归分析中 regstats怎么用? [推广有奖]

  • 5关注
  • 2粉丝

硕士生

42%

还不是VIP/贵宾

-

威望
0
论坛币
190 个
通用积分
1.6644
学术水平
4 点
热心指数
4 点
信用等级
4 点
经验
2881 点
帖子
108
精华
0
在线时间
107 小时
注册时间
2012-4-4
最后登录
2021-4-25

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
A = xlsread('D:\test.xls')
x = unique(A(:,1));
n = length(x);
for i =1:n   %回归次数
      B=A(A(:,1)==x(i),4);
      C=A(A(:,1)==x(i),2);
      D=A(A(:,1)==x(i),3);
      X=ones(length(B),1);
     [b,bint,r,rint,stats] = regress(B,[X,C,D]);


end
怎么写入求t值得代码?搜了下要用regstats,但是不会用。。。提示用法如下,但是运行后不对。。
    whichstats = {'tstat'};
     stats = regstats(B,[C,D],'linear',whichstats);
     stats= stats.tstat
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:MATLAB stats matla atlab 回归分析 回归分析 matlab

沙发
chelsea4 发表于 2014-11-23 17:43:59 |只看作者 |坛友微信交流群
MATLAB统计工具箱中提供了regstats函数,也可用来作多重线性或广义线性回归分析,它的调用方式如下:
regstats(y,X,model)
stats = regstats(…)
stats = regstats(y,X,model,whichstats)
(1)regstats(y,X,model)
作多重线性回归分析。输入参数X为自变量观测值矩阵(或设计矩阵),它是的矩阵。默认情况下,regstats函数自动在X第1列元素的左边加入一列1,不需要用户自己添加。输入参数y为因变量的观测值向量,是的列向量。可选的输入参数model是一个字符串,用来控制回归模型的类型,其可用的字符串如表1-2所示。
表1-2 regstats函数支持的model参数
model参数的参数值
说 明
'linear'
带有常数项的线性模型(默认情况)
'interaction'
带有常数项、线性项和交叉项的模型
'quadratic'
带有常数项、线性项、交叉项和平方项的模型
'purequadratic'
带有常数项、线性项和平方项的模型
在这种调用方式下,regstats函数会生成一个交互式图形用户界面(GUI),界面上带有回归诊断统计量列表,包括系数的估计值、因变量的预测值、残差、判定系数、调整的判定系数、F检验和t检验的相关结果等,共23个可选项。通过这个界面,用户可以很方便地将回归分析的各种结果导入MATLAB工作空间。
(2)stats = regstats(…)
返回一个结构体变量stats,它有24个字段,包括了回归分析的所有诊断统计量。这种调用方式不生成图形用户界面,stats的后23个字段分别与图形用户界面上的23个选项相对应。
(3)stats = regstats(y,X,model,whichstats)
仅返回由whichstats参数指定的统计量。whichstats可以是形如 'leverage' 的单个字符串,也可以是形如{'leverage' 'standres' 'studres'}的字符串元胞数组。若whichstats是字符串 'all' ,则返回所有统计量。
注意:当需要计算F统计量的观测值时,模型中应包含常数项。若模型中不包含常数项,regstats函数输出的F统计量的观测值是不正确的。在不考虑常数项的情况下,计算出的判定系数2R的值可能是负的,说明所用模型不适合用户的数据。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加JingGuanBbs
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-10-6 18:05