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