Please send reproducible code when you report problems.
I am unable to reproduce your problem. It seems to work just fine.
dat - read.table(e:/temp/crrdata.txt, sep=\t, header=TRUE)
names(dat)
[1] x1 x2 fstatus ftime
dim(dat)
[1] 198 4
cov=cbind(dat$x1,dat$x2)
z - crr(dat$ftime, dat$fstatus, cov)
z
convergence: TRUE
coefficients:
[1] 0.360900 -0.002484
standard errors:
[1] 0.23180 0.03413
two-sided p-values:
[1] 0.12 0.94
So, what is your problem?
Also, make sure that you understand what you are doing with competing risks
analysis. It seems like you do not have any censoring, i.e. all the subjects
are failing from one of three causes. Read the documentation of ?crr carefully.
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins University
Ph. (410) 502-2619
email: rvarad...@jhmi.edu
- Original Message -
From: kende jan kende...@yahoo.fr
Date: Tuesday, October 27, 2009 3:31 am
Subject: [R] Error in solve.default peforming Competing risk regression
To: R forum r-help@r-project.org
Dear all,
I am trying to use the crr function in the cmprsk
package version 2.2 to analyse 198 observations.I have receive the
error in solve.default.
Can anyone give me some
insights into where the problem is?
Thanks
here is my script :
cov=cbind(x1,x2)
z-crr(ftime,fstatus,cov))
and data file:
x1x2fstatusftime
0.021263
0.031113
0.031523
0.03349
0.031278
0.031190
0.041472
0.041357
0.041219
0.0429349
0.051207
0.051166
0.06355
0.06217
0.071264
0.07267
0.07337
0.082190
0.08372
0.08236
0.08196
0.092136
0.097681273
0.11167
0.1276
0.13107
0.1371
0.11116
0.11262
0.111256
0.111385
0.111266
0.111174
0.113334
0.123270
0.12149
0.121733
0.131271
0.132160
0.131105
0.14256
0.14249
0.151266
0.15146
0.151112
0.151246
0.161371
0.161140
0.161279
0.161174
0.1617180
0.17173
0.17354
0.171320
0.171213
0.181215
0.181282
0.181263
0.18399
0.21266
0.2159
0.2097611659
0.211111
0.21162
0.221386
0.241249
0.24279
0.261199
0.26195
0.261270
0.261189
0.26706321
0.27183
0.27145
0.291154
0.291221
0.291174
0.29241112
0.31208
0.3231
0.32182
0.332152
1.341595
0.371182
0.381231
0.381282
1.38160
0.421181
0.441103
0.453251
0.452176
0.471234
0.49326
1.53142
0.5551152
0.561229
1.571179
0.58189
0.62187361
0.66332
1.69114
1.69171
0.7594190
1.9462391
0.951138
1.969191
11.13129
01.171145
11.26146
11.361187
01.36174
11.48169
01.61102
11.62174
11.92103
11.91152
12.11159
12.161238
12.201140
02.231124
12.42161
12.451283
12.47142
12.582291
12.591182
12.73328
12.731459
13.0071179
13.04152
13.167385
13.2451559
13.251320
13.2331482
13.275182
13.37141982
13.3652142
13.411184
13.4833187
13.5324131
13.58181374
13.629149
13.63681297
13.7193
13.781364
13.884156
13.88166
13.891160
13.9361205
14.042161
14.38299
14.3921311
14.411257
14.48197
14.53176
14.510194
14.5721