Stack Overflow
Stack Overflow

Questions tagged [glmulti]

21 questions
3
votes
1 answer

glmulti syntax for mixed effects logistic regression in lme4

I am trying to compare the AIC (or AICc) values of various mixed effects logistic (or for some outcomes gamma) models containing all the variables in my dataset, a simplified version of which can be downloaded here:…
Mel
  • 700
  • 6
  • 31
2
votes
1 answer

metafor coupled with glmulti: How to exclude undesirable interactions

I am fitting the following: rma.glmulti.ran <- function (formula, data, random, ...) { rma.mv(as.formula(paste(deparse(formula))), Variance, random= ~1 | Experiment, data = msa, method="REML", ...) } msa_res <- glmulti(MSA ~…
Gabriele Midolo
  • 231
  • 1
  • 9
2
votes
2 answers

fitting quasi family using glmulti?

I used the glmulti function in the glmulti package to obtain the best glm model for poisson error distributed data. No problems there. Once I had obtained the best model, I used the Chi-square test to obtain p-values and test statistics for each of…
Candice-Lee Lyons
  • 21
  • 3
1
vote
0 answers

glmulti wrapper with offset term

I am using glmulti to assess the relative importance of several predictors on counts. As I am using glm.nb I am using the wrapper. I would like to add an offset term, but I am not sure where it should go or, if it is even possible. I tried adding it…
Kaddy Kuchens
  • 11
  • 2
1
vote
0 answers

Syntax for glmer function for use with glmulti?

Using glmer, I can run a logistic regression mixed model just fine. But when I try to do the same using glmulti, I get errors (described below). I think the problem is with the function I am specifying for use in glmulti. I want a function that…
Emily
  • 470
  • 5
  • 16
1
vote
1 answer

Does the glmulti function (from the gmulti package) need a set.seed value?

I am using glmulti to select a set of candidate generalized linear models and my variable importance values and 'best' model keep changing each time I run the model. I am struggling to understand why this is, does glmulti need a set.seed value to…
mrpargeter
  • 339
  • 3
  • 12
1
vote
0 answers

What interaction effect glmulti generates in A*B*C*D and A+B+C+D?

I am using glmulti to run hierarchical linear models and select the best model. I have 4 predictors (A, B, C, D) to the DV, and my goal is to run all main effect models plus all combination of interaction effects (i.e., A:B, A:C, A:D). How do the…
Lumos
  • 1,303
  • 2
  • 17
  • 32
1
vote
1 answer

parallel computing for glmulti in R in Windows

as per glmulti package document, chunks are arguments for using multi CPUs. when using exhaustive screening. But, even when I put 4 in both chunk and chunks and method='h' with family='binomial', R only use a single core. the function I…
Hemant Rupani
  • 145
  • 10
0
votes
1 answer

Using glmulti with GEE models

I can use glmulti with lmer models, but am struggling to make it work for a GEE. Here is an example showing an lmer version which works, and a GEE version I'm stuck on. library(tidyverse) library(lme4) library(geepack) library(glmulti) data <-…
Jeff
  • 57
  • 6
0
votes
1 answer

What code should I use for my GLMM, glmer, lmer?

I Want to use the dependent variable "Herps" (count data) to understand what response variables are important influencing herpetofauna species richness. The response variables include Total_Saplings, Total_Understorey_Vegetation,…
Hefin
  • 1
  • 1
0
votes
0 answers

Model averaging for rma.mv objects not working using AICcmodavg package

I have a set of candidate models that I want to average so that I get one estimate and variance for each predictor, weighted by the weight of the model. I've used the following code for glm objects but it's not working for rma.mv models (from…
Lauren Diaz
  • 75
  • 5
0
votes
1 answer

How to exclude interaction term from glmulti?

I am fitting the following: offspring.survival_STAB <-glmulti(ST~ STAB + Age + BS+ BSP + Sex, data=ST, exclude="BS:BSP", level=2, …
Ron Flatau
  • 1
0
votes
0 answers

Quit model fitting after no improvement

I am using glmulti to fit a logistic regression model. I have 6 predictor variables, with one predictor variable having 84 levels. Consequently, it is taking an exceptionally long time to fit the model. I thought I structured the code to stop…
user17047272
0
votes
0 answers

Problems using glmulti package in RStudio to predict best model

So I am trying to use the "glmulti" package find the best combination of variables (or best model) for my response variable Bio_class. Bio_class is a categorical variable that contains 5 outcomes (4 species and 1 absence option). I have 56 numerical…
SGilk
  • 1
  • 2
0
votes
0 answers

How to extract the quasi-BIC (QBIC) from a Quasipoisson model?

I would like to use GLMULTI for feature-selection by fitting a Quasipoisson model using quasi-BIC (QBIC) as information criterion. I use the following code: glmulti(y ~ ., data=data, level=2, fitfunction=glm, family = quasipoisson,…
Tea Tree
  • 882
  • 11
  • 26
1
2