I have a dataset containing information about a bunch of wireless devices. Specifically, the number of other wireless devices each device encounters. Based on other researcher's prior work in the same area, I have reason to believe my distribution may potentially be well fitted by the not so well known biPareto distribution, whose CCDF is defined as:
$1 - F(x) = (x/k)^{-\alpha}(\frac{x+kb}{k + kb})^{\alpha-\beta}$, for $ x > k$
and simply
$1 - F(x) = 1$, for $x \le k$.
Assume I already know $k = 1$. This leaves a function to fit with the three parameters $\alpha$,$b$,$\beta$, with the restriction that each parameter is $> 0$.
What I am unsure about is how one would go about determining the most suitable fitting parameter values, ideally in R.
I am quite the novice when it comes to fitting data, so aside from the syntactical considerations in R, I am also interested in any advice regarding what specific methods of fitting would be most appropriate for this particular function. For example:
- Does this function qualify for linear regression?
- Is least squares appropriate/applicable here? What about minimum mean squared error (MMSE)?
- How sensitive is any given fitting method to any "guiding" initial parameter values I give it (iterative methods?), or is the result determined analytically (always produces the same answer in the end?)
If it's at all relevant, I'm ultimately looking to use the Kolmogorov-Smirnov test on the data to see how well it matches the parameterized biPareto distribution.
