Heaviside機能を使用するにはどうしたらいいですか？

Question

は、私は非移植患者対移植の間に生存を分析するR.

library(survival) 
head(jasa) 

にスタンフォード心臓移植データを処理しています。しかし、私はHeaviside関数を使いたいです。

> transplant(t) : 0 if did not receive transplant by time, t i.e, if t<wait-time 
    1 if received transplant prior to time, i.e, 
    if t> wait-time 
> 
> msscore(t) 0 if t< wait-time 
      mscore if t>wait-time 

> age(t)  0 if t< wait-time 
      AGE if t>wait time 

これらのHeaviside関数を使用して生存プロセスを評価したいと思います。私は多くの文書を検索しましたが、私が知る限り、RのHeaviside関数を作ることはできません。 fBasicsでifelse(t>wait-time,1,0)とヘビサイド機能に加えて

Answer 1

おそらくlibrary(sos); findFn("Heaviside")で見つけた（ `ifelse（）`で

H <- function(x) as.numeric(x>0) 
mscore*H(t-(wait-time))