similar to: R indentation

I noticed the same thing as I reported below, when I upgraded from R-2.0.0-alpha-20040918 to today's beta version: 'R CMD check' failed as before; but after cleaning 'src' from .o and .so files, it worked again. So the two-day-old '.so' file is obviously different from today's. Has compiler directives changed in the two last days? Or what is happening? G?ran On

I use the recommendations in "R coding standards", i.e., I put ;;; C (add-hook 'c-mode-hook (lambda () (c-set-style "bsd"))) ;;; ESS (add-hook 'ess-mode-hook (lambda () (ess-set-style 'C++ 'quiet) (add-hook 'local-write-file-hooks (lambda

I need to order a long vector of integers with rather few unique values. This is very slow: > x <- sample(rep(c(1:10), 50000)) > system.time(ord <- order(x)) [1] 189.18 0.09 190.48 0.00 0.00 But with no ties > y <- sample(500000) > system.time(ord1 <- order(y)) [1] 1.18 0.00 1.18 0.00 0.00 it is very fast! This gave me the following idea: Since I don't care

Is it possible to let a function be known under two names without having two identical copies of the function body? --- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics fax: +46 90 786 6614 Ume? University http://www.stat.umu.se/egna/gb/ SE-90187 Ume?, Sweden e-mail: gb at stat.umu.se

I wanted to add some text to a plot and got (R-1.6.0, Linux): > text(x = c(1, 4), y = 5, labels = x) Error in text.default(x = c(1766, 1895), y = 5, labels = x) : Object "x" not found With the default value of 'labels': > text(x = c(1, 2), y = 5, labels = seq(along = x)) Error in seq(along = x) : Object "x" not found A scoping bug? :) But >

When reading an spss file, I get the following message: Warning message: ../totmorH.sav: Unrecognized record type 7, subtype 13 encountered in system file. What does it mean? Do I have to worry? The result looks as expected. I do not have (easy) access to SPSS, so I cannot check exactly that I got what I should have. Thanks, G??ran And, I'm using R-2.0.1 on Debian testing (from source).

I will be giving a course in survival analysis using R (of course!) for people who know nothing about the subject (including R), but know basic statistics. I'm looking for a suitable course book. Therneau & Grambsch (2000) is an excellent book, but too much for this course. I need somthing more elementary. I have a vague memory saying that such books exist, but I cannot find any for the

I find 'How do I convert factors to numeric?' in the FAQ (7.12), but not the other way around. Trivial maybe, but > codes(factor(c(2, 10)) [1] 2 1 > codes(factor(c(2, 9))) [1] 1 2 How do I get the levels attached to the codes in numeric order? G?ran --- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics fax: +46 90 786 6614 Ume? University

I have a data frame containing children, with variables 'year' = birth year, and 'm.id' = mother's id number. Let's assume that all the births of each mother is represented in the data frame. Now I want to create a subset of this data frame containing all children, whose mother's first birth was in the year 1816 or later. This seems to work: mid <-

A few days ago I uploaded to CRAN a new package called 'eha', which stands for 'Event History Analysis'. Its main focus is on proportional hazards modeling in survival analysis, and in that respect eha can be regarded as a complement and an extension to the 'survival' package. In fact eha requires survival. Eha contains three functions for proportional hazards

A few days ago I uploaded to CRAN a new package called 'eha', which stands for 'Event History Analysis'. Its main focus is on proportional hazards modeling in survival analysis, and in that respect eha can be regarded as a complement and an extension to the 'survival' package. In fact eha requires survival. Eha contains three functions for proportional hazards

Dear R-helpers, In a data set I got from a medical doctor there are six treatment groups and (about) 5 bivariate responses in each group. Using 'manova', it is easy to see significant differences in treatment effects, but the doctor is more interested in the correlation between the two responses (within groups). I'm willing to assume a common value over groups, and one way of

Now I try to add some C and Fortan code to my package, so the NAMESPACE file is useDynLib(eha) importFrom(survival, Surv) export(mlreg.fit, risksets) but I get ..... * checking R files for library.dynam ... OK * checking S3 generic/method consistency ... WARNING Error in .try_quietly({ : Error in library(package, lib.loc = lib.loc, character.only = TRUE, verbose = FALSE) :

I have been bitten by what is clearly described on the help page for 'sample', namely sampling from a population of size one. I agree that it is convenient to have an exception if 'length(x) == 1', but my suggestion is to enforce the exception only if 'x' is numeric. In any case, if x is not numeric and of length 1, all you get is an error message. I would like

Reading about 'R_max_col' in "Writing R extensions", Version 2.1.0,(2005-01-03), I find: "Given the nr by nc matrix matrix in row ("FORTRAN") order, ..." Looks like a contradiction to me, since FORTRAN stores matrices columnwise. So is this a documentation bug? -- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics

This is a question that I should have sent to 'ess-help', but I take my chances... In 'NEWS' (1.9.0) I read the good news o Underscore '_' is now allowed in syntactically valid names, and make.names() no longer changes underscores. Very old code ... but when I try it in emacs (ESS 5.1.24), my '_' are changed to ' <- '! What can I do?

I have two versions of R on my Debian testing: R-1.9.0 (precompiled) and R-1.9.1-alpha of 2004-06-07 (built on my system). In the latter, I get > system("emacs sim.R &") > emacs: standard input is not a tty while with R-1.9.0 it works as expected. I have absolutely no idea about what's going on, so I could use some help. Anyone? Thanks, G?ran -- G?ran Brostr?m

I am calling 'vmmin' several times from a C function (which is called via .C). It works very well, except for memory consumption. The cause is that vmmin allocates memory via R_alloc, and this memory is not freed as vmmin exits. Instead all the allocated memory is freed on return of the .C call. In one application, I have 2000 functions of 500 variables each to minimize. In each call to

I'm learning to use Sweave (it's really great!) and everything works fine, except that occasionally I get, after having sweaved many times, ---------------------------------------------------------- > Sweave("report.Snw") Writing to file report.tex Processing code chunks ... 1 : term hide Error: chunk 1 Error in file(file, encoding = encoding) : All connections

Guido, > thanks for signaling this. The --mno-cygwin was a my typos error > but I have to correct the readme.packages (--no-cygwin is a survival > from the about two year ago version of cygwin). The flag is set correctly in MkRules, though. > Anyway, note that in readme.packages we suggest to edit 'MkRules' > not MakeDll and to set at the beginning of that file