library(Bhat) #Control list for Plkhci, with a range (low,upp) on parameters plist=list(label=parnames,est=MLE, low=c(0,0,0,0,0),upp=c(1,100,10,100,10)) #Create a 5 by 2 matrix to contain the LRCIs LRCI.table=matrix(0,nrow=5,ncol=2, dimnames=list(parnames,c("LowerBound","UpperBound"))) source("Plkhci.R") #Defines function Plkhci for(i in 1:5) LRCI.table[i,]=Plkhci(plist,nllhood,parnames[i],prob=0.95, eps=0.000001,y=waiting) round(LRCI.table,4) #Print to 4 decimal places