Student's t-distribution is named in honor of William Sealy Gosset (1876-1937), who first determined it in 1908. Gosset was one of the best Oxford graduates in chemistry and mathematics in his generation. In 1899, he took up a job as a brewer at Arthur Guinness Son & Co, Ltd in Dublin, Ireland. Working for the Guinness brewery, he was interested in quality control based on small samples in various stages of the production process. Since Guinness prohibited its employees from publishing any papers to prevent disclosure of confidential information, Gosset had published his work under the pseudonym "Student", and his identity was not known for some time after the publication of his most famous achievements, so the distribution was named Student's or t-distribution, leaving his name less well known than his important results in statistics (Lovric 2011).
The t-distribution curve is, like the normal distribution curve, symmetric (bell shaped) about the mean. However, the t-distribution curve is flatter than the standard normal distribution curve. Consequently, the t-distribution curve has a lower height and a wider spread than the standard normal distribution.
The t-distribution has only one parameter, called the degrees of freedom $(df)$. The shape of a particular t-distribution curve depends on the number of degrees of freedom $(df)$. The number of degrees of freedom for a t-distribution is equal to the sample size minus one, that is,
$$df = n - 1$$As the sample size, $n$, and thus $df$ increases, the t-distribution approaches the standard normal distribution. The units of a t-distribution are denoted by t. The mean of the t-distribution is equal to $0$, and its standard deviation is $\sqrt{df/(df-2)}$ (Mann 2012).
x = np.arange(-4, 4, 0.1)
df = [1, 2, 8, 30]
color = ["red", "blue", "darkgreen", "gold"]
handle = ["df=1", "df=3", "df=8", "df=30"]
plt.figure(figsize=(10, 5))
plt.title("Comparison of t-Distributions", fontsize=16)
plt.plot(
x, stats.norm.pdf(x, 0, 1), color="black", linestyle="--", label="normal"
) ## add the normal distribution for comparison
for d, c, h in zip(df, color, handle):
y = stats.t.pdf(x, df=d, loc=0, scale=1)
plt.plot(x, y, color=c, label=h)
plt.axvline(
x=0,
color="black",
linestyle="dashed",
)
plt.text(-0.2, -0.05, ("$\mu = 0$"), fontsize=14)
plt.legend(fontsize=14)
plt.axis("off")
plt.show()
Citation
The E-Learning project SOGA-Py was developed at the Department of Earth Sciences by Annette Rudolph, Joachim Krois and Kai Hartmann. You can reach us via mail by soga[at]zedat.fu-berlin.de.
Please cite as follow: Rudolph, A., Krois, J., Hartmann, K. (2023): Statistics and Geodata Analysis using Python (SOGA-Py). Department of Earth Sciences, Freie Universitaet Berlin.