similar to: R "old versions" archive anywhere?

Hi again! Following up my previous posting below (to which no response as yet), I have located a report which situates this type of question in a longitudinal modelling context. http://www4.stat.ncsu.edu/~dzhang2/paper/glm.ps Generalized Linear Models with Longitudinal Covariates Daowen Zhang & Xihong Lin (This work seems to originally date from around 1999). They consider an outcome Y,

Hi Folks, Quick question. I have the following line in an R code file which runs fine on Linux: if(PNG) png(GraphName,width=12,height=15,units="cm",res=200) I learn that, when the same code was run on a Windows machine, there was the following error: Error in png(GraphName,width=12,height=15,units="cm",res=200): unused argument(s) (units = "cm") Sorry to

Hi Folks, I'm a bit puzzled by the following (example): N<-factor(sample(c(1,2,3),1000,replace=TRUE)) unique(N) # [1] 3 2 1 # Levels: 1 2 3 So far so good. Now: contrasts(N)<-contr.treatment(3, base=1, contrasts=FALSE) contrasts(N) # 1 2 # 1 1 0 # 2 0 1 # 3 0 0 whereas: contr.treatment(3, base=1, contrasts=FALSE) # 1 2 3 # 1 1 0 0 # 2 0 1 0 # 3 0 0 1 contr.treatment(3, base=1,

Hi Folks, Compare print(1234567890,digits=4) # [1] 1.235e+09 print(1234567890,digits=5) # [1] 1234567890 Granted that digits: a non-null value for 'digits' specifies the minimum number of significant digits to be printed in values. how does R decide to switch from the "1.235e+09" (rounded to 4 digits, i.e. the minumum, in "e" notation) to

Hi Folks, I'm wondering if there's a compact way to achieve the following. The "dream" is that, by analogy with rep(c(0,1),times=c(3,4)) # [1] 0 0 0 1 1 1 1 one could write rep(c(0,1),times=c(3,4,5,6)) which would produce # [1] 0 0 0 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 in effect "recycling" x through 'times'. The objective is to produce a vector of

Thanks much Ted. I actually had just tried what you suggest here before you posted, and resolved the problem. Thanks also for the other tips. I wrote x = as.vector(c(1:12)) because I thought that the mode of x might be the problem, the error message pointing to .Random.seed notwithstanding. On a related note, I did a brief test a couple weeks back where I ran a million random samples of 3 from

Greetings all! I have been searching the Site for a function, say "subsets", such that for instance subsets(10) would return a (say) matrix of indices to the 2^10 subsets of N items -- perhaps in the form of 2^10 rows each of which is 10 entries each either TRUE or FALSE. Or 1 or 0. Or ... I can of course write my own, using good old looping technology or similar, but it would be

Hi Folks, This must be easy but I've not managed to locate the solution! Basically: I'm using sink() to save successively obtained results, e.g. I construct a set of regression coefficients etc. with names rows and columne (I want to see the names in the output) as an object (say "Object"), and then I emit it with print(Object) So far so good. But I also want to print a

Hi Folks, I don't often grumble, but this time I've found myself inconvenienced by a posting stored on R-help archives: https://stat.ethz.ch/pipermail/r-help/2009-September/211095.html This was Karin Groothuis-Oudshoorn & Stef van Buuren's message on 10 September about the new version of MICE. This has been sent by software which inserted no line-breaks. As a result, each

Hi Folks, I have been using contour() to produce some contour plots (of a spatially-smooted density produced by kde2d()), with very satisfactory results. I now want access to the coordinates of the points on the contours, and it would seem that contour() does not return a value, so there is nothing from which these could be extracted. However, apparently contourLines() does, and it seems to be

Greetings Moderators! I moderated the message below just now (one of two identical test messages). I rejected it, with a covering note to the author (John Munroe <munddr at gmail.com>, who does not appear to be subscribed to the list) that I was doing so because approving it would serve no useful purpose, and pointing out that it had been held because the message headers matched a filter

Hi Folks, Has anyone got a compiled binary of the package 'plm' and/or of the package 'kinship' for Debian Etch? I'm asking because I'd like to install 'plm', but have failed because it depends on 'kinship', and 'kinship' refuses to compile on my Debian Etch. The underlying reason is the Debian refuses to allow me to install the development

Greetings, All! I want to draw contour lines in red, using contour(), but also have the contour labels (for the level-values) in black so that they will stand out against a coloured background already generated using filled.contour() (the background shades from green at low levels of "risk" to red at high levels). In any case, contour labels in red are already somewhat inconspicuous

Hi Folks, My Debian 4.0 Etch (installed Sep 2007, R version 2.7.1 (2008-06-23) installed at the same time from a Debian repo) has just flashed up that I can update as follows, using the update manager: r-base-core New version: 2.8.0-1~etch-cran.0 r-base-html New version: 2.8.0-1~etch-cran.0 r-base-latex New version: 2.8.0-1~etch-cran.0 r-cran-boot New version: 1.2.34-1~etch-cran.0

Hi Folks, Maybe I've missed it already being available somehow, but if the following isn't available I'd like to suggest it. If you're happy to let plot() choose its own limits, then of course plot(x,y) will do it. If you know what limits you want, then plot(x,y,xlim=c(x0,x1),ylim(y0,y1) will do it. But sometimes one would like to a) make sure that (e.g.) the y-axis has a

A follow-up to my previous query. Say one of the results returned by contourLines() is C.W <- contourLines(....) Then C.W is a list of (in this case 28) lists, each of which is a list with components $level (a single number), $x (the x coords of the points on a contour at that level) and $y (the y coords). I can of course get these individually with, for the 5th one for instance,

Hi Folks, This tip is probably lurking somewhere already, but I've just discovered it the hard way, so it is probably worth passing on for the benefit of those who might otherwise hack their way along the same path. Say (for example) you want to do a logistic regression of a binary response Y on variables X1, X2, X3, X4: GLM <- glm(Y ~ X1 + X2 + X3 + X4) Say there are 1000 cases in the

Hi Folks, I'm comparing some output from R with output from SPSS. The coefficients of the independent variables (which are all factors, each at 2 levels) are identical. However, R's Intercept (using default contr.treatment) differs from SPSS's 'constant'. It seems that the contrasts were set in SPSS using /CONTRAST (varname)=Simple(1) I can get R's Intercept to match

Greetings! I want to follow the evolution of a Nelder-Mead function minimisation (a function of 2 variables). Hence each simplex will have 3 vertices. Therefore I would like to have a function which can output the coordinates of the 3 vertices after each new simplex is generated. However, there seems to be no way (which I can detect) of extracting this information from optim() (the

Hi Folks. May I suggest a small addition to ?factor -- this is not explicit about how to force the ordering of levels in an ordered factor. One can guess it, but I think it would be a good idea to spell it out. For example, one can do: X <- c("Alf","Bert","Chris","Dave") F <-