search

Nosferican / Econometrics.jl

Econometrics in Julia
https://nosferican.github.io/Econometrics.jl/dev
ISC License
69 stars 19 forks source link

Ordinal Logistic Regression Wrong When Missing Levels #11

Closed Nosferican closed 5 years ago

Nosferican commented 5 years ago

Describe the bug Wrong estimates.

To Reproduce

using Econometrics, RDatasets
data = dataset("Ecdat", "Kakadu") |>
    (data -> select(data, [:RecParks, :Sex, :Age, :Schooling]))
data.RecParks = convert(Vector{Int}, data.RecParks)
data.RecParks = levels!(categorical(data.RecParks, ordered = true), collect(1:5))
data.RecParks = convert(Vector{Int}, data.RecParks) .- 1
data = data[data.RecParks .≠ 0,:]
data.RecParks = levels!(categorical(data.RecParks, ordered = true), collect(1:4))
model = fit(EconometricModel,
        @formula(RecParks ~ Age + Sex + Schooling),
        data)
Nosferican commented 5 years ago

Closed by #19.