lineage.m

%% lineage
% gets list of taxa to which a taxon belongs

%%
function classification = lineage(taxon)
% created 2015/09/18 by Bernd Brandt

%% Syntax
% classification = <../lineage.m *lineage*> (taxon) 

%% Description
% gets all taxa in the add_my_pet collection to which a particular taxon belongs
%
% Input:
%
% * character string with name of taxon
%
% Output:
% 
% * cell string of ordered taxa to which that taxon belongs, starting with Animalia

%% Remarks
% The root is Animalia. 
% The classification follows that of Wikipedia

%% Example of use
% classification  = lineage('Gorilla_gorilla')

  WD = pwd;                        % store current path
  taxa = which('lineage.pl');      % locate DEBtool_M/taxa/
  taxa = taxa(1:end - 10);         % path to DEBtool_M/taxa/
  cd(taxa)                         % goto taxa

  try
    classification = textscan(perl('lineage.pl', taxon), '%s'); 
    classification = classification{1};
  catch
    disp(['Warning from lineage: Name ', taxon, ' is not recognized as taxon'])
  end
  
  cd(WD)                           % goto original path

end