I have large sample of data that is approximately from a Pareto distribution with unknown parameters. Unfortunately the distribution is sufficiently heavy tailed that just taking the sample mean is not a good estimate. What is a good way to estimate the mean and variance in this case?