this is for holding javascript data
Mikhail Kuzmin Deleted File
almost 8 years ago
Commit id: 8f8934d5b643d47e82a149ab6749e5e2e5e94afd
deletions | additions
diff --git a/begin_lstlisting_language_R_N__.tex b/begin_lstlisting_language_R_N__.tex
deleted file mode 100644
index c840656..0000000
--- a/begin_lstlisting_language_R_N__.tex
+++ /dev/null
...
\begin{lstlisting}[language=R]
N <- 100
a<-0.5
b<-0.95
sigmaPsi<-1
sigmaEta<-50
k<-1:N
x<-k
z<-k
x[1]<-0
z[1]<- x[1]+ rnorm(1,mean = 1, sd = sigmaEta)
for (t in 1:(N-1)) {
x[t+1] <- a*x[t] + t + rnorm(1, mean =1 , sd = sigmaPsi)
z[t+1] <- b*x[t+1] + rnorm(1, mean = 1, sd = sigmaEta)
}
xOpt<-k
eOpt<-k
K<-k
y<-k
xOpt[1]<-z[1]
eOpt[1]=sigmaEta
for (t in 1:(N-1)) {
eOpt[t+1]=sqrt((sigmaEta*sigmaEta)*(eOpt[t]*eOpt[t]+sigmaPsi*sigmaPsi)/(sigmaEta*sigmaEta+eOpt[t]*eOpt[t] +sigmaPsi*sigmaPsi))
K[t+1]=eOpt[t+1]*eOpt[t+1]/(sigmaEta*sigmaEta)
xOpt[t+1]=(a*xOpt[t]+t)*(1-K[t+1])+K[t+1]*y[t+1]
}
plot(k,xOpt,type ="l",col="red")
lines(k,xOpt,col="red")
lines(k,z,col="green")
lines(k,x,col="blue")
\end{lstlisting}