Multinomial Random Effects
This code solves multinomial mixed effects models.
devtools::install_github('bdilday/mnre')
fit a binomial model using a formula interface
library(mnre)ev = mnre_simulate_multinomial_data_factors(nfct=2, K_class = 2, nlev=50, nobs=20000)mnre_mod = mnre_fit(y ~ 1 + (1|fct01) + (1|fct02), data=ev$fr, verbose=0)mnre_mod$theta[1] 0.9926043 0.9537682
fit a binomial model using lower-levelnd_min_fun
library(mnre)ev = mnre_simulate_multinomial_data_factors(nfct=2, K_class = 2, nlev=50, nobs=20000)nf = mnre::nd_min_fun(ev)ans = optim(c(1,1), nf, method = "L-BFGS", lower=1e-8)print(ans$par)[1] 0.9926043 0.9537682
Compare tolme4
glmer_mod <- glmer(ev$frm, data=ev$fr, family='binomial', nAGQ=0)glmer_mod@theta[1] 0.9925264 0.9537615
Fit a multinomial model
ev = mnre_simulate_multinomial_data_factors(nfct=1, K_class = 3, nlev=50, nobs=20000)ev$verbose = 0nf = mnre::nd_min_fun(ev)ans = optim(c(1, 1), nf, method = "L-BFGS", lower=1e-8)print(ans$par)[1] 0.9351018 1.0758567