データ解析のための統計モデリング入門をJuliaでやっていく(2.1)

2.1

using RData
using StatsBase

using Plots
gr()

resource_dir = normpath(joinpath(@__FILE__, "..", "..", "resources"))
data_path = joinpath(resource_dir, "data.RData")

objs = load(data_path)
data = get(objs, "data", [])
println(data)

# データ数
println("length", length(data))

# 標本平均、最小値、最大値、四分位数
println("stats", summarystats(data))

# ヒストグラム
histogram(data;bins=range(0, maximum(data), step=1.0))
png("sample.png")

# 標本分散
println("var", StatsBase.var(data))

# 標準偏差
println("std", StatsBase.std(data))

println("std(sqrt)", sqrt(StatsBase.var(data)))