@Dawoon did you take a look at this post? NIBR data analysis - #51 by vijay
an excerpt from there
avg_pk_plots = summary_observations_vs_time(ncapop,
color = :black, linewidth = 2, whiskerwidth = 8,
paginate = true,
limit = 1,
axis = (xlabel = "Time (hours)",
ylabel = "Concentration (μg/mL)",
yscale=log10, ytickformat=x -> string.(round.(x; digits=1)),
yticks = [0.1, 1, 10],
ygridwidth = 3,
yminorgridcolor = :darkgrey,
yminorticksvisible = true,
yminorgridvisible = true,
yminorticks = IntervalsBetween(10),
xminorticksvisible = true,
xminorgridvisible = true,
xminorticks = IntervalsBetween(5),
limits = (nothing, nothing, nothing, 30),
spinewidth = 2),
#facet = ( combinelabels = true,),
figure = (resolution = (800,800),
fontsize = 22),)