Сделал следующий скрипт для проверки гипотезы нормальности распределения beta по методу Монте-Карло
Код:
nulldata 1000
set seed 1
fi=0.99
series b = 0
series y = 0
loop j=1..1000
genr u = 2*uniform()-1
y[1]=u[1]
loop i=2..120
y[i]=fi*y[i-1]+u[i]
endloop
sum1=0
loop i=2..120
sum1=sum1+y[i]*y[i-1]
endloop
sum2=0
loop i=2..120
sum2=sum2+(y[i-1])^2
endloop
b[j]=sum1/sum2
endloop
Добавлено через 17 минут
всё работает