similar to: why the survival function estimate using package 'mstate' & package 'cmprsk' vary from sas and LTA (From WHO).

Dear all, in order to estimate transition-specific probabilities in a multi-state model i applied the 'probtrans()' function from the 'mstate'-package. Now, i am at loss with the following message (see attached example): Warning message: In probtrans(msf.0, predt = 0) : Negative diagonal elements of (I+dA); the estimate may not be meaningful. I am not very familiar with matrix

Hi, I was trying to get hyperlinks using xtable, but couldn't get the hyperlinks to function properly. For example, if I use ## Try to link NY times website to every figure in column 4 mat <- matrix(1:43,6,5) mat[,5] <- "http://nytimes.com" for(i in 1:nrow(mat)){ strr <- paste('<a href="', mat[i,5],'">', mat[i,4], '</a>',

I am new to the packge twitteR when using this package I came across several problems, can anybody help me on those problems? Many Thanks rm(list=ls()) library(bitops) library(RCurl) library(rjson) library(twitteR) library(digest) library(ROAuth) requestURL <- "https://api.twitter.com/oauth/request_token" accessURL =

> Date: Fri, 9 Sep 2011 01:07:50 -0700 (PDT) > From: komkon555 <komkon555@freenet.de> > Subject: Re: Re : Re : Re : [Xen-devel] Re: Patches for > VGA-Passthrough XEN 4.2 unstable > To: xen-devel@lists.xensource.com > Message-ID: <1315555670684-4785425.post@n5.nabble.com> > Content-Type: text/plain; charset=us-ascii > > > Weidong Han wrote: >

I am trying to use the crr function in the cmprsk package to analyze a large patient dataset (45000 +), The model has 100 + covariates and 5 competing risks. I am finding that R seems to get bogged down and even if I let it run for several hours I don't get anything back. Am I expecting too much, or are there ways to speed up the process? Any help is appreciated. Best, Spencer

Hi all, I am aware that crr lacks the "friendly" command structure of functions such as cph. All is clear to me about including covariates until I want to include a stratification term in the competing risk framework (no nice strat command). I am still a bit of a novice in R - I am looking for an example to help me with this, but can't seem to find one. Any advice appreciated (no

Hi, I need to fit model using crr, however my covariate is categorical with 3 levels. I use crr(time,status,agesplit,failcode=1,cencode=0) where agesplit is defined as <20,21-29,>30 years, so it takes 0, 1 or 2 for each patient. I hoped to get estimated coefficients for the levels 1 and 2 w.r.t level 0 as in coxph. But, I didn't. Could someone please help me to use crr in this

Hi, To use cuminc() from cmprsk package, if a subject has 2 events (both the event of interest and the event of competing risk), should I create 2 observations for this subject in the dataset, one for each event with different fstatus (1 and 2), or just 1 observation with whatever event that happened first? My analysis objective is calculate cumulative incidence for the event of interest.

Dear R-community, I would like to plot the numbers at risk for the different causes of failure at specific timepoints below a cumulative incidence function plot (plot.cuminc-function, cmprsk-package). For a Kaplan-Meier plot I know this is possible with the n.risk-argument in the survplot-function (rms-package), but to my knowledge no such readily-available functions are available for competing

Dear R-help list members, I am trying to create my own axes when plotting a cumulative incidence curve using the plot.cuminc function in the CMPRSK library. The default x-axis places tick marks and labels at 0, 20, 40, 60, and 80 (my data has an upper limit of 96), whereas I want them at my own specified locations. Here is my example code: library(cmprsk) attach(MYDATA) MYCUMINC <-

Dear R users, I d like to assess the effect of "treatment" covariate on a disease relapse risk with the package cmprsk. However, the effect of this covariate on survival is time-dependent (assessed with cox.zph): no significant effect during the first year of follow-up, then after 1 year a favorable effect is observed on survival (step function might be the correct way to say that

Dear Prof Gray and everyone, As our package developers discussed about incompatibility between Design and survival packages, I faced another problem with cmprsk- a survival dependent packacge. The problem is exactly similar to what happened to the Design package that when I just started running cuminc function, R was suddenly closed. These incidents suggest that maybe many other survival

Hello again, I m trying to use timereg package as you suggested (R2.7.1 on XP Pro). here is my script based on the example from timereg for a fine & gray model in which relt = time to event, rels = status 0/1/2 2=competing, 1=event of interest, 0=censored random = covariate I want to test library(timereg) rel<-read.csv("relapse2.csv", header = TRUE, sep = ",",

Dear R users, I d like to assess the effect of "treatment" covariate on a disease relapse risk with the package cmprsk. However, the effect of this covariate on survival is time-dependent (assessed with cox.zph): no significant effect during the first year of follow-up, then after 1 year a favorable effect is observed on survival (step function might be the correct way to say that ?).

I'm only doing one frame using speex_encode_int greatly simplifies my code I'm not sure why the sample I was working off of was converting the shorts to floats then calling the other encode/decode methods. Based off of your suggestions I tried the following but I get the same result. virtual Enigma::u8* Encode(Enigma::u8* inputBuffer,size_t inputSize, size_t& outputSize) {

Hi all I do not understand why I am getting the following error message. Can anybody help me with this? Thanks in advance. install.packages("cmprsk") library(cmprsk) result1 <-crr(ftime, fstatus, cov1, failcode=1, cencode=0 ) one.pout1 = predict(result1,cov1,X=cbind(1,one.z1,one.z2)) predict.crr(result1,cov1,X=cbind(1,one.z1,one.z2)) Error: could not find function

Alright noted, I changed me code so that the state is created in the constructor and destroyed in the destructor of the object. However I'm still getting the same issue although I'm sure that would have bit me sooner or later. The new code is as follows. virtual Enigma::u8* Encode(Enigma::u8* inputBuffer,size_t inputSize, size_t& outputSize) { short *in=(short*)inputBuffer;

Hi to all, I''d like to cut some log in /var/log/messages, as of netbios and ping entries. There are some particular rules in shorewall 1.4.5? I''ve tried with "run_iptables -A common -p udp --sport 138 -mstate --state NEW -j DROP" but it contiunes to send to log every netbios attempt. Also I don''t want to disable ping from loc to net, and from fw to net. Thanks

----- Mensagem Original -----=20 Eu boicotaria estas empresas s=F3 pelo fato de financiarem os enjoativos e = rancentos Cassetas, mas me surpreendi ao descobrir que tantos produtos que = eu costumo comprar s=E3o transg=EAnicos. Fui conferir no site do greenpeace= e =E9 verdade!!!!!!!!!!!!! Boicote neles! ----- Original Message -----=20 C A M P A N H A=20 N A T A L S E M B A I X A R I A

I found an issue with the data() command this evening when working on the survival package. 1. I have a lot of data sets in the package, almost all used in at least one vignette, help file, or test.? As a space saving measure, I have bundled many of them together, i.e., the file data/cancer.rda contains 19 data sets, many of them small. The resulting file (using xz compression) is quite a bit