search for: quined

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

predict.lm has been broken by recent changes to the patched branch. It fails for all singular fits. An example: library(MASS) data(quine) quine.hi <- aov(log(Days + 2.5) ~ .^4, quine) quine.nxt <- update(quine.hi, . ~ . - Eth:Sex:Age:Lrn) predict(quine.nxt) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 NA NA NA NA NA NA NA NA NA NA NA NA NA

samba@samba.org Giulio Orsero <giulioo@pobox.com> Mike Quin <m.r.quin@stir.ac.uk> Matthew Halliday <matthewh@fesa.co.uk> Carsten =3D?iso-8859-1?Q?Nordstr=3DF8m?=3D Jensen <north@get2net.= "Carlos Vinueza M." <carlos@ecuadorexplorer.com> Eric Dahnke <edahnke@istreetlabs.com> Giulio Orsero <giulioo@pobox.com> Giulio Orsero

...ilter.login_url CAS::Filter.validate_url CAS::Filter.server_name -- Ola Bini (http://ola-bini.blogspot.com) JvYAML, RbYAML, JRuby and Jatha contributor System Developer, Karolinska Institutet (http://www.ki.se) OLogix Consulting (http://www.ologix.com) "Yields falsehood when quined" yields falsehood when quined.

I get this error ... NoMethodError in StoreController#index undefined method `salable_items'' for Product:Class and here is my store_controller.rb ... class StoreController < ApplicationController def index @products = Product.salable_items end def self.salable_items find(:all, :conditions => "date_available <= now()", :order => "date_available

...is and make it available as a plugin, or something like that? Regards -- Ola Bini (http://ola-bini.blogspot.com) JvYAML, RbYAML, JRuby and Jatha contributor System Developer, Karolinska Institutet (http://www.ki.se) OLogix Consulting (http://www.ologix.com) "Yields falsehood when quined" yields falsehood when quined.

Hello all I've just spent a few joyless hours trying to get EBImage to install in R. I'm running Ubuntu (Hardy Heron), fully updated (including R and Imagemagick). EBImage installation seems to work, but when I 'library(EBImage)' I get the following: - - - - Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared library '/home/qilin/R/i486-pc-linux-

Your understanding isn't similar to mine. Mine says robust/resistant methods are for data with heavy tails, not heteroscedasticity. The common ways to approach heteroscedasticity are transformation and weighting. The first is easy and usually quite effective for dose-response data. The second is not much harder. Both can be done in R with nls(). Andy From: Quin Wills > > I am

i am trying to create a system so that different users have different priviliges on my application. i was going to modify my authorize function so that i could pass the required role to it and have it check the current user for that role. i am using before_filter and would like to have just one function instead of writing one for each role like so: authorize(role) instead of:

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

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

...ic configuration to a plugin? Would kinda come in handy, actually... -- Ola Bini (http://ola-bini.blogspot.com) JvYAML, RbYAML, JRuby and Jatha contributor System Developer, Karolinska Institutet (http://www.ki.se) OLogix Consulting (http://www.ologix.com) "Yields falsehood when quined" yields falsehood when quined.

Hi I'm running the latest R on a presumably up to date Linux server. 'Doing something silly I'm sure, but can't see why my saved PLMset objects come out all wrong. To use an example: Setting up an example PLMset (I have the same problem no matter what example I use) > library(affyPLM) > data(Dilution) # affybatch object > Dilution = updateObject(Dilution)

Hi All I have a requirement to consume a 3rd party web service from my Rails application. I am doing this in my action require ''soap/wsdlDriver'' factory = SOAP::WSDLDriverFactory.new(TRANSIDIOM_WSDL_URL) soap = factory.create_rpc_driver soap.wiredump_file_base="#{RAILS_ROOT}/log/transidiom.log" param = %(<Request

Hi I could only find some discussion on this wrt lattice graphics (which I'm not using). Apologies if I'm missing something obvious. I'd like to produce 3 rug plots under a kernel density plot for a population. The population is subdivided into 3 subpopulations, which I'd like the rug plots to highlight. Naturally, when I do 3 rug plots, they all plot over each other.

...placecompleter.html Usage is provided in the README and RDoc. Regards -- Ola Bini (http://ola-bini.blogspot.com) JvYAML, RbYAML, JRuby and Jatha contributor System Developer, Karolinska Institutet (http://www.ki.se) OLogix Consulting (http://www.ologix.com) "Yields falsehood when quined" yields falsehood when quined.

Hello all Simple question for the gurus. I'm trying to interactively label curves on a single plot. The labcurve() function from Hmisc seems like the way to do this (?). I just can't seem to get it to work. Toy example: > x <- 1:10 > y1 <- x^2 > y2 <- 2*x > plot(x,y1) > lines(x,y2) > labcurve(labels=c("curve1", "curve2"),

...a.bini at gmail.com Anda Abramovici - anda.abramovici at gmail.com -- Ola Bini (http://ola-bini.blogspot.com) JRuby Core Developer Developer, ThoughtWorks Studios (http://studios.thoughtworks.com) Practical JRuby on Rails (http://apress.com/book/view/9781590598818) "Yields falsehood when quined" yields falsehood when quined.