There is a function to find the likelihood of given GARCH(1,1) model parameters where sigma2 is expressed as the following:
```
for t in range(1, len(returns)):
sigma2 = gamma * parkinson[0] + alpha * returns[t-1]**2 + beta * parkinson[t]
```
Why isn't the third term parkinson[t-1]?