I would like to estimate the distribution of a very large population of known size but unknown mean and variance. I cannot assume anything about the underlying distribution. The values of observations are always discrete. I cannot sample the entire population but I would like to estimate the probability density function. I would also like some level of assurance of the correctness of the estimated distribution. What is the best way to go about this?

$\begingroup$ What do you mean by the statement "The values of observations are always discrete"? $\endgroup$ – Ashok Jun 20 '11 at 14:39

$\begingroup$ @Ashok: In this case, I am certain that they are nonnegative, nonzero integers. In general what I meant was that I believe that the distribution is lumped in some way rather than being continuously distributed over the entire range. $\endgroup$ – Misha Jun 20 '11 at 17:04
I can't comment so this will have to be in the form of an answer  my apologies for this.
First, I assume that what you want to do is to estimate the distribution of some property of the population, no? One way to go about it would be to construct the empirical distribution function: $\mathbb{F}(x) = \frac{1}{n}\sum _{i=1} ^{n} I(X_i \leq x) $ ($n$ the number of observations). Once you have done this, you can use various test to compare this to any distribution function that you'd like (i.e., one that you suspect produced the data).
If you are after a direct estimation of the density fucntion, and not any "known" distribution that it looks similar to, you can use a kernel density estimation: http://en.wikipedia.org/wiki/Kernel_density_estimation
There are tools in e.g. Matlab that will make this task fairly straightforward.
I would say that what method to use, and there are surely others, depends on that you want to do with the end result. I've only done this sort of thing in easy settings and there are probably more sophisticated ways to solve your problem. However, the above could perhaps set you off to a good start

$\begingroup$ Now that I read it again  is the problem to investigate the distribution of the population over a set of distinct "places" (such as nodes in a lattice or similar)? $\endgroup$ – Pierre Jun 20 '11 at 11:06

$\begingroup$ @Pierre Nyquist: It's not that I want to see the distribution over a set of distinct "places." You are right that I want to directly estimate the density function. I was actually wondering if Kernel density estimation might be a possibility but how do I know how large my sample size should be? If I use that approach, what can I say about the correctness of my estimation? $\endgroup$ – Misha Jun 20 '11 at 18:00