search for: glmlist

For the vcdExtra package, I'm exploring methods to generate and fit collections of glm models, and handling lists of such model objects, of class "glmlist". The simplest example is fitting all k-way models, from k=0 (null model) to the model with the highest-order interaction. I'm having trouble writing a function, Kway (below) to do what is done in the example below > factors <- expand.grid(A=1:3, B=1:2, C=1:3) > Freq <-...

...s behaves inconsistently in R. In R anovalist.lm does a sequential ANOVA using pairwise F tests, ignoring all the other objects, so the larger of the two models provides the denominator. In S anova.lmlist uses the denominator from the largest model (smallest residual df) in the set, as does anova.glmlist in both. I suggest that R's anovalist.lm is wrong (that might be a matter of opinion and can certainly be argued), but anovalist.lm and anova.glmlist should agree. ?anova says in R When given a sequence of objects, `anova' tests the models against one another in the order speci...

In R, we currently have the functions anovalist.lm and anova.glmlist S / S-plus has anova.lmlist anova.glmlist On the other hand, the [n]lme package (library) of Doug Bates and Jose Pinheiro has an "lmList" class and an anova.lmList(.) method for that. We are considering to use anovalist.lm and anovalist.glm instead of the S/S-plus names me...

...# Throws an error. print(anova(fit3,fit2)) # ``Works'' but gives negative degrees of freedom and sum of squares. Is this evidence of a ``bug''? Or am I being terribly young and naive to expect anova() to work at all in such circumstances? The error thrown is: Error in anova.glmlist(list(object, ...), test = test) : (list) object cannot be coerced to type 'double' cheers, Rolf Turner ###################################################################### Attention:\ This e-mail message is privileged and confid...{{dropped:9}}

...egate.default aggregate.ts agrep AIC AIC.default AIC.logLik alias .Alias alias.formula alias.lm alist all all.equal all.equal.character all.equal.default all.equal.factor all.equal.formula all.equal.language all.equal.list all.equal.numeric all.equal.POSIXct all.names all.vars anova anova.glm anova.glmlist anova.glm.null anovalist.lm anova.lm anova.lmlist anova.lm.null anova.mlm any aov aperm append apply approx approxfun apropos Arg args array arrows as.array as.call as.character as.character.default as.character.factor as.character.octmode as.character.POSIXt ... yinch zapsmall zip.file.extract fu...

...nsions of association plots (assoc) and sieve diagrams (sieve) for these models. See demo("mental-glm"), demo("vision-quasi") and others. o Some initial steps in collecting, handling and summarizing sets of related glm() and loglm() models provided by the functions glmlist(), loglmlist() summarize() methods, Kway() for fitting all k-way glm() models, etc. -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Street http://www.math.yorku.ca/SCS/friendl...

Dear R team, I don't get what I think I should get when using anova.glm() with family="gaussian" -- please ignore this and forgive me if this turns out to be another example of a fundamental misunderstanding on my part (a highly likely event!) For example: S <- as.factor(rep(c(rep("m",2),rep("f",2)),2)) A <-

...lhost). * Socket connections are now robust against spurious readability and spurious availability of an incoming connection. * blocking = FALSE is now respected also on the server side of a socket connection, allowing non-blocking read operations. * anova.glm() and anova.glmlist() computed incorrect score (Rao) tests in no-intercept cases. (Andr'e Gillibert, PR#17734) * summaryRprof() now should work correctly for the Rprof(*, memory.profiling=TRUE) case with small chunk size (and "tseries" or similar) thanks to a patch proposal by Benj...

...lhost). * Socket connections are now robust against spurious readability and spurious availability of an incoming connection. * blocking = FALSE is now respected also on the server side of a socket connection, allowing non-blocking read operations. * anova.glm() and anova.glmlist() computed incorrect score (Rao) tests in no-intercept cases. (Andr'e Gillibert, PR#17734) * summaryRprof() now should work correctly for the Rprof(*, memory.profiling=TRUE) case with small chunk size (and "tseries" or similar) thanks to a patch proposal by Benj...

...lhost). * Socket connections are now robust against spurious readability and spurious availability of an incoming connection. * blocking = FALSE is now respected also on the server side of a socket connection, allowing non-blocking read operations. * anova.glm() and anova.glmlist() computed incorrect score (Rao) tests in no-intercept cases. (Andr'e Gillibert, PR#17734) * summaryRprof() now should work correctly for the Rprof(*, memory.profiling=TRUE) case with small chunk size (and "tseries" or similar) thanks to a patch proposal by Benj...

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML