similar to: Nomogram

Dear list, I'm facing the following problem : A cox model with my sex variable interacting with several continuous variables : cph(S~sex*(x1+x2+x3)) And I'd like to make a nomogram. I know it's a bit tricky and one mights argue that nomogram is not a good a choice... I could use the parameter interact=list(sex=("male","female"),x1=c(a,b,c))... but with rcs or pol

Dear R-users, I am a novice R-user with some experience in using the RMS package for taking nomograms after various survival models. This time, I am trying to plot a nomogram after a Random Effects Cox, implemented by the "coxme" package. My questions are: 1. Is it possible to take a nomogram directly after the coxme survival function? 2. If not is there a way to take the linear

Hi R-help, ? I have fit a cox ph model to my data, but have beenreceiving an error when trying to fit a model to the nomogram. Here is the codeand corresponding error: ? ? >nomogramCF = nomogram(cph_age6_40avp4, +????????????????????lp.at= seq(-10,10,by =0.5),lp = TRUE, +?????????????????????? +??????????????????????funlabel="5year survival",

Inspired by the functions that Barry Rawlingson and Dave Forrest posted for searching Rwiki and R-help archive, I've made up a function that does the search on Prof. Baron's site (Thanks to Prof. Baron's help on setting up the query string!): RSiteSearch <- function(string, restrict="Rhelp", format="long", sortby="score",

Hi R users. I have a model of cox that already it is estimated (I have only the model estimated, I haven't data), how can I determine a nomogram with R? Is it posible to do nomograms in Design package? I think that's only when the model (Cox Regression in this case) is before estimated in R. Thanks

I've made a few changes in the indices in my R search site. http://finzi.psych.upenn.edu is the general URL, and RSiteSearch takes you to http://finzi.psych.upenn.edu/search.html, which is equivalent to http://finzi.psych.upenn.edu/nmz.html. The main change is that I have eliminated "Previous message" and "Next message" from each page. (I still have "Next in

The rms package includes the nomogram function, which generates a list object that can be passed to plot for graphical production of nomograms. I would like to remove the "linear predictor" line in the graph, which means (I suspect) removing it from the nomogram output object. I've looked at the nomogram output object, but it is not clear to me if or how it might be edited to

First, my site will be down December 27-28 because of a network upgrade at Penn. It will also be down at least one day before that, while I upgrade the operating system. (And another day some time in January because of a planned power outage.) Second, I have replaced the search engine in my R site: http://finzi.psych.upenn.edu/ I am now using Namazu instead of HtDig. The direct link to the

Hi all I have a question regarding differences in the way gregpexr works in R 2.3.0 and R 2.4.0. In R 2.3.0, this is what happens: > gregexpr(" [a-z] [a-z] ", " a b c d e f ", perl=T) [[1]] [1] 1 3 5 7 9 attr(,"match.length") [1] 5 5 5 5 5 ... while in R 2.4.0, this is what happens: > gregexpr(" [a-z] [a-z] ", " a b c d e f ", perl=T)

The CRAN host in Tampa, FL (cran.hostingzero.net) isn't responding, and hasn't responded in quite a while -- at least problems were reported more than a year ago (Oct 2007) http://finzi.psych.upenn.edu/R/Rhelp02a/archive/116706.html although someone apparently succeeded in July 2008 http://finzi.psych.upenn.edu/R/Rhelp02a/archive/137451.html Perhaps it's just flaky, and not

Here is a summary of responses to my original email (see my query at the bottom). Thank you to Achim Zeileis , Anders Nielsen, Pierre Kleiber and Dave Fournier who all helped out with advice. I hope that their responses will help some of you too. ***************************************** Check out glm.nb() from package MASS fits negative binomial GLMs.

There is both a technical and a theoretical element to my question... Should I be able to use the outputs which arise from the fastbw function as inputs to nomogram(). I seem to be failing at this, -- I obtain a subscript out of range error. That I can't do this may speak to technical failings, but I suspect it is because Prof Harrell thinks/knows it injudicious. However, I can't

Hi All, Sorry that I didn't provide enough information. I've been trying to import SAS xport files that contain multiple files using package foreign's read.xport. I first attempted this back in 2005 and had problems. Some of files that were present in the SAS xport file weren't being created in R. I submitted my problem to the community:

Hi, I used Dr. Harrell's rms package to make a nomogram. Below is my code for nomogram and calculate total points and probability *in original data set* used for building the nomogram. *My question is how I get the formula for calculating the survival probability for this nomogram. Then I can use this formula to do validation by using other data set. * f1 <- cph(Surv(retime,dfs) ~

Hello, I am constructing a nomogram using cph and nomogram commands in Dr. Harrell's Design/RMS package. The HR that I obtain for dichotomous and categorical variables are identical to those that I obtain using STATA stcox. However, the inter-quartile HR I obtain for continuous variables is obviously different, since STATA gives me HR for each unit (year, centimeter, etc) like coxph would

How can I get prediction intervals from a mixed-effects model? Consider the following example: library(nlme) fm3 <- lme(distance ~ age*Sex, data = Orthodont, random = ~ 1) df3.1 <- with(Orthodont, data.frame(age=seq(5, 20, 5), Subject=rep(Subject[1], 4), Sex=rep(Sex[1], 4))) predict(fm3, df3.1, interval='prediction') # M01 M01

Hi, I would like to emphasize ("zoom") the zone of a nomogram where the probability are > 0.01 (nomogram built with nomogram, Design). As a consequence, I don't need to draw the part of the "Total points" axis with score < 60 equivalent in my case to a linear predictor < 4.5 - As far as I know, this is not possible with the arguments of the function. - Changing

Dear All, I have the following function tstpar <- function(n = 200, want.pdf = FALSE, pdfFileName = NULL){ oldpar <- par(no.readonly = TRUE) on.exit(par(oldpar)) steps <- seq(from = 1, to = 8, by = 1) h <- 10; w <- 6 if(want.pdf){pdf(file = pdfFileName, onefile = TRUE, paper = "letter", width = w, height = h)} par(mfrow = c(4,2)) for(i

Dear colleagues, I hope someone can help me with my problem. I have fitted a cox model with the following syntax: # cox01def <-cph(Surv(TEVENT,EVENT) ~ ifelse(AGE>50, (AGE-50)^2,0) + BMI + # HDL+DIABETES +HISTCAR2 + log(CREAT)+ as.factor(ALBUMIN)+STENOSIS+IMT,data # = XC, x=T, y=T, surv=T) *1 Furthermore I have estimated my beta's also with a Lasso method - Coxpath ( from

Dear R-users, I have used the nomogram function from the rms package for a logistic regresison model made with lrm(). Everything works perfectly (r version 2.15.1 on a mac). My question is this: if my final model is not the one created by lrm, but I internally validated the model and 'shrunk' the regression coefficients and computed a new intercept, how can I build a nomogram using that