I looked into the code for the function and it doesn't seem to include that in the output unless a particular criterion (k<1
) is satisfied, where k
is defined within the function. That's not what the documentation says it should do. Your options are probably either to put a modified version of the function in your own workspace; what you need to do is to replace the last two lines
print(rval)
return(Confidence.Intervals)
with just
return(rval)
You won't get the confidence intervals you were getting, but you could probably add those to the output list if you wanted as well, with
rval$Confidence.Intervals <- Confidence.Intervals
return(rval)
It might also be worthwhile to contact the package maintainer, though it doesn't look like it's been updated since March 2012, for whatever that's worth.