データ解析のための統計モデリング入門を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)))