# Profile likelihood of N in binomial model

Let $(n_1,...,n_k)$ be a sample from a Binomial$(N,p)$ where both parameters are unknown. In many cases, the profile likelihood of $N$ is asymptotic in the sense that it never decays to $0$. An example of this is given here. The following code shows the profile likelihood of $N$ for this data set.

data = c(16, 18, 22, 25, 27)
prof = function(n) sum(dbinom(data,n,mean(data)/n,log=TRUE))
n = points = seq(25,500,1)
for(i in 1:length(n)) points[i] = exp(prof(n[i]))
m = max(points)
plot(n,points/m)


What can we do from a classical viewpoint in such scenario?

Shall we just say that a likelihood-confidence interval is $(N_0,\infty)$?

Note: Using the Bayesian approach, the posterior of $N$ is going to decay to $0$ for any prior, due to its integrability. The issue here is that the inferences are strongly influenced by the prior given that, in vague words, the likelihood has little information about $N$.