similar to: predict.lm is broken in 1.1.0-patched (2000-August-7) (PR#626)

anova() on three or more objects 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

Dear all, I would like to Boostrap the stepAIC() procedure from package MASS for variety of model objects, i.e., fn <- function(object, data, B = 2){ n <- nrow(data) res <- vector(mode = "list", length = B) index <- sample(n, n * B, replace = TRUE) dim(index) <- c(n, B) for (i in 1:B) { up.obj <- update(object, data = data[index[, i], ])

Dear all, I noticed the following while playing around with fitting log-linear models to contingency tables using R 1.8.1, but the problem also exists under R 1.9.0. A reproducible example uses the following contingency table: > library(MASS) > data(quine) > tmp <- with(quine, expand.grid(Eth=levels(Eth), Sex=levels(Sex), + Lrn=levels(Lrn), Age=levels(Age)))

Dear list, This is a question about overdispersion and the ML estimates of the parameters returned by the glm.poisson.disp (L. Scrucca) and glm.nb (Venables and Ripley) functions. Both appear to assume a negative binomial distribution for the response variable. Paul and Banerjee (1998) developed C(alpha) tests for "interaction and main effects, in an unbalanced two-way layout of counts

Hi Ming, On 11.11.2014 08:56, Ming Lei wrote: > On Tue, Nov 11, 2014 at 7:31 AM, Jens Axboe <axboe at kernel.dk> wrote: > > Known, I'm afraid, Ming is looking into it. Actually I had also tried to reproduce this bug, without success. But today I happened to know how to trigger the bug, by coincidence, during testing other things. Try to run xfstests/generic/034. You'll

Hi Ming, On 11.11.2014 08:56, Ming Lei wrote: > On Tue, Nov 11, 2014 at 7:31 AM, Jens Axboe <axboe at kernel.dk> wrote: > > Known, I'm afraid, Ming is looking into it. Actually I had also tried to reproduce this bug, without success. But today I happened to know how to trigger the bug, by coincidence, during testing other things. Try to run xfstests/generic/034. You'll

Hi, I have a dataset with the following properties: Y_i ~ N(mu_i, theta * (mu_i)^2) ln(mu_i) = B'Xi theta and beta's are the parameters here. I want to come up with a model to fit the data with the above property and test that model on the built in R dataset quine. Does nls() make sense in this case? Or is there any existing R package which can fit this model? -Shelly -- View

Does anybody know of aything that will help me do Quantitiative Comparative Analysis (QCA) and/or Fuzzy set analysis?? Or failing that Quine? ta rg Prof R Gott Durham Univesrity UK

Is there an R version of the very handy S-PLUS function spin? I've tried the help pages and apropos() without success. Malcolm Quine -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the

Dear R community, I would like to compare the degree of aggregation (or dispersion) of bacteria isolated from plant material. My data are discrete counts from leaf washes. While I do have xy coordinates for each plant, it is aggregation in the sense of the concentration of bacteria in high density patches that I am interested in. My attempt to analyze this was to fit negative binomial

Hi *! I am writing a C compiler to target the Lego Mindstorms NXT virtual machine. While I was implementing the optimizer and code generator I realized that it might be useful to use LLVM for this. The question I have is related to writing a LLVM target for this architecture. I have read most of the documentation available on the website and read through portions of the code (including

Hi all, I'm very new to R and I'm trying to construct a threaded binary tree using recursive functions. I'm very confused was wondering if anyone had any R sample code they would share. I've come across a lot of C++ code(nothing in R) and this is not helping. best, MK -- View this message in context:

Hi all, I have been previously been using scatterplot3d package to create some graphs but unfortunately it does not allow me to rotate the plot on all three axis. The cloud() function in the lattice package does allow me to do so. When I was using scatterplot3d I was using a script (Shown Below) to calculate the mean, quartiles and range limits for all three axis and I was representing that on the

Hi Folks, I'm playing with glm.nb() in MASS. Reference: the negative binomial distribution P(y) = (Gamma(theta+y)/(Gamma(theta)*y!))*(p^theta)*(1-p)^y y = 0,1,2,... in the notation of the MASS book (section 7.4), where p = theta/(mu + theta) so (1-p) = mu/(mu + theta) where mu is the expected value of Y. It seems from ?glm.nb that an initial value of theta is either supplied, or

during my seminar of advanced exploitation techniques (SEAT, [1]) i developed some methods to crack into system via DMA (e.g. via firewire). as part of this i developed a program that steals loaded ssh private keys from ssh-agents. i was astonished to find that the keys are not immediately removed from the agent when a timeout occurs, but only the next time the agent is queried via its socket. i

Hi I am trying to generate a design using Algdesign and came across terms of "exact design" and "approxiamte theory design", I did not find a reference to explain what they are, could some one shed some light about this on me? Another question is, I want to measure the main effects and at least two interactions, variables are factors, how do I ensure this in formula,

Dear list, I have a problem with a toy example: mtrx <- matrix(c(1,1,0,1,1,1,0,1,1,0,0,1), nrow=3) rownames(ma) <- letters[1:3] I would like to determine which is the minimum combination of rows that "covers" all columns with at least a 1. None of the rows covers all columns; all three rows clearly covers all columns, but there are simpler combinations (1st and the 3rd, or 2nd

Summary: loglin will fail if there is a global variable called c. I have fixed this for 0.64.2 and 0.65. On Sat, 22 May 1999, Cor en Aylin wrote: > Dear Brian, > > Using the mass library for R (the most recent port dated 16 may 1999) I Could you give version numbers in such reports please: at the moment there is an official version (5.3pl037) of MASS for R, plus a pre-release of MASS

On 7 Aug 2000, Peter Dalgaard BSA wrote: > ripley@stats.ox.ac.uk writes: > > > predict.lm has been broken by recent changes to the patched branch. > > > > It fails for all singular fits. An example: <snip> > Sometimes the fix for one bug uncovers another. We have > > > coef.aov > function (object, ...) > { > z <- object$coef >

On 7 Aug 2000, Peter Dalgaard BSA wrote: > thomas@biostat.washington.edu writes: > > > > > coef.aov > > > function (object, ...) > > > { > > > z <- object$coef > > > z[!is.na(z)] > > > } > > > > coef.lm > > > function (object, ...) > > > object$coefficients > > > > Ok, this