R/S02.fun.R

In TransModel: Fit Linear Transformation Models for Right Censored Data

S02.fun <-
function(ord.bz,Rt,r,n,ord.wt){
    Rt0 = c(0,Rt[1:(n-1)])
    s0 = numeric(n)
    
    for(i in 1:n){
          a = ord.wt*exp(-ord.bz)/(1+r*Rt0[i]*exp(-ord.bz))
          s0[i] = sum(a[i:n])
    }
    return(s0)  
}