Skip to content

Commit

Permalink
Formatting
Browse files Browse the repository at this point in the history
  • Loading branch information
aurorarossi committed Jan 17, 2024
1 parent 8bf04c9 commit fe0f3f5
Show file tree
Hide file tree
Showing 2 changed files with 18 additions and 16 deletions.
27 changes: 14 additions & 13 deletions plot/plots.jl
Original file line number Diff line number Diff line change
Expand Up @@ -66,14 +66,15 @@ function load_euclidean_data(len)
temp_corr_all = zeros(10, len)
# load data
for i in 1:10
@load "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/data_correct/euclidean_300_thre_02_005_098_v02_$(i).jld2" avEU clEU plEU tccEU
@load "data/euclidean_300_thre_02_005_098_v02_$(i).jld2" avEU clEU plEU tccEU
average_deg_all[i, :] = avEU
clustering_all[i, :] = clEU
path_len_all[i, :] = plEU
temp_corr_all[i, :] = tccEU
end
return average_deg_all, clustering_all, path_len_all, temp_corr_all
end

function load_euclideanBorder_data(len)
# initialize variables
average_deg_all = zeros(10, len)
Expand All @@ -82,7 +83,7 @@ function load_euclideanBorder_data(len)
temp_corr_all = zeros(10, len)
# load data
for i in 1:10
@load "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/data_correct/euclideanBorder_300_r_0_005_08_v01_$(i).jld2" avEU clEU plEU tccEU
@load "data/euclideanBorder_300_r_0_005_08_v01_$(i).jld2" avEU clEU plEU tccEU
average_deg_all[i, :] = avEU
clustering_all[i, :] = clEU
path_len_all[i, :] = plEU
Expand All @@ -99,7 +100,7 @@ function load_hyperbolic_data(len)
temp_corr_all = zeros(10, len)
# load data
for i in 1:10
@load "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/data_correct/hyperbolic_300_alpha_65_R_0_05_18_v_60_$(i).jld2" avHY clHY plHY tccHY
@load "data/hyperbolic_300_alpha_65_R_0_05_18_v_60_$(i).jld2" avHY clHY plHY tccHY
average_deg_all[i, :] = avHY
clustering_all[i, :] = clHY
path_len_all[i, :] = plHY
Expand All @@ -116,7 +117,7 @@ function load_hyperbolic_SB_data(len)
temp_corr_all = zeros(10, len)
# load data
for i in 1:10
@load "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/data_correct/hyperbolic_300_alpha_80_R_0_05_18_v_90_$(i).jld2" avHY clHY plHY tccHY
@load "data/hyperbolic_300_alpha_80_R_0_05_18_v_90_$(i).jld2" avHY clHY plHY tccHY
average_deg_all[i, :] = avHY
clustering_all[i, :] = clHY
path_len_all[i, :] = plHY
Expand Down Expand Up @@ -202,12 +203,12 @@ function myplot(measure, legendplace, yaxisscale, xaxisscale, variable, variable
orange = 5
violet = 4
grey = 9
if measure !="Temporal Path Length"
yticks = ([0, 0.25, 0.5,0.75,1], [L"0", L"0.25", L"0.5", L"0.75",L"1"])
if measure != "Temporal Path Length"
yticks = ([0, 0.25, 0.5, 0.75, 1], [L"0", L"0.25", L"0.5", L"0.75", L"1"])
else
yticks = ([10^0,10^1], [L"10^0", L"10^1"])
yticks = ([10^0, 10^1], [L"10^0", L"10^1"])
end
p = plot(av, variable[1, :], fillrange=variable[3, :], fillalpha=fillalpha, c=blue, label="", lw=0, palette=:Set1_9, legend=legendplace, legendfontsize=12, xlims=(0, 170), grid=false, dpi=1200, yaxis=yaxisscale, xaxis=xaxisscale, xticks = ([0, 50, 100, 150, 200], [L"0", L"50", L"100", L"150", L"200"]), yticks = yticks)
p = plot(av, variable[1, :], fillrange=variable[3, :], fillalpha=fillalpha, c=blue, label="", lw=0, palette=:Set1_9, legend=legendplace, legendfontsize=12, xlims=(0, 170), grid=false, dpi=1200, yaxis=yaxisscale, xaxis=xaxisscale, xticks=([0, 50, 100, 150, 200], [L"0", L"50", L"100", L"150", L"200"]), yticks=yticks)

plot!(avRE, variableRE[1, :], fillrange=variableRE[3, :], fillalpha=fillalpha, c=pink, label="", lw=0)
plot!(avRE, variableRE[2, :], label=L"\textrm{RTE}", lw=line_lw, c=pink)
Expand Down Expand Up @@ -252,11 +253,11 @@ function main()
# display(p3)
# display(p4)
# display(p5)
savefig(p1, "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/images/tsw.png")
savefig(p2, "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/images/tswsb.png")
savefig(p3, "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/images/tcc.png")
savefig(p4, "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/images/tclustering.png")
savefig(p5, "/user/aurossi/home/mri_networks/TemporalBrainNetworksCode/images/tpath.png")
savefig(p1, "images/tsw.png")
savefig(p2, "images/tswsb.png")
savefig(p3, "images/tcc.png")
savefig(p4, "images/tclustering.png")
savefig(p5, "images/tpath.png")
end

main()
7 changes: 4 additions & 3 deletions plot/simple_plot.jl
Original file line number Diff line number Diff line change
@@ -1,6 +1,7 @@
include("../optimization/findingParameters.jl")

function myplot()
subjects= readdlm("TemporalBrainNetworksCode/src/filtered-subjects-mod.txt", Int)
subjects = readdlm("src/filtered-subjects-mod.txt", Int)
Rrange = collect(0.5:0.5:14)
thresholds = append!(collect(0.2:0.05:0.90), collect(0.92:0.02:0.98))
average_deg_all, clustering_all, path_len_all, temp_clustering_all = load_fMRI_data(subjects, thresholds)
Expand All @@ -11,8 +12,8 @@ function myplot()
avEUb, clEUb, plEUb, tccEUb = load_euclideanBorder_data()
p = plot(average_deg_quantile, small_world_quantile[2, :], label="fMRI", xlabel="Average degree", ylabel="Small-worldness", legend=:topleft, xlims=(0, 170), ylims=(0, 0.6))
plot!(p, average_deg_HY[1, :, 1], small_world_HY[1, :, 1], label="HY")
plot!(p, avEU[1,:], (clEU ./ plEU)[1, :], label="EU")
plot!(p, avEUb[9,:], (clEUb ./ plEUb)[9, :], label="EU border")
plot!(p, avEU[1, :], (clEU./plEU)[1, :], label="EU")
plot!(p, avEUb[9, :], (clEUb./plEUb)[9, :], label="EU border")
end

myplot()

0 comments on commit fe0f3f5

Please sign in to comment.