Re: [R] bwplot change whiskers position to percentile 5 and P95

See also the panel.bpplot function in the Hmisc package. This gives you many options for extended box plots. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/bwplot-change-whiskers-position

[R] Frank Harrell's 2011 RMS Short Course-March 9,10 and 11 (fwd)

2011 Regression Modeling Strategies Short Course by Frank E. Harrell, Jr., Ph.D., Professor and Chair, Department of Biostatistics, Vanderbilt University School of Medicine. 2011 Dates: Wednesday, March 9 (8:00AM-11:00AM), Thursday, March 10 (8:00AM-4:00PM) and Friday, March 11, 2011

Re: [R] Error message in fit.mult.impute (Hmisc package)

I tried your code with the rms package (replacement for the Design package; see http://biostat.mc.vanderbilt.edu/Rrms) and it worked fine. Note that multiple imputation needs the outcome variable in the imputation model. Frank - Frank Harrell Department of Biostatistics, Vanderbilt

Re: [R] Multiple imputation for nominal data

The aregImpute function in the Hmisc package can do this through predictive mean matching and canonical variates (Fisher's optimum scoring algorithm). Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com

Re: [R] An example for predab.resample in r

. Note that .632 only seems to help when you are using an improper scoring rule (e.g., proportion classified correctly), not when using smoother scoring rules such as Brier score, ROC area, etc. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message

Re: [R] Bootstrap confidence intervals using bootcov from the rms package

Unfortunately, bootcov is not meant to operate on fit objects produced by fit.mult.impute. bootcov gets there too late in the process and does not know how to penalize for imputation. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context

Re: [R] leaps::regsubsets p-value

P-values have little meaning in this context. Nor do regression coefficient estimates and especially standard errors (they are biased). Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/leaps

Re: [R] Fit a distribution to an ECDF

Fitting curves to an ECDF will result in a fit that has the same precision as the ECDF if variances are calculated correctly. So why not stop with the ECDF as your estimator? Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http

Re: [R] Plotting number of patients at risk below survival curve

Use the rms package, a replacement for of Design Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Plotting-number-of-patients-at-risk-below-survival-curve-tp3048900p3050315.html Sent from the R help

Re: [R] calculating martingale residual on new data using

predictions (e.g., probability of survival at 2 years) are accurate. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/calculating-martingale-residual-on-new-data-using-predict-coxph-tp3050712p3052377.html Sent from

Re: [R] how do remove those predictor which have p value greater than 0.05 in GLM?

What would make you want to delete a variable because P 0.05? That will invalidate every aspect of statistical inference for the model. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/how-do-remove

Re: [R] xYplot error

If the x-axis variable is really a factor, xYplot will not handle it. You probably need a dot chart instead (see Hmisc's Dotplot). Note that it is unlikely that the confidence intervals are really symmetric. Frank On Tue, 27 Jul 2010, Kang Min wrote: Hi, I'm trying to plot a graph with

Re: [R] how to generate a random data from a empirical distribition

Easiest thing is to sample with replacement from the original data. This is the idea behind the bootstrap, which is sampling from the empirical CDF. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Rank ANCOVA

This has been shown to yield unreliable analyses. Use the more formal proportional odds ordinal logistic model. This is a generalization of the Wiloxon-Mann-Whitney-Kruskal-Wallis statistic. This is implemented in the rms package and elsewhere. Frank E Harrell Jr Professor and Chairman

Re: [R] xYplot error

, method=bands) That looks OK but I can't test it right now. Please continue to have a look, and if you still don't see the problem provide a tiny reproducible example with self-contained data I can access. Frank Thanks for your help. KM On Jul 27, 9:58 pm, Frank Harrell f.harr

Re: [R] how to generate a random data from a empirical distribition

This is true by definition. Read about the bootstrap which may give you some good background information. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Wed, 28 Jul 2010, xin wei wrote: hi,

Re: [R] Problems with normality req. for ANOVA

To add to David's note, the Kruskal-Wallis test is the nonparametric counterpart to one-way ANOVA. You can get a series of K-W tests for several grouping or continuous independent variables (but note these are SEPARATE analyses) using the Hmisc package's spearman2 function. The generalization

Re: [R] Problems with normality req. for ANOVA

In addition the poster did not tell us what is wrong with a nonparametric test. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Mon, 2 Aug 2010, Bert Gunter wrote: My sympathies, but I don't

Re: [R] Specifying interactions in rms package... error

Hi Rob, rms wants symmetry in the sense that the interactions need to use the same number and location of spline knots as the main effects. So if using the * notation omit the main effects (which are generated automatically) and live with the equal knots. Or use the restricted interaction

Re: [R] How to extract se(coef) from cph?

In an upcoming release of the rms package, all fit objects can be printed using LaTeX if putting LaTeX code directly to the console (this is optimized for Sweave). You will be able to say print(fit, latex=TRUE). Frank E Harrell Jr Professor and ChairmanSchool of Medicine

Re: [R] plot the dependent variable against one of the predictors with other predictors as constant

There are many ways to do this. Here is one. install.packages('rms') require(rms) dd - datadist(x, y); options(datadist='dd') f - ols(z ~ x + y) plot(Predict(f))# plot all partial effects plot(Predict(f, x)) # plot only the effect of x plot(Predict(f, y)) # plot only the effect of y f -

Re: [R] Logistic Regression in R (SAS -like output)

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Mon, 9 Aug 2010, Harsh wrote: Hello useRs, I have a problem at hand which I'd think is fairly common amongst groups were R is being adopted for

Re: [R] Logistic Regression in R (SAS -like output)

Note that stepwise variale selection based on AIC has all the problems of stepwise variable selection based on P-values. AIC is just a restatement of the P-Value. Frank Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics

Re: [R] Logistic Regression in R (SAS -like output)

Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Mon, 9 Aug 2010, Kingsford Jones wrote: On Mon, Aug 9, 2010 at 10:27 AM, Frank Harrell f.harr...@vanderbilt.edu wrote: Note that stepwise variale selection based on AIC has all

Re: [R] Identification of Outliners and Extraction of Samples

On Mon, 9 Aug 2010, Alexander Eggel wrote: Hello everybody, I need to know which samples (S1-S6) contain a value that is bigger than the median + five standard deviations of the column he is in. This is just an Why not the 70th percentile plus 6 times the difference in the 85th and 75th

Re: [R] Multiple imputation, especially in rms/Hmisc packages

On Mon, 9 Aug 2010, Mark Seeto wrote: Hello, I have a general question about combining imputations as well as a question specific to the rms and Hmisc packages. The situation is multiple regression on a data set where multiple imputation has been used to give M imputed data sets. I know how

Re: [R] Plotting confidence bands around regression line

Please give the prescription. The article is not available on our extensive online library. I wonder if the method can compete with the bootstrap. Frank Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt

Re: [R] Plotting confidence bands around regression line

Research Laboratory 3400 Spruce St. 7 Maloney Philadelphia, PA 19104 tel. (215) 662-3413 On 2010-08-10 12:29, Frank Harrell wrote: Please give the prescription. The article is not available on our extensive online library. I wonder if the method can compete with the bootstrap. Frank Frank E

Re: [R] extracting the standard error in lrm

On Wed, 11 Aug 2010, david dav wrote: Hi, I would like to extract the coefficients of a logistic regression (estimates and standard error as well) in lrm as in glm with summary(fit.glm)$coef Thanks David coef(fit) sqrt(diag(vcov(fit))) But these will not be very helpful except in the

Re: [R] Multiple imputation, especially in rms/Hmisc packages

tests to be exactly the same? Mark - you can see the code for this at the bottom of anova.rms. Compute W, divide by numerator d.f., then compute P-value using F with numerator and error d.f. Frank Thanks in advance; I really appreciate any help you can give. Mark Frank Harrell wrote

Re: [R] Plotting confidence bands around regression line

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Wed, 11 Aug 2010, Michal Figurski wrote: Peter, Frank, David and others, Thank you all for your ideas. I understand your lack of trust in PB

Re: [R] Plotting confidence bands around regression line

University On Wed, 11 Aug 2010, S Ellison wrote: Frank Harrell f.harr...@vanderbilt.edu 11/08/2010 17:02:03 This problem seems to cry out for one of the many available robust regression methods in R. Not sure that would be much more appropriate, although it would _appear_ to work. The PB

Re: [R] Growth Curves with lmer

Classification accuracy is an improper scoring rule, and one of the problems with it is that the proportional classified correctly can be quite good even if the model uses no predictors. [Hence omitting the intercept is also potentially problematic.] Frank E Harrell Jr Professor and

Re: [R] How to compare the effect of a variable across regression models?

David, In the Cox and many other regression models, the effect of a variable is context-dependent. There is an identifiability problem in what you are doing, as discussed by @ARTICLE{for95mod, author = {Ford, Ian and Norrie, John and Ahmadi, Susan}, year = 1995, title = {Model

Re: [R] val.prob in the Design package - Calibrated Brier Score

Please check the code. I hope that Brier is on the uncalibrated probabilities. Calibrated probabilities are from 1/(1+exp(-[a+b logit(uncalibrated probs)]) where a and b are maximum likelihood estimators (they will be 0 and 1 in training data). Frank Frank E Harrell Jr Professor and

Re: [R] Re : Re : How to compare the effect of a variable across regression models?

. _ De : Bert Gunter gunter.ber...@gene.com À : Biau David djmb...@yahoo.fr Cc : Frank Harrell f.harr...@vanderbilt.edu; r help list r-help@r-project.org Envoyé le : Ven 13 août 2010, 18h 22min 58s Objet : Re: [R] Re : How to compare the effect of a variable across regression models

Re: [R] Stepwise Regression + entry/exit significance level

The values of slentry and slstay that will avoid ruining the statistical properties of the result are slentry=1.0 and slstay=1.0. Frank Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Sat, 14 Aug

Re: [R] How to add lines to lattice plot produced by rms::bplot

Once you guys figure all this out, I'm glad to modify bplot to pass more arguments lattice if needed. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Fri, 13 Aug 2010, David Winsemius wrote:

Re: [R] confidence Intervals for predictions in GLS

install.packages('rms') require(rms) ?Gls ?plot.Predict Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Sat, 14 Aug 2010, Camilo Mora wrote: Hi everyone: Is there a function in R to calculate

Re: [R] How to add lines to lattice plot produced by rms::bplot

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Sat, 14 Aug 2010, David Winsemius wrote: On Aug 14, 2010, at 9:59 AM, Frank Harrell wrote: Once you guys figure all this out, I'm glad

Re: [R] HMisc/rms package questions

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Tue, 17 Aug 2010, Rob James wrote: 1) How does one capture the plots from the plsmo procedure? Simply inserting a routing call to a graphical

Re: [R] R reports

What do low level proc print and proc report have on Sweave or http://biostat.mc.vanderbilt.edu/wiki/pub/Main/StatReport/summary.pdf? If proc print and proc report are 4G, let's move back a generation. Frank E Harrell Jr Professor and ChairmanSchool of Medicine

Re: [R] ROCR predictions

At the heart of this you have a problem in incomplete conditioning. You are computing things like Prob(X x) when you know X=x. Working with a statistician who is well versed in probability models will undoubtedly help. Frank Frank E Harrell Jr Professor and ChairmanSchool of

Re: [R] logistic regression tree

It would be good to tell us of the frequency of observations in each category of Y, and the number of continuous X's. Recursive partitioning will require perhaps 50,000 observations in the less frequent Y category for its structure and predicted values to validate, depending on X and the

Re: [R] logistic regression tree

On Fri, 20 Aug 2010, Kay Cichini wrote: hello, my data-collection is not yet finished, but i though have started investigating possible analysis methods. below i give a very close simulation of my future data-set, however there might be more nominal explanatory variables - there will be no

Re: [R] R reports

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Sat, 21 Aug 2010, Donald Paul Winston wrote: Sweave and LaTex is way to much overhead to deal with. There should be a built in standard report()

Re: [R] R reports

Your notes are not well thought out. You'll find that r-help is a friendly place for new users that do not come in with an attitude. I once used SAS (for 23 years) and know it very well. I wrote the first SAS procedures for a graphics device, percentiles, logistic regression, and Cox

Re: [R] R reports

, Frank Harrell wrote: Your notes are not well thought out. You'll find that r-help is a friendly place for new users that do not come in with an attitude. I once used SAS (for 23 years) and know it very well. I wrote the first SAS procedures for a graphics device, percentiles, logistic

Re: [R] R reports

On Sat, 21 Aug 2010, Donald Paul Winston wrote: Good grief. Adding a report function is not going to make R less flexible. Don't you want to use a tool that's relevant to the rest of the world? That world is much bigger then your world. This is ridiculous. Looks like some people are

Re: [R] AUC

Samuel, Since the difference in AUCs has insufficient power and doesn't really take into account the pairing of predictions, I recommend the Hmisc package's rcorrp.cens function. Its method has good power and asks the question is one predictor more concordant than the other in the same

Re: [R] modify a nomogram

Update to the rms package which is the version now being actively supported. New features will not be added to Design. The nomogram function in rms separates the plotting into a plot method for easier understanding. You can control all axes - some experimentation can tell you if you can

Re: [R] Speeding up prediction of survival estimates when using `survifit'

Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Mon, 30 Aug 2010, Ravi Varadhan wrote: Hi, I fit a Cox PH model to estimate the cause-specific hazards (in a competing risks setting). Then , I

Re: [R] [Q] Goodness-of-fit test of a logistic regression model using rms package

don't feel bad. You just know that something somewhere is probably wrong with the model. I focus on directed tests such as allowing all continuous variables to have nonlinear effects or allowing selected interactions, and finding out how important the complex model terms are. Frank Harrell

Re: [R] Is there any package or function perform stepwise variable selection under GEE method?

The absence of stepwise methods works to your advantage, as these yield invalid statistical inference and inflated regression coefficients. Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt University On Thu,

Re: [R] Comparing COXPH models, one with age as a continuous variable, one with age as a three-level factor

On Thu, 2 Sep 2010, stephenb wrote: sorry to bump in late, but I am doing similar things now and was browsing. IMHO anova is not appropriate here. it applies when the richer model has p more variables than the simpler model. this is not the case here. the competing models use different

Re: [R] Linear Logistic Regression - Understanding the output (and possibly the test to use!)

the time to consult a statistician is before you have done any statistical analysis. Frank Harrell I have a couple of Kleinbaum's (et al) other texts and find them to be well written and reasoned, so I suspect the citation above would be as accessible as any. Thank you, that is useful

Re: [R] some questions about longitudinal study with baseline

Baseline should appear only as a baseline and should be removed from the set of longitudinal responses. This is often done with a merge( ) operation. Frank Frank E Harrell Jr Professor and ChairmanSchool of Medicine Department of Biostatistics Vanderbilt

[R] Where to find R-help options web page

. Nabble sent a link to turn off e-mail but the r-help mail service rejected nabble's command. Thanks Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Where-to-find-R-help-options-web-page-tp2535123p2535123.html

Re: [R] Where to find R-help options web page

(and to be able to post on it). Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Where-to-find-R-help-options-web-page-tp2535123p2535148.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Where to find R-help options web page

Aha. Thanks David. I somehow thought that that whole section was for administrators only. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Where-to-find-R-help-options-web-page-tp2535123p2535172.html

Re: [R] SF-8 (not 36) questionnaire scoring for R?

I know someone who has R code for SF-36 and perhaps SF-12. Aren't there copyright issues relating to SF-* even if it is reprogrammed? Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/SF-36

Re: [R] SF-8 (not 36) questionnaire scoring for R?

Yes the company behind that probably received federal funds for some of the research and has been very careful to minimize their contribution to the community. I didn't understand your parenthetical remark. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View

[R] [R-pkgs] New version of rms package on CRAN

CRAN has a significant update to rms. Windows and unix/linux versions are available and I expect the Mac version to be available soon. The most significant improvement is addition of latex=TRUE arguments to model fitting print methods, made especially for use with Sweave. Here is a summary

Re: [R] predict.lrm ( Design package)

You sent a private note about this which I just took the time to answer. Please send only one note, and please post my reply to you to r-help. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/predict

Re: [R] predict.lrm ( Design package)

Thanks Peter. I think you're right. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/predict-lrm-Design-package-tp2546894p2547146.html Sent from the R help mailing list archive at Nabble.com

Re: [R] predict.lrm ( Design package) poor performance?

the difficulty of the task and give a lower value. No cause for alarm. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/predict-lrm-Design-package-tp2546894p2550118.html Sent from the R help mailing list

Re: [R] Validation / Training - test data

Split sample validation is highly unstable with your sample size. The rms package can help with bootstrapping or cross-validation, assuming you have all modeling steps repreated for each resample. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View

Re: [R] Validation / Training - test data

It all depends on the ultimate use of the results. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Validation-Training-test-data-tp2718523p2719370.html Sent from the R help mailing list archive

Re: [R] Interpreting the example given by Frank Harrell in the predict.lrm {Design} help

- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Interpreting-the-example-given-by-Frank-Harrell-in-the-predict-lrm-Design-help-tp2883311p2891623.html Sent from the R help mailing list archive at Nabble.com

Re: [R] Interpreting the example given by Frank Harrell in the predict.lrm {Design} help

Why assign them at all? Is this a forced choice at gunpoint problem? Remember what probabilities mean. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Interpreting-the-example-given-by-Frank-Harrell

Re: [R] Interpreting the example given by Frank Harrell in the predict.lrm {Design} help

Well put Greg. The job of the statistician is to produce good estimates (probabilities in this case). Those cannot be translated into action without subject-specific utility functions. Classification during the analysis or publication stage is not necessary. Frank - Frank Harrell

Re: [R] Interpreting the example given by Frank Harrell in the predict.lrm {Design} help

You still seem to be hung up on making arbitrary classifications. Instead, look at tendencies using odds ratios or rank correlation measures. My book Regression Modeling Strategies covers this. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message

Re: [R] A problem about nomogram--thank you for you help

Please take the time to study the subject matter, and note that a nomogram is just a graphical method. It is not a statistical model or a process. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com

Re: [R] Interpreting the example given by Frank Harrell in the predict.lrm {Design} help

get predicted mean Y (see the Mean.lrm in the R rms package, a replacement for the Design package). Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Interpreting-the-example-given-by-Frank-Harrell

Re: [R] multiple moderated regression steps

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

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

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

, 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] GLS - Plotting Graphs with 95% conf interval

. 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] Named numeric vectors with the same value but different names return different results when used as thresholds for calculating true positives

. - 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] Sum weights of independent variables across models (AIC)

/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

]] __ 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] Design Survplot performance

-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

[R] Error in plotmath

=.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] How to make a nomogam and Calibration plot

! 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

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

.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

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] Goodness of fit of binary logistic model

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] Goodness of fit of binary logistic model

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

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] help annotating plot.xmean.ordinaly in rms package

-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

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] Revolutions Blog: July Roundup

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

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

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] fit.mult.impute() in Hmisc

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

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

/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()

://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] contrast package with interactions in gls model

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

: 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] contrast package with interactions in gls model

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] Labelling all variables at once (using Hmisc label)

-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

  1   2   3   4   5   >