similar to: tiny typo in ?prop.test: "if" for "is"

Dear Zongshan Li, I'm forwarding this back to r-help ... it's generally best to keep the conversation going in public, so that (a) other people can chime in with ideas and suggestions, (b) your question doesn't get lost if the original respondent doesn't have time to deal with it, (c) the answers are archived for future reference. Your example (I've attached it, hope it

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

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

with the following simple data frame dd = structure(list(z = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L ), .Label = c("a", "b"), class = "factor"), x = c(0.3, 0.2, 0.1, 0, 0, 0, 0.2, 0.3)), .Names = c("z", "x"), row.names = c(NA, -8L), class = "data.frame") I would like know if it's possible to use model.matrix() to construct the

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

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

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

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

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.

Executive summary: getRemotePkgDepends seems to use the wrong default string to (try to) download the CRAN package information. ======================== Someone asked me about the behavior of tools::pkgDepends() , and in investigating I discovered the following: ## tools:::getRemotePkgDepends ## undebug(tools:::getRemotePkgDepends) tools::pkgDepends("GDD") ## Error in

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

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

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

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 think this is probably known by someone, but I wanted to ask/comment: The 'adapt' package has been removed from CRAN because of an 'unclear' license. That makes sense, but it actually took a bit of digging for me to discover that, and if I had been a student I might not have figured it out. The package is simply missing from the CRAN compiled packages page; I did find

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

Dear list members, I want to perform in R the analysis "simultaneous confidence interval for multiple proportions", as illustrated in the article of Agresti et al. (2008) "Simultaneous confidence intervals for comparing binomial parameter", Biometrics 64, 1270-1275. If I am not wrong the R function implementing the Agresti et al. method is prop.test(). I ask an help because I

In the optim() documentation, the control parameter "maxit" says "There is no other stopping criteria." That should be "are", or "criterion" ... While I'm at it -- poking at the code a little bit, it looks as if the initial simplex is set from the initial point by displacing each parameter value by max(0.1,0.1*pmax(fabs(Bvec))), where Bvec is the