Как запустить svyquantile для подгрупп, которые имеют только значения n в R?

Question

Я пытался использовать следующую функцию для расчета предполагаемой медианы в г

myqt <- function (...) {

    a <- list( ... )

    this.qt <- if( is.null( a$quantile ) ) a[[3]] else a$quantile

    tt <- a[[2]]$variables

    # condition where svyquantile breaks
    if( sum( complete.cases( tt[ as.character( a[[1]] )[2] ] ) ) < 2 ){ 
        # re-create an empty svyquantile object of the same structure
        out <- matrix(rep(as.numeric(NA), length(this.qt)),nrow=1)
        dimnames(out) <- list( as.character(a[[1]])[2] , this.qt )

        # end the function
        return( out ) 

    # if svyquantile doesn't break, pass everything to it as normal
    } else svyquantile( ... )
}

Но вывод, который я получаю после использования этой функции, не имеет стандартных ошибок. Я хочу знать, как изменить эту функцию, чтобы я мог получить приблизительные значения медианы и ее стандартную ошибку. В случае медианы и стандартной ошибки в подгруппах, имеющих только значения n, следует возвращать NA.