similar to: Any package for best subset selection on random effects model?

Hi, I know leaps() computes the best subset selection for linear model, and the bestglm() computes the best subset selection for generalized linear model. Is there any package for best subset selection on random effects model, or mixed effects model? Thank you so much. Tao [[alternative HTML version deleted]]

Hi All: the package "MuMIn" can be used to select the model based on AIC or AICc. The code is as follows: data(Cement) lm1 <- lm(y ~ ., data = Cement) dd <- dredge(lm1,rank="AIC") print(dd) If I want to select the model by BIC, what code do I need to use? And when to select the best model based on AIC, what the differences between the function "dredge" in

Hi all, I'm having trouble with the model generating 'dredge' function in the MuMIn 'Multi-model Inference' package. Here's the script: globalmodel<- glm(TB~lat+protocol+tested+ streams+goats+hay+cattle+deer, family="binomial") chat<- deviance(globalmodel)/59 #There we 59 residual degrees of freedom in this global model. models<- dredge(globalmodel,

Dear list, I am trying to use MuMIn to compare all possible mixed models using the dredge function on binomial data but I am getting an error message that I cannot decode. This error only occurs when I use glmer. When I use an lmer analysis on a different response variable every works great. Example using a simplified glmer model global model: mod<- glmer(cbind(st$X2.REP.LIVE,

Hello, I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In the AIC list of all models, each model is differentiated by model number. Please kindly advise if it is possible to find the corresponding explanatory variable(s) for the model number. Q2 error message I tried to display sub-model with only

Hi All, I'm having trouble with the automatized model generation (dredge) function in the MuMIn package. I'm trying to use it to automatically generate subsets of models from a global cox proportional hazards model, and rank them based on AICc. These seems like it's possible, and the Mumin documentation says that coxph is supported. However, when I run the code (see below), it gives

The subject says it all really. Question 1. Here is some code created to illustrate my problem, can anyone spot where I'm going wrong? Question 2. The reason I'm following a manual specification of models relates to the fact that in reality I am using mgcv::gam, and I'm not aware that dredge is able to separate individual smooth terms out of say s(a,b). Hence an additional request,

I recently found the "bestglm" package while trolling CRAN for a function to save some keystrokes with simple (k<10) nested models. I am interested in using the best of several nested linear models which I have explored using lm(), glm() and now bestglm(). When I compare the AIC values for the 'best' candidate model I get the same AIC values with lm() and glm() but a

Hi All: Recently I am researching my dissertation about the quantile model selection by bayesian approach. I have the dependent variable(return) and 16 independent variables and I need to select the best variable for each quantile of return. And the method I used is the bayesian approach, which is based on calculating the posterior distibution of model identifier. In other words, I need to obtain

Hi! I am trying to know which habitat variables most affect bird counts in a radius of 100m. I obtained bird counts in 2751 spatial points, and measured percentage of 21 habitat variables in these points. I applied a mixed model using the "lmer" function to these data, but I do not know how to select the best model using AIC here. Is there a way to do this automatically with R?

Hello, I hope that you can bare with me. I am new to models, but I think I have a pretty godd understanding of how to run them now, including how to use AICc and Anova. The issue is that I have many factors that I wish to compare so doing each one at a time would take forever. I came across the MuMIn package and I was so excited, however I am getting an error message and i don't know why.

Hi all, I have been trying to run a bestglm in R for a while now and am struggling to get it to run. When I thought I had succeeded, the "output" it gave me was "NULL" and that's it. Below is my code: bestglmtest<-read.table("C:\\Documents and Settings\\clyons\\Desktop\\bestglmtest.txt",header=T,row.names=1) > names(bestglmtest) [1] "acc"

Dear List, I want to ask a AIC question based on package library(MuMIn) The relative importance of 16 explanatory variables are assessed using delta AIC in a generalized linear model. Please kindly advise if it is possible to show models with any two only certain variables. Thank you. Elaine I asked a similar question and got a great help for models with only one variable as below.

Hola a todos, He realizado un dredge (para obtener todos los modelos GAM posibles a parir de un full model), luego he seleccionado un confidence set (los modelos que no se diferencian en 2 en AIC) y he hecho un model averaging con ese confidence set. Ahora me gustaría aplicar ese modelo "average" ajustado sobre otro set de datos pero no se como especificar en R que use el mismo modelo

Hello, I'd like to sum the weights of each independent variable across linear models that have been evaluated using AIC. For example: > library(MuMIn) > data(Cement) > lm1 <- lm(y ~ ., data = Cement) > dd <- dredge(lm1, beta = TRUE, eval = TRUE, rank = "AICc") > get.models(dd, subset = delta <4) There are 5 models with a Delta AIC Score of

Hello R community, I am attempting to run multiple logistic regressions (multinomial, via package 'nnet'), with Automated Model Selection (dredge, package 'MuMIn'). The aim is to reduce the number of predictor variables by assessing relative performance of each variable, which can be done in a coarse fashion using the Automated Model Selection option in package 'MuMIn'

Hello! I've been running a looped AIC analysis using several modules including ape, nlme, and MuMIn, and during one particularly long analysis, R (ver 2.14.12) crashes several minutes into the routine with the simple message "R for windows GUI front-end has stopped working". I'm using a brand new laptop with Windows 7, i7 processor, 8GB RAM. I've tried it on both the

Hello I am not able to use the MuMIn package (version 1.7.7) for multimodel inference with a GLM Negative Binomial model (It does work when I use GLM Poisson). The GLM Negative Binomial gives the following error statement: Error in get.models(NBModel, subset = delta < 4) : object has no 'calls' attribute Here is the unsuccessful Negative Binomial code. > > BirdNegBin

Estimados amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si

Estimados amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si