Plots.jlを使って週毎の歩数のヒストグラムを表示してみた

ヒストグラムでの表示

ソースコード

using Plots

x = [
    49599, 38582, 56619, 48610, 50834, 49428,
    68770, 59745, 72272, 77897,
    93871, 74859, 77991, 49253,
    53963, 67515, 77303, 91320,
    51489, 74508, 74791, 74093, 61799,
    64005, 67184, 60391, 84259,
    49561, 59995, 68221, 68529, 58643
]

histogram(x, label="Weely step count")

savefig("histogram.png")

出力

平均、分散、歪度、尖度

ソースコード

x = [
    49599, 38582, 56619, 48610, 50834, 49428,
    68770, 59745, 72272, 77897,
    93871, 74859, 77991, 49253,
    53963, 67515, 77303, 91320,
    51489, 74508, 74791, 74093, 61799,
    64005, 67184, 60391, 84259,
    49561, 59995, 68221, 68529, 58643
]

μ = sum(x) / length(x)
σ = √sum((x .- μ) .^2) / length(x)

skewness = sum((x .- μ) ./σ .^3) / length(x)
kurtosis =  sum((x .- μ) ./σ .^4) / length(x) - 3

@show μ, σ
@show skewness, kurtosis

出力

yuu@penguin:~/src/eikou/0425$ julia statistic.jl 
(μ, σ) = (64871.84375, 2314.120320475039)
(skewness, kurtosis) = (-9.926167350636332e-24, -3.0)