lanxiaod 发表于 2009-8-3 11:50

if语句的代码与求出结果不符合,请问该如何解决和修改?

if语句的代码与求出结果不符合,请问该如何解决和修改?
x_farm有4个因子(列);4个对象(行),我想对x_farm进行因子剔除,通过计算相关系数和方差,要求两者(i,j)相关系数大于0.9,若i 的方差大于j 的方差,就剔除j 因子。我的代码如下
clear
x_farm=
cov=var(x_farm) %计算方差
C=corrcoef(x_farm)%计算相关系数
for i=1:length(x_farm)
for j=i:length(x_farm)
if C(j,i)>=0.9 && cov(i)>cov(j)
x_farm(:,j)=[]
else
x_farm(:,j)=x_farm(:,j)
end
end
end
结果如下:
x_farm =
167.6200 100.4584 100.1472 160.3753
171.3246 93.5399 97.7215 148.9536
159.2697 88.2294 92.2526 142.8844
177.3102 102.1516 102.4093 163.5383

cov =
56.9966 41.3726 19.0307 93.5436

C =
1.0000 0.7887 0.9010 0.7617
0.7887 1.0000 0.9740 0.9961
0.9010 0.9740 1.0000 0.9565
0.7617 0.9961 0.9565 1.0000
x_farm =
167.6200 100.4584 100.1472 160.3753
171.3246 93.5399 97.7215 148.9536
159.2697 88.2294 92.2526 142.8844
177.3102 102.1516 102.4093 163.5383

x_farm =
167.6200 100.4584 100.1472 160.3753
171.3246 93.5399 97.7215 148.9536
159.2697 88.2294 92.2526 142.8844
177.3102 102.1516 102.4093 163.5383

x_farm =
167.6200 100.4584 160.3753
171.3246 93.5399 148.9536
159.2697 88.2294 142.8844
177.3102 102.1516 163.5383
??? Attempted to access x_farm(:,4); index out of bounds because size(x_farm)=.
Error in ==> test at 10
x_farm(:,j)=x_farm(:,j)

nemia 发表于 2009-8-3 18:34

错误信息已经提示了,x_farm的第四列已经被剔出了,可是后面的运算还要访问第四列

friendchj 发表于 2009-8-5 14:45

回复 沙发 nemia 的帖子

可以考虑把x_farm赋给一个中间变量
页: [1]
查看完整版本: if语句的代码与求出结果不符合,请问该如何解决和修改?