Skip to main content

Feature Selection


1. Matlab using TreeBagger (it is actually like Random Forest)
==========================================

load ionosphere;

noBag = 5;
myBag  = TreeBagger( noBag  ,   X, Y, 'OOBPred','on' , 'oobvarimp' ,'on' );


//  increase in prediction error if the values of that variable are permuted across OOB observations.
//  The more increase in prediction Error ==> The more important the variable is
oobVarImp = myBag.OOBPermutedVarDeltaError


// re-substitution error
varImp = zeros( noBag, noFeature)
for i=1:noBag
   varimportance( myBag.Trees{i})
end

========== COMPLETE CODE==========

function fromRF

load ionosphere;
noBag = 5;
myBag  = TreeBagger( noBag  ,   X, Y, 'OOBPred','on');
varRanking = zeros( noBag  , size(X,2) ) ;

for i=1:noBag   
   [ val ,varRanking( i , :) ]= sort( varimportance( myBag.Trees{i}) ,'descend')   
end

// suppose finally taking top ranked 5 from all folds
topRank=5;
selectedFeat=[]
for i=1:noBag
   selectedFeat = union( selectedFeat , varRanking( i , 1:topRank) ); 
end

display('done');

end


Comments

Popular posts from this blog

MATLAB cross validation

// use built-in function samplesize = size( matrix , 1); c = cvpartition(samplesize,  'kfold' , k); % return the indexes on each fold ///// output in matlab console K-fold cross validation partition              N: 10    NumTestSets: 4      TrainSize: 8  7  7  8       TestSize: 2  3  3  2 ////////////////////// for i=1 : k    trainIdxs = find(training(c,i) ); %training(c,i);  // 1 means in train , 0 means in test    testInxs  = find(test(c,i)       ); % test(c,i);       // 1 means in test , 0 means in train    trainMatrix = matrix (  matrix(trainIdxs ), : );    testMatrix  = matrix (  matrix(testIdxs  ), : ); end //// now calculate performance %%  calculate performance of a partiti...

MATLAB confusion matrix

%  test_class  & predicted_class must be same dimension % 'order' - describes the order of label. Here labels are 'g' as positive and 'h' as negative [C,order] = confusionmat( test_class(1: noSampleTest), predicted_class, 'order', ['g' ;'h'] ) tp = C(1,1); fn = C(1,2); fp = C(2,1); tn = C(2,2); sensitivity = tp /( tp + fn ) specificity = tn /( fp + tn ) accuracy = (tp+tn) / (tp+fn+fp+tn) tpr = sensitivity fpr = 1-specificity precision = tp /( tp + fp ) fVal = (2*tpr*precision)/(tpr+precision)

Feature subset selection Using Genetic Algorithm in MATLAB

function callGeneticAlgo global mat global trainInd global testInd [trainInd,~,testInd] = dividerand(1420,0.7,0,0.3); global counter global errList counter = 1; errList = []; fileName=  '../features/alltopPNPDMF.feature' ; mat = load(fileName); [x,fval,exitflag,output,population,score] = gaFeaSelection(1588,100,10800); % param1 = #feature excludig label % param2 =  population size % param3 = sec to test (3 hour = 10800 sec) dlmwrite('selected.GA',x,'delimiter','\n'); display('Done'); end function [x,fval,exitflag,output,population,score] = gaFeaSelection (nvars,PopulationSize_Data,TimeLimit_Data) % This is an auto generated MATLAB file from Optimization Tool. % Start with the default options options = gaoptimset; % Modify options setting options = gaoptimset(options,'PopulationType', 'bitString'); options = gaoptimset(options,'PopulationSize', PopulationSize_Data); options = gaoptimset(options,'TimeLimit', T...