Questions tagged [mumin]

MuMIn (**Mu**lti**m**odel **In**ference) is an R package for model selection and model averaging based on information criteria (AICc and alike).

90 questions
1
vote
1 answer

Strange error with Dredge: MuMIn

I am wondering whether anyone could help me trouble shoot this error from a task with dredge package: Fixed term is "(Intercept)" |============================================================================ …
Joe
  • 347
  • 5
  • 20
1
vote
2 answers

Extract average model from MuMIn for latex output

I'm trying to extract two different averaged models from MuMIn for output to latex via texreg or stargazer. I'd like to have one table where I can compare two species' response to different sets of abiotic variables, that looks the same as one…
1
vote
1 answer

dredge (MuMIn package) error message: no applicable method for 'logLik'

I came across an error I haven't been able to find any information on, and I haven't been able to figure it out, yet. Here is an example: x1 <- rnorm(100,0,1); x2 <- rnorm(100,5,1); x3 <- rnorm(100,8,4) x4 <- rnorm(100,3,3); x5 <- rnorm(100,1,2);…
Jota
  • 17,281
  • 7
  • 63
  • 93
1
vote
2 answers

Rank model with QAIC using MuMIn

I need to rank candidate model using the QAIC criterion. This is what I've tried: library(MuMIn) model_global <- glm(vs ~ mpg + disp + wt, family = quasibinomial, mtcars) model_1 <- glm(vs ~ mpg, family = quasibinomial, mtcars) model_2 <- glm(vs…
luciano
  • 13,158
  • 36
  • 90
  • 130
1
vote
1 answer

Can no longer return beta coefficients using MuMIn package in R

When using the R package MuMIn (v1.9.5), and the command dredge, the beta = TRUE switch no longer functions as expected. # load Cement data set with MuMIn package data(Cement) # build global model fmd <- lm(y ~ X1 + X2 + X3 + X4, data = Cement) #…
Prophet60091
  • 589
  • 9
  • 23
0
votes
0 answers

Problem with dredge function in MuMIn "result is empty" and "missing values in object" despite dataset seeming to be complete

I get the following error message when I run the dredge function : Fixed term is "(Intercept)" Error in dredge(global.model = ChangeModel, evaluate = TRUE, rank = "AICc") : result is empty In addition: There were 50 or more warnings (use…
0
votes
0 answers

How to always keep two variables to interact with other variables simultaneously in model selection using dredge() in R?

The title may not be very explicit, I will use code to illustrate. My actual data got a four-way interaction, here I simplified it to three-way: colnames(iris) <- c("y", "a", "b", "c") m1 <- lm(y ~ a * b * (sin(c) + cos(c)), na.action =…
Juan_814
  • 51
  • 8
0
votes
0 answers

Error in emmeans for averaging object (mth$objs[[1]])(object, trms, xlev, grid, ...): Unable to match model terms

emmeans can't evaluate averaging object The problem is that emmeans gives an error when evaluating an averaging object form the MuMIn package, even when it says it should in this link I have been trying to debug this for a couple of days with no…
Derek Corcoran
  • 3,930
  • 2
  • 25
  • 54
0
votes
1 answer

How to interpret the output of MuMIn model.avg on GAMs

Say I have a series of GAMs that I would like to average together using MuMIn. How do I go about interpreting the results of the averaged smoothers? Why are there numbers after each smoother…
spops
  • 572
  • 1
  • 7
  • 25
0
votes
1 answer

Getting "+" sign in the results of MuMIn :: dredge

I am trying to MuMIn::dredge linear mixed-effect models lme4::lmer with categorical/continuous variables, the code is as follows: # Selection of variables of interest sig<-c("Age", "Sex", "BMI", "(1|HID)", "h_age", "h", "h_g", "smk_hs") # Model…
0
votes
1 answer

No difference in R package MuMIn model.avg output when using delta < 2 vs. 95% AICc weight?

I've developed a set of candidate models (several binomial GLMs with a null and two global models) and I am trying to use MuMIn to do model averaging on the top delta 2 AIC models. Here is a subset of the candidate models fit1 <- glm(formula = Y ~…
Tiff-D
  • 25
  • 5
0
votes
1 answer

Singular in MuMIn::dredge

I'm about to run a glm model to test my hypothesis. Once I ran the model I got a singularity warning massage: glmakde_mdDredge <- lmer(log(akdearea) ~ 1 + age_binary + sex + trimester_number + Year + trackNightnum + mean_max_displacement.s +…
0
votes
2 answers

Creating an AIC model selection table with model included

I have been using code as explained here:How to Create AIC Model Selection Table in R in LaTex format? However, I am running into an error with my most recent dataset and I think it has to do with having a larger number of model terms (I have 11…
Amanda Goldberg
  • 169
  • 1
  • 11
0
votes
1 answer

MuMin dredge gam with random effect

I am using the dredge function from the MuMin package for a gam with a random effect: The global model using bam from the package mgcv is below. The random effect iss(StationID, bs = 're') gam_global<-bam(Numberpertow ~…
user41509
  • 978
  • 1
  • 10
  • 31
0
votes
1 answer

Dredge on GAM returns NA results for logLik, AICc, delta and weight

I have a global moel with two smoothed terms and two random effects. mRI_qb <- gam(ri ~ SE_score + s(deg, k=7) + s(gs, k=7) + TL + species + sex + season + year + s(code, bs = 're') + s(station, bs = 're'), …
mikejwilliamson
  • 405
  • 1
  • 7
  • 17