## 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 EB estimator
## and its MSE by naive method
estEBnaive(data.dir = temp$direst, pcap = temp$pcap, param = temp1)
Run the code above in your browser using DataLab