similar to: Model ranking (AICc, BIC, QIC) with coxme regression

Dear R Users, May be this message should be directy send to Douglas Bates ... I just want to know if I can use the AIC value given in the output of an lmer model to classify my logistic models. I heard that the AIC value given in GLIMMIX output (SAS) is false because it come from a calculation based on pseudo-likelyhood. Is it the same for lmer ??? thanks, Arnaud Arnaud MOSNIER Biologiste

Hello, Is it possible to derive Pan's QIC (2001 Biometrics 57:120) from either a fitted gee() object in the gee package or from a geese() fit in the geepack package? If so, would anyone be kind enough to provide me with code to do so? I realize that QIC is part of the output from yags() but I would like to use one of the other functions. Thanks. Richard

Hi, I'm using 'geepack' to run Generalized Estimating Equations. I'm aware that I can use anova to compare two models, but would it be possible to test QIC on R? It seems that there were similar questions a couple of years ago, but the question has not been answered yet. I'd appreciate if someone could show me the code! Thank you, Sachi [[alternative HTML version

Dear all, I'm trying to compare two sets of variables, the first set is composed exclusively of numerical variables and the second regroups factors and numerical variables. I can't use a Manova because of this inclusion of numerical variables in the second set. The solution should be to perform a Mancova, but I didn't find any package that allow this type of test. I've already

Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are respected. What it means? Is it an error in the script? Commands and output: > permutest(rda.ind, perm=999, first=TRUE) Permutation test for rda  Call:

Dear all, I'm trying to get a Brownian bridge kernel (kernelbb) for each combination of two consecutive animal locations (see commands below) and put them, with a loop, inside a list. It works well at the beginning but after 42 runs, it appears the following warning : >Error in seq.default(yli[1], yli[2], by = diff(xg[1:2])) :  >  invalid (to - from)/by in seq(.) I looked at the

Dear R-helpers I have generated a suite of GLMs. To select the best model for each set, I am using the meta-analysis approach of de Luna and Skouras (Scand J Statist 30:113-128). Simply put, I am calculating AIC, AICc, BIC, etc., and then using whichever criterion minimizes APE (Accumulated Prediction Error from cross-validations on all model sets) to select models. My problem arises where I

I'm using the package geepack to fit GEE models. Does anyone know of methods for add1 and drop1 for a 'geeglm' model object, or perhaps a method for extractAIC based on the QIC of Pan 2001? I see there has been some mention of this on R-help a few years ago (RSiteSearch("QIC")). The package does provide an anova method for its model objects, and update() seems to work:

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

--On Monday, May 11, 2015 02:26:17 PM -0700 John R Pierce <pierce at hogranch.com> wrote: > never met a unix that didn't come with Perl already installed, or as a > base option SunOS-4 :) Didn't have emacs, either, nor an ANSI-C compiler. And the OS came on QIC-150 tape (ie: 150 MB total capacity). Not that that defeats the argument ... Devin

Hi, When doing a two-ways anova in R and comparing my same results with an SPSS output, I noticed that R calculated type I Sum of Squares. Is it possible to use Type III Sum of Squares? Thanks, S?bastien Plante Institut des Sciences de la Mer de Rimouski (ISMER) -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Besides this change on vfs objects, i delete winbindd_idmap.tdb and restarted winbind. Change smb.conf: vfs objects = acl_xattr full_audit full_audit:facility = LOCAL7 full_audit:priority = NOTICE full_audit:prefix = %u|%I|%m|%S full_audit:success = pread pwrite rename rmdir mkdir unlink chown full_audit:failure = none Here are my configurations: #

Rowland, now two members are working great. Past 24 hours and they didn´t lost uid/gid. root at fs1:/var/lib/samba# wbinfo -i edsontadeu edsontadeu:*:11106:10513::/home/DOMAIN/edsontadeu:/bin/false root at fs2:/var/lib/samba# wbinfo -i edsontadeu edsontadeu:*:11106:10513::/home/DOMAIN/edsontadeu:/bin/false Thanks! 2016-10-06 9:24 GMT-03:00 Rowland Penny <rpenny at samba.org>: >

Hello R-users: I am using "yags" for fitting GEE which is giving me the same result as "Proc GENMOD". Now I have couple of questions related to yags output. (By the way, someone told me to run the geeglm for the same analysis and I did run but did not get the same result as of genmod and don't know how to correct the geeglm codes so that all three will be same!)

Hi I get the following error when I try and get the AICc for a gls regression using qpcR: > AICc(gls1) Loading required package: nlme Error in n/(n - p - 1) : 'n' is missing My gls is like this: > gls1 Generalized least squares fit by REML Model: thercarnmax ~ therherbmax Data: NULL Log-restricted-likelihood: 2.328125 Coefficients: (Intercept) therherbmax 1.6441405

I have a case-control study that I'm analysing using the conditional logistic regression function clogit from the survival package. I would like to calculate the AICc of the models I fit using clogit. I have a variety of scripts that can calculate AICc for models with a logLik method, but clogit does not appear to use this method. Is there a way I can calculate AICc from clogit in R? Many

Hello, I am modelling in glmmADMB count data (I´m using a negative binomial distribution to avoid possitive overdispersion) with four fixed and one random effect. I´m also using MuMIn package to calculate the AICc and also to model averaging using the function dredge. What I do not understand is why dredge calculates a different value of the AICc and degrees of freedom than the function AICc

How can I extract K (number of parameters) from an AIC calculation, both to report K itself and to calculate AICc? I'm aware of the conversion from AIC -> AICc, where AICc = AIC + 2K(K+1)/(n-K-1), but not sure of how K is calculated or how to extract that value from either an AIC or logLik calculation. This is probably more of a basic statistics question than an R question, but I thank