similar to: design matrix construction question

the version 2 parser thinks I have unnecessary braces, but I can't find any. False positive or am I missing something? If a false positive, is there any way to work around the warning? * checking Rd files against version 2 parser ... WARNING Warning: ./man/dbetabinom.Rd:32-34: Unnecessary braces at ?{p(x) = % (C(N,x)*Beta(N-x+theta*(1-p),x+theta*p))/% Beta(theta*(1-p),theta*p)}?

Is there a reason that Reduce() doesn't take a "..." argument that would allow arbitrary extra arguments to be passed through to the function? Here is a little example of how this would be convenient: z <- list( data.frame(state=c("California"), cases=0), data.frame(state=c("California","Massachusetts"),

From: Vitalie S. <vitosmail <at> rambler.ru> Subject: Bug in cut.POSIXt Newsgroups: gmane.comp.lang.r.general Date: 2009-10-15 15:47:48 GMT (1 hour and 29 minutes ago) Hello Everyone, Before reporting decided to post here first: tt <- structure(c(1254238817, 1254238859, 1254238969, 1254239080), class = c("POSIXt",

I am having difficulty with evaluation/environment construction for a formula to be evaluated by model.matrix(). Basically, I want to construct a model matrix that first looks in "newdata" for the values of the model parameters, then in "object at data". Here's what I've tried: 1. model.matrix(~f,data=c(newdata,object at data)) -- fails because something (terms()?)

Consider the following situation: we have quantified algal concentrations for a variety of species using many samples at each of three years. It seems to make sense to generate a line plot (matplot-like), with each species plotted as a separate line, with the points connected to emphasize the temporal pattern. The problem: lots of overlapping error bars. The question: from both a

I have found the "simulate" method (incorporated in some packages) very handy. As far as I can tell the only class for which simulate is actually implemented in base R is lm ... this is actually a little dangerous for a naive user who might be tempted to try simulate(X) where X is a glm fit instead, because it defaults to simulate.lm (since glm inherits from the lm class), and the

Sweave does something clever with warnings, which I have so far been unable to figure out. There are a couple of threads on the list about this, but the best in here is a hack to redirect all the output and stick it back in. http://www.nabble.com/-R--Sweave-and-warning-messages-td7759353.html#a7759353

All R experts, How do I fit a dynamic Random effects model with a binary dependent variable in R Thanks JCM [[alternative HTML version deleted]]

The CRAN host in Tampa, FL (cran.hostingzero.net) isn't responding, and hasn't responded in quite a while -- at least problems were reported more than a year ago (Oct 2007) http://finzi.psych.upenn.edu/R/Rhelp02a/archive/116706.html although someone apparently succeeded in July 2008 http://finzi.psych.upenn.edu/R/Rhelp02a/archive/137451.html Perhaps it's just flaky, and not

Dear Daniel (and Duncan): back in January you reported that you were having trouble saving an rgl movie on Vista. I don't know if you ever got that resolved, but a student of mine is now having a similar problem, not specific to dynamic graphics (which would have to be a problem outside of rgl anyway). He may report more on the symptoms, but essentially rgl.snapshot produces black or funky

Does anyone have an address for a maintainer, or know what's going on? cheers Ben Bolker ben at bolker-lap2:~$ ping wiki.r-project.org PING econum.umh.ac.be (193.190.194.5) 56(84) bytes of data. ^C --- econum.umh.ac.be ping statistics --- 8 packets transmitted, 0 received, 100% packet loss, time 6999ms ben at bolker-lap2:~$ traceroute wiki.r-project.org traceroute to

[cc'ing to R-devel] It's moderately obscure to me too, but ... There are two sets of bessel code in the package, one in C++ from Chris Bond and the other in FORTRAN from Netlib (I think). The FORTRAN code is the one that's giving trouble, you might try just removing the FORTRAN code from the src directory and trying again.

I have an optimization question that I was hoping to get some suggestions on how best to go about sovling it. I would think there is probably a package that addresses this problem. This is an ordering optimzation problem. Best to describe it with a simple example. Say I have 100 "bins" each with a ball in it numbered from 1 to 100. Each bin can only hold one ball. This optimization is

Dear I am trying to find a 90% credible interval. I am using the following code. fit<-bugs( model.file=BUGScode, data=data, inits = list(geninits1,geninits2), parameters.to.save=keepers, n.chains=nchains, n.iter=runs, n.burnin=burn, n.thin=nthin, DIC= TRUE, bugs.directory="C:/Program Files/WINBUGS.14", \ ) But this is only giving 95%

I am currently attempting to hack the recently featured profileModels package so that it can handle models generated by the lme4 (mixed models) package. I'm getting really confused by different behavior of summary() before and after loading the lme4 package, and inside and outside the profileMethod() function. The basic behavior is that with lme4 loaded, and "obj" a fitted object

In the wake of discussions about documentation of R's random number generation (and whether there should be links to .Random.seed in other places), here is a proposed patch to Uniform.Rd from 2.8.1. (Watch out for line breaks.) [R-core may choose not to accept this, but at least it's easy if they want it.] cheers Ben Bolker *** Uniform.Rd.orig 2009-01-04 14:10:47.000000000 -0500

Another dumb ggplot2 question: a facet plot with free x scale, lots of different ranges, want to specify that the x axis ticks be sparse -- only 2 or 3. I was about to say that I could figure out how to manually set the ticks for each facet if necessary, but now that I've looked around some more (at scale_continuous, facet_grid, facet_wrap, the Book ...) I think I can't even do that.

I'm very excited by the new mle package now incorporated in stats4. If possible, I'd like to help develop it. In the past I wrote a similar package (mleprof, available from http://www.zoo.ufl.edu/bolker/R/src), and would like to see if there's anything that my package does that I could contribute (in particular, I'd like to make sure that the code is as robust as possible in

I'm probably missing something very basic here, but: I've written some C code that I load into R dynamically. In the course of this C code, I generate some multinomial random deviates. I initially used the publically available "randlib" library, which also implements its own random number generator and binomial deviates (which are used to generate the multinomial deviates).

Hi, I havw some porblem in R and i need your help. I have this vector and I want to change it to date. for example G=[05 12 2008] or g=[2] f=[3] y=[2208]