## 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)
Run the code above in your browser using DataLab