Displaying 5 results from an estimated 5 matches for "recidivate".
2008 Dec 28
1
cox regression warning/error messages
Hello,
I am hoping for some advice regarding warning/error messages I
received when running a Cox regression
# message 1 - obtained while creating a plot of residuals
> plot (NV.zph, main = "groupNUSM - UNFIT", var= 'groupNUSM')
Warning messages:
1: In approx(xx, xtime, seq(min(xx), max(xx), length.out = 17)[2 * :
collapsing to unique 'x' values
2: In
2008 Nov 29
0
Error in check(itp) : object does not represent a K sample problem with censored data
Hello,
I have two questions regarding a survival analysis I have been
working on. Below is the code to date.
The variables:
1) recidivism$intDaysUntilFVPO are the number of days before an
violent offence was committed - if no offence was committed than the
days between court hearing and end of data collection was recorded.
2) recidivism$intDaysUntilFNVPO are the number of days before a
2012 Nov 29
5
bootstrapped cox regression (rms package)
...ome success with the bootcov function in the rms
package, which at least generates confidence intervals similar to what is
observed in SPSS. However, the p-values associated with each predictor in
the model are not really close in many instances.
Here is the code I am using:
formula=Surv(months, recidivate) ~ fac1 + fac2 + fac3 + fac4 + fac5 + fac6
+ fac7 + fac8
fit=cph(formula, data=temp, x=T, y=T)
validate(fit, method="boot", B=9999, bw=F, type="residual", sls=0.05,
aics=0,force=NULL, estimates=TRUE, pr=FALSE)
out=bootcov(fit, B=9999, pr=F, coef.reps=T, loglik=F)
for (i in 1:8)...
2012 Nov 29
0
bootstrapped cox regression in rms package (non html!)
...ome success with the bootcov function in
the rms package, which at least generates confidence intervals similar
to what is observed in SPSS. However, the p-values associated with
each predictor in the model are not really close in many instances.
Here is the code I am using:
formula=Surv(months, recidivate) ~ fac1 + fac2 + fac3 + fac4 + fac5 +
fac6 + fac7 + fac8
fit=cph(formula, data=temp, x=T, y=T)
validate(fit, method="boot", B=9999, bw=F, type="residual", sls=0.05,
aics=0,force=NULL, estimates=TRUE, pr=FALSE)
out=bootcov(fit, B=9999, pr=F, coef.reps=T, loglik=F)
for (i in 1:8)...