similar to: Adjustment for multiple-comparison for log-rank test

Dear All, I would like to make a quadratic with a plateau model in R. Is there a package in R doing this? The bentcableAR package seems won't work. The link below describes what I am looking for in R exactly: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_nlin_sect033.htm -- Thanks so much! Orange help.ly2005@gmail.com [[alternative HTML version

Hi All, I am trying to shift from running mixed models in SAS using PROC MIXED to using lme4 package in R. In trying to match the coefficients of R output to that of SAS output, I came across this problem. The dataset I am using is this one: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect034.htm If I run the following code: proc mixed data=rc

Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a "weights=albumin" argument (I only chose albumin because it had no missing values, not because of any real relevance). Below are my code with the resulting error

thunderbird just crashed, and I've never seen this before. Any idea what happened? this was in the dovecot log: 2010-12-02 13:06:40 imap-login: Info: Maximum number of connections from user+IP exceeded (mail_max_userip_connections): user=<pbc>, method=PLAIN, rip=127.0.0.1, lip=127.0.0.1, secured 2010-12-02 13:06:40 imap-login: Info: Maximum number of connections from user+IP exceeded

Dear R-users: I can not reproduce figure 4 of *Statistical Tables and Plots using S and LATEX* by Prof. Frank Harrell with the following code: rm(list=ls()) library(Hmisc) getHdata(pbc) attach(pbc) age.groups <- cut2(age, c(45,60)) g <- function(y) apply(y, 2, quantile, c(.25,.5,.75)) y <- with(pbc, cbind(Chol=chol,Bili=bili)) # You can give new column names that are not legal S names

The data set in R is wrong. I've found mistakes on 2 lines in a quick look. I don't know if the data is incorrect in the Appendix of Fleming and Harrington as well (someone seems to have borrowed my copy), which is where the data set appears to have been taken from, given all the "-9" codes in it. (Note, Tom Fleming originally got the data from me, so I'm fairly

Dear All, I'm trying to teach myself latex along with the latex function in Hmisc and have hit a roadblock that I can't seem to get around. I'd greatly appreciate any pointers. I'm running R 2.6.0 on Windows XP and have Miktex 2.7 installed. I've reproduced the code below, taken from Frank Harrell's latexsummary introduction. My question relates to getting a pdf

Hello, I'm attempting to use the ratetable argument to survexp in the survival package. I use the example from the ?survexp help page below, and then slightly modify it to produce an error. library(survival) data(pbc) #fit a model without any factors pfit1 <- coxph(Surv(time, status > 0) ~ trt + log(bili) + log(protime) + age + platelet, data=pbc) #this works as expected

I would like to determine the probability of an event at a specific timepoint given the linear predictor of a given Cox model. For instance, assume that I fit the following model: data(pbc) fit <- coxph(Surv(time, status)~ age, data=pbc) To extract the value of the linear predictor for each patient in the dataset: prd <- predict(fit, newdata=pbc, type="lp")

Hi, When I run the following Cox proportional hazards model on the Mayo clinic's PBC data set (given in the "survival" package), the regression coefficients do not agree with the results presented in Table 4.6.3 (p. 195) of Fleming & Harrington's book. library(survival) data(pbc) ans.cox <- coxph(Surv(time, status) ~ log(bili) + log(alb) + age + log(protime) +

Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted model (e.g., to obtain a plot of the relative risk (or log r.r.) versus the predictors). When there is only 1 predictor in the model, termplot returns a "?". In this case, I have not been able to figure out how to create

R-helpers: I am using R 2.5 on Windows XP, packages all up to date. I have run into an issue with the MIcombine function of the mitools package that I hoped some of you might be able to help with. I will work through a reproducible example to demonstrate the issue. First, make a dataset from the pbc dataset in the survival package --------------- # Make a dataset library(survival) d <-

Hi there, Can anyone please help me because I am going to get crazy with the pbc data set. I just want to apply simple cox regression in the data set. I am a beginner in R but I don't think I am doing anything wrong. I have the book of Fleming and Harrington 1990. I perform cox regression by typing: out<- coxph(Surv(times/365,status)~log(bili)+log(proth)+edema+log(albumin)+age) out

R-helpers: I have a dataset that has 168 subjects and 12 variables. Some of the variables have missing data and I want to use the multiple imputation capabilities of the "mice" package to address the missing data. Given that mice only supports linear models and generalized linear models (via the lm.mids and glm.mids functions) and that I need to fit Cox models, I followed the previous

I have attempted to follow posting guidelines but I have failed to find out what I am doing wrong here. I am trying to use lines.survfit to plot a second curve onto a survival curve produced by plot.survfit. In my case this is to be a progression free survival curve superimposed upon an overall survival curve, but I will illustrate my problem using the example given in the help for

Hi: I have a few clarification questions about the elements returned by the coxph function used in conjuction with a frailty term. I create the following group variable: group <- NULL group[id<50] <- 1 group[id>=50 & id<100] <- 2 group[id>=100 & id<150] <- 3 group[id>=150 & id<200] <- 4 group[id>=200 & id<250] <- 5 group[id>=250

Dear R help, I have some readings in three dimensions (x, y, z) and an amplitude for each. I'd like to visualize the data using ggplot, using tile plots, as I have some additional point data I would like to eventually overlay on the tile plots. I would like to subset the data by sections, slices if you will, in the z dimension, and plot the data for that slice. I can do all of this, but am

Well, in my case, i have samba running as a server, with auth. to 2 samba pbc/bdc servers this server uses winbind to resolve usernames and groups. i get the error: [ date ] nsswitch/winbind_user.c:winbind_getpwnam(159) user 'collen' does not exist. [ date ] nsswitch/winbind_user.c:winbind_getpwnam(159) user 'COLLEN' does not exist. and the user in question can not acces the

Dear R-users, I'd like to announce the release of the new package JM (JM_0.1-0 available from CRAN) for the joint modelling of longitudinal and time-to-event data. The package has a single model-fitting function called jointModel(), which accepts as main arguments a linear mixed effects object fit returned by function lme() of package nlme, and a survival object fit returned by either

Dear R-users, I'd like to announce the release of the new package JM (JM_0.1-0 available from CRAN) for the joint modelling of longitudinal and time-to-event data. The package has a single model-fitting function called jointModel(), which accepts as main arguments a linear mixed effects object fit returned by function lme() of package nlme, and a survival object fit returned by either