Suppose we choose independently 25 numbers at random (uniform density) from the interval [0, 20]. We are interested in finding the normal densities that approximate the densities of their sum S25, their standardized sum S∗25, and their average A25.
We are interested in approximating the distributions of the sum S25, the standardized sum S∗25, and the average A25 of 25 independently chosen random numbers from the interval [0, 20] using normal densities.
The sum S25 is the total of 25 randomly chosen numbers from a uniform distribution between 0 and 20.
Therefore, the normal density approximating S25 is N(250,28.87).
The standardized sum S∗25 is the sum S25 minus its mean, divided by its standard deviation.
Hence, the normal density that approximates S∗25 is N(0,1).
The average A25 is the mean of 25 uniformly distributed random numbers between 0 and 20.
The normal density that approximates A25 is N(10,1.15).
By using the Central Limit Theorem, these normal densities provide a good approximation for the distributions of the sum, standardized sum, and average of the 25 randomly chosen numbers.
The sum S25 is the sum of 25 independent and identically distributed random variables from a uniform distribution over [0, 20].
E[U(0,20)]=0+202=10
E[S25]=25×E[U(0,20)]=25×10=250
Var(U(0,20))=(20−0)212=40012≈33.33
Var(S25)=25×Var(U(0,20))=25×33.33≈833.33
S25∼N(250,√833.33)
The standardized sum, S∗25, is obtained by subtracting the expected value (mean) and dividing by the standard deviation. The formula is given by:
S∗25=S25−E[S25]√Var(S25)
For S25, which is the sum of 25 independent random variables uniformly distributed from 0 to 20, we have:
Thus, when we standardize S25:
This is because the standardization process shifts the mean to 0 and scales the distribution to have a unit variance.
The average A25 is the sum S25 divided by 25.
E[A25]=E[S25]25=25025=10
Var(A25)=Var(S25)252=833.33625≈1.33
A25∼N(10,√1.33)
Below is the r code to deduce these normal approximations for S25, S∗25, and A25.
# Set seed for reproducibility
set.seed(123)
# Generate 25 random numbers from U(0, 20)
random_numbers <- runif(25, min = 0, max = 20)
# Sum S_25
S_25 <- sum(random_numbers)
# Standardized Sum S*_25
mean_S_25 <- mean(random_numbers) * 25
std_dev_S_25 <- sqrt(25 * (20^2 / 12))
S_star_25 <- (S_25 - mean_S_25) / std_dev_S_25
# Average A_25
A_25 <- mean(random_numbers)
# Display results
cat("Sum S_25:", S_25, "\n")
## Sum S_25: 297.7799
cat("Standardized Sum S*_25:", S_star_25, "\n")
## Standardized Sum S*_25: 0
cat("Average A_25:", A_25, "\n")
## Average A_25: 11.91119