similar to: Stepwise GLM selection by LRT?

Dear all, i am trying to redo the 'eyestudy' analysis presented on the site http://www.ats.ucla.edu/stat/stata/faq/relative_risk.htm with R (1.9.0), with special interest in the section on "relative risk estimation by poisson regression with robust error variance". so i guess rlm is the function to use. but what is its equivalent to the glm's argument "family"

Dear All, I'm analysing a negative binomial dataset from a population-based study. Many covariates were determined on household level, so all members of a household have the same value for those covariates. In STATA, there seems to be an option for 'clustered analysis' for neg-bin regression. Does an equivalent exist for R(MASS)'s glm.nb or a comparable function? Many thanks for

Dear R-list, maybe someone can help me with the following mixed-effects models problem, as I am unable to figure it out with the 'nlme-bible'. I would like to fit (in R, obviously) a so-called animal model (google e. g. "Heritability and genetic constraints of life-history" by Pettay et al.) to estimate the variance component that is due to genetic effects. The covariances of

Dear List, recently tried to reproduce the results of some custom model selection function after updating R, which unfortunately failed. However, I ultimately found the issue to be that testing with pchisq() in drop1() seems to have changed. In the below example, earlier versions (e.g. R 2.4.1) produce a missing P-value for the variable x, while newer versions (e.g. R 2.7.1) produce 0 (2.2e-16).

> Message: 24 > Date: Tue, 31 Jan 2006 18:22:52 +0000 > From: "Lutz Ph. Breitling" <lutz.breitling at gmail.com> > Subject: [R] Mixed-effects models / heterogeneous covariances > To: r-help at stat.math.ethz.ch > Message-ID: > <2e38a1c80601311022i2e1be92doa60b80b50b69eb0c at mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Dear

Dear R members, I would like to know if there is a stepwise function that gives LRT value instaed of AIC? thanks in advance, Luciana Pozzi CENPAT-CONICET Argentina [[alternative HTML version deleted]]

configure uses clock_gettime to check whether -lrt is needed - and don''t check other functions. With glibc 2.17 clock_gettime is part of libc, so use timer_gettime instead, which is in -lrt in old and new versions of glibc. Signed-off-by: Olaf Hering <olaf@aepfle.de> --- configure | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/configure b/configure index

Hi there, This is not a strictly related R question; I apologize for this. I'm interest in simulate the sampling distribution of the LRT, testing for trend among the levels of some categorical variable X in a regression model. (in practice this is achieved by assigning scores to the levels of X and fitting such numeric variable). To simulate the sampling distribution the steps are: 1.

On Wed, May 1, 2013 at 09:21:17AM -0700, Matt Wilson wrote: > On Mon, Apr 29, 2013 at 03:40:13PM +0200, Olaf Hering wrote: > > Ping > > Makes sense. > > Acked-by: Matt Wilson <msw@amazon.com> > > > On Wed, Jan 16, Olaf Hering wrote: > > > > > Ping > > > > > > On Wed, Dec 05, Olaf Hering wrote: > > > > > > >

Hi! I am resending the following message about problems with utmp handling. * In the meantime I had some request in private mail from people asking whether I have new information. * The problem is still persistant in 2.9p2. * My own new investigations show, that the problem only appears with protocol 2, not with protocol 1, I therefore only started to note it when protocol 2 became the

# HG changeset patch # User Lutz.Dube # Node ID ca2c8c262dbaa6b2521cec60af7cd688dabdf850 # Parent 8992134dcfd0b9e1e86f4111e68a8aa48bd33c3c Values of cpu_weight and cpu_cap are lost after xend restart For managed domains in state ''halted'' I always get default values for cpu_cap / cpu_weight after xend restart. This is because the names of parameters differ between a SXP file to

On my setup with glibc-2.8, benchmark_residual linkage fails with undefines references to clock_gettime(). Adding -lrt fixes that. The following is a small patch for it. Regards. -- O.S. diff --git a/configure.ac b/configure.ac index 993ac33..392485e 100644 --- a/configure.ac +++ b/configure.ac @@ -437,6 +437,11 @@ if test x$enable_stack_smash_protection = "xyes" ; then

My CentOS box with: $ uname -r 2.6.18-164.6.1.el5 gets the file order wrong with: ls -lrt * . Two files that were last modified in the same second are shown in order that is the reverse of the order in which they were actually last modified. Note also that the later file was first created after the earlier file was last modified. Is the precision of this only to the second? Mike.

On Wed, Sep 13, 2000 at 09:13:26AM +0000, Graham Murray wrote: > On upgrading to openssl 0.9.6-beta1, I find that openssh 2.2.0p1 fails > to connect. I did some more experiments and also saw the problems. They occur when using a 0.9.6-beta client to connect to 0.9.5a and 0.9.6-beta servers. They also occur when using a 0.9.5a client connecting to a 0.9.6-beta server. Connections fail with

Hi! I have experimented a bit with the latest OpenSSH from the CVS archive. I could realize some connections succesfully, but I experienced performance problem during the connection phase. It seems, that the client needs quite some computer time just after debug: Got SSH2_MSG_KEX_DH_GEX_GROUP. By inserting test-printouts, I verifyed that the dh_gen_key(dh); call seems to take that long. On a HP

Hello !! I'm recently having a debate with my PhD supervisor regarding how to write the result of a likelihood ratio test in an article I'm about to submit. I analysed my data using "lme" mixed modelling. To get some p-values for my fixed effect I used model simplification and the typical output R gives looks like this: model2 = update ( model1,~.-factor A) anova (model1,

Hi! The handling of the "status" information in bsd-waitpid.c and bsd-nextstep.c seems to be bit odd. Patch attached. Best regards, Lutz -- Lutz Jaenicke Lutz.Jaenicke at aet.TU-Cottbus.DE BTU Cottbus http://www.aet.TU-Cottbus.DE/personen/jaenicke/ Lehrstuhl Allgemeine Elektrotechnik Tel. +49 355 69-4129 Universitaetsplatz

Hello (Timo), yes, I know... maybe it's a bit far fetched, but I do have a use case. Would it be (reasonably easily) possible to implement IMAP ACL support for "imported" imapc namespaces/mailboxes? There is no local directory (like with the virtual plugin) yet... Greetings, Lutz -- Lutz Pre?ler <Lutz.Pressler at SerNet.DE> http://www.SerNet.DE/ SerNet Service

Hi! I am distributing 2.5.1p1 for production use on my system by now and prepare switching to protocol 2 as default protocol. I just noted, that ssh-agent can be used for protocol 1 and 2, but the keys kept in ssh-agent are not compared against keys in .ssh. Example: I have a DSA key in id_dsa which I load into ssh-agent on login. When connecting to an account accepting the key everything is

----- Forwarded message from SerNet Support Lutz Pre?ler <support at SerNet.DE> ----- From: SerNet Support Lutz Pre?ler <support at SerNet.DE> Subject: 2.1: imapc LIST problem / usage question To: dovecot at dovecot.org Date: Fri, 9 Sep 2011 12:00:46 +0200 Organization: SerNet Service Network GmbH Hello, 2.1.alpha1 (latest Debian auto build on amd64) As a test, namespace for gmail