I have a set of about 100000 numbers. Fitting a gaussian to my data I can visually see that the points follow a gaussian almost exactly. Using the normplot I see that my data again follow a gaussian except for a little noise on the tails.
Now, what I am looking for is a function that can give me a p-value that rejects the null hypothesis that my data aren't normal. Lilleforfs, and Jbtest have the null hypothesis that the data are normal. I can reject the null hopotheses if I subsample my data down to about 100 points.
What I really want is to reject the hypothesis that my data are not normally distributed, with some associated p-value. Is this possible?
edit: my data are integers in the range of 1 to 100.
Probably should have kept my notes from 3rd year stats.