Re: [R] FW: variable format

2007-09-07 Thread Frank E Harrell Jr
. __ R-help@stat.math.ethz.ch 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 E Harrell Jr Professor and Chair

Re: [R] Stratified weighted statistics

2007-08-29 Thread Frank E Harrell Jr
. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Excel

2007-08-29 Thread Frank E Harrell Jr
!!! Very well and very cogently expressed! cheers, Rolf Turner Yes! In addition I'd like to stress that the Excel model represents non-reproducible research. Frank Harrell -- Frank E Harrell Jr Professor and Chair School

Re: [R] validate (package Design): error message subscript out of bounds

2007-08-28 Thread Frank E Harrell Jr
will send a separate note with a new version of the function in question for you to source( ) to override the current function. Frank Best, Tore -Opprinnelig melding- Fra: Frank E Harrell Jr [mailto:[EMAIL PROTECTED] Sendt: 28. august 2007 03:17 Til: Wentzel-Larsen, Tore Kopi: r

Re: [R] how to include bar values in a barplot?

2007-08-27 Thread Frank E Harrell Jr
Donatas G. wrote: On Tuesday 07 August 2007 22:09:52 Donatas G. wrote: How do I include bar values in a barplot (or other R graphics, where this could be applicable)? To make sure I am clear I am attaching a barplot created with OpenOffice.org which has barplot values written on top of each

Re: [R] grouping scat1d/rug and plotting to 2 axes

2007-08-27 Thread Frank E Harrell Jr
the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] validate (package Design): error message subscript out of bounds

2007-08-27 Thread Frank E Harrell Jr
-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help

[R] Random number expert needed

2007-08-24 Thread Frank E Harrell Jr
with the agency. Thanks Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Ask for functions to obtain partial R-square (squared partial correlation coefficients)

2007-08-20 Thread Frank E Harrell Jr
. Xingwang Ye library(Design) # requires Hmisc f - ols(y ~ x1 + x2) p - plot(anova(f), what='partial R2') p The anova.Design function called above handles pooling related degrees of freedom and pools main effects with related interaction effects to get the total partial effect. Frank -- Frank E

Re: [R] Regulatory Compliance and Validation Issues

2007-08-20 Thread Frank E Harrell Jr
-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] Subject: Re: how to include bar values in a barplot?

2007-08-09 Thread Frank E Harrell Jr
, I disagree. Different elements in a display serve different purposes, according to the psychological aspects of visual preception. . . . -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Subject: Re: how to include bar values in a barplot?

2007-08-09 Thread Frank E Harrell Jr
. The Hmisc dotchart2 function makes this easy. Frank x - 1:5 names(x) - letters[1:5] bp - barplot(x) text(bp, x - .02 * diff(par(usr)[3:4]), x) On 8/9/07, Frank E Harrell Jr [EMAIL PROTECTED] wrote: [EMAIL PROTECTED] wrote: Greg, I'm going to join issue with your here! Not that I'll go

Re: [R] how to include bar values in a barplot?

2007-08-07 Thread Frank E Harrell Jr
in which the viewer adds some of the heights of the numbers to the heights of the bars. And in general dot charts are greatly preferred to bar plots, not the least because of their reduced ink:information ratio. -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] (Censboot, Z-score, Cox) How to use Z-score as the statistic within censboot?

2007-08-06 Thread Frank E Harrell Jr
guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] ROC curve in R

2007-08-04 Thread Frank E Harrell Jr
Dylan Beaudette wrote: On Thursday 26 July 2007 10:45, Frank E Harrell Jr wrote: Dylan Beaudette wrote: On Thursday 26 July 2007 06:01, Frank E Harrell Jr wrote: Note that even though the ROC curve as a whole is an interesting 'statistic' (its area is a linear translation of the Wilcoxon

Re: [R] proportional odds model

2007-08-02 Thread Frank E Harrell Jr
provides partial and score residual plots to check PO model assumptions. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch

Re: [R] logistic regression

2007-07-26 Thread Frank E Harrell Jr
commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https

Re: [R] ROC curve in R

2007-07-26 Thread Frank E Harrell Jr
through it. Thanks and Regards Rithesh M Mohan [[alternative HTML version deleted]] - Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] ROC curve in R

2007-07-26 Thread Frank E Harrell Jr
Dylan Beaudette wrote: On Thursday 26 July 2007 06:01, Frank E Harrell Jr wrote: Note that even though the ROC curve as a whole is an interesting 'statistic' (its area is a linear translation of the Wilcoxon-Mann-Whitney-Somers-Goodman-Kruskal rank correlation statistics), each individual

Re: [R] Constructing bar charts with standard error bars

2007-07-25 Thread Frank E Harrell Jr
commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Subscript out of bounds when using datadist() from Design library

2007-07-25 Thread Frank E Harrell Jr
://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 E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] Off-topic: simulate time dependent covariates

2007-07-22 Thread Frank E Harrell Jr
ratio treatment effect. The user specifies the hazard ratio as a function of time. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help

Re: [R] exclude1 in summary.formula from Hmisc

2007-07-20 Thread Frank E Harrell Jr
]] __ R-help@stat.math.ethz.ch 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 E Harrell Jr

Re: [R] exclude1 in summary.formula from Hmisc

2007-07-20 Thread Frank E Harrell Jr
are small are a bit overdone. Better might be to use the likelihood ratio tests, and many other tests are available. Frank Regards. David -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Package for .632 (and .632+) bootstrap and the cross-validation of ROC Parameters

2007-07-13 Thread Frank E Harrell Jr
, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0), CAT01 = c(4, 4, 2, 4, 4, 4, 4, 4, 4, 2, 1, 2, 2, 4, 4, 4, 1, 1, 2, 2, 1, 4, 1, 4, 1, 4, 2, 4, 1, 4, 2, 3, 1, 1, 3, 3, 4, 2, 4, 2, 1, 2, 2, 1, 1, please reply... Frank E Harrell Jr wrote: spime wrote: Hi

Re: [R] Package for .632 (and .632+) bootstrap and the cross-validation of ROC Parameters

2007-07-12 Thread Frank E Harrell Jr
seem to have advantages mainly if an improper, inefficient, discontinuous scoring rule such as the percent classified correct is used. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] variable and value labels

2007-07-11 Thread Frank E Harrell Jr
. Grateful for any help, Steve Powell In the Hmisc package see the functions spss.get, label, upData, and describe. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Incidence estimated from Kaplan-Meier

2007-07-05 Thread Frank E Harrell Jr
commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Formula syntax question

2007-07-03 Thread Frank E Harrell Jr
, e.g., form - as.formula(paste('y', paste(names(mydata)[continuous], collapse='+'), sep='~')) lrm(form, data=mydata) -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Harrell's C

2007-07-03 Thread Frank E Harrell Jr
: 917277112 [EMAIL PROTECTED] It's in the documention with the Hmisc package. Type ?rcorr.cens ?rcorrp.cens Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] logistic regression and dummy variable coding

2007-06-29 Thread Frank E Harrell Jr
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 E Harrell Jr Professor and Chair School of Medicine Department

Re: [R] Warning message in Design ...

2007-06-29 Thread Frank E Harrell Jr
@stat.math.ethz.ch 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 E Harrell Jr Professor and Chair School of Medicine

Re: [R] Marginal Effects of continuous variable in lrm model (Design package)

2007-06-26 Thread Frank E Harrell Jr
of the newly added variables) it doesn't work. I don't know what went wrong. Thank you very much. Thanks, Minyu Chen Please read the package's help files especially about the datadist function. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] help on the use of ldBand

2007-06-22 Thread Frank E Harrell Jr
/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 E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] BIC and Hosmer-Lemeshow statistic for logistic regression

2007-06-19 Thread Frank E Harrell Jr
spime wrote: I haven't find any helpful thread. How can i calculate BIC and Hosmer-Lemeshow statistic for a logistic regression model. I have used glm for logistic fit. See the Design package's lrm function and residuals.lrm for a better GOF test. -- Frank E Harrell Jr Professor

Re: [R] BIC and Hosmer-Lemeshow statistic for logistic regression

2007-06-19 Thread Frank E Harrell Jr
spime wrote: Is there any windows version of Design package??? Soon the new version will will make its way to Windows, probably in a day or two. Frank Frank E Harrell Jr wrote: spime wrote: I haven't find any helpful thread. How can i calculate BIC and Hosmer-Lemeshow

Re: [R] psm/survreg coefficient values ?

2007-06-18 Thread Frank E Harrell Jr
0.03789 -1.43 1.51e-01 Log(scale)-2.72870.03780 -72.18 0.00e+00 Scale= 0.0653 best, Spencer I have a case study using psm (survreg wrapper) in my book. Briefly, coefficients are on the log median survival time scale. Frank -- Frank E Harrell Jr Professor and Chair

Re: [R] error bars on survival curve

2007-06-17 Thread Frank E Harrell Jr
of R, you can have more control using the survplot function, with bars and bands options. Do ?survplot.survfit Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] selecting cut-off in Logistic regression using ROCR package

2007-06-16 Thread Frank E Harrell Jr
of categorization;PBC data} } Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] complex contrasts and logistic regression

2007-06-16 Thread Frank E Harrell Jr
?predict ?predict.lrm ?predict.Design. Also do ?gendata which will generate a data frame for getting predictors, with unspecified predictors set to reference values such as medians. Frank Nicholas -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] Design library installation problem

2007-06-13 Thread Frank E Harrell Jr
will be made available shortly. Please note that this would be a lot easier if we did not get a segfault when running R CMD CHECK --- a segfault that does not happen with the example code is run outside of CMD CHECK. Charles Thomas Dupont is working hard on this. Frank -- Frank E Harrell Jr

Re: [R] R vs. Splus in Pharma/Devices Industry

2007-06-11 Thread Frank E Harrell Jr
__ R-help@stat.math.ethz.ch 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 E

Re: [R] R is not a validated software package..

2007-06-08 Thread Frank E Harrell Jr
that such compliance is actually needed. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] R is not a validated software package..

2007-06-08 Thread Frank E Harrell Jr
-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt

Re: [R] R is not a validated software package..

2007-06-08 Thread Frank E Harrell Jr
South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Frank E Harrell Jr Sent: Friday, June 08, 2007 7:45 AM To: Giovanni Parrinello Cc: r-help@stat.math.ethz.ch Subject: Re: [R] R is not a validated software

Re: [R] Tools For Preparing Data For Analysis

2007-06-08 Thread Frank E Harrell Jr
Dale Steele wrote: For windows users, EpiData Entry http://www.epidata.dk/ is an excellent (free) tool for data entry and documentation.--Dale Note that EpiData seems to work well under linux using wine. Frank __ R-help@stat.math.ethz.ch mailing

Re: [R] Tools For Preparing Data For Analysis

2007-06-07 Thread Frank E Harrell Jr
. __ R-help@stat.math.ethz.ch 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 E

Re: [R] Harrell's C

2007-06-02 Thread Frank E Harrell Jr
of Hmisc not loading under Windows but Design is not currently available for Windows for R 2.5. We are working to fix that. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Problem with Weighted Variance in Hmisc

2007-06-01 Thread Frank E Harrell Jr
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 E Harrell Jr Professor and Chair School of Medicine

Re: [R] R's Spearman

2007-05-31 Thread Frank E Harrell Jr
, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] AIC for lrm(Hmisc/Design) model.

2007-05-29 Thread Frank E Harrell Jr
. For that you can do aic - function(fit) round(unname(fit$stats['Model L.R.'] - 2*fit$stats['d.f.']),2) f - lrm( ) aic(f) If unname doesn't exist in S-Plus as it does in R, you can remove that part. -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] normality tests

2007-05-25 Thread Frank E Harrell Jr
not use distribution-free methods? You should examine the power of the tests for n=20. You'll probably find it's not good enough to reach a reliable conclusion. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
Of Liaw, Andy Sent: Friday, May 25, 2007 12:04 PM To: [EMAIL PROTECTED]; Frank E Harrell Jr Cc: r-help Subject: Re: [R] normality tests [Broadcast] From: [EMAIL PROTECTED] On 25/05/07, Frank E Harrell Jr [EMAIL PROTECTED] wrote: [EMAIL PROTECTED] wrote: Hi all, apologies for seeking advice

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
in fairly quickly. Testing for normality prior to choosing a test statistic is generally not a good idea. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Liaw, Andy Sent: Friday, May 25, 2007 12:04 PM To: [EMAIL PROTECTED]; Frank E Harrell Jr Cc: r-help

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
Frank E Harrell Jr [EMAIL PROTECTED] To bilt.edu

Re: [R] clustered standarderrors using design package

2007-05-19 Thread Frank E Harrell Jr
-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch

Re: [R] How to analyse simple study: Placebo-controlled (2 groups) repeated measurements (ANOVA, ANCOA???)

2007-05-17 Thread Frank E Harrell Jr
will be appreciate! Thanks! Karl -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] Re : Bootstrap sampling for repeated measures

2007-05-15 Thread Frank E Harrell Jr
]] __ R-help@stat.math.ethz.ch 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 E Harrell Jr Professor

Re: [R] Nicely formatted summary table with mean, standard deviation or number and proportion

2007-05-13 Thread Frank E Harrell Jr
, stat.name=c(n, mean, sd))) treatment n mean sd 1 Drug 263 49.9 4.94 2 Placebo 237 50.1 4.97 # this works but table is rotated, and it count data has to be # treated separately. -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] Read SAS data into R

2007-05-11 Thread Frank E Harrell Jr
, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Follow-up about ordinal logit with mixtures: how about 'continuation ratio' strategy?

2007-05-10 Thread Frank E Harrell Jr
Paul Johnson wrote: This is a follow up to the message I posted 3 days ago about how to estimate mixed ordinal logit models. I hope you don't mind that I am just pasting in the code and comments from an R file for your feedback. Actual estimates are at the end of the post. . . . Paul, lrm

[R] Mantel-Haenszel relative risk with Greenland-Robins variance estimate

2007-05-08 Thread Frank E Harrell Jr
Does anyone know of an R function for computing the Greenland-Robins variance for Mantel-Haenszel relative risks? Thanks Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Mantel-Haenszel relative risk with Greenland-Robins variance estimate

2007-05-08 Thread Frank E Harrell Jr
factor other than occupational group. Thanks again Frank Frank E Harrell Jr

Re: [R] pseudo-R2 or GOF for regression trees?

2007-05-05 Thread Frank E Harrell Jr
;nonparametric regression;CART;simulation setup} } -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Error in if (!length(fname) || !any(fname == zname)) { :

2007-05-04 Thread Frank E Harrell Jr
(v1)+v2+v3+v4, data=h) -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Hmisc curve label size cex

2007-04-28 Thread Frank E Harrell Jr
( One= list( x,sin(x)), Two= list( x,cos(x)), Three=list( x,(x*x)), Four= list( x,exp(x)) ), pl=TRUE ) and got some nice output after reducing cex.* Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] Integrating R-programs into larger systems

2007-04-27 Thread Frank E Harrell Jr
but it will publish the model to a web server. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Not showing dvi with Hmisc latex()

2007-04-26 Thread Frank E Harrell Jr
Duncan Murdoch wrote: On 4/26/2007 9:20 PM, Gad Abraham wrote: Hi, I'm using latex() from Frank Harrell's Hmisc library to produce LaTeX files. By default, it calls xdvi and displays the dvi. How can I make xdvi not show? I couldn't find a clue in the extensive documentation. Unclass

Re: [R] Coercing data types for use in model.frame

2007-04-25 Thread Frank E Harrell Jr
. After model.frame is finished I'll change such variables to factors or matrices. Frank On Tue, 24 Apr 2007, Frank E Harrell Jr wrote: In the Hmisc package there is a new data class 'mChoice' for multiple choice variables. There are format and as.numeric methods (the latter creates

[R] Biostatistician Opportunities at Vanderbilt

2007-04-25 Thread Frank E Harrell Jr
The Department of Biostatistics at Vanderbilt University's School of Medicine has openings for biostatisticians at all levels. Details and application procedures may be found at http://biostat.mc.vanderbilt.edu/JobOpenings . For M.S. and B.S. biostatisticians we are especially interested in

[R] Coercing data types for use in model.frame

2007-04-24 Thread Frank E Harrell Jr
that model.frame will use automatically? I would use as.factor here. model.frame does not seem to use as.data.frame.foo for individual variables. Thanks Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] help comparing two median with R

2007-04-18 Thread Frank E Harrell Jr
Thomas Lumley wrote: On Tue, 17 Apr 2007, Frank E Harrell Jr wrote: The points that Thomas and Brian have made are certainly correct, if one is truly interested in testing for differences in medians or means. But the Wilcoxon test provides a valid test of x y more generally. The test

Re: [R] help comparing two median with R

2007-04-18 Thread Frank E Harrell Jr
Frank E Harrell Jr [EMAIL PROTECTED] To bilt.edu Thomas Lumley

Re: [R] help comparing two median with R

2007-04-17 Thread Frank E Harrell Jr
with the Hodges-Lehmann estimator: the median of all possible differences between an X and a Y. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-17 Thread Frank E Harrell Jr
://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 E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] Nonparametric Effect size indices

2007-04-13 Thread Frank E Harrell Jr
://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 E Harrell Jr Professor and Chair School of Medicine Department

Re: [R] sas.get problem

2007-04-11 Thread Frank E Harrell Jr
-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch

Re: [R] Reasons to Use R

2007-04-10 Thread Frank E Harrell Jr
save( ) and load( ). Frank --Todd -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-06 Thread Frank E Harrell Jr
-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Doing partial-f test for stepwise regression

2007-04-01 Thread Frank E Harrell Jr
and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing

Re: [R] what is the difference between survival analysis and (...)

2007-03-29 Thread Frank E Harrell Jr
Thomas Lumley wrote: On Wed, 28 Mar 2007, Frank E Harrell Jr wrote: Eric Elguero wrote: Hi everybody, recently I had to teach a course on Cox model, of which I am not a specialist, to an audience of medical epidemiologists. Not a good idea you might say.. anyway, someone in the audience

Re: [R] what is the difference between survival analysis and (...)

2007-03-28 Thread Frank E Harrell Jr
://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] weight factor in somers2 function

2007-03-28 Thread Frank E Harrell Jr
could be exactly included as weight factor? Thank you very much! Send me the smallest artificial example you can construct and I'll work on it. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt

Re: [R] how to get lsmeans?

2007-03-22 Thread Frank E Harrell Jr
to easy-to-explain predicted values. I would appreciate your thoughts. Thanks Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help

Re: [R] how to get lsmeans?

2007-03-22 Thread Frank E Harrell Jr
-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: Frank E Harrell Jr [mailto:[EMAIL PROTECTED] Sent: Thursday, March 22, 2007 8:41 AM To: John Fox Cc: 'Prof Brian Ripley'; 'r-help'; 'Chuck Cleland'; 'Liaw, Andy' Subject: Re: [R

Re: [R] how to get lsmeans?

2007-03-21 Thread Frank E Harrell Jr
, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] Stepwise Logistic Regression

2007-03-21 Thread Frank E Harrell Jr
almost the same result with the following code :-) candidates - setdiff(names(mydataframe), 'y') p - length(candidates) sample(candidates, sample(round(p/4):round(p/2), 1)) -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] R and clinical studies

2007-03-20 Thread Frank E Harrell Jr
__ R-help@stat.math.ethz.ch 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 E Harrell Jr Professor and Chair

Re: [R] R and clinical studies

2007-03-16 Thread Frank E Harrell Jr
analyses [and validating analysis file derivation]. Pivotal analyses can be re-done a variety of ways, in R or in separate programmable packages such as Stata. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt

Re: [R] ols Error : missing value where TRUE/FALSE needed

2007-03-14 Thread Frank E Harrell Jr
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] R and SAS proc format

2007-03-07 Thread Frank E Harrell Jr
that have PROC FORMAT value labels, an attribute 'sas.codes' keeps the original codes; these can be retrieved using sas.codes(variable name). This could be done outside the SAS import context also. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] R and SAS proc format

2007-03-06 Thread Frank E Harrell Jr
of times someone has sent me a SAS dataset and forgotten to send the PROC FORMAT value labels]. See the factor function for information about how R does this. Frank Harrell -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

Re: [R] R and SAS proc format

2007-03-06 Thread Frank E Harrell Jr
://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 E Harrell Jr Professor and Chair School of Medicine Department

Re: [R] bootstrap

2007-02-28 Thread Frank E Harrell Jr
/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 E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt

Re: [R] fitting of all possible models

2007-02-27 Thread Frank E Harrell Jr
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 E Harrell Jr Professor and Chair School of Medicine

Re: [R] rpart minimum sample size

2007-02-27 Thread Frank E Harrell Jr
, when relying on a single tree (hence the popularity of random forests, bagging, boosting). Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University __ R

Re: [R] fitting of all possible models

2007-02-27 Thread Frank E Harrell Jr
Indermaur Lukas wrote: Hi Frank I fitted a set of 12 candidate models and evaluated the importance of variables based on model averaged coefficients and SE (model weights =0.9). Variables in my models were not distributed in equal numbers across all models thus I was not able to evaluate

Re: [R] fitting of all possible models

2007-02-27 Thread Frank E Harrell Jr
Bert Gunter wrote: ... Below -- Bert Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Frank E Harrell Jr Sent: Tuesday, February 27, 2007 5:14 AM

Re: [R] bootcov and cph error

2007-02-27 Thread Frank E Harrell Jr
_ Dr. Scott Williams MD Peter MacCallum Cancer Centre Melbourne, Australia [EMAIL PROTECTED] -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University

Re: [R] Double-banger function names: preferences and suggestions

2007-02-25 Thread Frank E Harrell Jr
best) Any suggestions? Thanks, Hadley Personal preference is for 3. I wish I had used that throughout the Hmisc and Design packages for non-generics. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics

  1   2   3   4   5   6   7   >