similar to: Help with sample function

When plotting using the cex.axis argument to boxplot(), the size of the plotting symbols beyond the whiskers of the boxplot are being changes. Example: par(mfrow=c(2,1)) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=2) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=1) This is from the following line in bxp()

Greetings, How to I read in SPSS .sav files into R. Thank you. David ===================================== David Kaplan, PhD Professor of Education School of Education University of Delaware Newark DE 19716 Voice: 302-831-8696 Fax: 302-831-4110 email: mailto:dkaplan at udel.edu Homepage:

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

Hi, when I generated a survfit() object, I can get number of patients at risk at various time points by using summary(): fit<-survfit(Surv(time,status)~class,data=mtdata) summary(fit) class=1 time n.risk n.event survival std.err lower 95% CI upper 95% CI 9.9 78 1 0.987 0.0127 0.963 1 41.5 77 1 0.974 0.0179 0.940 1 54.0 76

Hi, Is it possible to add labels to variables in R (so as to have a better description of what the variables represent)? Thanks, Neil

Hello there, I am doing analysis on survival data. How do I pick out a probability of survival at a chosen landmark time(for example, 3 years, 4 years) from the result of "survfit"? or any other functions? As I know, the result of 'survfit' only have the probabilities for all event or cencor time. Thank you very much Lisa Wang Princess Margaret Hospital Toronto, Ca tel:

There are two ways to view weights. One is to treat them as case weights, i.e., a weight of 3 means that there were actually three identical observations in the primary data, which were collapsed to a single observation in the data frame to save space. This is the assumption of survfit. (Most readers of this list will be too young to remember when computer memory was so small that we had to

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung,

R Folk: Please forgive what I'm sure is a fairly na?ve question; I hope it's clear. A colleague and I have been doing a really simple one-off survival analysis, but this is an area with which we are not very familiar, we just happen to have gathered some data that needs this type of analysis. We've done quite a bit of reading, but answers escape us, even though the question below

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, using the lung dataset. The first problem is the warning: To me, it seems that there are too many dummies

Dear all - I have followed the thread the reply to which was lead by Thomas Lumley about using pweibull to generate fitted survival curves for survreg models. http://tolstoy.newcastle.edu.au/R/help/04/11/7766.html Using the lung data set, data(lung) lung.wbs <- survreg( Surv(time, status)~ 1, data=lung, dist='weibull') curve(pweibull(x, scale=exp(coef(lung.wbs)),

Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax:

Hi, I'm sure this should be simple but I can't figure it out! I want to get the median survival calculated by the survfit function and use the value rather than just be able to print it. Something like this: library(survival) data(lung) lung.byPS = survfit(Surv (time, status) ~ ph.ecog, data=lung) # lung.byPS Call: survfit(formula = Surv(time, status) ~ ph.ecog, data = lung) 1

Hi, I'm trying to do a simple survival analysis on some data, and I'm having the following problem (here's my code and the error message): out <- Surv(fup,event=status) Error in Surv(fup, event = status) : argument "time2" is missing, with no default >From reading the documentation, it seems that I should be able to simply write: Surv(time1, event) if my data is

First, here is your message as it appears on R-help. On 10/14/2012 05:00 AM, r-help-request@r-project.org 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

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]

Hello How do I use function 'MERGE" to combine the FILE A and FILE B below to make FILE C? Thank you FILE A 140 151 167 30.1 11.4 40 FILE B 140 167 5.7 30.3 FILE C 140 151 167 30.1 11.4 40 5.7 NA 30.3

I was looking at how the model.frame method for lm works and comparing it to my own for coxph. The big difference is that I try to retain xlevels and predvars information for a new model frame, and lm does not. I use a call to model.frame in predict.coxph, which is why I went that route, but never noted the difference till now (preparing for my course in Nashville). Could someone shed light

Hello, is there a R package that provides a log rank trend test for survival data in >=3 treatment groups? Or are there any comparable trend tests for survival data in R? Thanks a lot Markus -- Dipl. Inf. Markus Kreuz Universitaet Leipzig Institut fuer medizinische Informatik, Statistik und Epidemiologie (IMISE) Haertelstr. 16-18 D-04107 Leipzig Tel. +49 341 97 16 276 Fax. +49 341 97 16