similar to: optim/vmmin and R_alloc

Hello, I get a segfault when running glmmboot in my own package glmmML. Has happened many time before, but this time I get no hint of where in my C functions the error might be. I give the output below. Can this be an R bug? I suspect it has to do with repeated calls to 'vmmin' like this: for (...){ vmax = vmaxget(); vmmin(*p, b, &Fmin, bfun,

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

Dear listers, Am currently using MCMC approaches to estimate some parameters of my model. One parameter has to be updated using a tuned gamma distribution. So at each iteration I estimate the mean and variance of the density of the gamma approximation using "vmmin" (i also supply the gradient argument). For moderate replications the procedure works, but if I increase them R crashes.

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

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

> Date: Wed, 20 Mar 2002 14:31:03 +0100 (CET) > From: =?iso-8859-1?Q?G=F6ran_Brostr=F6m?= <gb at stat.umu.se> > Subject: [R] optim with gradient > > I want to maximise a function using 'optim' with a method that requires > the gradient, so I supply two functions, 'fun' for the function value > and 'd.fun' for its gradient. My question is: Since

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

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

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

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

I use Gnu emacs 21, Linux RH9 and R-1.7.0. I have succeeded to get the recommended indentation of 4 in C thru customization as described in 'R-exts', p.73, but I can't get it to work in R code. Can someone help, for instance by sending me the appropriate lines in '.emacs'? Thanks! --- G?ran Brostr?m tel: +46 90 786 5223 Department of Statistics

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

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'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

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

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