[R] sequential treatment of a vector for formula

[Wu Gong]

Thank you very much David. I'm sorry for this fault , hope it has not
confused Frostygoat.

I was clueless of recursive reference and I didn't meet any error when I
test the code. So I wonder if there are some useful tips to prevent making
this kind of faults:)

The revised code is followed.

### Kaplan-Meier estimate
### Compute surviving proportion of a cohort
## Original data is a vector of death days
example <- as.vector(c(8,2,1,5,6,7,7,0,8,10,13,8,11,11,11,2,7,1,5,6,8,6)) 

## Count the frequency of each death day
tb <- table(example)

## Output
## Interval(Start-End)  "Interval"
## At risk at start of interval "RSI"
## Deaths during interval "Deaths"
## At risk at end of interval "REI"
## Propoition surviving this interval "PSI" = REI/RSI
## Cumulative surviving at end of interval "CSI" = cumprod(PSI)
iC <- data.frame("Interval"=as.numeric(names(tb)),
	"RSI"=sum(tb)-cumsum(tb)+tb[],
	"Deaths"=tb[],
	"REI"=sum(tb)-cumsum(tb),
	"PSI"=(sum(tb)-cumsum(tb))/(sum(tb)-cumsum(tb)+tb[]),
	"CSI"=cumprod((sum(tb)-cumsum(tb))/(sum(tb)-cumsum(tb)+tb[])))

## Plot survival curve	
plot(iC$Interval,iC$CSI,"s")

-----
A R learner.
-- 
View this message in context: http://r.789695.n4.nabble.com/sequential-treatment-of-a-vector-for-formula-tp2232517p2233491.html
Sent from the R help mailing list archive at Nabble.com.