[RsR] library cox robust : weights?
Elisabetta Mattiolo
bett@to|o83 @end|ng |rom exc|te@|t
Tue Feb 7 22:12:34 CET 2012
I'm a italian student and I use R for my final degree thesis ... I want
study a Robustly Proportional Hazards (Bednarski's teory) for a dataset
of medical result.
I have exstimate with library coxrobust a coxr model (exponential
weights)
> s=Surv(survival,status)
site , cyto5.6 and zeb1 are my esplicative variables so the model is:
> exp=coxr(formula = s ~ site + cyto5.6 + zeb1, data = IHC, trunc =
0.95, f.w= "exp" , singular.ok = TRUE, model = FALSE)
with plot(exp) I can see 5 graphs
- the firs shows the standardizzed survival differences : one with Cox
model, and one green with Kaplan -Meir stimator
-other four show the same differences for four strata, defined by the
quartiles of the estimated linear predictor.
But my problem is the I want a graphic the robust exponential weight
(log trasformed) versus case number for the dataset..
If I ask R about weight of the model exp:
> exp$weights
NULL
If I write
> plot(exp$weights)
Error in plot.window(...) : 'xlim' devono essere finiti
Inoltre: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
3: In min(x) : no non-missing arguments to min; returning Inf
4: In max(x) : no non-missing arguments to max; returning -Inf
Can you help me?
Thanks and best regards
Elisabetta
<P><p><font face="Arial, Helvetica, sans-serif" size="2" style="font-size:13.5px">_______________________________________________________________<BR>1 richiesta di amicizia<br>
<a href="http://www.facebook.com/pages/Excite-Italia/101704826535618">Diventa fan di Excite su Facebook!</a></font>
[[alternative HTML version deleted]]
More information about the R-SIG-Robust
mailing list