Re: [R] making a plot
Thanks David and Adrés, it worked fine. Enrico. 2014-10-20 15:37 GMT-02:00 David Winsemius : > > On Oct 20, 2014, at 8:24 AM, Andrés Aragón wrote: > > > Enrico, > > > > This may help you: > > > > text(locator(1), "*", cex=1.5,adj=0.5 > > > > and > > > > text(locator(1), "º", cex=1.5,adj=0.5 > > Why not just use the values of x2 and y2 that were given to segments: > > > text( (ano+ranges)[1:3], 1:3, "*", cex=1.5,adj=0.5) > > text( (ano+ranges)[4:6], 4:6 , "º", cex=1.5,adj=0.5) > > > > > > > Draw your plot, then write the code, locate the cursor on your plot, put > > the symbols where you want itl and click. > > > > Regards, > > > > > > Andrés > > > > PS ?locator > > > > > > 2014-10-20 9:46 GMT-05:00 Enrico Colosimo : > > > >> Dear all, > >> > >> I am struggling to make a plot for my survival analysis > >> class. > >> > >> This is my script > >> > >> > >> labels<-c('1','2','3','4','5','6') > >> ano<-c(2001,2002,2003,2004,2006,2008) > >> ranges<-c(6,3,4,5,4,2) > >> dotchart(ano, labels=labels, xlab='ano', > >> ylab='Pacientes',pch=20,xlim=c(min(ano), max(ano+ranges))) > >> segments(ano,1:6,ano+ranges,1:6,pch=25,lty=1,lend=4) > >> > >> I need to put an asterix (failure) by the end of the three first lines > and > >> a small circle (censoring) > >> by the end of the last three. > >> > >> Someone can help me? > >> > >> Thanks, > >> Enrico. > >> > >>[[alternative HTML version deleted]] > > > -- > > David Winsemius > Alameda, CA, USA > > [[alternative HTML version deleted]] __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] making a plot
Dear all, I am struggling to make a plot for my survival analysis class. This is my script labels<-c('1','2','3','4','5','6') ano<-c(2001,2002,2003,2004,2006,2008) ranges<-c(6,3,4,5,4,2) dotchart(ano, labels=labels, xlab='ano', ylab='Pacientes',pch=20,xlim=c(min(ano), max(ano+ranges))) segments(ano,1:6,ano+ranges,1:6,pch=25,lty=1,lend=4) I need to put an asterix (failure) by the end of the three first lines and a small circle (censoring) by the end of the last three. Someone can help me? Thanks, Enrico. [[alternative HTML version deleted]] __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS
But it is very very unlikely to have a time with survival probability of 0. in a real data set. It would be necessary huge data set. A Monte Carlo simulation could put a little more light in this issue? 2012/4/16 Frank Harrell > Just generate some data where the estimated survival probability is 0. > at > a certain time. The log-log transformation blows up. > Frank > > Enrico Colosimo wrote > > > > What are the significant problems of the log-log transformations? > > Any papers published about it? > > Enrico. > > > > > > 2012/4/14 Frank Harrell <f.harrell@> > > > >> I used log-log in my book too until Terry Therneau alerted me to the > >> significant problems this creates. In the 2nd edition it will use log > >> S(t). > >> Frank > >> > >> Paul Miller wrote > >> > > >> > 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) > >> > > >> > Originally, it looked like the SAS results could be replicated using > >> > conf.type="plain". Applying this option to my actual data revealed > that > >> > this was not the case, however. > >> > > >> >>From Dr. Harrell, I learned that using conf.type="log-log" may not be > >> such > >> a good idea. Interestingly though, I've seen at least one instance where > >> experts in the R community use this option in their book. The book is > >> about > >> 10 years old. So maybe opinion about the use of this option has shifted > >> since then. > >> > > >> > Thanks, > >> > > >> > Paul > >> > > >> > [[alternative HTML version deleted]] > >> > > >> > > >> > __ > >> > R-help@ mailing list > >> > https://stat.ethz.ch/mailman/listinfo/r-help > >> > PLEASE do read the posting guide > >> > http://www.R-project.org/posting-guide.html > >> > and provide commented, minimal, self-contained, reproducible code. > >> > > >> > >> > >> - > >> Frank Harrell > >> Department of Biostatistics, Vanderbilt University > >> -- > >> View this message in context: > >> > http://r.789695.n4.nabble.com/Kaplan-Meier-analysis-95-CI-wider-in-R-than-in-SAS-tp4554559p4557695.html > >> Sent from the R help mailing list archive at Nabble.com. > >> > >> __ > >> R-help@ mailing list > >> https://stat.ethz.ch/mailman/listinfo/r-help > >> PLEASE do read the posting guide > >> http://www.R-project.org/posting-guide.html > >> and provide commented, minimal, self-contained, reproducible code. > >> > > > > [[alternative HTML version deleted]] > > > > __ > > R-help@ mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > > http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > > > - > Frank Harrell > Department of Biostatistics, Vanderbilt University > -- > View this message in context: > http://r.789695.n4.nabble.com/Kaplan-Meier-analysis-95-CI-wider-in-R-than-in-SAS-tp4554559p4561432.html > Sent from the R help mailing list archive at Nabble.com. > > __ > 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-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS
What are the significant problems of the log-log transformations? Any papers published about it? Enrico. 2012/4/14 Frank Harrell > I used log-log in my book too until Terry Therneau alerted me to the > significant problems this creates. In the 2nd edition it will use log > S(t). > Frank > > Paul Miller wrote > > > > 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) > > > > Originally, it looked like the SAS results could be replicated using > > conf.type="plain". Applying this option to my actual data revealed that > > this was not the case, however. > > > >>From Dr. Harrell, I learned that using conf.type="log-log" may not be > such > a good idea. Interestingly though, I've seen at least one instance where > experts in the R community use this option in their book. The book is about > 10 years old. So maybe opinion about the use of this option has shifted > since then. > > > > Thanks, > > > > Paul > > > > [[alternative HTML version deleted]] > > > > > > __ > > R-help@ mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > > http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > > > > > - > Frank Harrell > Department of Biostatistics, Vanderbilt University > -- > View this message in context: > http://r.789695.n4.nabble.com/Kaplan-Meier-analysis-95-CI-wider-in-R-than-in-SAS-tp4554559p4557695.html > Sent from the R help mailing list archive at Nabble.com. > > __ > 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-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] lme command
Hi, I am doing a longitudinal data set fit using lme. I used two forms of the lme command and I am getting two different outputs. FIRST out<-lme(Altura~Idade+Idade2+sexo+status+Idade:sexo+Idade:status+Idade2:sexo+Idade2:status, random=(list(ident=~Idade+Idade2))) SECOND out<-lme(Altura~Idade+Idade2+sexo+status+Idade:sexo+Idade:status+Idade2:sexo+Idade2:status, random= ~Idade+Idade2|ident,data=dados) I got weird results from the first one and could not understand the reason of it. All the results are exactly the same but the intercetp, and the two main terms sexo (gender) and status (treatment). That differences made a lot of difference in the final results. Anybody can tell me what is the differences between them? Thanks. Enrico. __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] glm output for binomial family
Hello, I am having some trouble running a very simple example. I am running a logistic regression entering the SAME data set in two different forms and getting different values for the deviance residual. Just look with this naive data set: # 1- Entering as a Bernoulli data set y<-c(1,0,1,1,0) x<-c(2,2,5,5,8) ajust1<-glm(y~x,family=binomial(link="logit")) ajust1 # Coefficients: (Intercept) x 1.3107 -0.2017 Degrees of Freedom: 4 Total (i.e. Null); 3 Residual Null Deviance: 6.73 Residual Deviance: 6.491 AIC: 10.49 # # 2- Entering as Binomial data set # ysim<-c(1,2,0) ynao<-c(1,0,1) x<-c(2,5,8) dados<-cbind(ysim,ynao,x) dados<-as.data.frame(dados) attach(dados) ajust2<-glm(as.matrix(dados[,c(1,2)])~x,family=binomial, data=dados) summary(ajust2) # Coefficients: (Intercept) x 1.3107 -0.2017 Degrees of Freedom: 2 Total (i.e. Null); 1 Residual Null Deviance: 3.958 Residual Deviance: 3.718 AIC: 9.104 = It seems that there is problem with the first fitting!!! Best Enrico Colosimo Dept Statistics, UFMG Brazil __ 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-guide.html and provide commented, minimal, self-contained, reproducible code.