search for: quine

...wo or more objects are used in the call, a similar table is produced showing the effects of the pairwise differences between the models, considered sequentially from first to last. neither of which are clear enough to disambiguate this. Example: library(MASS) data(quine) quine.hi <- aov(log(Days + 2.5) ~ .^4, quine) quine.nxt <- update(quine.hi, . ~ . - Eth:Sex:Age:Lrn) quine.lo <- aov(log(Days+2.5) ~ 1, quine) anova(quine.hi, quine.nxt, quine.lo) quine.hi1 <- glm(log(Days + 2.5) ~ .^4, data=quine) quine.nxt1 <- update(quine.hi1, . ~ . - Eth:Sex:Ag...

...LSE) } res } #################### library(MASS) # 'glm' objects x1 <- runif(100, -4, 4) x2 <- runif(100, -4, 4) y <- 1 + 2 * x1 + rnorm(100, sd = 3) dat <- data.frame(y, x1, x2) glmFit <- glm(y ~ x1 + x2, data = dat) fn(glmFit, data = dat) # 'aov' objects quine.hi <- aov(log(Days + 2.5) ~ .^4, quine) fn(quine.hi, data = quine) However, for "negbin" objects returned by glm.nb() the following problem occurs: quine.nb <- glm.nb(Days ~ .^4, data = quine) fn(quine.nb, data = quine) Any hints to overcome this are greatly appreciated. Than...

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 NA NA NA NA NA NA NA 21 22 2...

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:

...then the linear model, and so on. I hope I've understood right! However, this procedure if based on the assumption that the same value of theta applies across the board, regardless of the values of any covariates or factors. The MASS book (Section 7.4) illustrates the use of glm.nb() on the Quine data (for numbers of days absence from school in a school year by Australian children). There are four factors in addition to the "outcome" Days: Eth: Ethnicity (Aboriginal "A"/Non-Aboriginal "N") Lrn: Learning Ability (Slow learner"SL", Average Learner &qu...

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))) > n <- nrow(quine) > quine2 <- with(quine, + data.frame(tmp, + Count=as.vector(tapply(rep(1,n), list(Eth, Sex, Lrn, Age), sum)))) First fit a s...

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