作业帮MATLAB 多元回归 t检验
来源：学生作业帮 编辑：作业帮 分类：数学作业 时间：2024/05/27 05:08:02
MATLAB 多元回归 t检验
自变量是x1-x4,因变量y
改了一点程序向用还是用不了,求大神看下.
谁再能解释下程序的意思.
z=[x1,x2,x3,x41,x42];y=A(:,1); % example
[n,k]=size(z);
X=[ones(n,1),z];%构建结构阵X
B=X'*X; %求算信息阵A
C=inv(B); %求算信息阵的逆阵
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
自变量是x1-x4,因变量y
改了一点程序向用还是用不了,求大神看下.
谁再能解释下程序的意思.
z=[x1,x2,x3,x41,x42];y=A(:,1); % example
[n,k]=size(z);
X=[ones(n,1),z];%构建结构阵X
B=X'*X; %求算信息阵A
C=inv(B); %求算信息阵的逆阵
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
没问题,说说你哪里不明白?
再问： Warning: Matrix is singular to working precision.
Warning: Rank deficient, rank = 4, tol = 2.3232e-009.
SSy =
5.2308e+003
R2 =
NaN
为什么结果是这样？
再答： 你显示一下RSS，看你RSS有没有结果
再问： RSS =
NaN
SSy =
5.2308e+003
R2 =
NaN
我的数据跑回归是没有问题的，所以应该是t检验程序不对，请问您能看出来吗？
再答： 很明显，是你RSS 出问题了。补充一下，matlab做多元回归及检验有内置的命令，没有必要亲自写。
再问： MATLAB的T检验函数不是只能检验单一变量的吗？那多元回归改怎么办？
再答： 多元变量显著性检验用F统计量。建议好好看一下统计或者计量经济学的基础教材，欲速则不达，任何人都是从你现在的阶段过来的
再问： Warning: Matrix is singular to working precision.
Warning: Rank deficient, rank = 4, tol = 2.3232e-009.
SSy =
5.2308e+003
R2 =
NaN
为什么结果是这样？
再答： 你显示一下RSS，看你RSS有没有结果
再问： RSS =
NaN
SSy =
5.2308e+003
R2 =
NaN
我的数据跑回归是没有问题的，所以应该是t检验程序不对，请问您能看出来吗？
再答： 很明显，是你RSS 出问题了。补充一下，matlab做多元回归及检验有内置的命令，没有必要亲自写。
再问： MATLAB的T检验函数不是只能检验单一变量的吗？那多元回归改怎么办？
再答： 多元变量显著性检验用F统计量。建议好好看一下统计或者计量经济学的基础教材，欲速则不达，任何人都是从你现在的阶段过来的