I often have data, sample sizes typically around $N = 10000$, where I know the data follows a specific distribution -- either because it is simulated or I understand the physics processes that drive it.
What are the best methods of determining the distribution parameters?
The approach I used to take is to bin the data by some rule e.g. Freedman–Diaconis rule or the square-root-$N$ rule, and get the count number for each bin centre and then fit the function to these coordinates. This is very robust, and gives reasonable numbers but as I understand it, this isn't a good approach as this can result in information being "hidden" or lost -- and of course is sensitive to the number of bins and indeed the bin width. However this seems to be the most robust and reliable approach.
Other approaches I have tried include likelihood evaluations, often inbuilt into routines and functions in software such as Mathematica, give results but the more free parameters in a given distribution seems to make finding reasonable values more challenging.
What is the best approach to extract distribution parameters from a data set, when the distribution is known?