Re: [R] Survival Analysis with an Historical Control

Hi Chris, Thanks for your reply. Certainly no need to apologize for a delayed response. Appreciate your taking the time to answer my question. My concern was about the value of "16". My understanding is it is based on one expert's opinion. I wondered if this is in keeping with whatever the sta

[R] Using discrete choice data to estimate marginal willingness to pay (MWTP)

Hello Everyone,   Can anyone show me how to use discrete choice data to estimate marginal willingness to pay? I've searched for simple information or code on how to do this but haven't found much.   The closest I've come is a paper titled "Introduction to Choice Experiments using R." The paper s

[R] How good is R at making publication quality tables?

Hello Everyone,   I have just started learning R and am in the process of figuring out what it can and can't do. I must say I am very impressed with R so far and am amazed that something this good can actually be free.   Recently, I finished reading R for SAS and SPSS Users and have begun readin

Re: [R] How good is R at making publication quality tables?

Hello Everyone,   At last count, I had 14 responses to my question. Thanks very much to everyone who replied, for helping me to understand that R is in fact capable of making publication quality tables and is very flexible in rendering these outputs. I'm still working my way through some of the

[R] Reading SQL Server Tables using RODBC

Hello Everyone,   I'm still quite new to R and am having trouble reading SQL Server Tables using RODBC. I've looked though the RODBC documentation as well as material I found online and in the book "Data Manipulation with R," but I just can't quite seem to get things to work.   Right now, I have

[R] Reading SQL Server Tables using RODBC (Got it!)

Thanks to those who responded to my post about reading SQL server tables using RODBC. Someone in our IT department sent me a link that contained the following information:   SQL ODBC connection strings Standard Security:< br> "Driver={SQLServer};Server=Your_Server_Name;Database= Your_Database_Na

[R] Reading dates in R using SQL and otherwise (and some interesting behavior by the data editor)

Hello Everyone,   I am a newbie with about a month's worth of experience using R. I've just spent a little time learning how to work with date values. Mostly, this has involved converting text values into dates.   So far, I've managed to figure out how to do this in R proper using as.Date. How

Re: [R] Reading dates in R using SQL and otherwise (and someinteresting behavior by the data editor)

Hi Paul, Jim, Dennis, and others,   I got RODBC and SQL Server to create a date variable in R (see below). I'm pretty sure the commands I was issuing before were working but that I was stripping off the date class every time I submitted the command "fix(mydata)." So I kept getting the mistaken i

[R] Does "sink" stand for anything?

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Pa

[R] R and S-Plus: Two programs separated by a common language?

Hello Everyone,   My company purchased S-Plus before it was bought out by Tibco. My understanding is that we own version 7.0 outright.   So far, I've been learning in R but thought I might also try working in S-Plus. My understanding is that S-Plus has some useful extra features. Another potent

[R] SPLUS Seqtrial vs. R Packages for sequential clinical trials designs

Hello Everyone,   I’m a SAS user who has recently become interested in sequential clinical trials designs. I’ve discovered that the SAS based approaches for these designs are either too costly or are “experimental.” So now I’m looking for alternative software. Two programs that seem

[R] odfWeave: Trouble Getting the Package to Work

Hello All, Have recently started learning Sweave and Knitr. Am now trying to learn odfWeave as well. Things went pretty smoothly with Sweave and Knitr but I'm having some trouble with odfWeave. My understanding was that odfWeave should work in pretty much the same way as Sweave. With odfWeave,

Re: [R] odfWeave: Trouble Getting the Package to Work

expected '&' 1: paste(levels(iris$Species), collapse = &       ^ > > --- On Mon, 2/18/13, Max Kuhn wrote: From: Max Kuhn Subject: Re: [R] odfWeave: Trouble Getting the Package to Work To: "Paul Miller" Cc: "r-help@r-p

Re: [R] odfWeave: Trouble Getting the Package to Work

Hi Max,   It works! Updating the odfWeave and xml software did the trick. Thanks very much for your help.   Paul [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

[R] Converting matrix to data frame without losing an assigned dimname

Hello All, Would like to be able to convert a matrix to a dataframe without losing an assigned dimname. Here is an example that should illustrate what I'm talking about. tableData <- state.x77[c(7, 38, 20, 46), c(7, 1, 8)] names(dimnames(tableData)) <- c("State", "") tableData State F

Re: [R] Converting matrix to data frame without losing an assigned dimname

Hi David, Gabor, Arun, and Rolf,   Thanks for your replies to my question. Very interesting and very helpful. I particularly liked the one-line solutions proposed by David and Gabor. I had a suspicion that this conversion could be done using only a minimal amount of code.   For what I'm doing, i

[R] Table with n(%) formatting: A better way to do this?

Hello All, Am learning to create tables with n(%) formatting using R. Below is a working example. I think this is not bad but wondered if there are better ways of doing it. Although it can be quite humbling, seeing good code helps me assess my progress as an R programmer. Ultimately want to h

Re: [R] Table with n(%) formatting: A better way to do this?

Hello Dr. Murdoch, Thanks for your reply. Your tables package is interesting. Spent some time experimenting with it yesterday and it looks like something I can use. The n(%) format seems to be the standard for what I'm doing. So I can't directly use the output your package produces. It does pro

[R] odfWeave: Some questions about potential formatting options

Hello All, Learning to use the odfWeave package. I really like the package. It has good documentation, makes some very nice looking tables, and seems to have lots of options for customizing output. There are a few things I'd like to do that don't seem to be covered in the documentation though.

[R] Convert results from print(survfit(formula, ...)) into a matrix or data frame

Hello All, Below is some sample survival analysis code. I'd like to able to get the results from print(gehan.surv) into a matrix or data frame, so I can manipulate them and then create a table using odfWeave. Trouble is, I'm not quite sure how make such a conversion using the results from a pri

Re: [R] Convert results from print(survfit(formula, ...)) into a matrix or data frame

Hi Marc, Oh, yes. That is quite simple. Silly me for not recognizing this. Thanks very much for your help. Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] odfWeave: Some questions about potential formatting options

Hi Milan and Max, Thanks to each of you for your reply to my post. Thus far, I've managed to find answers to some of the questions I asked initially. I am now able to control the justification of the leftmost column in my tables, as well as to add borders to the top and bottom. I also downloade

Re: [R] odfWeave: Some questions about potential formatting options

uchet-Valat Subject: Re: [R] odfWeave: Some questions about potential formatting options To: "Paul Miller" Cc: "Max Kuhn" , r-help@r-project.org Received: Wednesday, April 17, 2013, 11:17 AM Le mardi 16 avril 2013 à 10:15 -0700, Paul Miller a écrit : > Hi Milan and Max, &

[R] Trouble Computing Type III SS in a Cox Regression using drop1 and Anova

Hello All, Am having some trouble computing Type III SS in a Cox Regression using either drop1 or Anova from the car package. Am hoping that people will take a look to see if they can tell what's going on. Here is my R code: cox3grp <- subset(survData, Treatment %in% c("DC", "DA", "DO"), c("PT

Re: [R] Trouble Computing Type III SS in a Cox Regression

Therneau wrote: > From: Terry Therneau > Subject: Re: Trouble Computing Type III SS in a Cox Regression > To: r-help@r-project.org, "Paul Miller" > Received: Wednesday, April 24, 2013, 5:55 PM > I should hope that there is trouble, > since "type III" is an u

Re: [R] Trouble Computing Type III SS in a Cox Regression

ning this is a > re-education issue, and I can't much help with that. Fortune nomination! cheers, Rolf --- On Thu, 4/25/13, Terry Therneau wrote: > From: Terry Therneau > Subject: Re: Trouble Computing Type III SS in a Cox Regression > To: "Paul Miller&qu

[R] Text mining? Text manipulation? Both? Predicting KRAS test results in cancer patients

Happy Friday Everyone,   Hope Friday afternoon doesn't turn out to be a terrible time to post a question. I've been doing a little data mining of patient text medical records as of late. I started out trying to predict whether or not cancer patients had received KRAS mutation testing and did qui

[R] Analysis of censored cost data

Hello All, Am intested in using R to analyze censored cost data. Looks like I'd need to estimate some means and associated variances. I've got a few papers describing this. These talk about the use of re-weighted estimators (Lin), inverse probability weighted estimators (Bang-Tsiatis), phase-ba

[R] Survival Analysis with an Historical Control

Hello All, I'm trying to figure out how to perform a survival analysis with an historical control. I've spent some time looking online and in my boooks but haven't found much showing how to do this. Was wondering if there is a R package that can do it, or if there are resources somewhere that s

Re: [R] Survival Analysis with an Historical Control

Hi Dr. Therneau, Thanks for your response. This is very helpful. My historical control value is 16 weeks. I've been having some trouble though determining how this value was obtained. Are you able to indicate how people normally go about determining a value for the historical control? Or do you

Re: [R] Survival Analysis with an Historical Control

Hi Chris, Thanks for pointing out the use of "View page source". Very helpful to know. Do you happen to know anything about how to perform the analysis itself? I haven't been able to find anything confirming that the approach described in my original email (below) is correct. Thanks, Paul

[R] SEM validation: Cross-Validation vs. Bootstrapping

Hello All, Recently, I was asked to help out with an SEM cross-validation analysis. Initially, the project was based on "sample-splitting" where half of cases were randomly assigned to a training sample and half to a testing sample. Attempts to replicate a model developed in the training sample

Re: [R] SEM validation: Cross-Validation vs. Bootstrapping

Hi Joshua, Thanks for your very helpful reply. I took a look at your mediation tutorial, MplusAutomation, and semutils. I may not need any of this for the current project, but it was very interesting and I could imagine using it in the future. I also looked at the Hastie, Tibshirani, & Friedman

Re: [R] Retain last grouping after a strsplit()

Hi Steven, Not sure if you want to understand regular expressions in general or just the solution to your particular problem. If it's the former and you'd be willing to read a book on the subject, I'd recommend "Mastering Regular Expressions" by Jeffrey Friedl. I'm about halfway through now, an

[R] Mixed Conditional Logit with nested data

Hello Everyone,   This is my first attempt to do something in R. As a precursor to a Willingness to Pay analysis, I want to conduct a Mixed Conditional Logit analysis but am unsure how to proceed because of some nesting within my data.   Below is some data and code that illustrate what I’m t

[R] Using McNemar's test to detect shifts in response from pre- to post-treatment

Hello Everyone,   I've been asked to check if there is a significant difference in the following:     Pre    Post Group A 15/19 14/19 Group B 14/19 10/19   My sense is that I need to perform McNemar's test on these data because responses are correlated within pa

[R] Sample size estimation for non-inferiority log-rank and Wilcoxon rank-sum tests

Hello Everyone,   I'm trying to conduct a couple of power analyses and was hoping someone might be able to help. I want to estimate the sample size that would be necessary to adequately power a couple of non-inferiority tests. The first would be a log-rank test and the second would be a Wilcoxon

[R] How many R packages are not free?

Hello Everyone,   Just finished reading A Handbook of Statistical Analyses using R by Everitt and Hothorn. I'll begin by saying that I quite liked the book. It's both little and mighty in the sense that it's very compact but contains a tremendous amount of useful material.   The last chapter of

Re: [R] How many R packages are not free?

Thanks to everyone who replied to my post. One person indicated that this question treads on "dangerous terrain." I had that sense too but really wanted to know how many packages might not be "free" in the sense that they could be used at no cost by everyone. I was a little disappointed to learn

[R] read.spss warning message (Unrecognized record type 7, subtype 18 encountered in system file)

Hello Everyone,   Trying to help someone recover the contents of an SPSS.sav file using read.spss. This seemed to work well but produced a warning message. My code and the warning are displayed below.   Spent some time looking for previous questions about this warning. Found a lot of questions p

[R] What is the SAS equivalent of this R glm() code?

    Hello Everyone,   I'm trying to use SAS to replicate some results obtained in R. I was wondering if anyone call tell me the SAS equivalent of the code that appears below.   fm.glm.x <- glm(resp ~ . - 1, data = as.data.frame(mm.x),  na.action = na.exclude, family = binomial(link = "probit"))

Re: [R] What is the SAS equivalent of this R glm() code?

ivalent of this R glm() code? To: "Paul Miller" Cc: r-help@r-project.org Received: Friday, December 3, 2010, 3:55 PM ummm... This is an R list. Shouldn't you be posting this query on a SAS list? -- Bert (and why would you want to do this anyway?! ... but we won't go there.)

Re: [R] What is the SAS equivalent of this R glm() code?

peter dalgaard Subject: Re: [R] What is the SAS equivalent of this R glm() code? To: "Paul Miller" Cc: r-help@r-project.org Received: Friday, December 3, 2010, 4:15 PM On Dec 3, 2010, at 22:42 , Paul Miller wrote: >  >  > Hello Everyone, >  > I'm trying to use SAS

[R] LaTeX, MiKTeX, LyX: A Guide for the Perplexed

Hello Everyone,   Been learning R over the past several months. Read several books and have learned a great deal about data manipulation, statistical analysis, and graphics.   Now I want to learn how to make nice looking documents and about "literate programming." My understanding is that R use

Re: [R] LaTeX, MiKTeX, LyX: A Guide for the Perplexed

Hello Everyone,   Thanks to all who replied to my earlier message. I got quite a few responses and several suggestions about how to create nice looking documents that incorporate output from R. It took some time to investigate these, and in some ways I'm more perplexed than when I started.   Af

[R] Code for The R Book

Hello Everyone,   Does anyone have the R code for this book? I contacted the author about this and was told that it used to be on the book's website, but that the publisher had asked that it be taken down. I'm hoping that someone will have downloaded the code when it was available, and that they

Re: [R] Code for The R Book

ked to do this and he sent me no code. So I was just hoping that someone would have downloaded the code previously and would be willing to share it with me.   Thanks,   Paul   --- On Thu, 12/16/10, Carlos Ortega wrote: From: Carlos Ortega Subject: Re: [R] Code for The R Book To: "Paul M

[R] Seeking feedback on my first attempt at R programming

Hello Everyone,   Below is my first attempt at R programming. The code replicates example 5.1 from Common Statistical Methods for Clinical Research with SAS Examples. I was hoping that people more experienced than myself would be willing to take a look and let me know what I did well and what co

Re: [R] Seeking feedback on my first attempt at R programming

Hi Gabor,   Thanks very much for your reply. I was hoping to ask a follow-up question. Can you explain why it's better to use transform than to use within? The main difference between the two functions seems to be that within allows one to use variables created in a step (e.g., A and B) to creat

[R] R Equivalent of SAS Datastep Line-Hold (@@) Specifier?

Hello Everyone,   Below is some SAS code that uses a "line hold specifier" to read multiple observations from each of several input lines of data. There are 3 patients per line in the in-stream data.   Is there a way in R to read this kind of data? I've looked in my books and online but haven't

Re: [R] R Equivalent of SAS Datastep Line-Hold (@@) Specifier?

1     109   M  47    15 6   A      1     111   F  69    70 Rather than use a textConnection, you would probably just use a file. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Joshua Wiley Sent: Friday, 30 July 2010 4:59 PM

[R] Setting up counting process data for survival analysis

Hello Everyone, I'm learning to do survival analysis in R using a time-varying covariate. I've managed to set up my data correctly using SAS and then to get the correct result using R. What I'd like help with is getting the data set up correctly in R without having to resort to using SAS. Bel

[R] Pairwise T-Tests and Dunnett's Test (possibly using multcomp)

Hello Everyone,   I've been learning to use R in my spare time over the past several months. I've read about 7-8 books on the subject. Lately I've been testing what I've learned by trying to replicate the analyses from some of my SAS books. This helps me make sure I know how to use R properly an

[R] Pairwise T-Tests and Dunnett's Test (possibly using multcomp)

Hello Everyone,   Figured out one part of the code. Setting the reference level for a factor is accomplished using the relevel funtion (pg. 383 of MASS; pg. 70 of Data Manipulation with R):   gad$dosegrp <- relevel(gad$dosegrp,3)   This works very well. Much better than using a format in SAS proc

Re: [R] Pairwise T-Tests and Dunnett's Test (possibly using multcomp)

Hi Josh,   Thanks for your reply. You're right about letting R's method dispatch system choose the method for summary and inserting more spaces in the code. I was just messing around with the code in the former case and forgot to change it back. As far as not having enough whitespace goes, I've

[R] Cox Proportional Hazards model with a time-varying covariate

Hello Everyone,   I'm learning how to perform various statistical analyses in R. I'm checking my understanding by replicating examples from my SAS books. Below is an attempt to replicate a Cox Proportional Hazards model with a time-varying covariate. I think I'm doing this correctly but am not c

Re: [R] Cox Proportional Hazards model with a time-varying covariate

Hi Terry and Kevin,   Thanks for your replies to my post. And thanks particularly to Terry for pointing out my confusion about the "exact partial likelihood" vs. the "exact marginal likelihood." I'm learning and so am likely to be confused about such things from time to time.   I had thought

[R] Test email: Please disregard

Having trouble posting. Thought it might have something to do with the particular message I'm sending. Sending this message to test that possibility. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

[R] Test email #2: Please ignore

Having trouble posting. This is a second test email to help determine if the problem has been resolved. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-g

[R] Multiple endpoint (possibly group sequential) sample size calculation

Hello everyone, I need to do a sample size calculation. The study two arms and two endpoints. The two arms are two different cancer drugs and the two endpoints reflect efficacy (based on progression free survival) and toxicity. Until now, I have been trying to understand this in terms of a on

Re: [R] Multiple endpoint (possibly group sequential) sample size calculation

Thanks for your reply Marc. I'll look into the information you sent. It ocurred to me that I could have done a better job of describing what I'm trying to do. It can be difficult to do that sometimes when your unfamiliar with a topic and are still trying to figure it out yourself. In investigat

[R] What's your favourite R book that doesn't have R in the title?

Hello All,   Recently, I discovered that there are appear to be some good R related books that don't have R in the title. For example, I'm currently reading Gellman and Hill's book Data Analysis using Regression and Multilevel/Hierarchical Models. The book uses R (and BUGS) but doesn't appear to

Re: [R] What's your favourite R book that doesn't have R in the title?

I can add one more title. The foundations of statistics: A simulation-based approach by Vasishth and Broe. Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/

[R] How can I get this function to work?

Hello All, Can anyone tell help me understand why the function below doesn't work and how I can fix it? Below are some sample data, some code that works on individual rows of the data, and my attempt to translate that code into a function. My hope is to get the function working and then to appl

Re: [R] How can I get this function to work?

Hello Bert and Sarah, Thank you for your replies. Helped me understand how people might perceive my question and why they might not respond. Spent some time learning about R's debugging tools this morning. Began to realize why my function didn't work. My second argument was the name of a vari

Re: [R] How can I get this function to work?

Hi Sarah, That I was making things too complicated doesn't surprise me. A skilled programmer makes everything look easy I think. And someone who is still learning does just the opposite. Am going to spend some time now looking through your tweaks. Thank you very much for your help. Paul

[R] Power analysis for Cox regression with a time-varying covariate

Hello All, Does anyone know where I can find information about how to do a power analysis for Cox regression with a time-varying covariate using R or some other readily available software? I've done some searching online but haven't found anything. Thanks, Paul _

Re: [R] Power analysis for Cox regression with a time-varying covariate

he actual analysis itself, I'll start out using the steps you've listed and see where that takes me.   Paul   --- On Fri, 7/13/12, Greg Snow <538...@gmail.com> wrote: From: Greg Snow <538...@gmail.com> Subject: Re: [R] Power analysis for Cox regression with a time-varying covari

Re: [R] Power analysis for Cox regression with a time-varying covariate

-varying covariate To: "Marc Schwartz" , "Greg Snow" <538...@gmail.com>, r-help@r-project.org, "Paul Miller" Received: Wednesday, July 18, 2012, 8:24 AM Marc gave the referencer for Schoenfeld's article.  It's actually quite simple. Sample size for a C

[R] A LaTeX question -- Hope people won't mind

Hello All, Hope people won't mind my posting a LaTeX question here. I know a lot of people who use R are also using LaTeX. I'm in a bit of a rush to complete a document and am having trouble with one aspect of the formatting. I'm creating a list of tables using: \listoftables I also have som

Re: [R] A LaTeX question -- Hope people won't mind

rom: Marc Schwartz > Subject: Re: [R] A LaTeX question -- Hope people won't mind > To: "Paul Miller" > Cc: r-help@r-project.org > Received: Monday, August 20, 2012, 3:27 PM > On Aug 20, 2012, at 3:13 PM, Paul > Miller > wrote: > > > Hello All, > >

[R] Sample Size Calculation using a Drop the Loser/Keep the Winner Design

Hello Everyone,   I need to perform a sample size calculation using a drop the loser/keep the winner design. I've been searching for examples of how to do this using R but haven't found much.   The most promising possibility thus far is an R function called DrpLsrNRst. This appears in the book

[R] Problem working with dates

Hello All, I've been reading books about R for awhile now and am in the process of replicating the SAS analyses from an old report. I want to be sure that I can do all the things I need to in R before using it in my daily work. So far, I've managed to read in all my data and have done some data

Re: [R] Problem working with dates

Hi Michael and David, Thank you both for your reply to my question. Problem solved. I'm finding that my level of success with R is a little uneven thus far. I'm sometimes surprised by the things I can do, but then am even more surprised by the simple things I struggle with. Appreciate your hel

[R] Coalesce function in BBmisc, emoa, and microbenchmark packages

Hello All, Need to coalesce some columns using R. Looked online to see how this is done. One approach appears to be to use ifelse. Also uncovered a coalesce function in the BBmisc, emoa, and microbenchmark packages. Trouble is I can't seem to get it to work in any of these packages. Or perhaps

Re: [R] Coalesce function in BBmisc, emoa, and microbenchmark packages

Hi Michael, So the coalesce functions don't do what I thought they did. I'll just use ifelse() then unless anyone makes it clear I should do something different. Thanks, Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/

Re: [R] Coalesce function in BBmisc, emoa, and microbenchmark packages

Hi Brian, This works very well. Still trying to develop some skill with R. So can't say I understand your function completely as yet, but will work on it. I had thought that your function might only work for two columns (because of the "function(x,y)" part), but the example below suggests it wi

[R] Reshaping data from long to wide without a "timevar"

Hello All, I was wondering if it's possible to reshape data from long to wide in R without using a "timevar". I've pasted some sample data below along with some code. The data are sorted by Subject and Drug. I want to transpose the Drug variable into multiple columns in alphabetical order. My

Re: [R] Reshaping data from long to wide without a "timevar"

Erlotinib, Paclitaxel" instead of "Gemcitabine, Paclitaxel, Erlotinib". That's what I mean when I say I want the columns in alphabetical order. Thanks, Paul --- On Tue, 3/20/12, R. Michael Weylandt wrote: > From: R. Michael Weylandt > Subject: Re: [R] Reshaping data

Re: [R] Reshaping data from long to wide without a "timevar"

Hello All, Tried some more Internet searches and came to the conclusion that one probably does need to create a "timevar" before reshaping from long to wide. Below is some code that creates the "timevar" and transposes the data. connection <- textConnection(" 005 1 Gemcitabine 005 2 Erlotinib

[R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

Hello All, Want very much to learn how to plot patient drug timelines. Trouble is I need to figure out how to do this today. So not much time for me to struggle with it. Hoping someone can just help me out a bit. Below are some sample data and code that produces what I think is the beginning o

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

anges? Or can anyone else show me how to make them? Thanks, Paul --- On Thu, 3/22/12, Bert Gunter wrote: > From: Bert Gunter > Subject: Re: [R] Plotting patient drug timelines using ggplot2 (or some other > means) -- Help!!! > To: "Paul Miller" > Cc: r-help@r-proje

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

le_key = 2, line = 1)") + opts(axis.text.x = theme_blank() ) --- On Thu, 3/22/12, R. Michael Weylandt wrote: > From: R. Michael Weylandt > Subject: Re: [R] Plotting patient drug timelines using ggplot2 (or some other > means) -- Help!!! > To: "Pau

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

et up correctly. The graphs themselves might be called something like plot11, plot12, plot21, and plot22. Thanks, Paul --- On Thu, 3/22/12, R. Michael Weylandt wrote: > From: R. Michael Weylandt > Subject: Re: [R] Plotting patient drug timelines using ggplot2 (or some oth

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

Hi Michael, This may be an ugly way of doing things but the "title = unique(paste( ..." bit in the code below automates the code for pattern, patient, and line. So now all I'd need is the looping part. If good ways of automating the code for pattern, patient, and line or of doing the looping a

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

Hi Michael, Added a little more to my code (see below). It now automatically sets the name of the file. It also does a better job of spacing the text for pattern and patient x line at the top of the graph. I really like the way this looks now. I just need to figure out how to loop through the

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

Petr suggested. I didn't think it would help but wanted to try just in case. No dice -- ggplot just doesn't seem to like lists. Thanks, Paul --- On Fri, 3/23/12, R. Michael Weylandt wrote: > From: R. Michael Weylandt > Subject: Re: [R] Plotting patient drug timelines using

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

s over the last day or so. Paul --- On Fri, 3/23/12, R. Michael Weylandt wrote: > From: R. Michael Weylandt > Subject: Re: [R] Plotting patient drug timelines using ggplot2 (or some other > means) -- Help!!! > To: "Paul Miller" > Cc: "Petr PIKAL" , r-he

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

insemius > Subject: Re: [R] Plotting patient drug timelines using ggplot2 (or some other > means) -- Help!!! > To: "Paul Miller" > Cc: "R. Michael Weylandt" , "Petr PIKAL" > , "Bert Gunter" , > r-help@r-project.org > Received: Friday, March

Re: [R] Plotting patient drug timelines using ggplot2 (or some other means) -- Help!!!

Hello All, Figured out how to get the gaps in the bars when a drug stops and the starts again (see below). Made the graph overlaid/superimposed/stacked and got the desired result. Not sure this is how an expert would do it. But it's simple and it works. Paul setwd("N:/Regimen Coding/0906/Pl

[R] Adding text for written comments to bottom of graphs

Hello All, Recently developed the code below for graphing patterns of chemotherapy administration. As someone just starting to use R in their work, I managed to figure out some parts of the code but needed help with others. setwd("N:/Regimen Coding/0906/Plots Test") getwd() TestData <- struct

Re: [R] Adding text for written comments to bottom of graphs

Hello Baptiste, What you've done is very interesting. Went through and tried to understand all the steps. Reminded me of studying languages in years gone by. Always found it easier to read and understand a sentence than to construct a sentence of one's own. This is particularly true when you're

Re: [R] Adding text for written comments to bottom of graphs

Hi Baptiste, Thanks for your help with this. Sorry for being slow to express my appreciation. I had intended to put some more time into tweaking the graphs before responding. Recently have been reading Hadley Wickham's ggplot2 book and have also located some materials on the knitr package. Didn

[R] Kaplan Meier analysis: 95% CI wider in R than in SAS

Hello All,   Am replicating in R an analysis I did earlier using SAS. See this as a test of whether I'm ready to start using R in my day-to-day work.   Just finished replicating a Kaplan Meier analysis. Everything seems to work out fine except for one thing. The 95% CI around my estimate for the

Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS

Hi Enrico, Not sure how SAS builds the CI but I can look into it. The SAS documentation does have a section on computational formulas at: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_lifetest_a000259.htm Although I can't provide my dataset, I can p

Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS

Hello Drs. Colosimo and Harrell,   Thank you for your replies to my question. From Dr. Colosimo, I was able to determine that the SAS results can be replicated by adding the option conf.type="log-log" to my code as in :   survobj <- survfit(survfrm, conf.type="log-log", data=Survival)   Originall

[R] Checking dates for entry errors

Hello Everyone,   I have a question about how best to check dates for entry errors. I recently discovered that R will read the incorrectly entered date "11/23/21931" without producing a warning or an error message at least under some circumstances.   > as.Date("11/23/21931", format = "%m/%d/%Y")

[R] Checking for invalid dates: Code works but needs improvement

Hello Everyone, Still new to R. Wrote some code that finds and prints invalid dates (see below). This code works but I suspect it's not very good. If someone could show me a better way, I'd greatly appreciate it. Here is some information about what I'm trying to accomplish. My sense is that th

Re: [R] R-help Digest, Vol 107, Issue 25

TestDates[DateNames][[i]][TestDates$Invalid==1], "\n") } Rui Barradas Message: 53 Date: Tue, 24 Jan 2012 08:54:49 -0800 (PST) From: Paul Miller To: r-help@r-project.org Subject: [R] Checking for invalid dates: Code works but needs improvement Message-ID: <1327424089.1149.y

[R] Checking for invalid dates: Code works but needs improvement

rint' it allows for several arguments and (very) simple formating. { cat("Error: Invalid date values in", DateNames[[i]], "\n", TestDates[DateNames][[i]][TestDates$Invalid==1], "\n") } Rui Barradas Message: 53 Date: Tue, 24 Jan 2012 08:54:49 -

  1   2   >