Module 3: Probability and Statistics Tasks
Module 3: Probability and Statistics Tasks
The method of moments involves equating the moments of the distribution to the sample moments. For a uniform distribution f(x) = (a+1)x^a on (0, 1), the first theoretical moment E(X) is 1/(a+2). Equating this to the sample mean gives an equation in terms of a, which can be solved for a. Further moments could be used for additional accuracy in estimation .
Generate 50 values of peak flows for assumptions of normal, log-normal, and exponential distributions from provided data. Calculate appropriate statistics (e.g., mean, variance) for each generated set. Analyze differences in these statistics between models, noting how assumptions about distribution shape impact results. Assess fit by comparing generated statistics to observed flow data using goodness-of-fit tests .
Generate 50 observations using λ = 0.7. The theoretical mean of an exponential distribution is 1/λ = 1.428. Simulate the sample and calculate the sample mean. Compare this sample mean to the theoretical mean to assess the accuracy of the simulation. Deviations are expected due to random sampling effects .
To estimate the parameter α for an exponential distribution, we use the cumulative distribution function (CDF). For an exponential distribution, the CDF is given by F(x) = 1 - e^{-αx}. Given that 60% of the values are less than 2.5, F(2.5) = 0.6. Hence, we have 1 - e^{-α*2.5} = 0.6. Solving this equation for α gives α = -ln(0.4)/2.5 ≈ 0.366 .
Generate independent samples of varying sizes (10, 20, 50, 100) with known mean (5) and standard deviation (2). Calculate the sample means and standard deviations for each size. Compare these calculated statistics to the expected theoretical values (5 and 2). Sampling variability affects small samples. Analyze whether larger sample sizes give statistics closer to theoretical values due to the Central Limit Theorem .
Chebyshev’s inequality states that for any random variable with mean μ and standard deviation σ, the probability that the value deviates from the mean by more than kσ is at most 1/k². For a mean of 675 units and a standard deviation of 220 units, the deviation of 450 units means k = 450/220 = 2.045. The maximum probability that the peak flow deviates more than 450 units from the mean is therefore at most 1/(2.045^2) ≈ 0.239 .
To obtain the maximum likelihood estimators (MLE) for parameters a and b of the given PDF f(x) = (x-a)/b * e^{-(x-a)/b}, where a < x < ∞ and b > 0, we must first construct the likelihood function L based on n observations xi. L = Product ( (xi - a)/b * e^{-(xi-a)/b} ). Taking the log of the likelihood function and differentiating with respect to a and b, and then setting these derivatives to zero provides the MLEs for the parameters. Detailed calculations are required to obtain specific expressions for a and b. .
The regression equation is calculated using the least squares method, which minimizes the sum of squares of the residuals. The formula is ŷ = b₀ + b₁x, where b₁ = covariance(x, y) / variance(x) and b₀ = mean_y - b₁*mean_x. Using the given rainfall and runoff data, the regression coefficients are calculated resulting in a regression equation of ŷ = 12.97 + 0.947x .
Using the linear regression equation obtained from historical data relating discharge and recharge: discharge = b₀ + b₁*recharge. With given data, calculate the regression coefficients b₀ and b₁. For recharge of 13 units, substituting into the regression equation gives a predicted discharge value using these coefficients .
To obtain the correlation coefficient, use the formula: r = Σ[(x_i - x̄)(y_i - ȳ)] / [ sqrt(Σ(x_i - x̄)²) * sqrt(Σ(y_i - ȳ)²) ], where x̄ and ȳ are the means of the x and y samples, respectively. Plugging in the data of yearly rainfall and runoff provided, the computations yield r ≈ 0.856, indicating a strong positive correlation .