*RMS Short Course 2018*
Frank E. Harrell, Jr., Ph.D., Professor
Department of Biostatistics, Vanderbilt University School of Medicine
fharrell.com @f2harrell
*May 15-18, 2018* With Optional R Workshop May 14
9:00am - 4:00pm
Alumni Hall
Vanderbilt University
Nashville Tennessee USA
See http://
Fixed 'maxiter' in the help file. Thanks.
Please give the original source of that dataset.
That dataset is a tiny sample of GUSTO-I and not large enough to fit this
model very reliably.
A nomogram using the full dataset (not publicly available to my knowledge)
is already available in http://bio
In recent versions of rms on CRAN there was a non-downward compatible
change. To get latex output for summary, anova, and print on fit objects
you leave off file="" (because we're usually using knitr anyway) and use
options(prType='latex')
anova(f) # LaTeX output
You can use options(prType='
A major new version of the rms package is now on CRAN. The most
user-visible changes are:
- interactive plotly graphic methods for model fits. The best example of
this is survplot for npsurv (Kaplan-Meier) estimates where the number of
risk pop up as you hover over the curves, and you can click
Hmisc 4.0-0 is now on CRAN. The package has undergone a massive update.
The most user-visable changes are;
- support for Rmarkdown html notebooks
- advanced html tables using the htmlTable package and summaryM function;
can copy and paste into word processors
- support for plotly interactive gr
This happens when you have not strat variables in the model.
--
Frank E Harrell Jr Professor and Chairman School of Medicine
Department of *Biostatistics* *Vanderbilt University*
On Thu, Jun 2, 2016 at 10:55 AM, Steve Lianoglou
wrote:
> Hello foks,
>
I use this function to update my installed R packages:
updatePac <- function (checkBuilt = FALSE, ...)
update.packages(repos = "http://cran.rstudio.com";, instlib =
"/usr/local/lib/R/site-library",
checkBuilt = checkBuilt, ...)
When I type updatePac() I get:
Warning: package 'MASS' in libr
This does seem to be a good situation for ordinal regression. The R rms
package's orm function allows for thousands of categories in Y. But it
doesn't handle censoring.
This discussion would be better for stats.stackexchange.com
Frank
--
---
Rolf I believe \textsf is the correct font to use for the symbol R but
not necessarily for the names of R variables and functions. I'd like
more discussion about that.
Frank
--
Frank E Harrell Jr Professor and Chair
Please note that Kirsten is cross-posting to stats.stackexchange.com
creating extra work for everyone.
--
Frank E Harrell Jr Professor and Chairman School of Medicine
Department of *Biostatistics* *Vanderbi
New versions of Hmisc and rms are available on CRAN. Changes are listed
below.
The most significant change to Hmisc is the addition of the ffCompress
function that creates an optimal ff package object for large data frames
by computing the maximum number of bits used by each numeric or logical
It's implemented in the R rms package's lrm and orm functions.
-----
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/Nagelkerke-Pseudo-R-squared-tp4710014p4710031.html
Sent from the R help mailing list
Thank you very much Terry. I'm still puzzled at why this worked a year
ago. What changed? I'd very much like to reverse the change by setting
an argument somewhere or manipulating the terms object.
I echo your sentiments about the general approach.
Frank
On 06/15/2015 09:05 AM, Therneau, T
Terry - your example didn't demonstrate the problem because the variable
that interacted with strata (zed) was not a factor variable.
But I had stated the problem incorrectly. It's not that there are too
many strata terms; there are too many non-strata terms when the variable
interacting with
For building design matrices for Cox proportional hazards models in the
cph function in the rms package I have always used this construct:
Terms <- terms(formula, specials=c("strat", "cluster", "strata"), data=data)
specials <- attr(Terms, 'specials')
stra<- specials$strat
Terms.ns <- Ter
Several updates have been made to the Hmisc package:
Changes in version 3.16-0 (2015-04-25)
* html.contents.data.frame: corrected html space character to add
semicolon
* ggplot.summaryP: added size of points according to denominators
* colorFacet: new function
* labelPlotmath: add
*RMS Short Course 2015*
Frank E. Harrell, Jr., Ph.D., Professor and Chair
Department of Biostatistics, Vanderbilt University School of Medicine
*March 3, 4, 5 & 6, 2015* With Optional R Workshop March 2
9:00am - 4:00pm
Student Life Center Board of Trust Room
Vanderbilt University
Nashville Tenness
Subject: Regression Modeling Strategies 4-Day Short Course March 2015
*RMS Short Course 2015*
Frank E. Harrell, Jr., Ph.D., Professor and Chair
Department of Biostatistics, Vanderbilt University School of Medicine
*March 3, 4, 5 & 6, 2015* With Optional R Workshop March 2
9:00am - 4:00pm
Student
classes created by survfit.formula
* logLik.Gls: added. Makes AIC(Gls object) work.
* NAMESPACE: several changes
See http://biostat.mc.vanderbilt.edu/rms and
https://github.com/harrelfe/rms for more information.
Frank Harrell
Vanderbilt University
handle Inf. Thanks: Benjamin Tyner
* DESCRIPTION, NAMESPACE: multiple function changes to work in R-devel
See http://biostat.mc.vanderbilt.edu/Hmisc,
https://github.com/harrelfe/Hmisc for more information.
Frank Harrell
Vanderbilt University
___
Greg I just re-copied the latest subplot and its help file from
TeachingDemos to Hmisc for the next release. Thanks for pointing this out.
Frank
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the post
When using cph in the rms package there is a function Mean that operates
on cph objects to produce an R function for computing the mean or
restricted mean life time.
Frank
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
Please try hard to create a simple reproducible example, then I'll work
on it.
[[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.
I can't think of an example where R does not work better than SAS except
for a few cases of mixed effects regression models and for processing
enormous datasets when the R user does not want to learn about the
latest R tools for large datasets. I quit using SAS in 1991 (in favor
of S-Plus and
The Hmisc package has had a number of updates/fixes:
Changes in version 3.14-4 (2014-04-13)
* rcspline.eval: stop new logic for ignoring outer values when there
are many ties when there are also many ties on interior values. Added
new logic to use interior unique values of x when the number
The first submission of the new greport package is now on CTAN. This
package facilitates graphical reporting of clinical trials an is highly
tied to LaTeX, Hmisc, and lattice. It creates hyperlinks and pop-up
tooltips in the resulting pdf report. Tables are greatly de-emphasized,
but hyperli
Changes are listed below.
NOTE: The next release of rms will replace the survfit.formula function
with a new function named npsurv ("nonparametric survival estimates") so
as to not conflict with the survival package. This is a
NON-DOWNWARD-COMPATIBLE change. To get Kaplan-Meier estimates in
Recent changes include the following.
Changes in version 3.14-2 (2014-02-26)
* latex.default: improved logic using new function in Misc: latexBuild
* latex.default: fixed bug with ctable=TRUE with no caption by
removing default label
* latex.default: improved formatting for insert.top
My yearly Regression Modeling Strategies course is expanded to 4 days
this year to be able relax the pace a bit. Details are below.
Questions welcomed.
-
*RMS Short Course 2014*
Frank E. Harrell, Jr., Ph.D., Professor and Chair
ith no
> result so far.
> Is there a way that i can see the progress of the read?
> Or is there another way to read the file with less computing time?
> I do not have access to SAS, the file was sent to me.
>
> Let me know what you guys think
> KR
> Hans
-
Frank H
Hmisc has had the following updates, with some of the new functions
demonstrated on http://biostat.mc.vanderbilt.edu/HmiscNew . Much of the
new work is motivated by the need to replace tables with graphics.
Changes in version 3.1
My yearly Regression Modeling Strategies course is expanded to 4 days
this year to be able relax the pace a bit. Details are below.
Questions welcomed.
-
*RMS Short Course 2014*
Frank E. Harrell, Jr., Ph.D., Professor and Chai
This is an unstable process. I suggest using the bootstrap to get a
confidence interval for the rank of each correlation coefficient among all
non-diagonal correlations.
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4
Thanks very much Rich and Duncan. latticeExtra's resizePanels function
was a perfect solution.
Frank
__
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-g
Is it possible to vary the size of the panels in lattice within one
page? The examine I have in mind is a 3 row by 2 column display where I
vary the y-axis scales in a dot plot and there are only a few levels on
the y-axis for one of the rows. I'd like to remove wasted space in that
row.
On
Thank you Bill, that worked perfectly.
Frank
__
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
I would like to do this:
f <- function(formula, data=NULL) {
gg <- sqrt
model.frame(formula, data=data)
}
x <- y <- 1:10
f(y ~ gg(x))
Error in eval(expr, envir, enclos) : could not find function "gg"
Is there a simple way to get access to gg from within the model.frame
invocation inside f?
The rms package has had several updates in version 4.1-0:
* Fixed orm.fit to not create penalty matrix if not needed
(penalties are not yet implemented anyway)
* Added yscale argument to plot.Predict
* Added Wald test simulation to orm help file
* Added example in help file for plot
A significant update to the Hmisc package is now available on CRAN for
all platforms. Hmisc source is now on github at
https://github.com/harrelfe/Hmisc and the full change log may be found
at https://github.com/harrelfe/Hmisc/commits/master
The most important updates are additions of new gra
You did not follow the posting guide, did you use pure ascii email, and
used illegal characters in your source code. This caused extra work.
Once I cleaned up your characters and made the example self-contained,
the labels worked fine for me. Here's the cleaned-up code:
library(rms)
x1 <- ru
Not certain about .por but this works with ordinary SPSS files:
require(Hmisc)
dat <- spss.get(...) # gets variable labels, etc.
contents(dat)
html(contents(dat), ...)
The last command produces a hyperlinked data dictionary, e.g., for each
variable the number of levels is given and you click o
This is dichotomania and is an inappropriate use of continuous
variables. Use an information-preserving approach such as rank
correlation. Also, this is not an R question.
Frank
---
Assume that I want to compare two methods, say skinfold measurement and
BMI, in how they classify subj
; cIndexCoxglmnet <- c(cIndexCoxglmnet,
>> 1-rcorr.cens(predict(glmnet.cox.meaningful,
>> reformat_testSet),Surv(reformat_testSet$time,reformat_testSet
>> $status))[1])
>>}
>> }
>>
>> #Get average C-Index
>> cIndex<- mean (unlist
require(rms)
?orm# ordinal regression model
For a case study see Handouts in
http://biostat.mc.vanderbilt.edu/CourseBios330
Since you have lost the original values, one part of the case study will
not apply: the use of Mean().
Frank
-
I have a dataset which has sever
Bill I found a workaround:
f <- ff(formula, lab)
f <- as.formula(gsub("`", "", as.character(deparse(f
Thanks for your elegant solution.
Frank
--
Thanks Bill. The problem is one of the results of convertName might be
'Heading("Age in Years")*age' (this is fo
uot;/","%in%","(",
":"))) {
for(i in seq_along(expr)[-1]) {
expr[[i]] <- Recall(expr[[i]], convertName = convertName)
}
} else if (is.name(expr)) {
expr <- as.name(convertName(expr))
}
expr
}
Bill Dunlap
Spotfi
quot;Female") * SBPz) *
Heading() * Gz + (AGEz + SBPz) * Heading() * TRIOz ~ Heading() *
COUNTRYz * Heading() * SEXz
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -Original Message-
> From: [hidden email] [mailto:[hidden email]] On Behalf
> Of William Dunlap
&
I really appreciate the excellent ideas from Bill Dunlap and Greg Snow.
Both suggestions almost work perfectly. Greg's recognizes expressions
such as sex=='female' but not ones such as age > 21, age < 21, a - b >
0, and possibly other legal R expressions. Bill's idea is similar to
what Dunca
I would like to be able to use gsub or gsubfn to process a formula and
to translate the variables but to ignore expressions in the formula.
Supposing that the R formula has already been transformed into a
character string and that the transformation is to convert variable
names to upper case an
ables (here on the x-axis) where the labels tell all. The last 2
lines in the code uses superpositioning for group but the col argument
affects the color of confidence bands; I didn't look into that further.
Frank Harrell
--
require(rms)
n = 30
group = factor(sample(c('a
Duncan,
I had read your excellent tables package vignette at
http://cran.r-project.org/web/packages/tables/vignettes/tables.pdf when
it first came out. It is extremely impressive. I'm glad to be reminded
to give it another look.
Is there a way to make the special symbols n and 1 refer to t
Our Hmisc package summary.formula function and its latex methods can
make some fairly advanced tables. But the tables have to be regular.
For example, all rows of the tables are based on the same data frame.
I'm thinking that what is needed is a ggplot2-like set of functions for
building a tab
[Sorry I can't quote past messages as I get mail on Nabble and have been
told I can't reply through Nabble].
Thanks Kennel for recommended a narrowing range for usr y-limits. That
does help quite a bit.
But I found a disappointing aspect of the graphics system: When you
change height= on th
Thanks Rich and Jim and apologies for omitting the line
x <- c(285, 43.75, 94, 150, 214, 375, 270, 350, 41.5, 210, 30, 37.6,
281, 101, 210)
But the fundamental problem remains that vertical spacing is not correct
unless I waste a lot of image space at the top.
Frank
--
Frank E Harrell Jr Pr
I can't get the points and "a" and "b" to render correctly (they are
squeezed into the axis) unless I make the height of the figure waste a
lot of space. I'd appreciate any ideas. The code is below. Thanks!
png('/tmp/z.png', width=480, height=100)
par(mar=c(3,.5,1,.5))
plot.new()
par(usr=c(-1
The rms ("Regression Modeling Strategies") package has undergone a
massive update. The entire list of updates is at the bottom of this
note. CRAN has the update for linux and will soon have it for Windows
and Mac - check http://cran.r-project.org/web/packages/rms/ for
availability. This rms
Thank you Marc
Frank
Marc Schwartz-3 wrote
> On Jul 10, 2013, at 1:29 PM, Frank Harrell <
> f.harrell@
> > wrote:
>
>> I have been confused about the appropriate e-mail address to use to make
>> announcements to r-help for major package update. In the past I
I have been confused about the appropriate e-mail address to use to make
announcements to r-help for major package update. In the past I've
submitted to r-packa...@lists.r-project.org without seeing the
announcement appear on r-help.
Thanks for any guidance.
Frank
--
Frank E Harrell Jr Profe
Mike,
Do something like:
require(rms)
dd <- datadist(mydatarame); options(datadist='dd')
f <- Rq(y ~ rcs(age,4)*sex, tau=.5) # use rq function in quantreg
summary(f) # inter-quartile-range differences in medians of y (b/c tau=.5)
plot(Predict(f, age, sex)) # show age effect on median as a co
You are using an informal shrinkage method. It is much better to use
penalized maximum likelihood estimation, built in to lrm.
If you really want to go the informal route, compute the linear
predictor from your final estimates and use ols( ) to predict that from
the component variables (you'l
Unless you have detailed simulations to back up the performance of this
method I would avoid it. It violates several statistical principles.
Frank
Hari wrote
> Hello R geeks,
>
> Waiting for an reply.
>
> Thanks,
> Hari
-
Frank Harrell
Department of Biostat
ase the summation portion 1/ (1-r)
> would change). Is there an easy way to do this?
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/rmean-in-survfit-tp4667668p4667708.html
Sent from
___
> 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
7;t know how to make a
> loop. Could you please help me with this or suggest an R package ? I would
> be very very grateful for help with this task !
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/how-to
f Environmental and Forest Sciences,
>
> University of Washington
>
> [[alternative HTML version deleted]]
>
> __
> R-help@
> mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the post
,
because its assumed discontinuities seldom exist in nature and most
relationships are not piecewise flat.
Frank
levanovd wrote
> Or even simpler (no need to specify labels):
>
> x<-runif(100,0,100)
> u <- cut(x, breaks = c(0, 3, 4.5, 6, 8, Inf), labels = FALSE)
-
Frank Har
tinfo/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
/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.
>
> __
> R-help@
> mailing list
> htt
Now I see it. S3method() wants two arguments so I need to create multiple
S3method() statements for each generic.
Frank
Frank Harrell wrote
> Right, I should have said import(Hmisc) instead of importFrom(Hmisc), but
> that does not explain the error message.
> Blaser Nello wrote
>>
-exts.pdf), importFrom
> needs to know what functions you are importing [e.g. importFrom(Hmisc,
> "latex") importFrom(stats, "anova")].
>
>
> -Original Message-
> From:
> r-help-bounces@
> [mailto:
> r-help-bounces@
> ] On Behalf Of Frank Har
quot;bad 'S3method' directive: ols, pphsm,
psm, rms, Rq, summary.rms, validate)")
ERROR: lazy loading failed for package ‘rms’
Any advice appreciated.
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789
What a terrific article. Thanks for sharing! The more we critically
examine how research is actually done the more frightened we become.
Frank
--
Frank E Harrell Jr Professor and Chairman School of Medicine
Department of Biostatistics Vanderbilt University
___
Does anyone know of R functions for doing composite quantile regression (Hou
and Yuan Ann Stat 36:1108, 2008)? The paper's authors do not talk about
software in their paper or on their web sites.
Thanks
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View
as it possible with significant
> t-values?
>
> cheers, Oleg
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/Model-Selection-based-on-individual-t-values-with-the-largest-possible-number-of-variables-in-reg
Public Health
> University College London
>
>
>
> __
> 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, minima
. Thanks in advance!
>
> Regards!!
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/Accuracy-some-classifiers-tp4661109p4661181.html
Sent from the R help mailing list archive at Nabble.com.
___
> 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 o
gt; sem.cnes <- sem(model.cnes, data=CNES)
> summary(sem.cnes)
>
> set.seed(12345) # for reproducibility
> system.time(boot.cnes <- bootSem(sem.cnes, R=5000))
> class(boot.cnes)
> boot.ci(boot.cnes)
>
> ------ snip --
>
> I hope this helps,
> Joh
mp; bca$fail) {
bca <- NULL
warning('could not obtain BCa bootstrap confidence interval')
} else {
bca <- bca$bca
m <- length(bca)
bca <- bca[c(m-1, m)]
}
list(np=np, bca=bca)
}
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
-
I'm not sure why you posted the original note. I quit using SAS in 1991 and
haven't needed it yet.
Frank
RogerJDeAngelis wrote
> Sorry about the double post. But I keep getting 'post' rejections, so I
> resubmitted about an hour later.
-
Frank Harrell
D
.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
Whoops - these 2 lines should have been omitted from the program:
n <- sprintf('%s (n%s=%g, n%s=%g)', v, nam[1],n[1], nam[2],n[2])
vn[var == v] <- n
Frank Harrell wrote
> I would like to have a lattice conditioning ( | var ) variable have
> expression() as value
at=1:2, labels=levels(trt)),
x=list(relation='free', limits=list(c(0,1),c(0,13)))),
ylab='Treatment', layout=c(1,2))
Error in unique.default(x) :
unimplemented type 'expression' in 'HashTableSetup'
-
Frank Harrell
Department of Biostatistic
Mar 1, 2013 at 7:16 PM, Frank Harrell <
>>
>>> f.harrell@
>>
>>> > wrote:
>>>> The lattice package uses special logic to allow for multiple
>>>> left-hand-side
>>>> variables in a formula, e.g. y1 + y2 ~ x. Is there an elegan
Achim this is perfect. I had not seen Formula before. Thanks for writing
it!
Frank
Achim Zeileis-4 wrote
> On Fri, 1 Mar 2013, Frank Harrell wrote:
>
>> Thank you Bill. A temporary re-arrangement of the formula will allow me
>> to
>> do the usual subset= na.actio
[5] "disp" "drat" "wt" "qsec"
> [9] "vs" "am" "gear" "carb"
>
> Bill Dunlap
> Spotfire, TIBCO Software
> wdunlap tibco.com
>
>
>> -Original Messag
Thanks for your reply Gabor. That doesn't handle a mixture of factor and
numeric variables on the left hand side.
Frank
Gabor Grothendieck wrote
> On Fri, Mar 1, 2013 at 7:16 PM, Frank Harrell <
> f.harrell@
> > wrote:
>> The lattice package uses special logic to
usual
invocation of model.frame( ) causes R to try to do arithmetic addition to
create a single dependent variable.
Thanks
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/Multiple-left-hand-side-variables
>>>> labels as 0= no complication,1=coronary heart disease,
>>>>>> 2=retinopathy,
>>>>>> 3=
>>>>>> nephropathy.
>>>>>>
>>>>>>
>>>>>> I want to select only the first complication that oc
FALSE)
>
> coefmat <- cbind(est, se, tval, pval)
> colnames(coefmat) <- c("Estimate", "Std. Error",
>"t value", "Pr(>|t|)")
>
> cat("\nCoefficients:\n")
> printCoefmat( coefmat, P
__
> 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.
-
*RMS Short Course 2013*
Frank E. Harrell, Jr., Ph.D., Professor and Chair
Department of Biostatistics, Vanderbilt University School of Medicine
*March 4, 5 & 6, 2013*
8:00am - 4:30pm
Student Life Center Board of Trust Room
Vanderbilt University
Nashville Tennessee USA
See http://biostat.mc.vande
_
> 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 Biosta
For lrm fits, predict(fit, type='mean') predicts the mean Y, not a
probability.
Frank
__
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 prov
81 wrote
> Dear R users,
>
> Please allow to me ask for your help.
> I am currently using Frank Harrell Jr package "rms" to model ordinal
> logistic regression with proportional odds. In order to assess model
> predictive ability, C concordance index is displayed and equals
e a tiny example and do the
calculations by hand.
Frank
blackscorpio81 wrote
> Dear R users,
>
> Please allow to me ask for your help.
> I am currently using Frank Harrell Jr package "rms" to model ordinal
> logistic regression with proportional odds. In order to ass
Chris,
I've fixed Surv in rms. The fix will be in the next release. For now you
can do source('http://biostat.mc.vanderbilt.edu/tmp/Surv.s') after issuing
require(rms).
Frank
Frank Harrell wrote
> Chris,
>
> Thanks for sending the specifics. It appears that I
Surv(left, right, type = "interval2") ~
>
> 1)
>
>
>
> Model LikelihoodDiscrimination
>
>Ratio Test Indexes
>
> Obs4LR chi2 0.00R2 0.000
>
> Events 6d.f. 0g
individual variables.
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
http://r.789695.n4.nabble.com/suggestions-about-import-SAS-results-to-R-tp465p4654554.html
Sent from the R help mailing list archive at Nabble.com
at there might be some issues with the newest SAS version
>>> but give it a try...
>>>
>>> ----------
>>> This message, and any attachments, is for the intended
>>> r...{{dropped:2}}
>>>
&
1 - 100 of 452 matches
Mail list logo