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