用regress函数求一元回归，出现warning,求大神解答

Alliris

所用程序如下：
load('mrso_bcc-csm1-1_2525.mat');
z1=reshape(z1,73,144,163)
for loni=1:144
    for lati=1:73
        Y=squeeze(z1(loni,lati,:));
        X=[ones(163,1) (1:163)'];
        [b,bint ,r,rint ,stats]=regress(Y,X);
        slope(loni,lati)=b(2);
        Rsqu(loni,lati) =stats(1);
    end   
end

出现错误提示：
In regress at 84
Warning: X is rank deficient to within machine precision.

不知道是什么原因造成的，求解答

风儿飘飘

用regress做一元线性回归，原始数据要增加一列为1的常数项，你的资料满足条件？

Alliris

风儿飘飘 发表于 2014-6-10 14:59
用regress做一元线性回归，原始数据要增加一列为1的常数项，你的资料满足条件？

就是Y需要增加一列为1的常数项吗？不太明白意思哎，在第三维增加？

风儿飘飘

Alliris 发表于 2014-6-10 15:56
就是Y需要增加一列为1的常数项吗？不太明白意思哎，在第三维增加？

实际上就是增加1列，放在第1列的位置

Alliris

风儿飘飘 发表于 2014-6-12 11:06
实际上就是增加1列，放在第1列的位置

那我这是三维数据的话，是在每一维都需要增加一列？

风儿飘飘

你的数据要说明白，三维分别代表什么？

Alliris

风儿飘飘 发表于 2014-6-12 14:36
你的数据要说明白，三维分别代表什么？

不好意思没说清楚，三维分别是经度，纬度，时间

风儿飘飘

Alliris 发表于 2014-6-12 14:46
不好意思没说清楚，三维分别是经度，纬度，时间

如果你的一元线性回归方程需要常数项时，自变量X才需要加入一列为1的常数项！

风儿飘飘

X=[ones(163,1) (1:163)'];是不是应该修改为X=[ones(163,1)， (1:163)'];

Alliris

风儿飘飘 发表于 2014-6-12 21:10
X=[ones(163,1) (1:163)'];是不是应该修改为X=[ones(163,1)， (1:163)'];

试了一下不是这个问题，还是非常感谢！