search for: sixten

...(129.116.226.162) The code for pnorm in R was adapted from a Fortran library published in the ACM TOMS journal. The published version had a typographical error, though, which was pointed out in a second article published three years after the original. The error was that a macro/variable named SIXTEN had a value of 1.6e0 instead of the correct 1.6e1 (that is, 16 not 1.6). The author of the code for R (Ross Ithaka?) must have worked from an outdated version of SPECFUN and not the corrected one that now appears in CALGO, STATLIB, et cetera. There is a macro in pnorm.c (in the nmath directory) w...

...ames in a data.frame, odd things happen: In R 1.8.1, ?? (å) doesn't work while the others do. In R 1.9.0, ?? (ö) doesn't work while the others do. Please find examples below. It would be nice if all three could be used in variable names. At least in Sweden :-) Thanks... Sixten. # # R 1.9.0: ?? is renamed to X. # > data.frame(a=1, ??=2, ??=3, ??=4) a ?? ?? X. 1 1 2 3 4 > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 9...

...or use with rweibull() since my error terms are assumed to be Weibull, but it does not make any sense. The mean survival time is easy to predict, but I would like to simulate individual survival times. I am probably missing something completely obvious. Any hints or advice are appreciated. Thanks Sixten > summary(mod1) Call: survreg(formula = Surv(tid, study$first.event.death) ~ regim + age + stadium2, data = study, dist = "weibull") Value Std. Error z p (Intercept) 11.6005 0.7539 15.387 2.01e-53 regimposto -0.1350 0.1558 -0.867 3.86e-01 age...

...age with the C source file removed will successfully build, and successfully load. Then of course there will be an error saying "C entry point ... not in load table" Therefore the DLL seems to be causing this problem. Any hints or advice would be highly appreciated. Thanks in advance, Sixten > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 2.1 year 2005 month 12 day 20 svn rev 3...

Seems like get() is what you are looking for; e.g., try dim(get("a")), names(get("a")), etc. HTH, Andy > From: Sixten Borg > > Hello, > > I have trouble converting a character string to a R object. > Let me describe this by an example; > > > dim(a) > [1] 270 14 > > dim("a") > NULL > > > names(a) > [1] "Var1" "Var2" "Var3&qu...

...ally the suggested work-around data.frame(..., check.names=F) which works very well. A mystery still to be solved is why two versions of R, running on the same machine on the same time, behaves differently. Please do not respond to this on the list. I very much welcome you not to respond at all. Sixten. >>> Prof Brian Ripley <ripley at stats.ox.ac.uk> 2004-06-24 10:20:55 >>> > Can we stop blaming R for things which are not its fault, especially as > that has already been pointed out twice this morning?

...m writing a script that imports several datasets. The file names of the datasets contain the '_' character which forces me to construct a valid dataset name for each file. Although I can do this by hand, I would like to know if there is any solution to my first approach. Thanks in advance, Sixten Sixten Borg IHE Box 2127 S-220 02 Lund Sweden tel: +46 46 32 91 07 fax:+46 46 12 16 04 www.ihe.se ******************************************************************* Note: The information contained in this message and attachments may be privileged and confidential and protected from disclosure...

...equire documentation of each function. The style that I see in the R help is that the functions are not enclosed like this in a list. I like the idea of having the functions collected in a single list, but I think the documentation might be messy. Any ideas, opinions, anyone? Thanks in advance, Sixten. Example: myreport <- report(filename="report.rtf") my.report$add.table(my.data.frame, "Table of ...") plot(runif(10)) my.report$add.picture("Plot of ...") or... r <- report(filename="report.rtf") r <- add.table(r, my.data.frame, "Table of.....

...rted to make sense (as you describe). My methodology with the transparency has a "stoneage" appearance, but it works quite well, if I may say so myself. I have written down some notes that describes what I did, but it's in Swedish so I guess it's not of much help. Kind regards, Sixten. >>> Ted Harding <Ted.Harding at nessie.mcc.ac.uk> 2005-03-20 17:09 >>> Hi Folks, I'm contemplating using locator() to digitise external graphics. To set context, I would be using X11 display on Linux. To pre-empt the obvious comment: I've found on the R site the...

...uld like to add one thing: Often, there already exists an R function that solves the problem at hand. Instead of writing your own function, search the help files [apropos(), help.search()]. What I some times find difficult though, is guessing what key words will produce relevant search hits... Mvh Sixten >>> <Lennart.Borgman at astrazeneca.com> 2004-02-17 15:36:12 >>> I have been lurking in this list a while and searching in the archives to find out how one learns to write fast R code. One solution seems to be to write part of the code not in R but in C. However after findi...

...timate (Prob{T > t} as function of t). Is there a more straightforward way of doing this? I have had difficulties using pweibull() with the coefficients reported by summary(<survreg object>). I am enclosing an outline of my code for reference if anyone is interested. Thanks in advance, Sixten --------------------------------------------------------------------- sure <- survreg(formula = Surv(time, dead) ~ age + group + sex, data = modb) nd <- data.frame( age=50, group=factor("A", levels=c("A", "B", "C")), sex=factor("M", l...

Not sure where to even start on this.... I'm hoping there's some debugging I can do... I have a loop that cycles through several different data sets (same structure, different info), performing randomForest growth and predictions... saving out the predictions for later study... I get about 5 hours in (9%... of the planned iterations.. yikes!) and R just freezes. This happens in