Description
Alex Wako
1a
Gaussian white noise are white noise that are IID normally distributed.
ρ(h) =
Gaussian white noise is stationary as all {wt} are normally distributed. Additionally ws and wt are independent when s is not equal to t since all Gaussian white noise are independent and identically distributed.
1b
set.seed(1)
# Let sigma squared equal to 10 w_100 = ts(rnorm(100, 0, 10)) sample_mean_100 = mean(w_100)
w¯ = 1.08887366914655
# Plotting the time series plot(w_100)
Time
This is the time series plot of the Gaussian white noise with n = 100, mean = 0, and a variance of 10.
# Plotting the acf of the time series acf(w_100, 20)
Series w_100
Lag
This is the auto-correlation function plot of the time series plot.
Theoretically, ρ(h) would equal 0 because all Gaussian white noises are independent, implying there are no correlation between wt and ws. Comparing ρ(h) and ρˆ(h), the graph for ρˆ(h) shows that aside from when there is no lag, the auto-correlation function is relatively close to 0. Because there is randomness in creating 100 samples of variables to create a normal distribution, there are points of lag where the auto-correlation function does not equal to 0.
1c
set.seed(2)
w_1000 = ts(rnorm(1000, 0, 10)) sample_mean_1000 = mean(w_1000)
w¯ = 0.619987355624838
# Plotting the time series plot(w_1000)
Time
This is the time series plot of the Gaussian white noise with n = 1000, mean = 0, and a variance = 10.
# Plotting the acf of the time series acf(w_1000, 20)
Series w_1000
Lag
This is the auto-correlation function plot of the time series plot.
Once again, the theoretical ρ(h) would equal 0 as the time series plot is still using Gaussian white noise.
Comparing ρ(h) to ρˆ(h), aside from when there is no lag, the acf is relatively close to 0, even closer than when n=100. Since there is still randomness in creating a normal distribution with n = 1000, there are points of lag where the acf is not equal to 0.
Comparing the ρˆ(h) when n = 100 and n = 1000, the range of acf is smaller for when n = 1000 as there are more more samples. The sample mean is also smaller. This is because of the Central Limit Theorem; with larger sample size, the distribution gets closer to its theoretical distribution (in this case with mean 0 and variance 10).
2
set.seed(3) w = rnorm(300, 0, 1)
s1 = c(10*cos((pi*1:100)/10), 5*cos((pi*101:200)/10), cos((pi*201:300)/10)) x1 = ts(s1 + w)
s2 = c(10*cos((pi*1:100)/10), 5*cos((pi*101:200)/5), cos((pi*201:300)/3)) x2 = ts(s2 + w)
plot(x1)
Time
The pattern in the first signal-plus-noise model has its amplitude decreasing as t increases in intervals of 100, but the frequency of oscillation of the noise stays the same throughout the time series plot. The noise in this plot removes the smoothness of the lines given by the signal and adds randomness (the spikes on the line).
plot(x2)
Time
The pattern in the second signal-plus-noise model has its amplitude decreasing as t increases in intervals of
100, and this time the frequency of oscillation of the noise increases in the same intervals. The noise in this plot adds the same features as the first plot but there is more randomness because the frequency of oscillation of the signal increases in intervals.
3
On paper
4
On paper
Reviews
There are no reviews yet.