R gamlss绘制平滑基线危险

在拟合下面的离散生存函数后，如何获得“平滑”基线风险图

burr<- function(xi) {
    linkfun <-function(mu) log(((1-mu)^(-xi)-1)/xi)
    linkinv <- function(eta) 1-(1+xi*exp(eta))^(-1/xi)
    mu.eta <- function(eta) exp(eta)*(1+xi*exp(eta))^(-1/xi-1) 
    valideta <- function(eta) TRUE
    link <- "burr"
    structure(list(linkfun = linkfun, linkinv = linkinv, mu.eta = mu.eta,
                   valideta = valideta, name = link), class = "link-gamlss") }

model00 <- gamlssNP(y ~ cs(FS$t,df=3)+v1+v2, 
             random=~1|id,family=BI(mu.link=burr(5)),  
             data=FS, K=1, mixture="gq", tol=0.5, 
             control = NP.control(EMcc=0.001, trace=FALSE), 
             g.control = gamlss.control(trace=FALSE))

---

Burr请您共享
dput（FS）
的输出，或者至少是
dput（FS[20，]）
的输出，好吗？在RHS中看到
数据
参数的列是很少见的。我也看不到审查机制。但我不是gamlss的用户。在任何情况下，您都应该在其中一个hte包中发布数据或引用示例数据集。请您共享
dput（FS）
的输出，或者至少是
dput（FS[20，]）
？在RHS中很少看到使用
$
访问的
数据
参数列。我也看不到审查机制。但我不是gamlss的用户。在任何情况下，您都应该在其中一个hte包中发布数据或引用示例数据集。