similar to: R-2.0.2 soon?

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of > Ted.Harding at nessie.mcc.ac.uk > Sent: Sunday, December 05, 2004 7:14 PM > To: r-help at stat.math.ethz.ch > Subject: Re: [R] What is the most useful way to detect > nonlinearity in lo > > > On 05-Dec-04 Peter Dalgaard wrote:

Hi, I finally managed to make it work (I just needed to have a X11 window open). Thank you very much for your help. Agnes -----Original Message----- From: Ted Harding [mailto:Ted.Harding at nessie.mcc.ac.uk] Sent: Friday, July 21, 2006 3:34 PM To: Paquet, Agnes Cc: r-help at stat.math.ethz.ch Subject: RE: [R] connection to X11 problem On 21-Jul-06 Paquet, Agnes wrote: > Dear List, >

Hi, Splus allows pasting a graphics object into the plotting window, which makes it possible to do what you describe below. Now I use R which doesn't seem to allow pasting the picture into the graphics window, so I copy the graph onto a transparency sheet, and stick it onto my screen using tape. The coordinates need to be converted to make sense (as you describe). My methodology with the

To estimate the covariance matrix of e you could use the sample covariance matrix of the residuals. If desired, use its cholesky decomposition to transform to make the error approximately uncorrelated, then refit (and back-transform the coefficient matrix). Stacking the columns of Y and replicating X won't do what you write; it forces each univariate regression to have the same coefficients.

Hi Folks, I don't think I'm up to implementing a function for drawing French Curves, as requested by "dream". However, the discussion about it took me back so vividly to the old days that, paranormally, I felt once again the urge to ascertain magnitudes as it used to be, and should be, done. I have therefore implemented, and hereby donate to the R community, a new function

Hi Folks, I'm looking for a neat procedure for the following: Given t0 such that pt(t0,df,ncp=0) = alpha (given) find ncp0 such that for given beta pt(t0,df,ncp=ncp0) = (1 - beta) (In other words, what's the ncp such that you get power (1-beta) to detect it, using a 1-sided test with size alpha when ncp = 0?) In the past I've done the groping by hand, but this time it needs

Hi Folks, A recent exchange on the 'octave' list led to the following paper being cited, which I had not met before: What Every Computer Scientist Should Know About Floating-Point Arithmetic, by David Goldberg, originally published in the March, 1991 issue of Computing Surveys. PDF and HTML versions are widely available on the web (see Google), e.g. at

Hi Folks, There has been past correspondence regarding Catherine Loader's Bell Labs (oops, Lucent) paper "Fast and Accurate Computation of Binomial Probabilities" which gives the algorithm on which R's dbinom() is based. The original URL given in the R documentation "?dbinom" is: http://cm.bell-labs.com/cm/ms/departments/sia/catherine/dbinom but this link is

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

Hi Folks, Given a series of n independent Bernoulli trials with outcomes Yi (i=1...n) and Prob[Yi = 1] = Pi, I want P = Prob[sum(Yi) = r] (r = 0,1,...,n) I can certainly find a way to do it: Let p be the vector c(P1,P2,...,Pn). The cases r=0 and r=n are trivial (and also are exceptions for the following routine). For a given value of r in (1:(n-1)), library(combinat) Set <- (1:n)

Hi Folks, While I was browsing in the R-help archives yesterday, I got curious about the time series of the sizes of the monthly archives in MB. This turned out to have an unexpected feature or two, which I leave to readers to explore for themselves. I'm now wondering at what point in time we might expect to be receiving 1000MB/month (30+MB/day). It's not that far away, it seems, but

On 12-Feb-06 benphalan at gmail.com wrote: > Full_Name: Ben Phalan > Version: 2.2.1 > OS: Win XP > Submission from: (NULL) (131.111.111.231) > > > I have noticed that floor returns the wrong number when there are more > than 15 > decimal places: > >> floor(6.999999999999999) > [1] 6 >> floor(6.9999999999999999) > [1] 7 > > There is a

> From: Robin Hankin > > [stuff about the CLT deleted] > > > > > So you can use R usefully to eveluate general statisical > > issues of this kind! > > > > absolutely! R is excellent for this sort of thing. I use it for > teaching stats all the time. > I'd say that without a tool like R you cannot learn statistics. > > > Consider

The Iraq Statistics thread is a very interesting topic, and no doubt a lot of us would like to chip in with our views and comments -- indeed we are likely to bring a more discriminating view to such discussion than might be the case on many other lists. So if this took on a life of its own then it could become extensive. At which point the R-help admin would quite rightly call time on us! Yet I

Hi Folks, This query arises out of a call for help from a friend who has been running R version 1.9 unbder Windows XP, and wanted advice on using install.packages() in order to install the Windows binary of the package "Epi". Since this environment is outside my personal experience, I could only help in a general way. I did manage to track his problem down to the fact that package

Hi Folks, I'm trying to locate the source code for a (typical) R graphics device, in order to study how it's done. The underlying reason is that I'm thinking of trying to create a graphics device for 'pic' (the diagram drawing component of [g]troff). I thought the xfig device would be a good place to look, since the format of an xfig file is similar in nature (though very

How about: nis.na <- complete.cases --- From: <Ted.Harding at nessie.mcc.ac.uk> Hi Folks, My recent response to Laura Quinn's query about matrix subsetting reminded of a question. I wrote: iDir <- ((Winds[,20]<45)|(Winds[,20]>315))&(!is.na(Winds[,20])) Now, I find "!is.na" a bit awkward to type, so I might prefer to type it as "nis.na".

Hi Folks, I'm trying to find a neat solution to an apparently simple problem, but one which turns out to be a bit more intricate and tricky than one might expect. Suppose I have numbers given to a large number of digits. For example 1234567021 where (though I don't know this beforehand) only the last 3 digits will be varying (and all 3 will vary). What I want is, give a vector x of

Hi Folks, R has known speed issues for recursive definitions. There was a thread "Extremely slow recursion in R?" in August 2006 (24 Aug, from Jason Liao), with some interesting comparisons between programming languages. I'm re-opening the topic, with an ulterior motive (stated at the end). The starting point was the question "How many distinct increasing sequences of length

Hi Folks, Can someone satisfy my curiosity (well, at least about the following query!) Reading the draft proofs of an article I've been asked to look through, I find the typesetter has set what would normally be entered as source("xyz.R") or source('xyz.R') as source(`xyz.R') i.e. it has come out with an opening backquote, then xyz.R, then a closing forward