prostoc.r<- function(n.oss,n.traiet) {# programma per generazione di processi proc.stoc<-matrix(0,n.oss,1) for (i in 1:n.traiet) {set.seed(i+100) proc.stoc<-as.matrix(cbind(proc.stoc,rnorm(n.oss))) } ts.plot(proc.stoc[,2:ncol(proc.stoc)],lty=seq(1,n.traiet,1),col=seq(1,n.traiet,1)) } prostoc.r(100,3) arma.sim<- function(n.oss) {# Simulazione di processi ARMA # MA(1) sim.ma1<-arima.sim(list(ma=c(-0.6)),n=n.oss) ts.plot(sim.ma1,ylab="MA(1) - theta=-0.6") windows() par(mfrow=c(2,1)) acf(sim.ma1) pacf(sim.ma1) par(mfrow=c(1,1)) # AR(1) sim.ar1<-arima.sim(list(ar=c(0.6)),n=n.oss) windows() ts.plot(sim.ar1,ylab="AR(1) - phi=0.6") windows() par(mfrow=c(2,1)) acf(sim.ar1) pacf(sim.ar1) par(mfrow=c(1,1)) } arma.sim(1000)