similar to: No subject

Dear listers I am a very newbie with graphs in R. I have a pulmonary function prediction equation in the form of PVC = 1.1 - 0.45*age in years + 0.011*height in cm. How can I draw the corresponding nomogram?. I read the help for the design.nomogram function but it is too difficult for me. Excuse my ignorance. Any direct help will be appreciated. Thanks in advance.

Dear lister Few months ago I posted a question about the drawing of a nomogram through R routines. Any news? thanks _________________________________________________________________ http://join.msn.com/?page=features/junkmail -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Dear Lister I need to perform Generalized estimating equations on some data. Has 'R' any function or routine to do it? (GEE is same as GLM but with correlated y's) Thanks in advance Mostafa _________________________________________________________________ Protect your PC - get McAfee.com VirusScan Online http://clinic.mcafee.com/clinic/ibuy/campaign.asp?cid=3963

Dear colleague: Actually that is what is done. When using the z-test between proportions in two different groups, or using chi-squared to test the null hypothesis of equal proportions of two or more groups, the null hypothesis is that H0: p1=p2=p3.....=p where this p is the pooled proportion, the proportion in the summed groups = n1+n2+n3+..../N1+N2+N3.... in the z-test between two

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

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

Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one without stratifying, however, I will need to stratify to meet PH assumptions. This is where I go wrong, but I'm not sure where. Non-Stratified Nomogram:

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

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

Hi All, I would like to use anonymous cvs, but it appears not to be working (again?). There was a discussion back in Jan-Feb about whether to continue supporting it, but it seemed that Tony Rossini got it working and the discussion left off there. Did someone decide to disable it, or is it just not working properly? Here's the details: $ cvs -d

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

Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11, surv=T,

I am a new user of R for Windows, enthusiast about the many functions of the Design and Hmisc libraries. I combined the results of a Cox regression model after multiple imputation (of missing values in some covariates). Now I got my vector of coefficients (and of standard errors). My question is: How could I use directly that vector to run programs such as 'nomogram', 'calibrate',

Hello, I am using Dr. Harrell's rms package to make a nomogram. I was able to make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: library(rms) f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11,

Hi R users, I'd like to draw a nomogram using a competing-risks regression (crr function in R), rather than a cox regression. However, the nomogram function provided in the Design package is not good for this purpose. Do you have any suggestion. I really appreciate your help Many thanks F.Abdollah, MD San-Raffele hospital Milan, Italy -- View this message in context:

Dear R users, I ran into this curious problem: > d <- rnorm(100) > d.density <- density(d) > sum( d.density$x * d.density$y) [1] 2.517502 Admittedly the method of computing the mass under the density curve at line 3 is crude. But 2.5 is pretty far from 1, the value it should be. I tried a few other dataset and got similar result. Am I missing something obvious? Or is the return

Hi all R users I did a logistic regression with my binary variable Y (0/1) and 2 explanatory variables. Now I try to draw my nomogram with predictive value. I visited the help of R but I have problem to understand well the example. When I use glm fonction, I have a problem, thus I use lrm. My code is: modele<-lrm(Y~L+P,data=donnee) fun<- function(x) plogis(x-modele$coef[1]+modele$coef[2])

Dear R users, I am a new R user and something stops me when I try to write a academic article. I want to make a nomogram to predict the risk of prostate cancer (PCa) using several factors which have been selected from the Logistic regression run under the SPSS. Always, a calibration plot is needed to validate the prediction accuracy of the nomogram. However, I tried many times and read a lot of

Hi R experts, I am trying to do a prognostic model validation study, using cancer survival data. There are 2 data sets - 1500 cases used to develop a nomogram, and another of 800 cases used as an independent validation cohort. I have validated the nomogram in the original data (easy with the Design tools), and then want to show that it also has good results with the independent data using 60