From c09d680851f6bbf0e20170e4a30be98b3c0485f5 Mon Sep 17 00:00:00 2001 From: Wataru Otsubo Date: Wed, 23 Oct 2024 19:47:28 +0900 Subject: [PATCH] update: skew_stats.jl --- examples/clock_skew_stats.svg | 1459 +++++++++++++++++---------------- examples/skew_stats.jl | 97 +-- 2 files changed, 807 insertions(+), 749 deletions(-) diff --git a/examples/clock_skew_stats.svg b/examples/clock_skew_stats.svg index 451bc2f..cd4a945 100644 --- a/examples/clock_skew_stats.svg +++ b/examples/clock_skew_stats.svg @@ -2,1290 +2,1347 @@ - + - + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + - - - - - - - - - - - - - - - - - - - - - - + - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - + - - + + - + - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - + - + - + - + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + - + - + - + - + - - - - + - - - - - - + - + + + + - - - - + - - + + - + + + + + + + + + + + + + + + + + + + + + - + - + - - - - - - - - - - - - - - - - - - - + - - + + - + + + + - + + + + + + + + + + + + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + - + + + + + + + + + + - - - - - - - - - - - - - - - - - - - + - - - - - - - - - - - - - + - + - - - - - - - + - - - - - - - + - - - - - - - + - + + + + + + + - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - - - - - - - - - - + + + + + + + + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + + + + - - + + - - + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - + - - - - - - - - - - - - - + + + + + + + + + + diff --git a/examples/skew_stats.jl b/examples/skew_stats.jl index 3e6c267..7ff3cbc 100644 --- a/examples/skew_stats.jl +++ b/examples/skew_stats.jl @@ -256,9 +256,52 @@ df_skew_stats = let df end +# ╔═╡ 893253c3-f0b2-401f-b892-b23291bcf5c1 +fig_skew_stats = let + fig, ax, sc = scatter( + df_skew_stats.mean_skew, + df_skew_stats.std_skew, + marker = :x, + color = (Makie.wong_colors()[1], 0.8), + axis = (title = "skew mean vs std", xlabel = "mean", ylabel = "std"), + ) + text!( + ax, + df_skew_stats.mean_skew, + df_skew_stats.std_skew, + text = string.(df_skew_stats.psbid), + color = (:gray, 0.5), + ) + fig +end + # ╔═╡ 6467dcaa-6bd6-45c7-8c08-b310a09b8b0b save("clock_skew_stats.svg", fig_skew_stats) +# ╔═╡ 19f85f1b-07d3-48c2-9ca6-a5a8eb1ab746 +let + fig, ax, sc = scatter( + df_skew_stats.mean_skew, + df_skew_stats.std_skew, + marker = :x, + color = (Makie.wong_colors()[1], 0.8), + axis = ( + title = "skew mean vs std", + xlabel = "mean", + ylabel = "std", + limits = ((-0.6, 1.2), (-0.01, 0.09)), + ), + ) + text!( + ax, + df_skew_stats.mean_skew, + df_skew_stats.std_skew, + text = string.(df_skew_stats.psbid), + color = (:gray, 0.5), + ) + fig +end + # ╔═╡ 79e2f5d8-4609-4e9f-949e-6dc1f88c0b19 df_skew_stats_abnormals = let df = filter([:mean_skew, :std_skew] => ((m, s) -> m > -5 && s > 1), df_skew_stats) @@ -324,11 +367,12 @@ md""" #### psbid: 462 - run: 105, 132 - psbid 460と同じく105がだめ -#### psbid: 545 (**問題の**) +#### psbid: 545 - run: 126, 132 - どちらも測定結果自体には変なところはない - どちらも1回だけ立ち上がりがある - 立ち上がりもそれほど長くない +- 126でQSPIpが失敗してるので結果は使わない、132を使う - 126が電源が不安定なときだったかもしれないが、記録がない #### psbid: 799 - run: 217, 236 @@ -373,7 +417,6 @@ df_skews_selected = combine( riseup = sdf.riseup |> first, ) elseif sdf.psbid[1] == 545 - @error "temp" df = filter(:runid => ==(132), sdf) @assert nrow(df) == 1 ( @@ -532,6 +575,9 @@ let fig end +# ╔═╡ 46b2a3cd-d2e6-4277-8b65-9c61f25f69e8 +3/57 + # ╔═╡ 13bb4978-b98d-44a3-a4b6-4241cadc609b let bins = range( @@ -660,52 +706,6 @@ let fig end -# ╔═╡ 19f85f1b-07d3-48c2-9ca6-a5a8eb1ab746 -fig_skew_stats = let - fig, ax, sc = scatter( - df_skew_stats.mean_skew, - df_skew_stats.std_skew, - marker = :x, - color = (Makie.wong_colors()[1], 0.8), - axis = ( - title = "skew mean vs std", - xlabel = "mean", - ylabel = "std", - limits = ((-0.6, 1.2), (-0.01, 0.09)), - ), - ) - text!( - ax, - df_skew_stats.mean_skew, - df_skew_stats.std_skew, - text = string.(df_skew_stats.psbid), - color = (:gray, 0.5), - ) - fig -end - -# ╔═╡ 893253c3-f0b2-401f-b892-b23291bcf5c1 -# ╠═╡ disabled = true -#=╠═╡ -fig_skew_stats = let - fig, ax, sc = scatter( - df_skew_stats.mean_skew, - df_skew_stats.std_skew, - marker = :x, - color = (Makie.wong_colors()[1], 0.8), - axis = (title = "skew mean vs std", xlabel = "mean", ylabel = "std"), - ) - text!( - ax, - df_skew_stats.mean_skew, - df_skew_stats.std_skew, - text = string.(df_skew_stats.psbid), - color = (:gray, 0.5), - ) - fig -end - ╠═╡ =# - # ╔═╡ Cell order: # ╟─f25e7e08-8a73-4cac-ac7c-d310725c558d # ╠═7c69d12c-80a5-11ef-2674-e155a3483342 @@ -771,6 +771,7 @@ end # ╠═13b4112f-96ab-41a1-8ef1-d940771f6ece # ╠═6a856a55-2611-41d0-a2c7-3645c066fc3c # ╠═4a4ef945-b312-44ed-ab62-ce01fc33f926 +# ╠═46b2a3cd-d2e6-4277-8b65-9c61f25f69e8 # ╠═13bb4978-b98d-44a3-a4b6-4241cadc609b # ╠═6c3c7669-d594-425d-bf05-9aa217c7656e # ╠═e9b3f541-d87d-4424-8b82-be1b7b3273d8