データ解析のための統計モデリング入門をJuliaでやっていく(3.5)
3.5
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_f = glm(@formula(y ~ f), data, Poisson(), LogLink()) println(fit_f) println(loglikelihood(fit_f)) fit_all = glm(@formula(y ~ x + f), data, Poisson(), LogLink()) println(fit_all) println(loglikelihood(fit_all))