symulacje

##########################################################################################
###                                                                Oszacowanie macierzy przejścia – proces uregulowany                                                       ###
##########################################################################################

N=1000; m=100; sigma=5
x=rnorm(N,m,sigma)
y=matrix(NA,N,1)
z=(x-m)/sigma
for (i in 1:N)
{

if (abs(z[i])<1) y[i]=1 else if (abs(z[i])<2)  y[i]=2 else if (abs(z[i])<3) y[i]=3 else  y[i]=4

}
P=matrix(0,4,4)
for (i in 2:N)
{

k=y[i-1]
s=y[i]
P[k,s]=P[k,s]+1

}
P