Re: [R] Mixed-effects model for pre-post randomization design

Thank you very much, Ben. Your answer has been very useful. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] help in anova

I suspect that your "treatment" levels are perfectly nested inside the "soilgroup" levels. If this were the case you can either use one of the factor in your analysis or the other, depending on what you want to analyze. Il giorno Wed, 11 Feb 2015 14:47:18 + Mir Salam ha scritto: > Dear all,

[R] Mixed-effects model for pre-post randomization design

DeaR userRs, I recently read this Liang-Zeger article: http://sankhya.isical.ac.in/search/62b1/fpaper7.html in which (among other things) they adopt a random intercept model for pre-post designed trials, using a conditional likelihood approach (I didn't think it possible with only two measuremen

Re: [R] Median expected survival

Hi, Lars. I don't understand well your question. why don't you simply type pred_leuk Call: survfit(formula = leuk.cox, newdata = leuk_new) records n.max n.start events median 0.95LCL 0.95UCL 1 4242 42 30 1 1 23 2 4242 42 30 7

Re: [R] survplot invert number at risk labels

Hi all, and thank you for your interested answers. The problem could be already solved. It is present on my old R 2.15.1 (x86_64-pc-linux-gnu), but after reproducing the error I got the expected result on a R 3.0.2 through a web interface. However, since I spent a considerable amount of time to un

[R] survplot invert number at risk labels

deaR user, I found an unexpected behaviour of the rms::survplot.survfit function, that is giving me inverted labels for the patient-at-risk rows. The problem is that, for some reason, the survival::summary.survfit function changes the order of two of the suvfit object's strata when called (on my d

[R] [OT] Advice for medium size data management

me to bother at all. I'm not asking for prepackaged solutions, rather for help in documentation seeking and links to useful documentation or other threads (for example: is it worthwhile using parallel computing?) Thank you to anyone for reading this email. Marco Barbàra. P.S.: I work on a Debia

[R] Why do not prevent users from creating a local copy of c() ?

Hi all, sometime in the past I accidentaly created a copy of the c() primitive function inside my ``laboratory'' workspace, which I normally use to experiment and learn. As a consequence, c() stopped to work correctly and started to return lists of symbols instead of vectors, and today I had to sp

Re: [R] help with a survplot

Il giorno Mon, 2 May 2011 07:30:10 -0700 (PDT) Frank Harrell ha scritto: > Please elaborate. It is simply that, generally speaking, i don't like adding numbers to a plot. I eventually realized that riskset cardinality may be a useful indication. However, do not discriminate between events and ce

Re: [R] help with a survplot

Thank you very much. Despite prof. Harrell's support (for whom I feel great esteem) I still remain doubtful about this feature. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://w

[R] Adjustment for multiple-comparison for log-rank test

DeaR experts, I was asked for a log-rank pairwise survival comparison. I've a straightforward way to do this using the SAS system: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#/documentation/cdl/en/statug/63033/HTML/default/statug_lifetest_sect019.htm What

[R] Modeling time varying effects in with cph: how to ?

t this is generally not recommended, isn't it? Any sort of comment will be appreciated. Thank you very much. Marco Barbàra. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Time dependent Cox model

eed time interaction... > > The cox.zph function is primarily graphical; I would respond to your > question with "is it good to look at scatterplots before fitting a > linear model?" My answer to this is emphatically yes. > > Terry Therneau > > Thank y

Re: [R] [Off topic?] Time dependent Cox model fitting and validation

this is self-reply only to insert my real name __ 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,

Re: [R] How to plot the sample size on top of the bars of a barchart plot?

Il giorno Wed, 19 May 2010 08:46:06 -0700 (PDT) Javier ha scritto: > > Dear users, > > Is there anyone how knows how to plot the sample size on top of the > bars of a barchart plot? Do you mean something like this? mp <- barplot( VADeaths[1,], ylim = c( 0 , 18 ) ) text( mp , VADeaths[1,] + 0

Re: [R] about chisq.test()

> > Is this result reliable? No. chisq.test is a non parametric goodness-of-fit test for *counts*. In my opinion you need a one sample t.test. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] problem with an old .RData file

to look inside this file to search for differences with working ones, but from within R I cannot save an ascii version of this file due to this problem, and i do not know what kind of compression was applied (gunzip says «not in gzip format»). Is there a way of inspecting this file? thank you

[R] Observed responses in 'augPred' data frame - Wrong order ?

CTYPE=it_IT;LC_NUMERIC=C;LC_TIME=it_IT;LC_COLLATE=it_IT;LC_MONETARY=C;LC_MESSAGES=it_IT;LC_PAPER=it_IT;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=it_IT;LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [