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

2011-11-29 Thread Frank Harrell
' ), 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] fitting a distribution to a ecdf plot

2011-11-28 Thread Frank Harrell
__ 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] SPSS - R

2011-11-26 Thread Frank Harrell
, 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] The contrast and Design libraries

2011-11-25 Thread Frank Harrell
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] Thank you

2011-11-24 Thread Frank Harrell
://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] Correlation matrix removing insignificant R values

2011-11-23 Thread Frank Harrell
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] Similar function for Redun() from Hmisc ?

2011-11-22 Thread Frank Harrell
= 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?

2011-11-22 Thread Frank Harrell
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] the standard error of the quantile

2011-11-19 Thread Frank Harrell
to estimate the standard error of the quantile in survival analysis. Can anyone help me,please? - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/the-standard-error-of-the-quantile-tp4086479p4086494.html Sent from

Re: [R] the standard error of the quantile

2011-11-19 Thread Frank Harrell
. 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] linear against nonlinear alternatives - quantile regression

2011-11-16 Thread Frank Harrell
, Frank Harrell lt;f.harrell@gt; wrote: Just to address a piece of this - in the case in which you are currently focusing on only one quantile, the rms package can help by fitting restricted cubic splines for covariate effects, and then run anova to test for nonlinearity (sometimes a dubious

Re: [R] Odds ratios from lrm plot

2011-11-11 Thread Frank Harrell
wrote: The code library(Design) f - lrm(y~x1+x2+x1*x2, data=data) plot(f) produces a plot of log odds vs x2 with 0.95 confidence intervals. How do I get a plot of odds ratios vs x2 instead? Thanks - Frank Harrell Department of Biostatistics, Vanderbilt University -- View

Re: [R] linear against nonlinear alternatives - quantile regression

2011-11-06 Thread Frank Harrell
of the nonlinearity would be imposed at all levels of tau. -- David On Nov 5, 2011, at 10:43 AM, Frank Harrell lt;f.harrell@gt; wrote: Just to address a piece of this - in the case in which you are currently focusing on only one quantile, the rms package can help by fitting

Re: [R] linear against nonlinear alternatives - quantile regression

2011-11-05 Thread Frank Harrell
]] __ 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

Re: [R] Select some, but not all, variables stepwise

2011-11-04 Thread Frank Harrell
. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Select-some-but-not-all-variables-stepwise-tp3990002p3990026.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Select some, but not all, variables stepwise

2011-11-04 Thread Frank Harrell
Why do any of that? Those procedures are statistically invalid. Frank AndreE wrote: Thanks Francisco - I've actually realized that ?step can do pretty much exactly what I want. Andre - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context

Re: [R] How to interpret Spearman Correlation

2011-11-02 Thread Frank Harrell
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/How-to-interpret-Spearman-Correlation

Re: [R] How to use IML with R and SAS

2011-10-31 Thread Frank Harrell
/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

Re: [R] ACE/AVAS with bootstrap

2011-10-28 Thread Frank Harrell
. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/ACE-AVAS-with-bootstrap-tp3948613p3950068.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] interpreting bootstrap corrected slope [rms package]

2011-10-24 Thread Frank Harrell
behaves in a non-ordinal fashion. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/interpreting-bootstrap-corrected-slope-rms-package-tp3928314p3932923.html Sent from the R help mailing list archive at Nabble.com

Re: [R] interpreting bootstrap corrected slope [rms package]

2011-10-23 Thread Frank Harrell
the model as is would be fine for future predictions? Thanks for your time, Adam - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/interpreting-bootstrap-corrected-slope-rms-package-tp3928314p3930552.html Sent

Re: [R] interpreting bootstrap corrected slope [rms package]

2011-10-23 Thread Frank Harrell
in the initial and final model? - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/interpreting-bootstrap-corrected-slope-rms-package-tp3928314p3931493.html Sent from the R help mailing list archive at Nabble.com

Re: [R] interpreting bootstrap corrected slope [rms package]

2011-10-22 Thread Frank Harrell
, minimal, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/interpreting-bootstrap-corrected-slope-rms-package-tp3928314p3928467.html Sent from the R help mailing list archive

Re: [R] Cross posting (was Restricted Cubic Splines within survfit.cph)

2011-10-12 Thread Frank Harrell
Terry I would just add that if someone contacts the maintainer and does not follow the advice of the maintainer, still has problems, and posts a message to the list, then some time is wasted. Frank On 10/12/2011 07:38 AM, Terry Therneau wrote: You wrote: It may be best to either write to

Re: [R] How to test if two C statistics are significantly different?

2011-10-11 Thread Frank Harrell
__ 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

Re: [R] restricted cubic spline within survfit.cph in the package rms

2011-10-11 Thread Frank Harrell
-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/restricted-cubic-spline-within-survfit-cph-in-the-package-rms

Re: [R] Generalized Additive Models: How to create publication-ready regression tables

2011-10-09 Thread Frank Harrell
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

Re: [R] About stepwise regression problem

2011-10-07 Thread Frank Harrell
(mintemp,bs=cr)+s(RH,bs=cr)+s(solar,bs=cr)+s(windspeed,bs=cr)+s(transport,bs=cr),family=gaussian (link=log),groupD,methods=REML) , is it generalized additive mixed model? 3. what the different if I use other criteria such as AIC or BIC? Anyway, thank all of you! - Frank Harrell

Re: [R] model selection using logistf package

2011-09-27 Thread Frank Harrell
the logistf package but as I stepwise simplify the model I don't know if the simplification is motivated or not... Can anyone explain to me how I should approach this problem? I have data with y=0 and 1 and three x variables for 328 observations. Thanks! - Frank Harrell Department

Re: [R] Is there a latex summary function in the quantreg package for just 1 tau?

2011-09-27 Thread Frank Harrell
]] __ 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

Re: [R] Standardized beta-coefficients in regression

2011-09-26 Thread Frank Harrell
with B.R. but unfortunately the life of an applied statistician is complex sometimes :-) - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Standardized-beta-coefficients-in-regression-tp791616p3843763.html Sent

Re: [R] How to turn a LaTeX Sweave file (Rnw) into .HTML/.odf/.docx? (under windows)

2011-09-23 Thread Frank Harrell
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/How-to-turn-a-LaTeX-Sweave-file-Rnw-into-HTML-odf-docx-under-windows-tp3835230p3836939.html

Re: [R] Could R run Cox proportional hazard model with multiple failure time data

2011-09-18 Thread Frank Harrell
You sent that question to me individually as well as posting to the group. Kindly pick the most appropriate place and send to only that. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Could-R-run-Cox

Re: [R] ROCR package question for evaluating two regression models

2011-09-03 Thread Frank Harrell
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/ROCR-package-question

Re: [R] How to retrieve bias-corrected probability from calibrate.rms

2011-09-01 Thread Frank Harrell
]] __ 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

Re: [R] How to retrieve bias-corrected probability from calibrate.rms

2011-09-01 Thread Frank Harrell
Shanghai, China* 2011/9/1 Frank Harrell lt;f.harr...@vanderbilt.edugt; cal - calibrate(fit, ...); note that cal is a matrix. colnames(cal) will tell you what to pick, in this case cal[,'calibrated.corrected']. Be sure to follow the posting guide. Frank yz wrote: Dear R users

Re: [R] Differences in SAS and R defaults

2011-08-29 Thread Frank Harrell
-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/Differences-in-SAS-and-R-defaults-tp3776102p3776621.html Sent from the R help mailing

Re: [R] Ordinal logistic regression p-values

2011-08-29 Thread Frank Harrell
, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Ordinal-logistic-regression-p-values-tp3777674p368.html Sent from the R help mailing list archive at Nabble.com

Re: [R] How to find the accuracy of the predicted glm model with family = binomial (link = logit)

2011-08-26 Thread Frank Harrell
/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/How-to-find-the-accuracy-of-the-predicted-glm-model-with-family-binomial-link

Re: [R] survplot() for cph(): Design vs rms

2011-08-25 Thread Frank Harrell
://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

Re: [R] How to do cross validation with glm?

2011-08-24 Thread Frank Harrell
-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/How-to-do-cross-validation-with-glm-tp3765994p3766108.html Sent from the R help

Re: [R] val.surv

2011-08-21 Thread Frank Harrell
, fun = fun, n = nrow(S), d = sum(S[, 2]), units = units), class = val.survh) return(w) But a much better option would be to report the error to Frank Harrell. I'm copying him since I think your .txt files probably reached the list as well as my mailbox

Re: [R] rank analysis - reinventing the wheel?

2011-08-21 Thread Frank Harrell
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] Labelling all variables at once (using Hmisc label)

2011-08-19 Thread Frank Harrell
-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/Labelling-all-variables-at-once-using-Hmisc-label-tp3745660p3756459.html Sent from the R

Re: [R] Hmisc::rcorr on a 'data.frame'?

2011-08-19 Thread Frank Harrell
commented, minimal, self-contained, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Hmisc-rcorr-on-a-data-frame-tp3755868p3756462.html Sent from the R help mailing list archive at Nabble.com

Re: [R] help Dxy and C-index calculation

2011-08-19 Thread Frank Harrell
Replace Design with rms (for general reasons not related to your question; Design is about to be obsolete). Negate Dxy. The linear predictor for the Cox model is relative log hazard. Higher hazard means shorter survival time. For other survival models the model is stated in terms of

Re: [R] contrast package with interactions in gls model

2011-08-18 Thread Frank Harrell
interaction with function contrast of contrast package. I thanks for your reply, Regards. Marylin Bejarano - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/contrast-package-with-interactions-in-gls-model

Re: [R] contrast package with interactions in gls model

2011-08-17 Thread Frank Harrell
://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/contrast-package-with-interactions-in-gls-model

Re: [R] Labelling all variables at once (using Hmisc label)

2011-08-17 Thread Frank Harrell
: 11  Date: Tue, 16 Aug 2011 04:22:07 -0700 (PDT)  From: Frank Harrell lt;f.harr...@vanderbilt.edugt;  To: r-help@r-project.org  Subject: Re: [R] Labelling all variables at once (using Hmisc label)  Message-ID: 1313493727519-3746928.p...@n4.nabble.com  Content-Type: text/plain; charset=UTF-8

Re: [R] Labelling all variables at once (using Hmisc label)

2011-08-16 Thread Frank Harrell
/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

Re: [R] calibration curve for cph()

2011-08-16 Thread Frank Harrell
://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] fit.mult.impute() in Hmisc

2011-08-13 Thread Frank Harrell
iterations: 5 - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/fit-mult-impute-in-Hmisc-tp3419037p3741881.html Sent from the R help mailing list archive at Nabble.com

[R] useR! 2012 Nashville TN: Your Input Needed

2011-08-10 Thread Frank Harrell
So far we have received over 70 responses to our survey. If you have NOT responded already and are likely to attend useR! 2012, please take the extremely short survey today. The link is below. More information about Nashville may be seen at http://biostat.mc.vanderbilt.edu/UseR-2012 In a day

Re: [R] Revolutions Blog: July Roundup

2011-08-09 Thread Frank Harrell
clue]. It should have used the normal distribution. Cheers, Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Revolutions-Blog-July-Roundup-tp3731426p3731597.html Sent from the R help mailing list archive

Re: [R] help annotating plot.xmean.ordinaly in rms package

2011-08-07 Thread Frank Harrell
-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/help-annotating-plot-xmean-ordinaly-in-rms-package-tp3724864p3724937.html Sent from

[R] useR!2012 Nashville TN: Your Input Needed

2011-08-07 Thread Frank Harrell
Dear R and S-Plus Community: The 2012 R User Conference - useR!2012 - will be held in Nashville Tennessee USA, June 12-15, 2012 on the campus of Vanderbilt University. We would like to begin estimating the number of attendees, their area of interest, and the number seeking hotel vs.

Re: [R] Goodness of fit of binary logistic model

2011-08-06 Thread Frank Harrell
exists for that case. Paul On Fri, Aug 5, 2011 at 7:58 PM, Frank Harrell lt;f.harr...@vanderbilt.edugt; wrote: Please provide the data or better the R code for simulating the data that shows the problem. Then we can look further into this. Frank - Frank Harrell Department

Re: [R] help with predict for cr model using rms package

2011-08-06 Thread Frank Harrell
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/help-with-predict

Re: [R] Goodness of fit of binary logistic model

2011-08-05 Thread Frank Harrell
Please provide the data or better the R code for simulating the data that shows the problem. Then we can look further into this. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Goodness-of-fit-of-binary

Re: [R] How to make a nomogam and Calibration plot

2011-08-03 Thread Frank Harrell
.789695.n4.nabble.com/file/n3714477/untitled.jpg http://r.789695.n4.nabble.com/file/n3714477/%E9%99%84%E4%BB%B62.jpg http://r.789695.n4.nabble.com/file/n3714477/Dataset.xls Dataset.xls - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http

Re: [R] How to make a nomogam and Calibration plot

2011-08-03 Thread Frank Harrell
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/How-to-make-a-nomogam-and-Calibration-plot

Re: [R] How to make a nomogam and Calibration plot

2011-08-01 Thread Frank Harrell
! Ping Tang http://r.789695.n4.nabble.com/file/n3710068/Dataset.xls Dataset.xls http://r.789695.n4.nabble.com/file/n3710068/R_help.doc R_help.doc - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/How

[R] Error in plotmath

2011-07-31 Thread Frank Harrell
=.5) The problem also occurs if I use instead ylab=expression(prob(X = x))) All is well if I remove = but I need =. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Error-in-plotmath-tp3708153p3708153.html

Re: [R] Design Survplot performance

2011-07-21 Thread Frank Harrell
-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/Design-Survplot-performance-tp3684580p3684626.html Sent from the R

Re: [R] Sum weights of independent variables across models (AIC)

2011-07-14 Thread Frank Harrell
/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/Sum-weights-of-independent-variables-across-models-AIC-tp3666306p3668513.html

Re: [R] validate survival with val.surv

2011-07-14 Thread Frank Harrell
]] __ 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

Re: [R] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

2011-07-12 Thread Frank Harrell
. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Named-numeric-vectors-with-the-same-value-but-different-names-return-different-results-when-used-as-s-tp3660833p3662030.html Sent from the R help mailing list archive

Re: [R] GLS - Plotting Graphs with 95% conf interval

2011-07-11 Thread Frank Harrell
. Thanks - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/GLS-Plotting-Graphs-with-95-conf-interval-tp3659814p3660953.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Speed Advice for R --- avoid data frames

2011-07-06 Thread Frank Harrell
, reproducible code. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Speed-Advice-for-R-avoid-data-frames-tp3640932p3648681.html Sent from the R help mailing list archive at Nabble.com

Re: [R] multiple moderated regression steps

2011-07-01 Thread Frank Harrell
at second covariance block or third covariance block. like this,for value of Covariance of IV*Mod1, IV*Mod2 coefficients which covariance block must be taken as correct??? thaks so much Note: excuse me for my English - Frank Harrell Department of Biostatistics, Vanderbilt

Re: [R] multiple moderated regression steps

2011-07-01 Thread Frank Harrell
for interactions. for plottin slopes unstandardised regression coefficient were used. - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/multiple-moderated-regression-steps-tp3637807p3638669.html Sent from the R

[R] Issue with dataset inclusion in CRAN packages

2011-06-26 Thread Frank Harrell
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

2011-06-26 Thread Frank Harrell
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] Competing-risks nomogram

2011-06-25 Thread Frank Harrell
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

Re: [R] Competing-risks nomogram

2011-06-24 Thread Frank Harrell
-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

2011-06-24 Thread Frank Harrell
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] Rms package - problems with fit.mult.impute

2011-06-23 Thread Frank Harrell
. - 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] Regularized logistic regression

2011-06-21 Thread Frank Harrell
. __ 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] Accessor functions in lattice graphics

2011-06-20 Thread Frank Harrell
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

[R] Accessor functions in lattice graphics

2011-06-19 Thread Frank Harrell
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

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

2011-06-16 Thread Frank Harrell
, 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...

2011-06-16 Thread Frank Harrell
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

Re: [R] Somers Dyx

2011-06-13 Thread Frank Harrell
]] __ 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] rcspline.plot query

2011-06-12 Thread Frank Harrell
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] rcspline.plot query

2011-06-11 Thread Frank Harrell
__ 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] Logistic Regression

2011-06-10 Thread Frank Harrell
/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' = ?

2011-06-10 Thread Frank Harrell
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] predict with model (rms package)

2011-06-08 Thread Frank Harrell
, 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

2011-06-07 Thread Frank Harrell
? - 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] Interpreting Quantile Regression

2011-06-04 Thread Frank Harrell
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

[R] New version of rms package on CRAN

2011-06-03 Thread Frank Harrell
. 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] Removal of elements from nomograms

2011-06-02 Thread Frank Harrell
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] Problem with % in an example when running R CMD check

2011-05-31 Thread Frank Harrell
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] Several Regression by combinations variables

2011-05-30 Thread Frank Harrell
, 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

2011-05-30 Thread Frank Harrell
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

2011-05-30 Thread Frank Harrell
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

Re: [R] Multinomial Logistical Model

2011-05-25 Thread Frank Harrell
+ 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] How to calculate confidence interval of C statistic by rcorr.cens

2011-05-22 Thread Frank Harrell
://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

2011-05-22 Thread Frank Harrell
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] changes in coxph in survival from older version?

2011-05-19 Thread Frank Harrell
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

<    1   2   3   4   5   >