0
$\begingroup$

I have a sequence of points which was obtained from an iterative algorithm, and I computed the order of convergence $p$ of the method using the formula

$$ p \approx \frac{\log({\rm err}(k+2))-\log({\rm err}(k+1))}{\log({\rm err}(k+1)) - \log({\rm err}(k))}. $$

The first computed ratio is much larger compared to the other ones. Without removing this value from the complete list of values, I obtained $p = 1.2615$ as the order of convergence of my algorithm. However, when I computed the order of convergence after removing the outlier (the first calculated ratio), I got $p = 1.0495$ (see figure).

My question is, which of the two is the correct order of convergence of my algorithm?

|cite|improve this question
$\endgroup$
  • $\begingroup$ I'd say the first one is an outlier. Is this data obtained from a one-time run or is each point an average of multiple measurements? In this case I would actually go for non-linear least-squares as a method of fitting to prevent outliers from polluting the results too much. $\endgroup$ – Nox Jan 28 at 14:08
  • $\begingroup$ Why not using a least squares approach to find the slope of your line? $\endgroup$ – nicoguaro♦ Jan 28 at 14:47
  • $\begingroup$ @Nox The data is a result of a one-time run and the above lines and slopes were obtained via the polyfit function in Matlab. So, which of the two should I take as the order of convergence for my algorithm? $\endgroup$ – Julienne Franz Jan 29 at 4:58
  • $\begingroup$ @nicoguaro The above lines and slopes were actually obtained via the polyfit function in Matlab. $\endgroup$ – Julienne Franz Jan 29 at 4:59
  • 1
    $\begingroup$ If possible re-run and check if the outlier is still there. Honestly I think it's an outlier and can be ignored, making the 1.0495 value more likely to be correct. $\endgroup$ – Nox Jan 29 at 14:02

Your Answer

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Browse other questions tagged or ask your own question.