Deriving the Cauchy density function
Let $X$ and $Y$ be independent standard Normal RVs then $\frac{X}{Y}$ has a Cauchy distribution. Note that if $T = \frac{X}{Y}$ then $1/T$ also has a Cauchy distribution by symmetry.
Deriving the pdf
Let's do this by finding the cdf and taking the derivative.
$$F_{\frac{X}{Y}}(t) = P \left( \frac{X}{Y} \le t \right) = P \left( \frac{X}{|Y|} \le t \right) $$where the absolute value can be added due to symmetry. We want the absolute value because $Y$ can be negative and we have to be careful when we multiply both sides because the inequality can switch. So
$$= P \left( X \le |Y|t \right) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^\infty e^{-y^2/2} \left( \int_{-\infty}^{|y|t} \frac{1}{\sqrt{2\pi}} e^{-x^2/2}dx \right) \, dy = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^\infty e^{-y^2/2} \, \Phi(|y|t) \, dy$$Note that $f(-y)=f(y)$ and so the function is symmetric and since we are now integrating over positive $y$ we can get rid of the absolute value. This gives us
$$= \frac{2}{\sqrt{2\pi}} \int_{0}^\infty e^{-y^2/2} \, \Phi(yt) \, dy$$The above is the cdf of a Cauchy RV so we can take the derivative with respect to $t$ to get the pdf.
$$f_{\frac{X}{Y}}(t) = \frac{d}{dt} F_{\frac{X}{Y}}(t) = \frac{2}{\sqrt{2\pi}} \int_{0}^\infty e^{-y^2/2} \, \frac{d}{dt} \Phi(yt) \, dy = \frac{2}{\sqrt{2\pi}} \int_{0}^\infty e^{-y^2/2} y f(yt) \, dy $$Above we interchanged the integral and $\frac{d}{dt}$ which you can do for "well behaved" functions (Leibniz integral rule). We also used chain rule for the derivative.
$$= \frac{2}{\sqrt{2\pi}} \frac{1}{\sqrt{2\pi}} \int_{0}^\infty e^{-y^2/2} y e^{-y^2t^2/2} \, dy = \frac{1}{\pi} \int_{0}^\infty e^{-y^2(1+t^2)/2} y \, dy $$To calculate the integral use substitution
$\quad u = y^2(1+t^2)/2$ $\quad du = y(1+t^2) \, dy$ $$= \frac{1}{\pi(1+t^2)} \int_{0}^\infty e^{-u} \, du = \frac{1}{\pi (1+t^2)} = f_{\frac{X}{Y}}(t)$$for all $t$ is the pdf of a Cauchy RV.
Articles
Personal notes I've written over the years.
- When does the Binomial become approximately Normal
- Gambler's ruin problem
- The t-distribution becomes Normal as n increases
- Marcus Aurelius on death
- Proof of the Central Limit Theorem
- Proof of the Strong Law of Large Numbers
- Deriving Multiple Linear Regression
- Safety stock formula derivation
- Derivation of the Normal Distribution
- Comparing means of Normal populations
- Concentrate like a Roman
- How to read a Regression summary in R
- Notes on Expected Value
- How to read an ANOVA summary in R
- The time I lost faith in Expected Value
- Notes on Weighted Linear Regression
- How information can update Conditional Probability
- Coupon collecting singeltons with equal probability
- Coupon collecting with n pulls and different probabilities
- Coupon collecting with different probabilities
- Coupon collecting with equal probability
- Adding Independent Normals Is Normal
- The value of fame during and after life
- Notes on the Beta Distribution
- Notes on the Gamma distribution
- Notes on Conditioning
- Notes on Independence
- A part of society
- Conditional Expectation and Prediction
- Notes on Covariance
- Deriving Simple Linear Regression
- Nature of the body
- Set Theory Basics
- Polynomial Regression
- The Negative Hyper Geometric RV
- Notes on the MVN
- Deriving the Cauchy density function
- Exponential and Geometric relationship
- Joint Distribution of Functions of RVs
- Order Statistics
- The Sample Mean and Sample Variance
- Probability that one RV is greater than another
- St Petersburg Paradox
- Drunk guy by a cliff
- The things that happen to us