How to fit histogram data with normfit using signed X-axis ?

2 次查看（过去 30 天）
Octopode
Octopode 2013-1-16
Hi,
I'm puzzled with normfit.
I have intensity of velocity data on a given X-direction, converted as a 1-D vector Vx. It has a size of S= 256 X 256 elements (65536) elements. They are signed values, of the order of 10^-2.
No problem with dfittool. I get sensible values for mean and standard deviation. But i need to automate the process, for batch processing. No GUI. Thus, thought i used normfit to fit a PDF. 1) Is this the right way to go ? Is normfit designed to fit a PDF with signed bins values ?
I thought i'd output the histogram data and built the PDF from there : So i did :
[nVx,Vxbins]=hist(Vx,nbins)
Then to obtain pdf i did :
pdf = nVx ./ (S*diff(Vxbins(1:2)))
A posteriori, i checked that actually sum(nVx) = S = 65336 (makes sense to me)
From there, not sure i should normalize to the max... does normfit require that ? My max value in the pdf is about 15.
From there, i do
[mu , sigma, blah...] = normfit(pdf)
Result is crap, i get something completely different from dfittool.
What am i missing ? Is normfit not designed for this situation ? If not, what is the right way to go ?
Thanks

请先登录，再回答此问题。

采纳的回答

Sarah Wait Zaranek
Sarah Wait Zaranek 2013-1-16
Probably the easiest way to handle this, if you like what the distribution fitting tool is doing, is to autogenerate code from the tool itself. You can do that via File -> Generate Code. That will allow you to do batch processing.
  1 个评论
Octopode
Octopode 2013-1-16
I wouldn't mind doing that on last resort. However, doing the "generate code" is to be done manually once, and applies only to my code. Batch processing would not be possible if the code/context changes as i would need to "re-generate" the diftool-related code. If i read whatever code is generated, i eventually will find what exactly dfittool is doing but normfit is meant to avoid that pain I guess...

请先登录，再进行评论。

更多回答（2 个）

owr
owr 2013-1-16
You should be calling normfit on the raw data vector itself, not your pdf approx:
[mu , sigma, blah...] = normfit(Vx)
  1 个评论
Octopode
Octopode 2013-1-16
Ok !! I thought i already tried that at first before going for the fit of the pdf. Obviously i didn't try properly. It works now, using Vx instead of the histogram output. That was really simple indeed, i feel a bit stupid. Thanks

请先登录，再进行评论。

Octopode
Octopode 2013-1-16
编辑：Octopode 2013-1-16
Thank you for your answer. I just tried generating the code. I found out that they are using fitdist in a way i didn't use before. They are not using directly normfit. So i assume normfit is more low level and assumes X to be positive (while my data are signed), although the help doesn't say anything about that. So, using fitdist is much easier and answers my question, although the use of normfit stays really obscure. (concerning the domain for X...)

类别

在 Help Center 和 File Exchange 中查找有关 Half-Normal Distribution 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by