<font face="arial,helvetica,sans-serif">Hi there,<br><br>thanks for the references. I managed to solve the issue using the reference Paul sent.<br><br>I was trying to fit a density profile to a distribution of star. The usual way of doing it is to build a histogram of dN/dR (number of points/radial bins). But that introduces biases due to the binning choice. <br>

<br>A way to avoid binning is to use a likelihood approach. So I followed this paper 2008ApJ...684.1075M. <br><br>If anyone is interested in the code I wrote for that please ask!<br><br>Cheers and thanks again. <br clear="all">

</font><div>-- Eduardo Balbinot<br></div>
<br><br><div class="gmail_quote">On Sat, Dec 22, 2012 at 9:14 PM, Paul Barrett <span dir="ltr">&lt;<a href="mailto:pebarrett@gmail.com" target="_blank">pebarrett@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


Eduardo,<br>
<br>
What type of data are you trying to fit; time series, spectral or image?<br>
<br>
You might want to take a look at what the gamma-ray community has<br>
done, since almost all of their analysis involves photon lists.  You<br>
could also try searching the astrostatistics website.<br>
<br>
Personally, I would use a Bayesian approach.  I have successfully used<br>
the Nested Sampling approach to fit FUSE and Chandra spectra. This<br>
algorithm succeeded where the traditional Levenberg-Marquardt<br>
algorithm did not.  The algorithm is described in &quot;Data Analysis: A<br>
Bayesian Approach&quot; by D.S. Sivia.  One of the best introductions to<br>
the subject in my opinion.<br>
<br>
 -- Paul<br>
<div><div><br>
On Sat, Dec 22, 2012 at 1:22 PM, Eduardo Balbinot<br>
&lt;<a href="mailto:eduardo.balbinot@gmail.com" target="_blank">eduardo.balbinot@gmail.com</a>&gt; wrote:<br>
&gt; A bit off-topic, but....<br>
&gt;<br>
&gt; Anyone could point to a good guide to fitting unbinned data using a maximum<br>
&gt; likelihood estimator?<br>
&gt;<br>
&gt; It would be a nice idea to have it on astropy with most used density<br>
&gt; profiles (Hubble, Plummer, King, etc...). If I manage to do it we could<br>
&gt; discuss including it on astropy.<br>
&gt;<br>
&gt; Thanks!<br>
&gt; -- Eduardo Balbinot<br>
&gt;<br>
</div></div>&gt; _______________________________________________<br>
&gt; AstroPy mailing list<br>
&gt; <a href="mailto:AstroPy@scipy.org" target="_blank">AstroPy@scipy.org</a><br>
&gt; <a href="http://mail.scipy.org/mailman/listinfo/astropy" target="_blank">http://mail.scipy.org/mailman/listinfo/astropy</a><br>
&gt;<br>
</blockquote></div><br>