similar to: Survival Plots by Strata

All, I had a question about the JGR console and whether or not I can manipulate the location where line wrapping occurs. I have searched 'JGR' in the R listserve archives and attempted to find console manipulation on the JGR website to no avail and could use some direction. TIA, Bret As an example, the below output wraps every 4th value, leaving about 2/3 of the console empty. >

Dear all, I am using Windows and R 2.9.2 for my analyses. I have a large dataset and I am particularly interested in looking at time to an event for a continuous variable. I would like to produce a plot of log(relative risk) or relative risk (also known as hazard ratio) against the continuous variable. I have spent a long time looking for advice on how to do this but my search has proved

R Users, I was wondering if someone might point me in the right direction. I am using a Cox model (survival package) to evaluate survival of pen-reared birds (time to event data collected daily) and I have been trying to determine how I can estimate a 'daily' survival rate and std error from the results of a Cox model? Using survfit (see input data below) I computed the predicted

Hi All, The following is a snippet of my code. It works fine but it is very slow. Is it possible to speed it up by using different data structure or better solution? For 40000 runs, it takes 8 minutes now. Thanks a lot fun_activation = function(s_k, s_hat, x_k, s_hat_len) { common = intersect(s_k, s_hat); if(length(common) != 0) { index = match(common, s_k);

Hello everybody. First of all, I would like to present myself. I'm a french student in public health and I like statistics though I'm not that good in mathematics (but I try to catch up). I've discovered R recently while trying to find a statistical program in order to avoid rebooting my computer under windows when I need to do some statistical work. And here is my first question.

Hi, I want to fit multiple cox models using the coxph() function. To do this, I use a for-loop and save the relevant results in a separate matrix. In the example below, only two models are fitted (my actual matrix has many more columns), one gives a warning message, while the other does not. Right now, I see all the warning message(s) after the for-loop is completed but have no idea which model

Hello, coxph does not use any information that are in the dataset between event times (or "death times") , since computation only occurs at event times. For instance, removing observations when there is no event at that time in the whole dataset does not change the results: > set.seed(1) > data <- as.data.frame(cbind(start=c(1:5,1:5,1:4),stop=c(2:6,2:6,2:5),status=c(rep(

I'm looking at memory efficiency for some of the survival code. The following fragment appears in coxph.fit coxfit <- .C("coxfit2", iter=as.integer(maxiter), as.integer(n), as.integer(nvar), stime, sstat, x= x[sorted,] , ... Does this make a second copy of x to pass to the routine (my

Hello, I'm trying to estimate a cox model with a frailty variable and time-dependent covariate (below there is the statement I use and the error message). It's seems to be impossible, because every time I add the time-dependent covariate the model doesn't converge. Instead, if I estimate the same model without the time-dependent covariate it's converge. I'd like knowing if

Dear Thomas & Dear List, the fitting function `coxph.fit' called by `coxph' may fail to estimate the regression coefficients when some values of the design matrix are very large. For example library(survival) ### load example data load(url("http://www.imbe.med.uni-erlangen.de/~hothorn/coxph_fit.Rda")) method <- "efron" ### copied from `coxph.fit' coxfit

dear all, I would like to know whether survdiff and survReg function in the survival package work for left-truncated and right-censored data. If not, what other functions can i use to make comparison between two survival curves with LTRC data. thanks for any help given sing yee

Hi, I looked though documentation and also checked last 9 months of archives, but could not find the corresponding information. We are having an issue with incremental backup with arkeia. Arkeia checks if ctime/mtime of the file are changed, if so, it backs up. I realized that when I ran rsync from a cammnd line, it changed ctime, and preserved atime of the file in a destination as

It looks like I need to elaborate further to get a feedback. I checked the rsync source code and it is using utime() to restore atime file attribute. It is fine to change ctime for that transferred file in this case. What we are having problem is that when rsync gets kicked off and transfers one file to the destination, this action changes ctime of "all" files in the same

I have a debug tree with shared objects enabled (to help streamline the development process), and I see those unit tests failing where they don't fail on a release/static build. Anyone seeing these, too? FAIL: LLVM-Unit :: Object/ObjectTests/StringTableBuilderTest.Basic (19422 of 20473) ******************** TEST 'LLVM-Unit :: Object/ObjectTests/StringTableBuilderTest.Basic' FAILED

Hi, I'm trying to convert a vector containing dates in character format ("dd/mm/yy"): mdy.date (from date package) seems to be able to do that, but it returns to me a vector containing julian dates... but negative! for example: 16/12/03 is converted into -20470 it is because R recognizes year ../03 as 1903, instead of 2003, but how can I do to solve this problem? (of course, I

Hi, We seem to be having an interesting issue with Asterisk whereby, it keeps routing calls coming in to the 'default' context.... regardless of what changes occur to h323.conf. <SNIP> [POP-A] type=user host=1.2.3.4 context=international </SNIP> == Starting H323/ip$1.2.3.4:12914/16313 at default,12126599878,1 failed so falling back to exten 's' == Starting

buen día colegas: resulta que hace unos años hice un programa para hacer regresiones lineales pero los datos venían en esta forma: tiempo degradacion elemento1 9204 14 10306 15 9860 16 11803 17 12211 18 11333 19 13490 126 12329 229 13325 230 13214 231 13645 232 14492 233 14570 234 14777 20 4187 33 5607 34 5803 35 5997 36 5529 37 5716 38 6164 325 15451 428 19062 429 17753 430 20465 431 20966 432

This one puzzles me a lot: [thba at vink layout]$ nc script/xw_functions.ample NEdit: getcwd() fails: No such file or directory NEdit: getcwd() fails: No such file or directory [thba at vink layout]$ ll script/xw_functions.ample -rw-rw-r-- 1 thba thba 16829 Oct 25 16:59 script/xw_functions.ample [thba at vink layout]$ pwd /home/thba/workarea/colibri/design/ana/layout [thba at vink layout]$

good day to you i have compiled evertything related to ip tables as modules, but still i get this error when trying to use /etc/shorewall/routes can someone tell me the spicific module i need to have? here are my info intranet linux # shorewall version 2.4.2 intranet linux # uname -a Linux intranet 2.6.12-gentoo-r10 #10 SMP Sun Sep 11 15:01:49 SAST 2005 i686 AMD Athlon(tm) XP 2400+

hail, I partitioned the disk this way: fdisk da0 ******* Working on device /dev/da0 ******* parameters extracted from in-core disklabel are: cylinders=12161 heads=255 sectors/track=63 (16065 blks/cyl) Figures below won't work with BIOS for partitions not in cyl 1 parameters to be used for BIOS calculations are: cylinders=12161 heads=255 sectors/track=63 (16065 blks/cyl) Media sector size