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) 错误信息已经提示了,x_farm的第四列已经被剔出了,可是后面的运算还要访问第四列
回复 沙发 nemia 的帖子
可以考虑把x_farm赋给一个中间变量
页:
[1]