Introduction to the package YPPE

library(YPPE)
#> Loading required package: survival

mle <- yppe(Surv(time, status)~trt, data=gastric, n_int = 10, approach = "mle", init = 0)
summary(mle)
#> Call:
#> yppe(formula = Surv(time, status) ~ trt, data = gastric, n_int = 10, 
#>     approach = "mle", init = 0)
#> 
#> Short-term coefficients:
#>     Estimate  StdErr z.value  p.value   
#> trt  1.77113 0.61843  2.8639 0.004185 **
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#> 
#> Long-term coefficients:
#>     Estimate   StdErr z.value   p.value    
#> trt -0.98230  0.29576 -3.3212 0.0008962 ***
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#> 
#> --- 
#> loglik = -582.5366   AIC = 1189.073
coef(mle)
#>     psi[1]     phi[1] 
#>  1.7711332 -0.9823035
vcov(mle)
#>            psi[1]      phi[1]
#> psi[1]  0.3824616 -0.10323293
#> phi[1] -0.1032329  0.08747657