R/vectorRao.R

In saebnocov: Small Area Estimation using Empirical Bayes without Auxiliary Variable

#' @title Vector g in Newton Raphson Method by J.N.K.Rao
#'
#' @param alpha An alpha estimate value on iterating process
#' @param beta A beta estimate value on iterating process
#' @param ni The number of sample in each area
#' @param yi The number of "success" value in each area
#'
#' @return This function returns a value of vector g.
#' @export
#'
#' @import descr
#' @import dplyr
#' @importFrom stats aggregate
#' @importFrom stats weighted.mean
#' @importFrom rlang .data
#'
#' @examples
#' ## load dataset with no weight value
#' data(dataEB)
#' temp = pcapdir(dataEB[,-c(3)])
#'
#' ## estimates alpha and beta parameter
#' ## in EB estimate with Moment method by J.N.K.Rao
#' temp1 = alphabetaEB(data.dir = temp$direst ,pcap = temp$pcap,
#'                       method = "rao", opt = "moment",
#'                       maxiter = 100,tol = 0.00001)
#'
#' ##calculates vector g
#' vectorRao(alpha = temp1$alpha_cap, beta = temp1$beta_cap,
#'  ni = temp$direst$ni, yi = temp$direst$yi)
#'
vectorRao <- function(alpha,beta,ni,yi){
  g1 <- sapply(1:length(ni),function(x) sum(1/(alpha+0:yi[x]-1))-sum(1/(alpha+beta+0:ni[x]-1)))
  ni.yi <- ni-yi
  g2 <- sapply(1:length(ni),function(x) sum(1/(beta+0:ni.yi[x]-1))-sum(1/(alpha+beta+0:ni[x]-1)))
  return(cbind(g1=sum(g1),g2=sum(g2)))
}