[R] Error in solve.default peforming Competing risk regression

2009-10-27 Thread kende jan

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.5721568
14.6196
14.664159
14.825136
15.151459
15.1852111
15.333158
15.32177
15.348338
15.473277
15.5441453
15.7331363
15.866132
15.9157
15.9411870
16128
16.09153
16.302123
16.6671196
16.647245
16.7271132
16.304339
16.746161
16.786150
17125
17.044   1252
17.164242
17.6903134
17.802180
18.3981133
18.666132
18.687341
18.794129
18.887193
19.172348
19.7961275
111.361230
111.56147
112.12175
112.142168
112.411318
116.50167
117.211510


  
[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Error in solve.default peforming Competing risk regression

2009-10-27 Thread Ravi Varadhan

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