WTSI Mutational Signature Framework
23 Oct 2012
19 Jul 2013)
Framework for deciphering signatures of mutational processes from a set of mutational catalogues
function [genomes mutationTypesToRemoveSet] = removeWeak(genomesOriginal, removeWeakMutationTypes)
% Ludmil B. Alexandrov
% Cancer Genome Project
% Wellcome Trust Sanger Institute
% This software and its documentation are copyright 2012 by the
% Wellcome Trust Sanger Institute/Genome Research Limited. All rights are reserved.
% This software is supplied without any warranty or guaranteed support whatsoever.
% Neither the Wellcome Trust Sanger Institute nor Genome Research Limited
% is responsible for its use, misuse, or functionality.
[sortedTotalMutationTypes indexTotalMutationTypes] = sort(sum(genomesOriginal, 2), 'ascend');
totalMutationTypes = sum(sum(genomesOriginal));
totalMutatationTypesToRemove = sum(cumsum( sortedTotalMutationTypes )/ totalMutationTypes < removeWeakMutationTypes);
mutationTypesToRemoveSet = indexTotalMutationTypes(1:totalMutatationTypesToRemove);
genomes = genomesOriginal;
genomes(mutationTypesToRemoveSet, :) = ;