How long is a distribution considered normal?

Question

I have a dataset of metric distances ($n=5800$) and plotted those as a histogram. My initial thought was that this distribution looks normal. But after performing a Shapiro Wilk test and plotting the distribution like below I am not that sure anymore. In the Shapiro Wilk test the p value was way to low.

My question is: Can this distribution still be considered as a Gaussian distribution? The mean of my data is -0.027 and standard deviation is 0.72.

Also a qq plot doesn't look like a normal distribution to me. But I am also not sure how to define this distribution.

When limiting my data range to -1.4 and 1.4 the qq plot looks like the following. How can this be interpreted?

Edit: Did a new plot with different software this time. Now the plot doesn't look normally distributed to me anymore.

Answer 1

Your first plot is a good example that histograms or density estimates (if using too much smoothing) is a poor way of identifying distributions. I first made a histogram with many small bins, to my eyes that do not look normal (not shown). The details can be seen much clearer if showing the histogram on a log scale:

with the best-fitting normal density superposed. R Code is here:

    myhist <- hist(ydata, prob=TRUE, breaks="FD", plot=FALSE)   
    plot(myhist$mids, myhist$density, log="y", type="h", lwd=5, 
         lend=2, main="histogram: density on log scale")
    plot( function(x) dnorm(x,mean(ydata), sd(ydata)), 
          from=min(myhist$mids), to=max(myhist$mids), col="red", 
          lwd=3, add=TRUE)

Could be even more interesting if showing error bars on the histogram, see Confidence interval for the height of a histogram bar.

Another helpful tool for distribution identification is the Cullen and Frey graph implemented in fitdistrplus (R):

This is quite interesting. The big blue point is the data shown on a squared skewness versus kurtosis plane, the yellow points are bootstrapped versions. All points are way away from the parts of the plane with named distributions, none of them have high enough kurtosis. The suggestion in another answer of looking into t distribution could be good. R code for last plot:

    library(fitdistrplus)
    descdist(ydata, boot=50)  

Answer 2

Do you have any alternative candidate distributions? When fitting distributions, I've always considered multiple options, relevant to the studied system. This distribution $\textrm{ad oculum}$ looks, especially after looking at qq plots, like the Student's t to me, yet this is just an unverified observation.

Answer 3

Considering $n=5800$, which is a huge sample, I'd not accept those qqplots. Your second qqplot suggests your smaller observations are too small to be from a normal distribution, and your larger observations are too big to be from a normal distribution. Your tail probabilities are too heavy for a normal distribution.

I don't know what the distribution is. Your data is probably some sort of ratio of random variables. Perhaps the distribution is Cauchy? or t?