Processing math: 100%

Exercise 4, Page 362

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 S25, and their average A25.

Normal Density Approximations

We are interested in approximating the distributions of the sum S25, the standardized sum S25, and the average A25 of 25 independently chosen random numbers from the interval [0, 20] using normal densities.

Sum S25

The sum S25 is the total of 25 randomly chosen numbers from a uniform distribution between 0 and 20.

Calculation:

  • Mean (μS25): The mean of a single uniform random number between 0 and 20 is 0+202=10. For 25 such independent numbers, μS25=25×10=250.
  • Variance (σ2S25): The variance for a single uniform number between 0 and 20 is (200)212. Thus, for 25 numbers, σ2S25=25×40012833.33.
  • Standard Deviation (σS25): 833.3328.87.

Therefore, the normal density approximating S25 is N(250,28.87).

Standardized Sum S25

The standardized sum S25 is the sum S25 minus its mean, divided by its standard deviation.

Calculation:

  • Mean (μS25): 0
  • Standard Deviation (σS25): 1

Hence, the normal density that approximates S25 is N(0,1).

Average A25

The average A25 is the mean of 25 uniformly distributed random numbers between 0 and 20.

Calculation:

  • Mean (μA25): The mean is the same as that of a single number, 10.
  • Variance (σ2A25): 40012×251.3333.
  • Standard Deviation (σA25): 1.33331.15.

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.

Sum S25

The sum S25 is the sum of 25 independent and identically distributed random variables from a uniform distribution over [0, 20].

Calculations

  • The expected value (mean) of a single uniform random variable U(0,20) is given by:

E[U(0,20)]=0+202=10

  • For the sum of 25 such random variables:

E[S25]=25×E[U(0,20)]=25×10=250

  • The variance of a single uniform random variable is:

Var(U(0,20))=(200)212=4001233.33

  • For the sum of 25 such random variables:

Var(S25)=25×Var(U(0,20))=25×33.33833.33

  • Therefore, the normal density approximating S25 is:

S25N(250,833.33)

Calculation of Standardized Sum S25

The standardized sum, S25, is obtained by subtracting the expected value (mean) and dividing by the standard deviation. The formula is given by:

S25=S25E[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.

Average A25

The average A25 is the sum S25 divided by 25.

Calculations

  • The expected value is:

E[A25]=E[S25]25=25025=10

  • The variance is:

Var(A25)=Var(S25)252=833.336251.33

  • So, the normal density approximating A25 is:

A25N(10,1.33)

Below is the r code to deduce these normal approximations for S25, S25, 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