% reconstruct an experimental setting

% Input
% M is the sensing matrix
% qMeasurement are the measuremented rare allele frequencies, i.e. # reads_from_the_rare_allele divided by the  # total_number_of_reads_from_that_locus

% Since a full end-to-end experimente of such data does not exist we provide a simulated example of 176 individuals with 2 carriers.


clear
load simulatedExpData % the data contains M and qMeasurement

M0 = normalizeRows(M,2)/2;
tau = 0.005*max(abs(M0'*qMeasurement));
fractionalOutput = applyGPSR(qMeasurement,M0,tau);
discreteOutput = modifySolution(fractionalOutput,M0,qMeasurement);

disp(['carriers are: ',num2str(find(discreteOutput'))])