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

2011-08-06 Thread Frank Harrell
ether the model > fits. > > Since the model is basically a full model for a 2x2 table, it is not > surprising to me that "goodness of fit" tests behave poorly. In fact, I > would conjecture that no sensible g.o.f. test exists for that case. > >> >> Paul >

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

2011-08-06 Thread Frank Harrell
> >> p.cr <- predict(ord.cr, all.data2.stand, type='mean', codes=TRUE) > Error in model.frame.default(Terms, newdata, na.action = na.action, ...) : > variable lengths differ (found for 'mean.ova.energy') > In addition: Warning message: > 'newdata

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

2011-08-07 Thread Frank Harrell
ect.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, Vand

[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. lower-

Re: [R] Revolutions Blog: July Roundup

2011-08-09 Thread Frank Harrell
sed is any 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 li

[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 or

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

2011-08-13 Thread Frank Harrell
binomial family taken to be 1) Null deviance: 80.508 on 80 degrees of freedom Residual deviance: 57.965 on 77 degrees of freedom AIC: 65.97 Number of Fisher Scoring iterations: 5 - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in c

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

2011-08-16 Thread Frank Harrell
arah Goslee > http://www.functionaldiversity.org > [[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > ht

Re: [R] calibration curve for cph()

2011-08-16 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,

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

2011-08-17 Thread Frank Harrell
uld be grateful if somebody tell me I'm doing wrong with my contrast > function script. > > Thanks in advance, > > > Marylin Bejarano > PHd candidate in Ecology Institute of Mexico's National Autonomous > University > > [[alternative HTML version d

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

2011-08-17 Thread Frank Harrell
quot;Warning: data and labels are not of > same length\n") > return(x) > } > > Thanks > >> Message: 11 >> Date: Tue, 16 Aug 2011 04:22:07 -0700 (PDT) >> From: > Frank Harrell <f.harr...@vanderbilt.edu> >> To: > r-help@r-project.org >

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

2011-08-18 Thread Frank Harrell
rying to contrast > 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.

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

2011-08-19 Thread Frank Harrell
osting 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/Labelling-all-variables-

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

2011-08-19 Thread Frank Harrell
st.net/~garbl/stylemanual/e.htm#e-mail > > __ > 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, minim

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 surv

Re: [R] val.surv

2011-08-21 Thread Frank Harrell
d the fitted values but is failing in the situation where > the nrows of the newdata does not match that of the fit. > > At this point one option might be trying to generate a structure that > matches the val.surv output by going through the code and building it > up bit by

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

2011-08-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 Departm

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

2011-08-24 Thread Frank Harrell
g/posting-guide.html >> > and provide commented, minimal, self-contained, >> reproducible code. >> > >> >> -- Brian D. Ripley,          >>         rip...@stats.ox.ac.uk >> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/ >> University of Oxford,          >>  

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

2011-08-25 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. &g

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

2011-08-26 Thread Frank Harrell
ovide commented, minimal, self-contained, reproducible code. >> > > __ > 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.htm

Re: [R] Differences in SAS and R defaults

2011-08-29 Thread Frank Harrell
elp > 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/D

Re: [R] Ordinal logistic regression p-values

2011-08-29 Thread Frank Harrell
o 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/Ordinal-logi

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

2011-09-01 Thread Frank Harrell
anghai, China* > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html &

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

2011-09-01 Thread Frank Harrell
gt; *Yao Zhu* > *Department of Urology > Fudan University Shanghai Cancer Center > Shanghai, China* > > > 2011/9/1 Frank Harrell <f.harr...@vanderbilt.edu> > >> cal <- calibrate(fit, ...); note that cal is a matrix. colnames(cal) >> will >> tell you wha

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

2011-09-03 Thread Frank Harrell
with that? > Thanks a lot,Andra > > > > [[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide >

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] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Frank Harrell
>> Robert W. Baer, Ph.D. >> Professor of Physiology >> Kirksville College of Osteopathic Medicine >> A. T. Still University of Health Sciences >> 800 W. Jefferson St. >> Kirksville, MO 63501 >> 660-626-2322 >> FAX 660-626-2965 >> >> _

Re: [R] Odds Ratios in rms package

2012-06-20 Thread Frank Harrell
obtained, especially in the > presence of interactions? > > Best regards, > Sebastian > > __ > 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 >

Re: [R] Odds Ratios in rms package

2012-06-21 Thread Frank Harrell
elp >> PLEASE do read the posting guide >> http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > -- > Peter Dalgaard, Professor > Center for Statistics, Copenhagen Business School > Solbjerg Plads 3, 2000

Re: [R] rms package-superposition prediction curve of ols and data points

2012-06-26 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. > > David Winse

Re: [R] Simple mean trajectory (ordinal variable)

2012-06-29 Thread Frank Harrell
gt; >> __ >> 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. > > David Winse

Re: [R] Hosmer-Lemeshow test for Cox model

2012-07-05 Thread Frank Harrell
0) { > cutyhat = cut(yhat, > breaks = quantile(yhat, probs=seq(0, >1, 1/g)), include.lowest=TRUE) > obs = xtabs(cbind(1 - y, y) ~ cutyhat) > expect = xtabs(cbind(1 - yhat, yhat) ~ cutyhat) > chisq = sum((obs - expect)^2/expect) > P = 1 - pchisq(chisq, g -

Re: [R] Plotting the probability curve from a logit model with 10 predictors

2012-07-09 Thread Frank Harrell
t;>> and provide commented, minimal, self-contained, reproducible code. >>>> >>> >>> >>> >>> -- >>> >>> Bert Gunter >>> Genentech Nonclinical Biostatistics >>> >>> Internal Contact Info: >>>

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

2011-06-16 Thread Frank Harrell
(.5,p+.25), type='n', axes=FALSE, 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),

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

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

[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

Re: [R] Accessor functions in lattice graphics

2011-06-20 Thread Frank Harrell
> 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. >> > > > > -- > "Men by nature long t

Re: [R] Regularized logistic regression

2011-06-21 Thread Frank Harrell
ting guide >> http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. >> > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEAS

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

2011-06-23 Thread Frank Harrell
d]] > > __ > 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, repr

Re: [R] Competing-risks nomogram

2011-06-24 Thread Frank Harrell
> Many thanks > > F.Abdollah, MD > San-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

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

2011-06-24 Thread Frank Harrell
iour here. > > Many thanks > > Eleni Rapsomaniki > > Research Associate/Statistician > Strangeways Research Laboratory > Department of Public Health and Primary Care > > University of Cambridge >   > > __ > R-help@r-project

Re: [R] Competing-risks nomogram

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

[R] Issue with dataset inclusion in CRAN packages

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

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 Mil

Re: [R] multiple moderated regression steps

2011-07-01 Thread Frank Harrell
uestion: In three way spreadshhet it is necessary to enter > variance of IV*Moderator coefficient. ın matrıx which covariance block is > correct to enter this value? should ı look at second covariance block or > third covariance block. > like this,for value of Covariance of IV*Mod1, IV*Mod

Re: [R] multiple moderated regression steps

2011-07-01 Thread Frank Harrell
in the coefficient table) and slope test were done 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-mod

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

2011-07-06 Thread Frank Harrell
t;>> data frame into a matrix, do the operations, and then copy them back. >>>> not ideal, either. >>>> >>>> In my opinion, this is an R design flow. Data frames are the >>>> fundamental unit of much statistical analysis,

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

2011-07-11 Thread Frank Harrell
mial, I am getting the same line of fit regardless of what is > on my X axis (depsite different variables have positive and negative > relationships). > > any imput on any of these problems would be appreciated. > > Thanks > - Frank Harrell Department of Biostatistics,

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

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 of Biostatisti

Re: [R] validate survival with val.surv

2011-07-14 Thread Frank Harrell
China* > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and

Re: [R] Design Survplot performance

2011-07-21 Thread Frank Harrell
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, Vand

[R] Error in plotmath

2011-07-31 Thread Frank Harrell
quot;, X <= x, "]"))) lines(ecdf(x), cex=.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

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

2011-08-01 Thread Frank Harrell
e his/her time > to help for my questions. > > Warmest regards! > > 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, Vanderbil

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

2011-08-03 Thread Frank Harrell
> probability of cancer using R? > > 2. How to make the Calibration plot (similar to the below figure 2) which > used to validate the prediction accuracy of the nomogram using R? And how > to calculate the concordance index (C-index) ? > > > http://r.789695.n4.nabble.

Re: [R] Hire person to convert R code to SAS

2012-07-23 Thread Frank Harrell
g to look at this, or could point > me to someone who would, it would be very much appreciated. > > Thanks in advance! > > -will > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/H

Re: [R] stepwise in svyglm???

2012-07-24 Thread Frank Harrell
HTML version deleted]] > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible c

Re: [R] stepwise in svyglm???

2012-07-25 Thread Frank Harrell
is it valid? ;-) > > -- Bert > > On Tue, Jul 24, 2012 at 3:00 PM, Frank Harrell <f.harrell@> wrote: >> Stepwise variable selection is invalid in this context. >> Frank >> >> Diana Marcela Martinez Ruiz wrote >>> >>> Hello, >>> &g

Re: [R] optimal cut off with Epi package

2012-07-26 Thread Frank Harrell
ain="ROC with Epi package", MI=TRUE, MX=TRUE, PV=TRUE) > > > > ### REPORT: lr.eta=0.491 > > > [[alternative HTML version deleted]] > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/

Re: [R] Hmisc improveProb() and PredictABEL reclassification () function and continuous NRI

2012-07-30 Thread Frank Harrell
ailing 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] sensitivity and specificity in svyglm??

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

Re: [R] error: in catg (xi, name=nam, label=lab): "LO2" has <2 category levels

2012-09-07 Thread Frank Harrell
st predictors are (dichotomous) nominal variables as is the problematic > "LO2". Does anyone know what the problem is and how I can correct it? > > Kind regards, > > Tobias > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message i

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

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

2011-11-16 Thread Frank Harrell
t; > -- > David > > On Nov 5, 2011, at 10:43 AM, Frank Harrell <f.harrell@> 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 spli

Re: [R] the standard error of the quantile

2011-11-19 Thread Frank Harrell
I need 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-tp4086479

Re: [R] the standard error of the quantile

2011-11-19 Thread Frank Harrell
r help. > 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) >

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

2011-11-22 Thread Frank Harrell
> V <- redun(~., data = data.frame, r2 = 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 Harr

Re: [R] Bezier curve in R?

2011-11-22 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 Biostati

Re: [R] Correlation matrix removing insignificant R values

2011-11-23 Thread Frank Harrell
in > the instances were the R value 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/Correlat

Re: [R] Thank you

2011-11-24 Thread Frank Harrell
development.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm > > __ > 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, s

Re: [R] The contrast and Design libraries

2011-11-25 Thread Frank Harrell
; CF10 3AX > Tel: 02920 875885 > E-mail: lelloj@.ac > [[alternative HTML version deleted]] > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-proj

Re: [R] SPSS -> R

2011-11-26 Thread Frank Harrell
> 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.nabb

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

2011-11-28 Thread Frank Harrell
that. > > Thank you. > > Kind regards > Maria > > __ > 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 prov

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

2011-11-29 Thread Frank Harrell
lity', > '10-Year Survival 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

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

2011-12-01 Thread Frank Harrell
bility', > '10-Year Survival 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) > > How can I show failure probability instead of survival probability? > > I would very much apprecia

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

2011-12-18 Thread Frank Harrell
*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 Pl

Re: [R] Stepwise in lme

2011-12-22 Thread Frank Harrell
ted a singularity which then stopped > lme from working > > Hmm, perhaps I should abandon such abominations ... 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 Harr

[R] Call for Abstracts for useR! 2012

2011-12-24 Thread 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 course offered on June 11. Participants are encouraged to submit a one-page abstract for oral or poster presentation at the conferenc

Re: [R] R 2.14.0 Design library

2012-01-04 Thread Frank Harrell
or some time, with (AFAIK) the > same functions. > > http://cran.r-project.org/web/packages/rms/index.html > > > Cheers > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide &

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

2012-01-09 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. > - Fra

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

2012-01-09 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. >> > > __ > R-help@ mailing lis

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

2012-01-09 Thread Frank Harrell
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

2012-01-10 Thread Frank Harrell
[in above command i am getting Error in Match(Y = Y, Tr = Tr, X = ps0) : > length(Y) != length(Tr)] > > Matchbalance(fserv_cd~ A4,match.out=rr1,nboots=1000,data=Passport) > > Please suggest me right approach for Propensity score matching using > classification trees(if possible code me) . >

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

2011-09-23 Thread Frank Harrell
> PLEASE do read the posting guide >> http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. >> > > > > -- > Joshua Wiley > Ph.D. Student, Health Psychology > Programmer Analyst II, ATS Statistical Con

Re: [R] Standardized beta-coefficients in regression

2011-09-26 Thread Frank Harrell
lm)$coefficients); > coef.table <- cbind(coef.table, std.coef); > print(coef.table); > > I do agree with B.R. but unfortunately the life of an applied statistician > is complex sometimes :-) > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this mes

Re: [R] model selection using logistf package

2011-09-27 Thread Frank Harrell
ration warning message so I tried using 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 > o

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 Depart

Re: [R] About stepwise regression problem

2011-10-07 Thread Frank Harrell
s="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? &g

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

2011-10-09 Thread Frank Harrell
version deleted]] > > ______ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproduc

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

2011-10-11 Thread Frank Harrell
glieder: Prof. Dr. Guido Sauter (Vertreter des Vorsitzenden), > Dr. Alexander Kirstein, Joachim Prölß, Prof. Dr. Dr. Uwe Koch-Gromus > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide

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

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

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] interpreting bootstrap corrected slope [rms package]

2011-10-22 Thread Frank Harrell
z.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:

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

2011-10-23 Thread Frank Harrell
r the middle Y > category have a wide range then keeping 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.78969

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

2011-10-23 Thread Frank Harrell
be included > 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 arc

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

2011-10-24 Thread Frank Harrell
dictors with which the > response variable 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.h

Re: [R] ACE/AVAS with bootstrap

2011-10-28 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/ACE-AVAS-with-bootstrap-tp3948613p3950068.html Sent from the R help mailing list archive at Nabble.com.

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

2011-10-31 Thread Frank Harrell
; > -- > Joshua Wiley > Ph.D. Student, Health Psychology > Programmer Analyst II, ATS Statistical Consulting Group > University of California, Los Angeles > https://joshuawiley.com/ > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help &g

Re: [R] How to interpret Spearman Correlation

2011-11-02 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. >> >> >> > >> > >> &g

<    1   2   3   4   5   >