[R] bug submitting bug? [Fwd: Undelivered Mail Returned to Sender]

socket/lmtp] said: 452 4.2.2 Over quota (in reply to RCPT TO command) (B) Full_Name: Chris Andrews Version: 2.7.1 OS: Windows Submission from: (NULL) (72.89.97.235) # The following code shows that the legend is misplaced when a title is included AND xlim specifies a decreasing x axis. matpl

[R] Brief Mention of R in New York Times (July 29, 2009)

In an article about the IBM purchase of SPSS, R gets a brief mention as a "price pressure" on software suppliers. http://www.nytimes.com/2009/07/29/technology/companies/29ibm.html?_r=1&hpw "The broad consolidation wave in business intelligence software, analysts say, will bring increasing pric

Re: [R] Question on survival

lau pel wrote > Hi, > I'm going crazy trying to plot a quite simple graph. > i need to plot estimated hazard rate from a cox model. > supposing the model i like this: > coxPhMod=coxph(Surv(TIME, EV) ~ AGE+A+B+strata(C) data=data) > with 4 level for C. > how can i obtain a graph with 4 estimated (be

Re: [R] Problems with coxph and survfit in a stratified model with interactions

rm wrote > I’m trying to set up proportional hazard model that is stratified with > respect to covariate 1 and has an interaction between covariate 1 and > another variable, covariate 2. Both variables are categorical. In the > following, I try to illustrate the two problems that I’ve encountered,

Re: [R] Error in Surv(time, status) : Time variable is not numeric

R is finding the R function 'time' rather than your variable time in your dataset. Perhaps adding 'data=mydata' to your function call will resolve your issue. Chris Braem M wrote: > > Dear, > > I want to analyze two-level survival data using a shared frailty model, > for which I want to use

Re: [R] understanding the output from survival analysis

Ullrika Sahlin ekol.lu.se> writes: > > Why do I get different sign of the coefficients of covariates when I run the > semi-parametric proportional hazard model (coxph) compared to the parametric > proportional hazard model (survreg)? > > Anyone with experience in extracting information from su

Re: [R] plotting density for truncated distribution

Another option mydata <- rnorm(10) mydata <- mydata[mydata>0] plot(density(c(mydata, -mydata), from=0)) If you want the area under the curve to be one, you'll need to double the density estimate dx <- density(c(mydata, -mydata), from=0) dx$y <- dx$y * 2 plot(dx) Chris Jeroen Ooms wrote: