https://github.com/tansey/smoothfdr
Tip revision: c5b693d0a66e83c9387433b33c0eab481bd4a763 authored by Wesley Tansey on 08 May 2020, 15:42:20 UTC
Fixed bug in easy that created too large a support for the alternative distribution
Fixed bug in easy that created too large a support for the alternative distribution
Tip revision: c5b693d
hmrfRun.m
function []=hmrfRun(expdir, L)
% L = str2num(Lstr);
input_info = [3, 130, 130, 0.1];
reglocafile = strcat(expdir, 'bufferregions.csv');
pvaluefile = strcat(expdir, 'bufferpvalues.csv');
signalfile = strcat(expdir, 'buffer_signals_bh_hmrf.csv');
zvaluefile = strcat(expdir, 'bufferdata.csv');
output_estimate = strcat(expdir, 'buffer_estimate_result.csv');
output_LIS = strcat(expdir, 'buffer_lis.csv');
output_process = strcat(expdir, 'buffer_process_output.csv');
result=fdrBH95(input_info,reglocafile,pvaluefile,signalfile);
initials = initialvalue(signalfile, zvaluefile, reglocafile, L);
input_info = [3, 130, 130, L, 1000, 5000, 1e6, 5000];
computelis(input_info, initials, zvaluefile, reglocafile, output_estimate, output_LIS, output_process);
end