erdos_renyi_dirichlet_ks_test.m

function [alphas_ERdirichlet_ks_test, vals_ERdirichlet_ks_test] = erdos_renyi_dirichlet_ks_test(social, structural, size_)

%ERDOS RENYI POWER GRAPH
sz = 200;

%Run for size_ iterations     
vals_erDirichlet = zeros(1,sz);
alphas_erDirichlet = zeros(1,sz);

for j=1:size_
    social = gen_ER_Dirichlet(social);

    %Normalize the social graph by the maximum value in the structural graph
    social1 = social .* max(max(structural)) / max(max(social)); 

    for i=1:sz
        %Create thresholds
        delta = 0;
        alpha = max(max(structural)) / (i+150);
        alphas_erDirichlet(i) = alpha+alphas_erDirichlet(i);
        less_idx = find(structural <=(alpha - delta));
        more_idx = find(structural > (alpha + delta));
        %Create social values
        weak_social = social1(less_idx);
        strong_social = social1(more_idx);
        %Create PMFs
        v1 = linspace(min(min(social1)), max(max(social1)), 20);
        v2 = linspace(min(min(social1)), max(max(social1)), 20);
        pmf_weak = (hist(weak_social, v1));
        pmf_weak = pmf_weak ./ length(weak_social(:));
        pmf_strong = (hist(strong_social, v2));
        pmf_strong = pmf_strong ./ length(strong_social(:));
        %Create CDF for weak edges
        cdf_weak_social = cumsum(pmf_weak);
        cdf_strong_social = cumsum(pmf_strong);
    
        %%Ks test
        vals_erDirichlet(i) = vals_erDirichlet(i) + max(abs(cdf_strong_social - cdf_weak_social));
    end
end
alphas_ERdirichlet_ks_test = alphas_erDirichlet./size_;
vals_ERdirichlet_ks_test = vals_erDirichlet./size_;