similar to: How to generate a prediction equation for a stratified survival model that was fitted by cph() in Design package

Dear all, I am totally confused by the following R output, but don't have a clue for it. > a <- 1 - 0.2 > a == 0.8 [1] TRUE > a <- 1 - 0.8 > a == 0.2 [1] FALSE > a <- 1 - 0.5 > a == 0.5 [1] TRUE > a <- 1 - 0.6 > a == 0.4 [1] TRUE > a <- 1 - 0.9 > a == 0.1 [1] FALSE My R version is Windows XP R version 2.8.1 (2008-12-22).

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 hoping for assistance in regards to examining the contribution of stratified variables in a cox regression. A previous post by Terry Therneau noted that "That is the point of a strata; you are declaring a variable to NOT be proportional hazards, and thus there is no single "hazard ratio" that describes it". Given this purpose of stratification, in the

Using R version 1.8.1 for Windows, I obtain an error message using the following code. The data frame was constructed in the counting process style, where V1 is the start time, V2 is the stop time, and V3 is the censoring indicator. There are no zero-length time intervals. Variable V4 is the stratification factor (gender: F,M). S<-Surv(V1,V2,V3) fit<-survfit(S ~ V4,data=test.dat)

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,

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,

Dear All, I was writing a small wrapper to bootstrap a classification algorithm, but if we generate the indices in the "usual way" as: bootindex <- sample(index, N, replace = TRUE) there is a non-zero probability that all the samples belong to only one class, thus leading to problems in the fitting (or that some classes will end up with only one sample, which will be a problem

Rers: What is the preferred library/function for doing stratified random sampling from a dataset, given I want to control the number of samples (rather than the proportion of samples) per strata? Thanks! --j -- Jonathan A. Greenberg, PhD Postdoctoral Scholar Center for Spatial Technologies and Remote Sensing (CSTARS) University of California, Davis One Shields Avenue The Barn, Room

Dear all, I used fit.mult.impute in Dr. Harrell's Design package to fit a cox ph regression model on five imputed datasets, where all missing predictors were filled by multiple imputation using R package Mice. Are there any functions able to do bootstrapping or cross-validation for the aggregated model? I tried function 'validate' and 'calibrate' in Design package, but

Dear R-help, I am having a problem with the interpretation of result from validate.cph in the Design package. My purpose is to fit a cox model and validate the Somer's Dxy. I used the hypothetical data given in the help manual with modification to the cox model fit. My research problem is very similar to this example. This is the model without stratification: > library(Design) > f1

Hello, I'm fitting poisson regression to mortality data and wish to stratify by age. Is there any way to perform this stratification and use the glm function in R? Thanks, Hannah Murdoch

Dear Experts, I am trying to do a time-stratified case-crossover analysis on air pollution data and number of myocardial infarctions. In order to avoid model selection bias, I started with a simple simulation. I'm still not sure if my simulation is right. But the results I get from the "ts-case-crossover" are much more variable than those from a glm. Is this: a. Due to the simple

Hi All, I'm using R's randomForest package (and it's quite awesome!) but I'd really like to do some stratified sampling with a regression problem. However, it appears that the package was designed to only accommodate stratified sampling for classification purposes (see https://stat.ethz.ch/pipermail/r-help/2006-November/117477.html). As Andy suggests in the link just

This is a follow-up to a query that was posted regarding some problems that emerge when running anova analyses for cox models, posted by Mathias Gondan: Matthias Gondan wrote: >* Dear List,*>**>* I have tried a stratified Cox Regression, it is working fine, except for*>* the "Anova"-Tests:*>**>* Here the commands (should work out of the box):*>**>*

Hi all: For stratified count data,how to perform regression analysis? My data: age case oc count 1 1 1 21 1 1 2 26 1 2 1 17 1 2 2 59 2 1 1 18 2 1 2 88 2 2 1 7 2 2 2 95 age: 1:<40y 2:>40y case: 1:patient 2:health oc: 1:use drug 2:not use drug My purpose: Anaysis whether case and

Hi, Has anyone done (or know of) any nice R activities that help introductory students ( and teachers :) ) better understand the concepts of simple vs stratified vs cluster sampling? Any links? David -- View this message in context: http://r.789695.n4.nabble.com/SRS-Stratified-and-Cluster-sampling-tp4655099.html Sent from the R help mailing list archive at Nabble.com.

Dear List, I have tried a stratified Cox Regression, it is working fine, except for the "Anova"-Tests: Here the commands (should work out of the box): library(survival) d = colon[colon$etype==2, ] m = coxph(Surv(time, status) ~ strata(sex) + rx, data=d) summary(m) # Printout ok anova(m, test='Chisq') This is the output of the anova command: > Analysis of Deviance Table

Dear R-List,   I would like to have a large number of stratified random subsamples drawn from my dataframe and automatically test for correlation differences in every subsample.   Let this be my dataframe   df<-data.frame(group=c(rep(1,5),rep(2,5),rep(3,5)),a=c(3,4,5,6,3,4,5,4,5,4,1,2,1,2,1),b=c(1,2,3,4,5,3,4,3,4,5,6,5,6,2,3),c=c(2,2,3,3,5,1,1,6,6,5,6,1,1,2,1))   Then I would like to have n

Dear R community, I have data on beta&standard error (for the main effect of variable x), stratified by sex for my dataset. I wish to calculate the sex-interaction effect (as beta&se) from these two stratified datasets. Is there a package to do this? If not, any advice how to do it manually? Thank you very much and best regards, Georg. ************************ Georg Ehret, JHU,

I?m trying to set up proportional hazard model that is stratified with respect to covariate 1 and has an interaction between covariate 1 and another variable, covariate 2. Both variables are categorical. In the following, I try to illustrate the two problems that I?ve encountered, using the lung dataset. The first problem is the warning: To me, it seems that there are too many dummies