数学建模社区-数学中国

标题: matlab做回归分析时，能做t检验吗？ [打印本页]

作者: Reader_Founder 时间: 2013-8-17 21:04
标题: matlab做回归分析时，能做t检验吗？
求高人指教，谢谢~
在线等

作者: cyshumo 时间: 2013-8-17 22:42
进行多元线性回归统计数F, t 测验的小程序：
clear,clc
x=rand(50,10);y=rand(50,1); % example
[n,k]=size(x);
X=[ones(n,1),x];%构建结构阵X，
A=X'*X; %求算信息阵A，
C=inv(A); %求算信息阵的逆阵，
b=X\y, % 求算回归统计数向量，其中第一行为回归截距a,
RSS=y'*y-b'*X'*y, %求算离回归平方和，
MSe=RSS/(n-k-1),%求算离回归方差，
Up=b.*b./diag(C);%求算偏回归平方和，其中第一行是a与0差异的偏平方和，
F=Up/MSe,%F测验，其中第一行为a与0差异的F值，
sb=sqrt(MSe*diag(C)); %求算回归统计数标准误，
t=b./sb, % 回归统计数的 t 测验，其中第一行为a与0差异的t测验值。
[t, t.^2, F],%验证t^2=F
SSy=var(y)*(n-1)
R2=(SSy-RSS)/SSy
顺便说一下，你的ttest(x,m)的 t 测验指的是单个样本（平均数）与 m 之间差异显著性的 t 测验，而非多元线性回归系数的 t 测验。
这是一段编程，不知有没有用。

作者: Reader_Founder 时间: 2013-8-17 23:44

cyshumo 发表于 2013-8-17 22:42
进行多元线性回归统计数F, t 测验的小程序：
clear,clc
x=rand(50,10);y=rand(50,1); % example

真心感谢啊~

作者: Reader_Founder 时间: 2013-8-17 23:45

lvyanlong 发表于 2013-8-17 22:07
可以的，http://pan.baidu.com/share/link? ... 7&uk=3090271520

可是访问出错了

作者: Reader_Founder 时间: 2013-8-17 23:49

cyshumo 发表于 2013-8-17 22:42
进行多元线性回归统计数F, t 测验的小程序：
clear,clc
x=rand(50,10);y=rand(50,1); % example

那怎么做多元回归系数的显著性检验呢

作者: lvyanlong 时间: 2013-8-18 08:07
用h=ttest(x,m,0.05)语句进行t检验，结果h值为1，拒绝原假设，没通过检验。

作者: magic2728 时间: 2013-8-18 12:38
t检验可以用在回归分析的检验中，和是不是matlab无关。。。

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)