Timeline for Fitting data to an Normal Inverse Gaussian distribution
Current License: CC BY-SA 3.0
25 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Apr 13, 2017 at 12:55 | history | edited | CommunityBot | replaced http://mathematica.stackexchange.com/ with https://mathematica.stackexchange.com/ | |
| Jul 1, 2016 at 7:14 | vote | accept | Xavier_B | ||
| Jun 7, 2016 at 18:07 | answer | added | Frank | timeline score: 4 | |
| Feb 23, 2016 at 16:41 | comment | added | Xavier_B | @J.M. I put all the data in the csv file for anomalies in temperature in France over 30 years on a 4-week rolling basis. Numbers can b epositive or negative. They should allow to evaluate the best fitting NIG. | |
| Feb 23, 2016 at 16:36 | comment | added | Xavier_B | Excuse my ignorance, how do i put a csv file on pastebin? OK got it, link to pastebin is: pastebin.ca/3380053 | |
| Feb 23, 2016 at 15:59 | comment | added | J. M.'s missing motivation | I would say that's a bit too much. Why not reduce/downsample your data a bit so that it is more manageable as an example. In any case, if you can convert your spreadsheet to a *.csv file, you can put it on pastebin and link to it here. | |
| Feb 23, 2016 at 15:54 | comment | added | Xavier_B | @J.M. how can I send you an excel file with the 30 years' worth of weekly temperature anomalies in France? These anomalies are supposed to follow a NIG distribution. that 1675 numbers. | |
| Feb 22, 2016 at 3:07 | comment | added | J. M.'s missing motivation | Do you have any sample data that can be used by potential answerers? | |
| Feb 22, 2016 at 3:06 | history | edited | J. M.'s missing motivation | CC BY-SA 3.0 | added 20 characters in body |
| Dec 25, 2015 at 20:28 | comment | added | Michael E2 | I think they're all probably closely related: DistributionFitTest, EstimatedDistribution, FindDistributionParameters; maybe FindDistribution. I can't think of anything else. Maybe someone with a better statistics will come along. Things can get a little slow on SE this time of the year. | |
| Dec 25, 2015 at 20:15 | comment | added | Xavier_B | @MichaelE2 You are actually correct: the NIG is from the hyperbolic family of functions, and in the case I am using, is a PDF. As you mentioned it could be used as one of the Hyperparabolic pdf for DistributionFitTest purposes, I tried it, but, as I said got parameters for which the hyperbolic function does not match the data. But maybe, I am not using the correct function to test whether a function describes a dataset. Do you know of some other Mathematica test I could use? | |
| Dec 25, 2015 at 15:56 | comment | added | Michael E2 | Of course, I should have mentioned: The "constant" means that to get the functions to be identical, the parameters or a function would have to be adjusted. Unless Mathematica is wrong, it means NIG2 is not a PDF, except for values of parameters that make the constant equal to 1; since you don't claim it is a PDF, maybe that's unimportant. | |
| Dec 25, 2015 at 15:49 | comment | added | Michael E2 | I meant PDF[HyperbolicDistribution[3/2, α/δ, β/δ, δ, μ], x]/NIG2[x, {α, β, μ, δ}] // FullSimplify[#, δ > 0]& simplifies to a "constant" depending only on α and β. Since I've no data to test and I'm not that familiar with your problem, I can't really say much more. I thought the built-in distribution might be used more easily than a user-constructed one. One thing I don't understand is that the linked question deals with time-series data & processes, but I don't see the connection with your problem, at least in the details presented. (I'm only vaguely familiar with time series.) | |
| Dec 25, 2015 at 9:25 | comment | added | Xavier_B | @Michael E2: I tried plotting the PDF[HyperbolicDistribution] using the parameters you mention, and the values for them as defined above, but the functions don't overlap: they merely have approximately the same mean and they don't fit the data. Are you sure about the parameter definition? When I use the HyperbolicDistribution as the distribution definition for the DistributionFitTest, I get the following parameters : {3/2,7.71259,-1.04967,18.447,2.38682}, and the tails for this distribution are too thin compared to the data (I don't know how to insert a picture of the graph). | |
| Dec 25, 2015 at 9:05 | comment | added | Xavier_B | @Rom38 There are a number of fittedDistribution tests that I can use in Mathematica, but only with PDF known to Mathematica. The NIG is not one of them. If the NIG is the correct PDF to use, then, the parameters can been estimated using the literature. However, I now wish to make sure that this NIG is effectively the a good fit for the data. | |
| Dec 24, 2015 at 20:12 | comment | added | Michael E2 | I don't know if it's any help, but your NIG2 is proportional to PDF[HyperbolicDistribution[3/2, α/δ, β/δ, δ, μ], x] (as a function of x). | |
| Dec 24, 2015 at 15:21 | comment | added | george2079 | Your calculation has a complex result with an imaginary part that is pretty small, but not so small to just use Chop. This is likely the source of your error. Is it supposed to be real? Do you just want to take the real part? Another thing, you might set K1[w_] = 2 BesselK[1, w] and verify the conditional is satisfied after you get your fit. | |
| Dec 24, 2015 at 15:10 | history | edited | george2079 | CC BY-SA 3.0 | deleted 35 characters in body |
| Dec 24, 2015 at 6:27 | comment | added | Rom38 | I did't understand your question. What do you want to do with this function? Just use it to fit your data with calculation of the optimal alfa, beta etc? | |
| Dec 24, 2015 at 2:56 | history | edited | Xavier_B | CC BY-SA 3.0 | added 279 characters in body |
| Dec 23, 2015 at 18:59 | comment | added | george2079 | Can you edit the question to add the definition of K1 (at a glance I don't see a "third kind" Bessel function directly available..) | |
| Dec 23, 2015 at 18:39 | comment | added | Xavier_B | Sorry, forgot to say that K1 is the Bessel Function of the third kind: K1[w_] = 1/2 \!( *SubsuperscriptBox[([Integral]), (0), ([Infinity])](Exp[(- *FractionBox[(w ((t + t^(((-1)))))), (2)])] [DifferentialD]t)); | |
| Dec 23, 2015 at 17:08 | history | edited | m_goldberg | CC BY-SA 3.0 | Improved formatting |
| Dec 23, 2015 at 15:44 | comment | added | george2079 | whats K1? (15) | |
| Dec 23, 2015 at 15:24 | history | asked | Xavier_B | CC BY-SA 3.0 |