Re: [R] the standard error of the quantile

. But I need to estimator the standard error of the quantile in survival analysis, because my data is censored. For example~ T-c(84,240,261,332,348,437,521,565) S-c(0,1,1,0,1,0,1,0) ##0 for censoed; 1 for event G-rep(1,8) ori_s.surv-survfit(Surv(T,S)~G) - Frank Harrell Department

Re: [R] Similar function for Redun() from Hmisc ?

= 0.8) It takes enormously high time for execution, is there anything wrong in the script? Suggest any other similar function available for dropping redundant variables. Thanks in advance! ~A - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message

Re: [R] Bezier curve in R?

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/Bezier-curve-in-R

Re: [R] Correlation matrix removing insignificant R values

is significant (for 10 observations it should be above 0.64) Can someone help with this? - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Correlation-matrix-removing-insignificant-R-values-tp4099412p4099719

Re: [R] Thank you

://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

Re: [R] The contrast and Design libraries

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/The-contrast-and-Design-libraries

Re: [R] SPSS - R

, minimal, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/SPSS-R-tp4110995p4111006.html Sent from the R help mailing list archive at Nabble.com

Re: [R] fitting a distribution to a ecdf plot

__ 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

Re: [R] Nomogram with stratified cph in Design package-- failure probability

' ), lp=F, fun.at=c(at.surv, at.surv),label.every=1, force.label=FALSE, cex.axis=0.8, verbose=TRUE, cex.var=0.8) Thanks Min - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Nomogram-with-stratified-cph

Re: [R] Nomogram with stratified cph in rms package, how to get failure probability

any assistance in this matter. Thank you Very much. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Nomogram-with-stratified-cph-in-rms-package-how-to-get-failure-probability-tp4123249p4129173.html Sent from

[R] Regression Modeling Strategies 3-Day Short Course March 2012

*RMS Short Course 2012* Frank E. Harrell, Jr., Ph.D., Professor and Chair Department of Biostatistics, Vanderbilt University School of Medicine *March 7, 8 9, 2012* 8:00am - 4:30pm Student Life Center Board of Trust Room Vanderbilt University http://biostat.mc.vanderbilt.edu/RMSShortCourse

Re: [R] Stepwise in lme

... does anyone have some advice on how I might get such a crude method working or must I turn back and fight with stepAIC? Thank you! - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Stepwise-in-lme

[R] Call for Abstracts for useR! 2012

The international R User Conference useR! 2012 will be held June 12-15 at Vanderbilt University in Nashville, Tennessee, USA, with a pre-conference full-day course offered on June 11. Participants are encouraged to submit a one-page abstract for oral or poster presentation at the

Re: [R] R 2.14.0 Design library

-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

Re: [R] Joint confidence interval for fractional polynomial terms

://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/Joint-confidence-interval-for-fractional-polynomial

Re: [R] Joint confidence interval for fractional polynomial terms

, reproducible code. __ 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

[R] Call for Tutorial Proposals and Abstracts for useR! 2012

for contributed presentations and posters. Thanks! Frank Harrell The international R User Conference useR! 2012 will be held June 12-15 at Vanderbilt University in Nashville, Tennessee, USA, with a pre-conference full-day

Re: [R] Propensity score matching in R using Classification tree method

=rr1,nboots=1000,data=Passport) Please suggest me right approach for Propensity score matching using classification trees(if possible code me) . ps0-propensity score values - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r

Re: [R] kruskal wallis post hoc?

-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/kruskal

Re: [R] general question on Spotfire

__ 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

Re: [R] kruskal wallis post hoc?

-Whitney U test. Frank Peter Dalgaard-2 wrote On Jan 12, 2012, at 14:11 , Frank Harrell wrote: The Kruskal-Wallis test is a special case of the proportional odds ordinal logistic model. Eh? Can you elaborate on that? I would expect that at best it is equivalent to some _test_ in a polr

[R] Regression Modeling Strategies 3-Day Short Course March 2012

This course covers a variety of regression modeling and model validation methods as well as the R rms package. Please email interest to eve.a.ander...@vanderbilt.edu - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com

Re: [R] general question on Spotfire

functions in S-Plus ever happened. Frank Uwe Ligges-3 wrote On 12.01.2012 17:38, Frank Harrell wrote: As a slight aside, Tibco/Spotfire originally planned to provide a capability to load R packages into S-Plus. This always seemed to me to be a hard thing to do, and if my understanding

Re: [R] read in Stata and SPSS with value labels/formats

@ 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

Re: [R] How do I compare 47 GLM models with 1 to 5 interactions and unique combinations?

(Shells, TotalEggs-Shells) ~ HTL:Veg:TotalEggs:Aeventexhumed, data.to.analyze, family=binomial) Please advise, thanks! J - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/How-do-I-compare-47-GLM-models

Re: [R] How do I use the cut function to assign specific cut points?

figuring the code for 20 and =30? Please help. Thank you. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/How-do-I-use-the-cut-function-to-assign-specific-cut-points-tp4329788p4330380.html Sent from the R help

Re: [R] How do I compare 47 GLM models with 1 to 5 interactions and unique combinations?

not be pretty but works for me. An AIC difference equal or less than 2 is a tie, anything higher is evidence for the model with the less AIC (Sakamoto et al., Akaike Information Criterion Statistics, KTK Scientific Publishers, Tokyo). I wouldn't go quite as far as Frank Harrell did

Re: [R] pls. help -Logisitic Regression-Urgent

have 124 variables(mix of categorical and continuous) how would i select(R code) attribute - ie. independent variable from those 124 variable(Variable selection) so that i can build a model using glm pls. give me codes in R- on how to do that - Frank Harrell Department of Biostatistics

Re: [R] How do I compare 47 GLM models with 1 to 5 interactions and unique combinations?

. Model selection without shrinkage (penalization) seems to offer benefits but this is largely a mirage. Frank Rubén Roa wrote -Mensaje original- De: Bert Gunter [mailto:gunter.berton@] Enviado el: jueves, 26 de enero de 2012 21:20 Para: Rubén Roa CC: Ben Bolker; Frank Harrell

Re: [R] How do I compare 47 GLM models with 1 to 5 interactions and unique combinations?

Ruben you are mistaken on every single point. But I see it's not worth continuing this discussion. Frank Rubén Roa wrote -Mensaje original- De: r-help-bounces@ [mailto:r-help-bounces@] En nombre de Frank Harrell Enviado el: viernes, 27 de enero de 2012 14:28 Para: r-help@ Asunto

[R] useR! 2012: Earlybird Registration for International R Users Conference, Nashville TN 12-15 2012

the conference that you may also want to take advantage of. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/useR-2012-Earlybird-Registration-for-International-R-Users-Conference-Nashville-TN-12-15-2012

Re: [R] Post hoc test for lm() or glm() ?

-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/Post-hoc-test-for-lm-or-glm-tp4352761p4352799.html Sent from the R

Re: [R] glht (multicomparisons) with a binomial response variable

? Any advice, please? - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/glht-multicomparisons-with-a-binomial-response-variable-tp4360898p4361171.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Outlier removal techniques

, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Outlier-removal-techniques-tp4372652p4373592.html Sent from the R help mailing list archive at Nabble.com

Re: [R] question about val.surv in R

/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/question-about-val-surv-in-R-tp350p3507015.html Sent from the R help

Re: [R] Hosmer-Lemeshow 'goodness of fit'

SD 118.5308396 118.3771115 0.1435944 Z P 1.0705717 0.2843620 - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r

Re: [R] recommendation on B for validate.lrm () ?

. Is this how the procedure works in R? 2. Is the reason you recommend repeating k-fold 100 times because the partitioning is random, ie not 1st 10th, 2nd 10, et cetera so you might obtain slightly different results? - Frank Harrell Department of Biostatistics, Vanderbilt University

Re: [R] Hosmer-Lemeshow 'goodness of fit'

- Department of Internal Medicine University of Virginia, Charlottesville, Virginia - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Hosmer-Lemeshow-goodness-of-fit-tp3508127p3509102.html Sent from the R help

Re: [R] Draw a nomogram after glm

/nomogram_2.jpeg Please can I do to have linear predictor between 0 to 1. Thanks very much for your Help Komine - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Draw-a-nomogram-after-glm-tp3498144p3511758.html

Re: [R] Total effect of X on Y under presence of interaction effects

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/Total-effect-of-X-on-Y-under-presence

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

Using anything other than deviance (or likelihood) as the objective function will result in a suboptimal model. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Can-ROC-be-used-as-a-metric-for-optimal

Re: [R] Can ROC be used as a metric for optimal model selection for randomForest?

, 2011 at 8:11 AM, Frank Harrell lt;f.harr...@vanderbilt.edugt; wrote: Using anything other than deviance (or likelihood) as the objective function will result in a suboptimal model. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message

Re: [R] Question on approximations of full logistic regression model

__ 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. - Frank Harrell Department of Biostatistics, Vanderbilt

Re: [R] Question on approximations of full logistic regression model

structure. I am wondering which CIs are better. I guess bootstrapping one, but is it right? I would appreciate your help in advance. -- KH (11/05/16 12:25), Frank Harrell wrote: I think you are doing this correctly except for one thing. The validation and other inferential calculations

Re: [R] changes in coxph in survival from older version?

://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

Re: [R] changes in coxph in survival from older version?

narrow down the variables based on univarate analysis and build a multivariate model from that. Many thanks in advance. ...Tao - Original Message From: Frank Harrell lt;f.harr...@vanderbilt.edugt; To: r-help@r-project.org Sent: Mon, May 16, 2011 11:25:20 AM Subject: Re

Re: [R] changes in coxph in survival from older version?

-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/changes-in-coxph-in-survival-from-older-version

Re: [R] can not use plot.Predict {rms} reproduce figure 7.8 from Regression Modeling Strategies (http://biostat.mc.vanderbilt.edu/wiki/pub/Main/RmS/course2.pdf)

://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

Re: [R] logistic regression lrm() output

://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

Re: [R] Covariable Logistic Regression In R

/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

Re: [R] Covariable Logistic Regression In R

it. Thank you, Ana On Wed, May 18, 2011 at 3:14 PM, Frank Harrell lt;f.harr...@vanderbilt.edugt;wrote: First of all, use the correct terminology for the statistical method. It is not productive to identify 'significant' genes unless you use quite complex methods. This kind of endeavor would

Re: [R] changes in coxph in survival from older version?

to validate an existing model? Anything different? Many thanks! ...Tao - Original Message From: Frank Harrell lt;f.harr...@vanderbilt.edugt; To: r-help@r-project.org Sent: Tue, May 17, 2011 10:51:02 AM Subject: Re: [R] changes in coxph in survival from older version

Re: [R] How to calculate confidence interval of C statistic by rcorr.cens

://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

Re: [R] How to calculate confidence interval of C statistic by rcorr.cens

whether this is enough improvement though. -- Kohkichi (11/05/22 23:27), Frank Harrell wrote: S.D. is the standard deviation (standard error) of Dxy. It already includes the effective sample size in its computation so the sqrt(n) terms is not needed. The help file for rcorr.cens has

Re: [R] Multinomial Logistical Model

+ mldata$dr0 + mldata$dr1 + mldata$yr_tch_exp + mldata$yr_tch_exp_s + mldata$ismgdr + mldata$ismgyr_t_A + mldata$pair) Error: inherits(object, formula) is not TRUE - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4

Re: [R] Several Regression by combinations variables

, minimal, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Several-Regression-by-combinations-variables-tp3560123p3560662.html Sent from the R help mailing list archive

[R] Syntax for lattice scales argument with |a*b

I have not been able to figure out the format for specifying limits, at, labels in scales - list(relation='free', limits=..., at=..., labels=...) when there is more than one paneling variable, e.g., xyplot(y ~ x | a*b). Thanks for any guidance. Frank - Frank Harrell Department

Re: [R] Syntax for lattice scales argument with |a*b

Thanks very much David and Bert. That did it. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Syntax-for-lattice-scales-argument-with-a-b-tp3560973p3561131.html Sent from the R help mailing list archive

[R] Problem with % in an example when running R CMD check

of an .Rd file is: f - ols(y ~ pol(x1,2) + pol(x2,2) + pol(x1,2) %ia% pol(x2,2), subset=group==g) Putting \ before the % gets around the problem but I didn't have to do that before. Any ideas? Thanks Frank - Frank Harrell Department of Biostatistics, Vanderbilt University

Re: [R] Removal of elements from nomograms

this output in the first place. As ever, thanks in advance are likely due to Frank Harrell, without whom many things would be much more difficult. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] New version of rms package on CRAN

. More details are at http://biostat.mc.vanderbilt.edu/Rrms - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/New-version-of-rms-package-on-CRAN-tp3570655p3570655.html Sent from the R help mailing list archive

Re: [R] Interpreting Quantile Regression

This is not really an R question, and it indicates that you have a good deal of studying to do about quantile regression before you rely on it. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Interpreting

Re: [R] Logistic Regression

? - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Logistic-Regression-tp3578962p3579392.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] predict with model (rms package)

, minimal, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/predict-with-model-rms-package-tp3581229p3582758.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Logistic Regression

/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

Re: [R] 'PROC CONTENTS' = ?

See the contents function in the Hmisc package. karena wrote: Is there any function in R which does the same thing as what 'PROC CONTENTS' does in SAS? thanks, karena - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r

Re: [R] rcspline.plot query

__ 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. - Frank Harrell Department of Biostatistics

Re: [R] rcspline.plot query

than 6 non missing observations. George - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/rcspline-plot-query-tp3590233p3591951.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Somers Dyx

]] __ 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. - Frank Harrell Department of Biostatistics, Vanderbilt

[R] lines(..., lwd=3) inaccuracy

, xlab='',ylab='') text(rep(.5,p), 1:p, v, adj=1) maxabsr - max(abs(r[row(r) != col(r)])) for(i in 1:p) { for(j in 1:p) { if(i = j) next lines(c(i,i),c(j,j+r[i,j]/maxabsr/2), lwd=3) lines(c(i-.2,i+.2),c(j,j), col=gray(.7)) } text(i, i, v[i], srt=-45, adj=0) } - Frank Harrell

Re: [R] When models and anova(model) disagree...

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/When-models-and-anova-model-disagree

[R] Accessor functions in lattice graphics

I know about the current.row, current.column, and panel.number functions that are useful within panel functions written for lattice. Are there easy ways to obtain the names of the conditioning variables (those appearing after |) and their values for the current panel? Thanks Frank - Frank

Re: [R] Accessor functions in lattice graphics

with writing one based on details like thoseyou have elucidated. But this is really not good programming practice, imo. Contrary and more informed views welcome. Cheers, Bert On Mon, Jun 20, 2011 at 7:17 AM, David Winsemius lt;dwinsem...@comcast.netgt; wrote: Frank Harrell wrote: I know

Re: [R] Regularized logistic regression

. __ 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. - Frank Harrell Department of Biostatistics, Vanderbilt

Re: [R] Rms package - problems with fit.mult.impute

. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Rms-package-problems-with-fit-mult-impute-tp3619814p3621045.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Competing-risks nomogram

-Raffele hospital Milan, Italy - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Competing-risks-nomogram-tp3621907p3622291.html Sent from the R help mailing list archive at Nabble.com

Re: [R] UnoC function in survAUC for censoring-adjusted C-index

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

Re: [R] Competing-risks nomogram

is risky. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Competing-risks-nomogram-tp3621907p3625011.html Sent from the R help mailing list archive at Nabble.com

[R] Issue with dataset inclusion in CRAN packages

the package much more compact. Thanks for listening. I'm writing to r-help because this may applied to other R packages as well. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Issue-with-dataset-inclusion

Re: [R] Issue with dataset inclusion in CRAN packages

I was wrong about this. The dataset is small. Most of the space is taken up by a nice tutorial on rpart.plot. Still I would favor linking to datasets rather than duplicating part of them. Thanks Frank Frank Harrell wrote: I was glad to see the new rpart.plot package by Stephen Milborrow. I

Re: [R] Does psm::Surv handle interval2 data?

-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Does-psm-Surv-handle-interval2-data-tp4655472p4655475.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Does psm::Surv handle interval2 data?

Chris, I've fixed Surv in rms. The fix will be in the next release. For now you can do source('http://biostat.mc.vanderbilt.edu/tmp/Surv.s') after issuing require(rms). Frank Frank Harrell wrote Chris, Thanks for sending the specifics. It appears that I've let Surv in rms fall behind

Re: [R] Difference between R and SAS in Corcordance index in ordinal logistic regression

example and do the calculations by hand. Frank blackscorpio81 wrote Dear R users, Please allow to me ask for your help. I am currently using Frank Harrell Jr package rms to model ordinal logistic regression with proportional odds. In order to assess model predictive ability, C concordance

Re: [R] Difference between R and SAS in Corcordance index in ordinal logistic regression

, Please allow to me ask for your help. I am currently using Frank Harrell Jr package rms to model ordinal logistic regression with proportional odds. In order to assess model predictive ability, C concordance index is displayed and equals to 0.963. This is the code I used with the data attached

Re: [R] Difference between R and SAS in Corcordance index in ordinal logistic regression

For lrm fits, predict(fit, type='mean') predicts the mean Y, not a probability. Frank __ 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

Re: [R] Relative Risk in logistic regression

-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/Relative-Risk-in-logistic-regression-tp4657040p4657163.html Sent from the R help

[R] Regression Modeling Strategies 3-Day Short Course March 2013

*RMS Short Course 2013* Frank E. Harrell, Jr., Ph.D., Professor and Chair Department of Biostatistics, Vanderbilt University School of Medicine *March 4, 5 6, 2013* 8:00am - 4:30pm Student Life Center Board of Trust Room Vanderbilt University Nashville Tennessee USA See

Re: [R] Nomogram after Cox Random Effect (frailty) model

/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

Re: [R] Plotting survival curves after multiple imputation

. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Plotting-survival-curves-after-multiple-imputation-tp4658552p4659505.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Selecting First Incidence from Longitudinal Data

/xgsu00/ __ 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

[R] Multiple left hand side variables in a formula

invocation of model.frame( ) causes R to try to do arithmetic addition to create a single dependent variable. Thanks Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Multiple-left-hand-side-variables

Re: [R] Multiple left hand side variables in a formula

Thanks for your reply Gabor. That doesn't handle a mixture of factor and numeric variables on the left hand side. Frank Gabor Grothendieck wrote On Fri, Mar 1, 2013 at 7:16 PM, Frank Harrell lt; f.harrell@ gt; wrote: The lattice package uses special logic to allow for multiple left-hand

Re: [R] Multiple left hand side variables in a formula

@ [mailto: r-help-bounces@ ] On Behalf Of Frank Harrell Sent: Friday, March 01, 2013 4:17 PM To: r-help@ Subject: [R] Multiple left hand side variables in a formula The lattice package uses special logic to allow for multiple left-hand-side variables in a formula, e.g. y1 + y2 ~ x

Re: [R] Multiple left hand side variables in a formula

Achim this is perfect. I had not seen Formula before. Thanks for writing it! Frank Achim Zeileis-4 wrote On Fri, 1 Mar 2013, Frank Harrell wrote: Thank you Bill. A temporary re-arrangement of the formula will allow me to do the usual subset= na.action= processing afterwards. Nice idea

Re: [R] Multiple left hand side variables in a formula

, Frank Harrell lt; f.harrell@ gt; wrote: The lattice package uses special logic to allow for multiple left-hand-side variables in a formula, e.g. y1 + y2 ~ x. Is there an elegant way to do this outside of lattice? I'm trying to implement a data summarization function that logically takes

[R] Expressions in lattice conditional variables

) : unimplemented type 'expression' in 'HashTableSetup' - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Expressions-in-lattice-conditional-variables-tp4660089.html Sent from the R help mailing list archive

Re: [R] Expressions in lattice conditional variables

Whoops - these 2 lines should have been omitted from the program: n - sprintf('%s (n%s=%g, n%s=%g)', v, nam[1],n[1], nam[2],n[2]) vn[var == v] - n Frank Harrell wrote I would like to have a lattice conditioning ( | var ) variable have expression() as values because I want panel labels

Re: [R] xYplot help

]] __ 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

Re: [R] SAS and R complement each other

I'm not sure why you posted the original note. I quit using SAS in 1991 and haven't needed it yet. Frank RogerJDeAngelis wrote Sorry about the double post. But I keep getting 'post' rejections, so I resubmitted about an hour later. - Frank Harrell Department of Biostatistics

Re: [R] Multiple rms summary plots in a single device

, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Multiple-rms-summary-plots-in-a-single-device-tp4631330p4631423.html Sent from the R help mailing list archive at Nabble.com

<    1   2   3   4   5   >