# Does a best fit make sense in this case?

I have collected data from 60 candidates, 30 female and 30 male, measuring their pain tolerance, or rather time able to withstand pain. I have plotted the times against the $n$th candidate: Would it make sense to compute curves of best fit? My research focus is on comparing the tolerance of males and females. I do not think it would make sense for this graph since there is no relation between the candidate number and time.

• Are you really sure that the sampling was random? – Tim Dec 4 '14 at 13:49

The main thing this graph does is make me suspicious of your randomization.

If your goal is to compare male and female tolerance for pain graphically, you could use either overlaid density plots or parallel boxplots or given that you have 30 men and 30 women, parallel strip plots.

If you want to compare them statistically then a t-test is a place to start - this compares the means; but you might also want to compare quantiles, in which case quantile regression might be useful (although your N is a bit small).

• The t-test and LinReg will be done later on. Thank you for confirming my doubts, didn't think a best fit would work in this case but had to make sure. – Omer Cooper Dec 4 '14 at 14:08

(This answer takes an optimistic view that "randomised $n$th volunteer" is a strange way of referring to values sorted in magnitude.)

If you've ordered values by magnitude, as appears, this is a standard plot in statistical graphics often known as a quantile plot. Some would say it's a cumulative distribution plot with axes interchanged. Quantile plots rarely appear in introductory courses or texts, but in my view they are neglected; just as box plots and stem-and-leaf plots are often oversold, quantile plots are often undersold (as said, to the point of not appearing at all).

The horizontal scale is proportional to cumulative probability, i.e. fraction $\le$ the given value. Fitting a best fit curve could make perfect sense if you used a theoretical quantile function (a.k.a. inverse cumulative distribution function) of interest.

It's more common to assess fit with a particular distribution (e.g. normal, exponential, gamma) through a quantile-quantile plot of some kind.

In this case, given a focus on the difference between males and females, then a quantile-quantile plot comparing the two groups directly is indicated as a better way to address your research problem.

Here it would be an excellent way to show that the distributions are similar for low thresholds and different for higher thresholds, a fact which makes conventional focus on just comparing means and medians of limited value: the two distributions do not just differ additively.

We don't have your data (if you posted them, we could play), but here is a quantile-quantile plot for ozone concentrations at Yonkers and Stamford. For equal numbers in two groups, the recipe could hardly be simpler: plot

 lowest on group 1 vs lowest on group 2
2nd lowest on group 1 vs lowest on group 2
...
highest on group 1 vs highest on group 2


(The idea doesn't depend on groups being equally numerous: for groups with unequal numbers, we need to interpolate between values in the larger group to get values comparable in cumulative probability with those in the smaller group.) The data come from http://www.stat.purdue.edu/~wsc/visualizing.datatables/ozone In this particular case, they point to logarithmic transformation as the way forward.

More on quantile-quantile plots, with references: Benefits of using QQ-plots over histograms

No, I don't think it makes sense to fit a curve to this data. What would you use the best fit for? You can't really describe or predict anything meaningful from it.

It would make more sense to compare male/female pain tolerances using the mean or median of the two samples. I'd start exploring by making a histogram of the two data sets and see where that leads you.

• Thank you for confirming, I will not include a best fit. Histogram does sound like a much better choice here. – Omer Cooper Dec 4 '14 at 14:04