similar to: Ranking submodels by AIC (more general question)

I know this a newbie question, but I've only just started using AIC for model comparison and after a bunch of different keyword searches I've failed to find a page laying out what the differences are between the AIC scores assigned by AIC() and mle.aic() using default settings. I started by using mle.aic() to find the best submodels, but then I wanted to also be able to make comparisons

Folks, I'm using stepAIC(MASS) to do some automated, exploratory, model selection for binomial and Poisson glm models in R 1.3. Because I wanted to experiment with the small-sample correction AICc, I dug around in the code for the functions glm.fit stepAIC dropterm.glm addterm.glm extractAIC.glm and came across something I just don't understand. stepAIC() passes dropterm.glm() a

Hello, I have a question pertaining to the stepwise regression which I am trying to perform. I have a data set in which I have 14 predictor variables accompanying my response variable. I am not sure what the difference is between the function "mle.stepwise" found in the wle package and the functions "step" or "stepAIC"? When would one use

Hi all -- I'm seeing the following commands reliably produce a crash in RGui, version 2.0.1, for both my home and office machine: > rm(list = ls(all = TRUE)); > load("dataset.R"); > library("wle"); > data.wle = wle.lm(abortion ~ year * lib.con + age + gender + + urbanism + census + income + church.att + children + educ + + religion.imp, data =

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 Sacha, On second thought, perhaps this is more the direction that you want ... X2 = cbind(X_train,y_train) colnames(X2)[3] = "y" regr2<-randomForest(y~x1+x2, data=X2,maxnodes=10, ntree=10) regr regr2 #Make prediction predictions= predict(regr, X_test) predictions2= predict(regr2, X_test) HTH, Eric On Tue, May 9, 2023 at 6:40?AM Eric Berger <ericjberger at gmail.com>

Subject: Eliminate repeated components from a vector X-Mailer: VM 7.00 under 21.4 (patch 6) "Common Lisp" XEmacs Lucid Reply-To: fjmolina at lbl.gov FCC: /home/f/.xemacs/mail/sent Does anyone know how I can eliminate repeated elements from a vector?

Dear [R]-list, I'm looking for a classic equivalent of the wle.gamma function (library wle) that estimate robustly the shape and the scale parameters of gamma data. I have a vector of iid gamma rv : >data=rgamma(100,shape=10,scale=3) and a vector of their weights: >weights=c(rep(.5/70,70),rep(.25/20,20),rep(.25/10,10)) and want to estimate the scale and shape of the gamma

I'm trying to develop a linear model for crop productivity based on variables published as part of the SSURGO database released by the USDA. My default is to just run lm() with continuous predictor variables as numeric, and discrete predictor variables as factors, but some of the discrete variables are ordinal (e.g. drainage class, which ranges from excessively drained to excessively poorly

I would like to try a likelihood ratio test in place of waldtest. Ideally I'd like to provide two glm models, the second a submodel of the first, in the style of lrt (http://www.pik-potsdam.de/~hrust/tools/farismahelp/lrt.html). [lrt takes farimsa objects] Does anyone know of such a likelihood ratio test? Chris Elsaesser, PhD Principal Scientist, Machine Learning SPADAC Inc. 7921

Hi, all When I am using mle.cv(wle), I find a interesting problem: I can't do leave-one-out cross-validation with mle.cv(wle). I will illustrate the problem as following: > xx=matrix(rnorm(20*3),ncol=3) > bb=c(1,2,0) > yy=xx%*%bb+rnorm(20,0,0.001)+0 > summary(mle.cv(yy~xx,split=nrow(xx)-1,monte.carlo=2*nrow(xx),verbose=T), num.max=1)[[1]] mle.cv: dimension of the split subsample

Buenas tardes, Estoy intentando ajustar distribuciones utilizando un vector de ponderación en los datos (Weighted Likelihood). ¿Existen paquetes en R que resuelven esto? He mirado ya el paquete "wle" pero no me permite introducir los pesos mediante los cuales ponderar los datos. En un primer momento, se me ha ocurrido realizar lo siguiente: repetir cada elemento del vector datos

I haven't seen this cryptic warning before: > update.views('Robust') Warning message: In update.views("Robust") : The following packages are not available: covRobust, distr, FRB, MASS, mblm, multinomRob, mvoutlier, quantreg, RandVar, rgam, RobAStBase, robfilter, RobLox, RobRex, robust, RobustAFT, robustbase, ROptEst, ROptRegTS, rrcov, sandwich, wle >

I would consider this is a question for a statistics forum such as stats.stackexchange.com, not R-help, which is about R programming. They do sometimes intersect, as here, but I think you need to *understand what you're doing* before you write the R code to do it. Obviously, IMO. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and

======================================================================= Simple CRAN packages which do not compile without modifications (all others do) ======================================================================= -- akima /usr/bin/ld: multiple definitions of symbol _idlc_ -- fracdiff /usr/bin/ld: multiple definitions of symbol _gammfd_ (and others) -- odesolve --

Dear R-helpers, In the case of two categorical factors, say a and b, once I have fixed the constrasts, the model matrix is set according to these contrasts with "lm", and the t-tests for the significance of the parameters provided by "summary" indeed concern the comparison of the model with each submodel obtained by removing the corresponding column of the model matrix.

Dear R-users, I would like to announce the release of the new package JMbayes available from CRAN (http://CRAN.R-project.org/package=JMbayes). This package fits shared parameter models for the joint modeling of normal longitudinal responses and event times under a Bayesian approach using JAGS, WinBUGS or OpenBUGS. The package has a single model-fitting function called jointModelBayes(),

Dear R-users, I would like to announce the release of the new package JMbayes available from CRAN (http://CRAN.R-project.org/package=JMbayes). This package fits shared parameter models for the joint modeling of normal longitudinal responses and event times under a Bayesian approach using JAGS, WinBUGS or OpenBUGS. The package has a single model-fitting function called jointModelBayes(),

I have re-built all the compiled packages for R 1.2.0 for Windows, and the set available in CRAN/bin/windows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5

I have re-built all the compiled packages for R 1.2.0 for Windows, and the set available in CRAN/bin/windows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5