Fix notebook

RunningAnalysis.jl

@@ -13,6 +13,8 @@ begin
 	using Dates
 	using Plots
 	using Statistics
+	using StatsBase
+	using StatsPlots
 end
 
 # ╔═╡ 9fee4ffb-83c0-46df-b84a-ef902a5e99ab
@@ -42,7 +44,6 @@ begin
 	# Start with some average time in each stage. Timestamp should be in a more usable
 	# format.
 	metrics.ParsedTimestamp = DateTime.(metrics.Timestamp, dateformat"y-m-d H:M:S")
-	# Equivalent: metrics[:ParsedTimestamp], metrics[:, :ParsedTimestamp]
 	describe(metrics)
 end
 
@@ -151,15 +152,24 @@ md"""
 begin
 	@pyimport garmin_fit_sdk
 	
-	fit_files = readdir("fit_files"; join=true)
+	fit_files = readdir("./fit_files"; join=true)
 
 	function decode_fit_file(fp::AbstractString)::Tuple{Any,Any}
 		stream = garmin_fit_sdk.Stream.from_file(fp)
 		decoder = garmin_fit_sdk.Decoder(stream)
+		try
 			decoder.read()
+		finally
+			stream.close()
+		end
 	end
 
-	decoded_fit_files = [dc[1] for dc in decode_fit_file.(fit_files)]
+	decoded = []
+	for file ∈ fit_files
+		push!(decoded, decode_fit_file(file))
+	end
+
+	decoded_fit_files = [dc[1] for dc in decoded]
 end
 
 # ╔═╡ 065c7811-0203-44fd-bb88-49946d674fd8
@@ -177,17 +187,193 @@ md"""
 
 # ╔═╡ 39a46471-663e-4395-b4bd-2380932d7189
 begin
-	decoded_fit_files[1]["record_mesgs"]
+	function string_keys(d::Dict{Any,Any})::Dict{String,Any}
-
+		Dict([k => v for (k, v) in d if k isa AbstractString])
-	function extract_record_and_timestamp(fit_file::Dict{Any, Any})::Tuple{DateTime, Vector{Dict{Any, Any}}}
-		ts = DateTime(fit_file["activity_mesgs"][1]["timestamp"])
-		records = filter(x -> x isa AbstractString)
-		ts, fit_file["record_mesgs"]
 	end
 	
-	extract_record_and_timestamp.(decoded_fit_files)
+	function extract_record_and_timestamp(fit_file::Dict{Any, Any})
+		ts = DateTime(fit_file["activity_mesgs"][1]["timestamp"])
+		messages = string_keys.(fit_file["record_mesgs"])
+
+		try
+			ts, DataFrame(messages)
+		catch
+			nothing
+		end
+	end
+
+	decoded_files = extract_record_and_timestamp.(decoded_fit_files)
 end
 
+
+# ╔═╡ ed9e2e1f-cebe-4a0b-969a-93d344578a28
+md"""
+
+## Mapping my data to DataFrame
+
+- Extract the data needed
+- Write function that pulls data out of PyDict
+
+"""
+
+# ╔═╡ b7854b76-eb0e-4d78-9bd6-9b0258a725ba
+describe(decoded_files[1][2])
+
+# ╔═╡ 3e1d6c7b-41b5-435a-88f0-b0f3ab442662
+begin
+	# A generic graph_run function
+	function graph_run(decoded_file::Tuple{DateTime, DataFrame})
+		label = "Run on $(decoded_file[1])"
+		df = decoded_file[2]
+		l = @layout [a b; c d]
+		p = plot(df."timestamp", df."heart_rate"; title="Heart Rate", ylabel="BPM", xticks=nothing)
+		p2 = plot(df.timestamp, df.altitude; title="Altitude", xticks=nothing)
+		p3 = scatter(df.timestamp, df.cadence, title="Cadence", xticks=nothing)
+		p4 = plot(df.timestamp, df.speed; title="Speed", xticks=nothing)
+		plot(p, p2, p3, p4, layout = l)
+	end
+end
+
+# ╔═╡ 2f4cdc08-27a1-4e92-8709-8021cae771c6
+graph_run(decoded_files[14])
+
+# ╔═╡ a9d1b364-65d5-44bc-965f-61a5f2b47a36
+md"""
+## Graphing Runs
+
+- Show speed, cadence, heartrate, and altitude
+- Trail runs appear to have varied cadence, road runs are more stable
+"""
+
+# ╔═╡ 0c5d6cdb-7255-4712-92bc-eb3f3182d7eb
+graph_run(decoded_files[19])
+
+# ╔═╡ e4461dae-f648-4c34-a963-17dfb3f0220a
+md"""
+## Trends across Training
+
+- What was my average heartrate for each run?
+- Was my average heartrate effected by sleep?
+"""
+
+# ╔═╡ d1588e94-a450-4d8c-9811-69873deed252
+begin
+	function calc_avgs(data::Union{Nothing, Tuple{DateTime, DataFrame}})
+		if isnothing(data)
+			return nothing
+		end
+		try
+			return (timestamp=data[1], mean_hr=mean(data[2].heart_rate), mean_cadence=mean(data[2].cadence), mean_alt=mean(data[2].altitude))
+		catch
+			return nothing
+		end
+	end
+	
+	avg_heartrates = DataFrame(filter(f -> !isnothing(f), calc_avgs.(decoded_files)))
+
+	ps = scatter(avg_heartrates.timestamp, avg_heartrates.mean_hr, title="Avg Heartrate over Time", label="hr", xlabel="Timestamp", ylabel="Beats per minute", ylims=(115, 160))
+	hspan!(ps, [145, maximum(avg_heartrates.mean_hr)]; color=:red, alpha=0.2, label="endurance run")
+	hspan!(ps, [minimum(avg_heartrates.mean_hr), 145]; color=:green, alpha=0.2, label="recovery run")
+end
+
+# ╔═╡ 0b6bce6e-22e5-4a8d-95f4-4708deb07cdd
+md"""
+## Determining if Heartrate was affected by sleep
+
+- Will need to use the workouts to determine which days were endurance runs
+- Find sleep data for same day
+- Stitch together to see results
+"""
+
+# ╔═╡ 4561f42b-2658-4c74-b642-6b2f92cd3598
+begin
+	# Copy the existing dataframe.
+	hr_sleep_df = deepcopy(avg_heartrates)
+
+	# Get the days that we have an endurance run
+	my_endurance_days = workouts[:, :WorkoutDay]
+
+	my_endurance_selector = [d in my_endurance_days for d in Date.(hr_sleep_df.timestamp)]
+
+	my_endurance_days_df = hr_sleep_df[my_endurance_selector, :]
+	my_endurance_days_df.Date = Date.(my_endurance_days_df.timestamp)
+	# sleep_hours = metrics[metrics.Type .== "Sleep Hours", [:Value, :ParsedTimestamp]]
+	sleep_on_endurance_days = [d in my_endurance_days for d in Date.(sleep_hours.ParsedTimestamp)]
+	sleep_on_endurance_days_df = sleep_hours[sleep_on_endurance_days, :]
+	sleep_on_endurance_days_df.Date = Date.(sleep_on_endurance_days_df.ParsedTimestamp)
+
+	sleep_hr_df = leftjoin(my_endurance_days_df, sleep_on_endurance_days_df; on = :Date)
+	ns = scatter(sleep_hr_df.mean_hr, sleep_hr_df.Value; xlims=(120, 160), title="Heartrate over Sleep", xlabel="BPM", ylabel="Hours Slept")
+	vspan!(ns, [140, maximum(sleep_hr_df.mean_hr)]; color=:red, alpha=0.2, label="endurance run")
+	vspan!(ns, [minimum(sleep_hr_df.mean_hr), 140]; color=:green, alpha=0.2, label="recovery run")
+end
+
+# ╔═╡ 27a59d46-f838-4d8f-9c66-a7edc1eea8dd
+md"""
+## Findings
+
+- No discernible relationship
+- Appears to have clusters, left half is recovery, right half is endurance
+
+## Next Question
+
+- Does Cadence Affect HR?
+"""
+
+# ╔═╡ 504ed9c8-ff54-4a25-84b9-d832fc435486
+begin
+	cadence_hr_df = deepcopy(avg_heartrates)
+	# cadence_hr_df.Date = Date.(cadence_hr_df.timestamp)
+	scatter(cadence_hr_df.mean_hr, cadence_hr_df.mean_cadence; xlabel="HR", ylabel="Cadence", title="HR vs. Cadence")
+end
+
+# ╔═╡ 873bcfa3-a989-426a-8402-6b81a8bbb46a
+md"""
+## Findings
+
+- No relationship, because I always ran with the same cadence!
+
+## Next Question
+
+- Does altitude affect HR?
+"""
+
+# ╔═╡ b10b4931-e60f-4791-8dae-d6d58cd71e66
+begin
+	altitude_hr_df = deepcopy(avg_heartrates)
+	scatter(altitude_hr_df.mean_hr, altitude_hr_df.mean_alt; xlabel="HR", ylabel="Mean Altitude", title="HR vs. Altitude")
+end
+
+# ╔═╡ b14e7de2-b91f-40d9-8aa8-237d82f25286
+md"""
+# Findings
+
+- Altitude does appear to affect heart rate.
+- Majority of runs were not at altitude, show normal distribution
+- Runs at higher altitudes appear to cluster at higher HR
+
+# Next Questions
+
+- Do my data follow the normal distribution?
+"""
+
+# ╔═╡ 1fd645ee-0680-4de8-b1eb-1c94c2069ae9
+begin
+	d = fit(UnitRangeTransform, avg_heartrates.mean_hr)
+	hr_normalized = StatsBase.transform(d, avg_heartrates.mean_hr)
+
+	avg_heartrates.hr_normalized = hr_normalized
+	@df avg_heartrates histogram(:mean_hr; bins=25, title="HR Histogram")
+end
+
+# ╔═╡ c3286275-3cff-4771-9ae6-78a960ba0484
+md"""
+## Findings
+
+- Kind of normally distributed, but have two clusters
+- Recovery vs. Endurance again!
+"""
+
 # ╔═╡ 00000000-0000-0000-0000-000000000001
 PLUTO_PROJECT_TOML_CONTENTS = """
 [deps]
@@ -198,6 +384,8 @@ Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 PyCall = "438e738f-606a-5dbb-bf0a-cddfbfd45ab0"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
+StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
+StatsPlots = "f3b207a7-027a-5e70-b257-86293d7955fd"
 
 [compat]
 CSV = "~0.10.9"
@@ -205,6 +393,8 @@ Conda = "~1.8.0"
 DataFrames = "~1.5.0"
 Plots = "~1.38.10"
 PyCall = "~1.95.1"
+StatsBase = "~0.33.21"
+StatsPlots = "~0.15.4"
 """
 
 # ╔═╡ 00000000-0000-0000-0000-000000000002
@@ -213,15 +403,45 @@ PLUTO_MANIFEST_TOML_CONTENTS = """
 
 julia_version = "1.8.5"
 manifest_format = "2.0"
-project_hash = "68d73a09dbba44f41a25a749e93b301504e20a2a"
+project_hash = "d3e66daec1ac6ba734bfea3d7957f5868c5983b5"
+
+[[deps.AbstractFFTs]]
+deps = ["ChainRulesCore", "LinearAlgebra"]
+git-tree-sha1 = "16b6dbc4cf7caee4e1e75c49485ec67b667098a0"
+uuid = "621f4979-c628-5d54-868e-fcf4e3e8185c"
+version = "1.3.1"
+
+[[deps.Adapt]]
+deps = ["LinearAlgebra", "Requires"]
+git-tree-sha1 = "cc37d689f599e8df4f464b2fa3870ff7db7492ef"
+uuid = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
+version = "3.6.1"
 
 [[deps.ArgTools]]
 uuid = "0dad84c5-d112-42e6-8d28-ef12dabb789f"
 version = "1.1.1"
 
+[[deps.Arpack]]
+deps = ["Arpack_jll", "Libdl", "LinearAlgebra", "Logging"]
+git-tree-sha1 = "9b9b347613394885fd1c8c7729bfc60528faa436"
+uuid = "7d9fca2a-8960-54d3-9f78-7d1dccf2cb97"
+version = "0.5.4"
+
+[[deps.Arpack_jll]]
+deps = ["Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "OpenBLAS_jll", "Pkg"]
+git-tree-sha1 = "5ba6c757e8feccf03a1554dfaf3e26b3cfc7fd5e"
+uuid = "68821587-b530-5797-8361-c406ea357684"
+version = "3.5.1+1"
+
 [[deps.Artifacts]]
 uuid = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"
 
+[[deps.AxisAlgorithms]]
+deps = ["LinearAlgebra", "Random", "SparseArrays", "WoodburyMatrices"]
+git-tree-sha1 = "66771c8d21c8ff5e3a93379480a2307ac36863f7"
+uuid = "13072b0f-2c55-5437-9ae7-d433b7a33950"
+version = "1.0.1"
+
 [[deps.Base64]]
 uuid = "2a0f44e3-6c83-55bd-87e4-b1978d98bd5f"
 
@@ -248,6 +468,12 @@ git-tree-sha1 = "4b859a208b2397a7a623a03449e4636bdb17bcf2"
 uuid = "83423d85-b0ee-5818-9007-b63ccbeb887a"
 version = "1.16.1+1"
 
+[[deps.Calculus]]
+deps = ["LinearAlgebra"]
+git-tree-sha1 = "f641eb0a4f00c343bbc32346e1217b86f3ce9dad"
+uuid = "49dc2e85-a5d0-5ad3-a950-438e2897f1b9"
+version = "0.5.1"
+
 [[deps.ChainRulesCore]]
 deps = ["Compat", "LinearAlgebra", "SparseArrays"]
 git-tree-sha1 = "c6d890a52d2c4d55d326439580c3b8d0875a77d9"
@@ -260,6 +486,12 @@ git-tree-sha1 = "f84967c4497e0e1955f9a582c232b02847c5f589"
 uuid = "9e997f8a-9a97-42d5-a9f1-ce6bfc15e2c0"
 version = "0.1.7"
 
+[[deps.Clustering]]
+deps = ["Distances", "LinearAlgebra", "NearestNeighbors", "Printf", "Random", "SparseArrays", "Statistics", "StatsBase"]
+git-tree-sha1 = "7ebbd653f74504447f1c33b91cd706a69a1b189f"
+uuid = "aaaa29a8-35af-508c-8bc3-b662a17a0fe5"
+version = "0.14.4"
+
 [[deps.CodecZlib]]
 deps = ["TranscodingStreams", "Zlib_jll"]
 git-tree-sha1 = "9c209fb7536406834aa938fb149964b985de6c83"
@@ -345,6 +577,12 @@ git-tree-sha1 = "bfc1187b79289637fa0ef6d4436ebdfe6905cbd6"
 uuid = "e2d170a0-9d28-54be-80f0-106bbe20a464"
 version = "1.0.0"
 
+[[deps.DataValues]]
+deps = ["DataValueInterfaces", "Dates"]
+git-tree-sha1 = "d88a19299eba280a6d062e135a43f00323ae70bf"
+uuid = "e7dc6d0d-1eca-5fa6-8ad6-5aecde8b7ea5"
+version = "0.4.13"
+
 [[deps.Dates]]
 deps = ["Printf"]
 uuid = "ade2ca70-3891-5945-98fb-dc099432e06a"
@@ -353,6 +591,28 @@ uuid = "ade2ca70-3891-5945-98fb-dc099432e06a"
 deps = ["Mmap"]
 uuid = "8bb1440f-4735-579b-a4ab-409b98df4dab"
 
+[[deps.DensityInterface]]
+deps = ["InverseFunctions", "Test"]
+git-tree-sha1 = "80c3e8639e3353e5d2912fb3a1916b8455e2494b"
+uuid = "b429d917-457f-4dbc-8f4c-0cc954292b1d"
+version = "0.4.0"
+
+[[deps.Distances]]
+deps = ["LinearAlgebra", "SparseArrays", "Statistics", "StatsAPI"]
+git-tree-sha1 = "49eba9ad9f7ead780bfb7ee319f962c811c6d3b2"
+uuid = "b4f34e82-e78d-54a5-968a-f98e89d6e8f7"
+version = "0.10.8"
+
+[[deps.Distributed]]
+deps = ["Random", "Serialization", "Sockets"]
+uuid = "8ba89e20-285c-5b6f-9357-94700520ee1b"
+
+[[deps.Distributions]]
+deps = ["ChainRulesCore", "DensityInterface", "FillArrays", "LinearAlgebra", "PDMats", "Printf", "QuadGK", "Random", "SparseArrays", "SpecialFunctions", "Statistics", "StatsBase", "StatsFuns", "Test"]
+git-tree-sha1 = "180538ef4e3aa02b01413055a7a9e8b6047663e1"
+uuid = "31c24e10-a181-5473-b8eb-7969acd0382f"
+version = "0.25.88"
+
 [[deps.DocStringExtensions]]
 deps = ["LibGit2"]
 git-tree-sha1 = "2fb1e02f2b635d0845df5d7c167fec4dd739b00d"
@@ -364,6 +624,12 @@ deps = ["ArgTools", "FileWatching", "LibCURL", "NetworkOptions"]
 uuid = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
 version = "1.6.0"
 
+[[deps.DualNumbers]]
+deps = ["Calculus", "NaNMath", "SpecialFunctions"]
+git-tree-sha1 = "5837a837389fccf076445fce071c8ddaea35a566"
+uuid = "fa6b7ba4-c1ee-5f82-b5fc-ecf0adba8f74"
+version = "0.6.8"
+
 [[deps.Expat_jll]]
 deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
 git-tree-sha1 = "bad72f730e9e91c08d9427d5e8db95478a3c323d"
@@ -382,6 +648,18 @@ git-tree-sha1 = "74faea50c1d007c85837327f6775bea60b5492dd"
 uuid = "b22a6f82-2f65-5046-a5b2-351ab43fb4e5"
 version = "4.4.2+2"
 
+[[deps.FFTW]]
+deps = ["AbstractFFTs", "FFTW_jll", "LinearAlgebra", "MKL_jll", "Preferences", "Reexport"]
+git-tree-sha1 = "f9818144ce7c8c41edf5c4c179c684d92aa4d9fe"
+uuid = "7a1cc6ca-52ef-59f5-83cd-3a7055c09341"
+version = "1.6.0"
+
+[[deps.FFTW_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "c6033cc3892d0ef5bb9cd29b7f2f0331ea5184ea"
+uuid = "f5851436-0d7a-5f13-b9de-f02708fd171a"
+version = "3.3.10+0"
+
 [[deps.FilePathsBase]]
 deps = ["Compat", "Dates", "Mmap", "Printf", "Test", "UUIDs"]
 git-tree-sha1 = "e27c4ebe80e8699540f2d6c805cc12203b614f12"
@@ -391,6 +669,12 @@ version = "0.9.20"
 [[deps.FileWatching]]
 uuid = "7b1f6079-737a-58dc-b8bc-7a2ca5c1b5ee"
 
+[[deps.FillArrays]]
+deps = ["LinearAlgebra", "Random", "SparseArrays", "Statistics"]
+git-tree-sha1 = "fc86b4fd3eff76c3ce4f5e96e2fdfa6282722885"
+uuid = "1a297f60-69ca-5386-bcde-b61e274b549b"
+version = "1.0.0"
+
 [[deps.FixedPointNumbers]]
 deps = ["Statistics"]
 git-tree-sha1 = "335bfdceacc84c5cdf16aadc768aa5ddfc5383cc"
@@ -478,16 +762,34 @@ git-tree-sha1 = "129acf094d168394e80ee1dc4bc06ec835e510a3"
 uuid = "2e76f6c2-a576-52d4-95c1-20adfe4de566"
 version = "2.8.1+1"
 
+[[deps.HypergeometricFunctions]]
+deps = ["DualNumbers", "LinearAlgebra", "OpenLibm_jll", "SpecialFunctions"]
+git-tree-sha1 = "432b5b03176f8182bd6841fbfc42c718506a2d5f"
+uuid = "34004b35-14d8-5ef3-9330-4cdb6864b03a"
+version = "0.3.15"
+
 [[deps.InlineStrings]]
 deps = ["Parsers"]
 git-tree-sha1 = "9cc2baf75c6d09f9da536ddf58eb2f29dedaf461"
 uuid = "842dd82b-1e85-43dc-bf29-5d0ee9dffc48"
 version = "1.4.0"
 
+[[deps.IntelOpenMP_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "0cb9352ef2e01574eeebdb102948a58740dcaf83"
+uuid = "1d5cc7b8-4909-519e-a0f8-d0f5ad9712d0"
+version = "2023.1.0+0"
+
 [[deps.InteractiveUtils]]
 deps = ["Markdown"]
 uuid = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 
+[[deps.Interpolations]]
+deps = ["Adapt", "AxisAlgorithms", "ChainRulesCore", "LinearAlgebra", "OffsetArrays", "Random", "Ratios", "Requires", "SharedArrays", "SparseArrays", "StaticArrays", "WoodburyMatrices"]
+git-tree-sha1 = "721ec2cf720536ad005cb38f50dbba7b02419a15"
+uuid = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
+version = "0.14.7"
+
 [[deps.InverseFunctions]]
 deps = ["Test"]
 git-tree-sha1 = "6667aadd1cdee2c6cd068128b3d226ebc4fb0c67"
@@ -533,6 +835,12 @@ git-tree-sha1 = "6f2675ef130a300a112286de91973805fcc5ffbc"
 uuid = "aacddb02-875f-59d6-b918-886e6ef4fbf8"
 version = "2.1.91+0"
 
+[[deps.KernelDensity]]
+deps = ["Distributions", "DocStringExtensions", "FFTW", "Interpolations", "StatsBase"]
+git-tree-sha1 = "4a9513ad756e712177bd342ba6c022b515ed8d76"
+uuid = "5ab0869b-81aa-558d-bb23-cbf5423bbe9b"
+version = "0.6.6"
+
 [[deps.LAME_jll]]
 deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
 git-tree-sha1 = "f6250b16881adf048549549fba48b1161acdac8c"
@@ -562,6 +870,10 @@ git-tree-sha1 = "8c57307b5d9bb3be1ff2da469063628631d4d51e"
 uuid = "23fbe1c1-3f47-55db-b15f-69d7ec21a316"
 version = "0.15.21"
 
+[[deps.LazyArtifacts]]
+deps = ["Artifacts", "Pkg"]
+uuid = "4af54fe1-eca0-43a8-85a7-787d91b784e3"
+
 [[deps.LibCURL]]
 deps = ["LibCURL_jll", "MozillaCACerts_jll"]
 uuid = "b27032c2-a3e7-50c8-80cd-2d36dbcbfd21"
@@ -651,6 +963,12 @@ git-tree-sha1 = "cedb76b37bc5a6c702ade66be44f831fa23c681e"
 uuid = "e6f89c97-d47a-5376-807f-9c37f3926c36"
 version = "1.0.0"
 
+[[deps.MKL_jll]]
+deps = ["Artifacts", "IntelOpenMP_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "Pkg"]
+git-tree-sha1 = "2ce8695e1e699b68702c03402672a69f54b8aca9"
+uuid = "856f044c-d86e-5d09-b602-aeab76dc8ba7"
+version = "2022.2.0+0"
+
 [[deps.MacroTools]]
 deps = ["Markdown", "Random"]
 git-tree-sha1 = "42324d08725e200c23d4dfb549e0d5d89dede2d2"
@@ -690,16 +1008,39 @@ uuid = "a63ad114-7e13-5084-954f-fe012c677804"
 uuid = "14a3606d-f60d-562e-9121-12d972cd8159"
 version = "2022.2.1"
 
+[[deps.MultivariateStats]]
+deps = ["Arpack", "LinearAlgebra", "SparseArrays", "Statistics", "StatsAPI", "StatsBase"]
+git-tree-sha1 = "91a48569383df24f0fd2baf789df2aade3d0ad80"
+uuid = "6f286f6a-111f-5878-ab1e-185364afe411"
+version = "0.10.1"
+
 [[deps.NaNMath]]
 deps = ["OpenLibm_jll"]
 git-tree-sha1 = "0877504529a3e5c3343c6f8b4c0381e57e4387e4"
 uuid = "77ba4419-2d1f-58cd-9bb1-8ffee604a2e3"
 version = "1.0.2"
 
+[[deps.NearestNeighbors]]
+deps = ["Distances", "StaticArrays"]
+git-tree-sha1 = "2c3726ceb3388917602169bed973dbc97f1b51a8"
+uuid = "b8a86587-4115-5ab1-83bc-aa920d37bbce"
+version = "0.4.13"
+
 [[deps.NetworkOptions]]
 uuid = "ca575930-c2e3-43a9-ace4-1e988b2c1908"
 version = "1.2.0"
 
+[[deps.Observables]]
+git-tree-sha1 = "6862738f9796b3edc1c09d0890afce4eca9e7e93"
+uuid = "510215fc-4207-5dde-b226-833fc4488ee2"
+version = "0.5.4"
+
+[[deps.OffsetArrays]]
+deps = ["Adapt"]
+git-tree-sha1 = "82d7c9e310fe55aa54996e6f7f94674e2a38fcb4"
+uuid = "6fe1bfb0-de20-5000-8ca7-80f57d26f881"
+version = "1.12.9"
+
 [[deps.Ogg_jll]]
 deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
 git-tree-sha1 = "887579a3eb005446d514ab7aeac5d1d027658b8f"
@@ -750,6 +1091,12 @@ deps = ["Artifacts", "Libdl"]
 uuid = "efcefdf7-47ab-520b-bdef-62a2eaa19f15"
 version = "10.40.0+0"
 
+[[deps.PDMats]]
+deps = ["LinearAlgebra", "SparseArrays", "SuiteSparse"]
+git-tree-sha1 = "67eae2738d63117a196f497d7db789821bce61d1"
+uuid = "90014a1f-27ba-587c-ab20-58faa44d9150"
+version = "0.11.17"
+
 [[deps.Parsers]]
 deps = ["Dates", "SnoopPrecompile"]
 git-tree-sha1 = "478ac6c952fddd4399e71d4779797c538d0ff2bf"
@@ -830,6 +1177,12 @@ git-tree-sha1 = "0c03844e2231e12fda4d0086fd7cbe4098ee8dc5"
 uuid = "ea2cea3b-5b76-57ae-a6ef-0a8af62496e1"
 version = "5.15.3+2"
 
+[[deps.QuadGK]]
+deps = ["DataStructures", "LinearAlgebra"]
+git-tree-sha1 = "6ec7ac8412e83d57e313393220879ede1740f9ee"
+uuid = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
+version = "2.8.2"
+
 [[deps.REPL]]
 deps = ["InteractiveUtils", "Markdown", "Sockets", "Unicode"]
 uuid = "3fa0cd96-eef1-5676-8a61-b3b8758bbffb"
@@ -838,6 +1191,12 @@ uuid = "3fa0cd96-eef1-5676-8a61-b3b8758bbffb"
 deps = ["SHA", "Serialization"]
 uuid = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 
+[[deps.Ratios]]
+deps = ["Requires"]
+git-tree-sha1 = "6d7bb727e76147ba18eed998700998e17b8e4911"
+uuid = "c84ed2f1-dad5-54f0-aa8e-dbefe2724439"
+version = "0.4.4"
+
 [[deps.RecipesBase]]
 deps = ["PrecompileTools"]
 git-tree-sha1 = "5c3d09cc4f31f5fc6af001c250bf1278733100ff"
@@ -867,6 +1226,18 @@ git-tree-sha1 = "838a3a4188e2ded87a4f9f184b4b0d78a1e91cb7"
 uuid = "ae029012-a4dd-5104-9daa-d747884805df"
 version = "1.3.0"
 
+[[deps.Rmath]]
+deps = ["Random", "Rmath_jll"]
+git-tree-sha1 = "f65dcb5fa46aee0cf9ed6274ccbd597adc49aa7b"
+uuid = "79098fc4-a85e-5d69-aa6a-4863f24498fa"
+version = "0.7.1"
+
+[[deps.Rmath_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "6ed52fdd3382cf21947b15e8870ac0ddbff736da"
+uuid = "f50d1b31-88e8-58de-be2c-1cc44531875f"
+version = "0.4.0+0"
+
 [[deps.SHA]]
 uuid = "ea8e919c-243c-51af-8825-aaa63cd721ce"
 version = "0.7.0"
@@ -886,6 +1257,10 @@ version = "1.3.18"
 [[deps.Serialization]]
 uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"
 
+[[deps.SharedArrays]]
+deps = ["Distributed", "Mmap", "Random", "Serialization"]
+uuid = "1a1011a3-84de-559e-8e89-a11a2f7dc383"
+
 [[deps.Showoff]]
 deps = ["Dates", "Grisu"]
 git-tree-sha1 = "91eddf657aca81df9ae6ceb20b959ae5653ad1de"
@@ -922,6 +1297,17 @@ git-tree-sha1 = "ef28127915f4229c971eb43f3fc075dd3fe91880"
 uuid = "276daf66-3868-5448-9aa4-cd146d93841b"
 version = "2.2.0"
 
+[[deps.StaticArrays]]
+deps = ["LinearAlgebra", "Random", "StaticArraysCore", "Statistics"]
+git-tree-sha1 = "c262c8e978048c2b095be1672c9bee55b4619521"
+uuid = "90137ffa-7385-5640-81b9-e52037218182"
+version = "1.5.24"
+
+[[deps.StaticArraysCore]]
+git-tree-sha1 = "6b7ba252635a5eff6a0b0664a41ee140a1c9e72a"
+uuid = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
+version = "1.4.0"
+
 [[deps.Statistics]]
 deps = ["LinearAlgebra", "SparseArrays"]
 uuid = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
@@ -938,16 +1324,38 @@ git-tree-sha1 = "d1bf48bfcc554a3761a133fe3a9bb01488e06916"
 uuid = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
 version = "0.33.21"
 
+[[deps.StatsFuns]]
+deps = ["ChainRulesCore", "HypergeometricFunctions", "InverseFunctions", "IrrationalConstants", "LogExpFunctions", "Reexport", "Rmath", "SpecialFunctions"]
+git-tree-sha1 = "f625d686d5a88bcd2b15cd81f18f98186fdc0c9a"
+uuid = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
+version = "1.3.0"
+
+[[deps.StatsPlots]]
+deps = ["AbstractFFTs", "Clustering", "DataStructures", "DataValues", "Distributions", "Interpolations", "KernelDensity", "LinearAlgebra", "MultivariateStats", "NaNMath", "Observables", "Plots", "RecipesBase", "RecipesPipeline", "Reexport", "StatsBase", "TableOperations", "Tables", "Widgets"]
+git-tree-sha1 = "e0d5bc26226ab1b7648278169858adcfbd861780"
+uuid = "f3b207a7-027a-5e70-b257-86293d7955fd"
+version = "0.15.4"
+
 [[deps.StringManipulation]]
 git-tree-sha1 = "46da2434b41f41ac3594ee9816ce5541c6096123"
 uuid = "892a3eda-7b42-436c-8928-eab12a02cf0e"
 version = "0.3.0"
 
+[[deps.SuiteSparse]]
+deps = ["Libdl", "LinearAlgebra", "Serialization", "SparseArrays"]
+uuid = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
+
 [[deps.TOML]]
 deps = ["Dates"]
 uuid = "fa267f1f-6049-4f14-aa54-33bafae1ed76"
 version = "1.0.0"
 
+[[deps.TableOperations]]
+deps = ["SentinelArrays", "Tables", "Test"]
+git-tree-sha1 = "e383c87cf2a1dc41fa30c093b2a19877c83e1bc1"
+uuid = "ab02a1b2-a7df-11e8-156e-fb1833f50b87"
+version = "1.2.0"
+
 [[deps.TableTraits]]
 deps = ["IteratorInterfaceExtensions"]
 git-tree-sha1 = "c06b2f539df1c6efa794486abfb6ed2022561a39"
@@ -1027,6 +1435,18 @@ git-tree-sha1 = "b1be2855ed9ed8eac54e5caff2afcdb442d52c23"
 uuid = "ea10d353-3f73-51f8-a26c-33c1cb351aa5"
 version = "1.4.2"
 
+[[deps.Widgets]]
+deps = ["Colors", "Dates", "Observables", "OrderedCollections"]
+git-tree-sha1 = "fcdae142c1cfc7d89de2d11e08721d0f2f86c98a"
+uuid = "cc8bc4a8-27d6-5769-a93b-9d913e69aa62"
+version = "0.6.6"
+
+[[deps.WoodburyMatrices]]
+deps = ["LinearAlgebra", "SparseArrays"]
+git-tree-sha1 = "de67fa59e33ad156a590055375a30b23c40299d3"
+uuid = "efce3f68-66dc-5838-9240-27a6d6f5f9b6"
+version = "0.5.5"
+
 [[deps.WorkerUtilities]]
 git-tree-sha1 = "cd1659ba0d57b71a464a29e64dbc67cfe83d54e7"
 uuid = "76eceee3-57b5-4d4a-8e66-0e911cebbf60"
@@ -1268,5 +1688,22 @@ version = "1.4.1+0"
 # ╠═7cb0daaf-6701-48c8-88ee-49a78d46d0a8
 # ╠═065c7811-0203-44fd-bb88-49946d674fd8
 # ╠═39a46471-663e-4395-b4bd-2380932d7189
+# ╠═ed9e2e1f-cebe-4a0b-969a-93d344578a28
+# ╠═b7854b76-eb0e-4d78-9bd6-9b0258a725ba
+# ╠═3e1d6c7b-41b5-435a-88f0-b0f3ab442662
+# ╠═2f4cdc08-27a1-4e92-8709-8021cae771c6
+# ╠═a9d1b364-65d5-44bc-965f-61a5f2b47a36
+# ╠═0c5d6cdb-7255-4712-92bc-eb3f3182d7eb
+# ╠═e4461dae-f648-4c34-a963-17dfb3f0220a
+# ╠═d1588e94-a450-4d8c-9811-69873deed252
+# ╠═0b6bce6e-22e5-4a8d-95f4-4708deb07cdd
+# ╠═4561f42b-2658-4c74-b642-6b2f92cd3598
+# ╠═27a59d46-f838-4d8f-9c66-a7edc1eea8dd
+# ╠═504ed9c8-ff54-4a25-84b9-d832fc435486
+# ╠═873bcfa3-a989-426a-8402-6b81a8bbb46a
+# ╠═b10b4931-e60f-4791-8dae-d6d58cd71e66
+# ╠═b14e7de2-b91f-40d9-8aa8-237d82f25286
+# ╠═1fd645ee-0680-4de8-b1eb-1c94c2069ae9
+# ╠═c3286275-3cff-4771-9ae6-78a960ba0484
 # ╟─00000000-0000-0000-0000-000000000001
 # ╟─00000000-0000-0000-0000-000000000002