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

2011-11-04 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/Select-some-but-not-all-variabl

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

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

2011-11-05 Thread Frank Harrell
> [[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-project.org/posting-guide.html > and provid

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

2011-11-06 Thread Frank Harrell
Rq() accepts a vector argument for tau. I >>> seem to remember that Koencker's rq() does.. Normally I would >>> consult the help page, but the power is still out here in Central >>> Connecticut and I am corresponding with a less capable device. I am >>>

Re: [R] Computing line= for mtext

2012-03-01 Thread Frank Harrell
his was provided for me > by David Winsemius in response to my query on this list in October 2011 > https://stat.ethz.ch/pipermail/r-help/2011-October/292806.html > > The email thread also includes comments by Deepayan Sarkar and Paul > Murrell. > > Rich > > On Wed, Fe

Re: [R] Computing line= for mtext

2012-03-02 Thread Frank Harrell
ank Richard M. Heiberger wrote > > Frank, > > this is it. It uses Peter's idea. > > plot(1:10) > axis(side=2, 1:10, las=1, line=-31.5, lwd=0) > axis(side=4, 1:10, las=1, labels=FALSE) > > Rich > > > On Thu, Mar 1, 2012 at 6:52 PM, Fran

Re: [R] Computing line= for mtext

2012-03-02 Thread Frank Harrell
help, Frank ilai-2 wrote > > On Fri, Mar 2, 2012 at 1:17 PM, Frank Harrell <f.harrell@> wrote: >> Hi Rich and Peter, >> >> What I am trying to do is the right-justify a vector of numbers to the >> right >> of the y-axis so that the leftmost digit of all of

Re: [R] Computing line= for mtext

2012-03-03 Thread Frank Harrell
text function instead of mtext (remembering to > set par(xpd=...)), then use the grconvertX and grconvertY functions to > find the location to plot at (possibly adding in the results from > strwidth or stheight). > > On Thu, Mar 1, 2012 at 4:52 PM, Frank Harrell <f.harrell@> wr

Re: [R] Computing line= for mtext

2012-03-03 Thread Frank Harrell
Whoops. I spoke too soon. That last approach does not survive shrinking the graphics window. I'll go back to the axis and par approach. Frank Frank Harrell wrote > > Thanks very much for the ideas Baptiste and Greg. I think this is a way > to go: > > # Right justifies

Re: [R] 632 estimator using boot package

2012-03-05 Thread Frank Harrell
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: htt

Re: [R] Nagelkerke R2

2012-03-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 Biostatisti

Re: [R] 632 estimator using boot package

2012-03-05 Thread Frank Harrell
aybe direclty using Boot package may not be a good idea for > assessing the overfitting in the regression. > > Jim > > > --- On Mon, 5/3/12, Frank Harrell <f.harrell@> wrote: > >> From: Frank Harrell <f.harrell@> >> Subject: Re: [R] 632 estimator usi

Re: [R] Summary values from Glm function (rms package)

2012-03-22 Thread Frank Harrell
original data install of using the gamma model). See the Hmisc function areg.boot and associated functions for example code. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Summary-values-from-Glm

Re: [R] survplot function

2012-03-27 Thread Frank Harrell
b: www.uk-sh.de > > __ > 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, reproducibl

[R] Deadlines Approaching: R User Conference Registration, Hotels, Student Scholarships

2012-04-04 Thread Frank Harrell
DEADLINE FAST APPROACHING – 8th Annual International R User Conference useR! 2012, Nashville, Tennessee USA Registration Deadlines: Early Registration: Passed Regular Registration: Mar 1- May 12 Late Registration: May 13 – June 4 On-Site Registration: June 12 – June 15 Please note: Nashville is

Re: [R] How do Sweave users collaborate with Word users?

2012-04-09 Thread Frank Harrell
gt;>> and provide commented, minimal, self-contained, reproducible code. >>>> >>>[[alternative HTML version deleted]] >>> >>> __________ >>> R-help@ mailing list >>> https://stat.ethz.ch/mai

Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS

2012-04-13 Thread Frank Harrell
 lty = c(1,2)) > > legend(0.75,0.19, >  legend = c("gD2","PBO"), >  lty = c(1,2), title = "Treatment") > > summary(surv.by.vac) > print(surv.by.vac) >   > > __ > R-help@ mailing list >

Re: [R] Kaplan Meier analysis: 95% CI wider in R than in SAS

2012-04-14 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] Kaplan Meier analysis: 95% CI wider in R than in SAS

2012-04-16 Thread Frank Harrell
where >> experts in the R community use this option in their book. The book is >> about >> 10 years old. So maybe opinion about the use of this option has shifted >> since then. >> > >> > Thanks, >> > >> > Paul >> > >> > [[alternative HTML version deleted]] >>

Re: [R] kruskal wallis post hoc?

2012-01-12 Thread Frank Harrell
g guide http://www.R-project.org/posting-guide.html >>>and provide commented, minimal, self-contained, reproducible code. >>> >>> >> >> >> > [[alternative HTML version deleted]] > > > __ > R

Re: [R] general question on Spotfire

2012-01-12 Thread Frank Harrell
ase, and I would have flipped back to straight R for > that anyway. > > Terry Therneau > > __ > R-help@ mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-g

Re: [R] kruskal wallis post hoc?

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

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

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

Re: [R] general question on Spotfire

2012-01-15 Thread Frank Harrell
e needed R functions in S-Plus ever happened. Frank Uwe Ligges-3 wrote > > On 12.01.2012 17:38, Frank Harrell wrote: >> As a slight aside, Tibco/Spotfire originally planned to provide a >> capability >> to load R packages into S-Plus. This always seemed to me to be a

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

2012-01-19 Thread Frank Harrell
ore > > > Any thoughts or suggestions? Thanks a lot! > > Jun Xu, PhD > Assistant Professor > Department of Sociology > Ball State University > > __ > R-help@ mailing list > https

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

2012-01-25 Thread Frank Harrell
l21.glm <- glm(cbind(Shells, TotalEggs-Shells) ~ HTL:Veg:TotalEggs, > data.to.analyze, family = binomial) > Model37.glm <- glm(cbind(Shells, TotalEggs-Shells) ~ > HTL:Veg:TotalEggs:Aeventexhumed, data.to.analyze, family=binomial) > > Please advise, thanks! > J > ---

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

2012-01-26 Thread Frank Harrell
lt;20, 20-25, 25-30 and >= 30. I am having > difficulty figuring the code for <20 and >=30? Please help. Thank you. > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/How-do-I-use-the-cut-function

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

2012-01-26 Thread Frank Harrell
not be pretty but works for me. > > >> An AIC difference equal or less than 2 is a tie, anything higher is >> evidence for the model with the less AIC (Sakamoto et al., Akaike >> Information Criterion Statistics, KTK Scientific Publishers, Tokyo). > > >

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

2012-01-27 Thread Frank Harrell
; My question here is - i have 124 variables(mix of categorical and > continuous) > how would i select(R code) attribute - ie. independent variable from those > 124 variable(Variable selection) so that i can build a model using glm > > pls. give me codes in R- on how to do that &g

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

2012-01-27 Thread Frank Harrell
odel > selection, Jim Lindsey being influent too. Wanting that your models say > some something about nature is another strong conditioning factor. Put > this two together and it becomes hard to do model-averaging. And it has > nothing to do with religion (yuck!). > >

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

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

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

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

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

2012-02-02 Thread Frank Harrell
-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.na

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

2012-02-06 Thread Frank Harrell
1. > july - june == 01.3182 0.7899 1.669 0.2870 > > I'm afraid "Tukey" it is not the correct contrast, but, What should i use > for a binary response variable? > Any advice, please? > - Frank Harrell Department of Biostatistics, Vande

Re: [R] Outlier removal techniques

2012-02-09 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 --

Re: [R] Matrix algebra in R to compute coefficients of a linear regression.

2012-02-18 Thread Frank Harrell
orrect results. > A better way to calculate the beta's is to calculate them directly without > explicitly calculating an invers > > # Better > # Compute betas, from t(X)X beta = t(Z) y > # without computing inverse explicitly > betaAlt <- solve(t(data) %*% data, num) > betaAlt > > Even better is th

Re: [R] Reporting Kaplan-Meier / Cox-Proportional Hazard Standard Error, km.coxph.plot, survfit.object

2012-02-20 Thread Frank Harrell
her publications. > > What is your approach? > > Thanks for your help, > PCJ > > __ > 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 p

Re: [R] how to How to compare the concordance index between two nomograms with R

2012-02-20 Thread Frank Harrell
by Cox > regression. My question is how to compare the concordance index between > these two nomograms with R? > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/how-to-How-to-compare-the-concordance-index-betwee

[R] useR! 2012 Early Bird Registration Ends Soon

2012-02-20 Thread Frank Harrell
ssions are welcomed. See http://biostat.mc.vanderbilt.edu/UseR-2012 for details. Frank ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/useR-2012-Early-Bird-Registration-Ends-Soon-tp4405954p4405954.html Sent from

Re: [R] stepwise selection for conditional logistic regression

2012-02-22 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.nabbl

Re: [R] How would you program an Adverse Events statistical table using R code?

2012-02-25 Thread Frank Harrell
s >> >> __ >> 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-contai

[R] Computing line= for mtext

2012-02-29 Thread Frank Harrell
'cin')[1] mtext(s, side=4, las=1, at=5, adj=1, line=w-.5, cex=1) mtext(s, side=4, las=1, at=7, adj=1, line=2*(w-.5), cex=2) Thanks Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Computing-line

Re: [R] Create a function "automatically" from lm formula and coefficients?

2012-02-29 Thread Frank Harrell
require(rms) f <- ols(y ~ rcs(age,4)*sex + blood.pressure + rcs(height,5)) latex(f) # typeset algebraic form of model Function(f) # create an R function that computes Y hat Frank ----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: h

Re: [R] Re : Create a function "automatically" from lm formula and coefficients?

2012-03-01 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. > -

Re: [R] Gls function in rms package

2012-04-18 Thread Frank Harrell
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/Gls-function-in-rms-package-tp4569521p4569829.html Sent fr

Re: [R] Splitting data into test and train (80:20) kepping attributes similar

2012-04-26 Thread Frank Harrell
for this to happen. Consider resampling instead. Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Splitting-data-into-test-and-train-80-20-kepping-attributes-similar-tp4583928p4589554.html Sent from the R

Re: [R] Gls function in rms package

2012-04-30 Thread Frank Harrell
x, data=d, correlation = corARMA(p=1)) #This works >>> >>> I would rather use Gls than gls so that I can represent a variable >>> with a spline using rcs. I'm using version 3-5.0 of rms in R 2.15.0. >>> >>> Thanks for any help you can give. >>

Re: [R] question on jitter in plot.Predict in rms

2012-05-01 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

[R] DEADLINE FAST APPROACHING for useR! 2012

2012-05-03 Thread Frank Harrell
DEADLINE FAST APPROACHING – 8th Annual International R User Conference useR! 2012, Nashville, Tennessee USA Registration Deadlines: Early Registration: Passed Regular Registration: Mar 1- May 12 Late Registration: May 13 – June 4 On-Site Registration: June 12 – June 15 Please note: Nashville is

[R] May 12 DEADLINE FAST APPROACHING for useR! 2012

2012-05-07 Thread Frank Harrell
8th Annual International R User Conference useR! 2012, Nashville, Tennessee USA ===> Regular registration rates and most blocks of discounted ===> hotel rooms are available only 5 more days Registration Deadlines: Early Registration: Passed Regular Registration: Mar 1- May 12 Late Registration:

[R] useR! 2012: Call for Late-Breaking Posters; REGULAR REGISTRATION ENDS 12May

2012-05-11 Thread Frank Harrell
*** Call for Late-breaking Posters *** Abstracts may be submitted for posters presenting recent developments and late-breaking applications of R, on topics as indicated in the earlier call for abstracts: http://biostat.mc.vanderbilt.edu/UseR-2012#Call_for_Abstracts_and_Tutorial Late-breaking

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

2012-05-25 Thread Frank Harrell
rms_3.5-0Hmisc_3.9-3 survival_2.36-14 > data.table_1.8.0 > > loaded via a namespace (and not attached): > [1] cluster_1.14.2 grid_2.15.0lattice_0.20-6 tools_2.15.0 > > __ > R-help@ mailing list > https://stat.ethz

<    1   2   3   4   5