part_index.m

数据挖掘算法中空间同位规则挖掘的matlab源码

function[p] = part_index(sample_size_2, sample, thrd)

p = zeros(45, 5);

index = 1;
for itype = 0:9
	for jtype = itype+1:9
    		ab = find(sample_size_2(:,2) == itype & sample_size_2(:,6) == jtype & sample_size_2(:,9)<thrd);
    		nor = nnz(unique(sample_size_2(ab,1)));
    		denorm = nnz(unique(find(sample(:,2) == itype)));
                p1 = nor/denorm;
                   
                nor = nnz(unique(sample_size_2(ab,5)));
                denorm = nnz(unique(find(sample(:,2) == jtype)));
                p2 = nor/denorm;
    
    		p (index, 2) = itype;
    		p (index, 3) = jtype;
    		p (index, 4) = p1;
    		p (index, 5) = p2;
    
    		if (p(index,4)>p(index,5))
    			p(index,1) = p(index,5);
    		else
    			p(index,1) = p(index,4);
    		end;
    		index = index +1;
	end
end;