similar to: GEE models

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

Thank you very much Christophe by your advice. However, my situation was the opposite to yours, because I'm behind a proxy but my sudo account couldn't see it, and when looked for the proxy configuration of the sudo account with : $ sudo printenv | grep proxy the result was empty. Finally, thanks to all the advice received in this thread, I could fix the problem by editing the

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

Hi, I'm actually trying to rank a set of candidate models with an information criterion (AICc, QIC, BIC). The problem I have is that I use mixed-effect cox regression only available with the package {coxme} (see the example below). #Model1 >spring.cox <- coxme (Surv(start, stop, Real_rand) ~ strata(Paired)+R4+R3+R2+(R3|Individual), spring) I've already found some explications in

Hello, I'm having a look at section 13.1 from SIP Connect v2 doc (see [1]). It refers to RFC6442 which gives the following example (sorry for its length): INVITE sips:bob at biloxi.example.com SIP/2.0 Via: SIPS/2.0/TLS pc33.atlanta.example.com;branch=z9hG4bK74bf9 Max-Forwards: 70 To: Bob <sips:bob at biloxi.example.com> From: Alice <sips:alice at

I realize that matrix indexing has been addressed in various flavors, but I''m stumped and didn''t find anything in the archives. It''s not clear if it is an igraph issue or a more general problem. Thanks in advance for your patience. I am using igraph to read a gml file (http://www-personal.umich.edu/~mejn/netdata/football.zip ). The gml file contains vertex attributes

Dear Christophe, Yes. It was helpful indeed. The four steps you outlined successfully helped me to install the newest version of R. I have been fiddling with this for days and I am quite grateful to all who contributed. I formally posted to r-help@r-project.org before somebody suggested posting to r-sig-debian@r-project.org. Best regards Ogbos 2011/5/31 Christophe Bonenfant

Dear Christophe, Yes. It was helpful indeed. The four steps you outlined successfully helped me to install the newest version of R. I have been fiddling with this for days and I am quite grateful to all who contributed. I formally posted to r-help@r-project.org before somebody suggested posting to r-sig-debian@r-project.org. Best regards Ogbos 2011/5/31 Christophe Bonenfant

Dear list, Thanks for your time. I was surprised. While I was still waiting for your response, I entered R and discovered that it has been updated!!!! I hurried logged out and went home since I had overstayed battling with it. I wonder how the update happened with the numerous errors that were reported. I wish I will be able to find out, otherwise, I am off with my analysis with the new version :)

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:

Server on the other sites, and get a response. However, no client has resolve the names of the Servers. I thought I could get around this by adding the following line into my smb.conf: name resolve order = hosts lmhosts wins bcast (though I'm sure there are better ways of achieving this!!) But this doesn't solve it! My Samba Server includes a WIN Server and a DHCP Server, which issues

I have compiled R version 1.3.1 on a Linux system. No errors or warnings are reported. However, "make check" fails thus: ===== make[5]: Leaving directory `/proj/informatics-build/GML/R/R-1.3.1/src/library' running code in base-Ex.R ... ../../bin/R --vanilla < base-Ex.R > base-Ex.Rout make[4]: *** [base-Ex.Rout] Error 139 ===== The base-Ex.Rout being produced at this

Hello, I am exporting a graph in "gml" format from igraph 0.5.5. When I open the file in Cytoscape v.2.7 I can't map the color attribute that I assign to the graph when I export it from igraph. Could anyone help me in this? Thanks. -- Warm Regards, Sandeep Amberkar BioQuant,BQ26, Im Neuenheimer Feld 267, D-69120,Heidelberg Tel: +49-6221-5451354 [[alternative HTML version

Hello, I have a problem with the princomp method, it seems stupid but I don't know how to handle it. I have a dataset with some regular data and some outliers. I want to calculate a PCA on the regular data and get the scores for all data, including the outliers. Is this possible on R? Thank you for helping!!! -- View this message in context:

Hi, I'm an experienced R user on windows, but new to linux. This may be an issue with how R (or Ubuntu) was install on my computer, but maybe not. I just installed Ubuntu 10.04, I have administrator privileges and am using R 2.11. Here goes: I'd like to change the location where R installs packages not included in the base installation (using install.packages()). Right now it is

Hi there, How may I smooth spline two vectors with the smoothing parameter selected by generalized maximum likelihood (GML) .? Thanks a lot. -- View this message in context: http://r.789695.n4.nabble.com/Smoothing-spline-with-smoothing-parameters-selected-by-generalized-maximum-likelihood-tp4210679p4210679.html Sent from the R help mailing list archive at Nabble.com.

--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

Dear R-list. I'm doing af logistic analyses using gml. The model explaines variations in Adverse events infections (0 og 1) using age as explanatory variable. model2d<-glm(formula=AEorSAEInfecBac~Age,family=binomial("logit"),data=emrisk) I want to get predictions with 95% confidence limits for age 30 and age 60. I've been reading the "google" and "search