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

3.4

using DataFrames
using CSV

using GLM

using StatsPlots


resource_dir = normpath(joinpath(@__FILE__, "..", "..", "resources"))
data_path = joinpath(resource_dir, "data3a.csv")

data = DataFrame(CSV.File(data_path))
fit = glm(@formula(y ~ x), data, Poisson(), LogLink())
println(fit)

println(loglikelihood(fit))

xx = DataFrame(x=range(minimum(data.x), maximum(data.x), length=100))
yy = predict(fit, xx)

@df data scatter(:x, :y, group=:f)
plot!(xx.x, yy)
png("predict.png")