similar to: Simulation from a model fitted by survreg.

Hello all, I upgraded from R 1.8.1 to 1.9.0 (Windows XP), and spotted an odd thing. The last three letters in the Swedish alphabet are ??, ?? and ??. (In case they don't show correctly: they are a with a ring, a with two dots, and o with two dots (HTML: å ä ö). When I use these as variable names in a data.frame, odd things happen: In R 1.8.1, ?? (å)

Full_Name: James Michael Rath Version: all (I think) OS: doesn't matter Submission from: (NULL) (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

Summary: The locale setting in the operating system seems to be involved in what confused me a little bit. Thank you all for your help, especially 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

Dear list members, This is a question is about building an R package under windows and cygwin. Please bear with me. I have a package in R that compiles well on my stationary computer (WINDOWS NT, R 1.8.1 and prior), and the resulting (package_version).zip file works well there. It contains an R script file and a C file. After installing tools to build the package on my laptop, the package

In "An Introduction to R" (See R help menu), there is an example of a function 'open.account' that makes use of the lexical scope in R. I have a set of functions that can be used to output R tables and graphics into a single report document. (I am aware that several tools can do this already). While reorganizing my code, I realize that I can collect my functions in a list, in

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

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" "Var4" "Var5" "Var6" "Var7" "Var8" "Var9" [10] "Var10" "Var11" "Var12"

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]

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

Hello all, I want to extract the survival function estimate from a model fitted by survreg(). Using predict.survreg(..., type="quantile", p=seq(0,1,0.001)), gives the quantiles, which I managed to turn around into a survival function estimate (Prob{T > t} as function of t). Is there a more straightforward way of doing this? I have had difficulties using pweibull() with the

Hej Lennart, I would 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

Hi Friends I am new here and have a problem Year Market Winner BID 1 1990 ABC Apple 0.1260 2 1990 ABC Apple 0.1395 3 1990 EFG Pear 0.1350 4 1991 EFG Apple 0.1113 5 1991 EFG Orange 0.1094 For each year and separately for the two

Hi! When checking validity of a model for a large number of experimental data I thought it to be interesting to check the information provided by the summary method programmatically. Still I could not find out which method to use to get to those data. Example (not my real world data, but to show the point): [BEGIN] > DNase1 <- subset(DNase, Run == 1) > fm1DNase1 <- nls(density ~

Hello everybody! I have an ARIMA model for a time series. This model was obtained through an auto.arima function. The resulting model is a ARIMA(2,1,4)(2,0,1)[12] with drift (my time series has monthly data). Then I perform a 12-step ahead forecast to the cited model... so far so good... but when I look the plot of my forecast I see that the result is really far from the behavior of my time

Hi all, I try to represent a multiple curve graphic where the x-axis is the temperature and the different y-axes are the different X (X22,X43,X44...) some X corresponds to the same molecule (22 and 44 are for CO2 for instance) so I use the same colour for them. I wanna mix the linetype with the colour to be able to visually see the difference between X43 and X45 The best I have done up to now

Hello R users, I am new to R and am wondering if anyone can help me out with the following issue: I wrote a function to build ts models using different inputs, but when R displays the call for a model, I cannot tell which variables it is using because it shows the arguments instead of the real variables passed to the function. (e.g Call: lm(formula = dyn(dep ~ lag(dep, -1) + indep)) --->

I'm fitting a GLMM to some questionnaire data. The structure is J individuals, nested within I areas, all of whom answer the same K (ordinal) questions. The model I'm using is based on so-called continuation ratios, so that it can be fitted using the lme4 package. The lmer function fits the model just fine, but using mcmcsamp to judge the variability of the parameter estimates produces

Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one................. When plotting dendrograms created by hclust you can "identify" clusters by clicking on the graphics and returning a list of what is contained in each cluster. However I'd like to be able to "zoom in" on specific clusters and plot

Hi, I am having a hard time getting what the summary.lm-output for an ANOVA / ANCOVA means. Examples I find always seem to deal with simpler cases than what I meet in my data. My main problem is understanding the output when getting significant INTERACTION TERMS (what never occurs in examples :(). The following is the output after summary.lm(ancova) where "week" is continuous,

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................