# analizza le sequenze causali

# function pausa
pausa <- function() {
  cat ("\n >> Guarda il plot e dai enter per continuare\n")
  scan()
}

# leggi le sequenze (il file 10sequenze.dat e' in formato R,
# quindi basta eseguirlo e la variabile s e' importata nello script)
source("10sequenze.R")

str(s)

nr.series  <- length(s[1,])   # nr di serie
len.series <- length(s[,1])   # lunghezza di ciascuna serie

# media e sigma di ciascuna serie
media=sigma=rep(0,nr.series)
cat(" medie e sigma: \n")
for (i in 1:nr.series){
 media[i] = mean(s[,i])
 sigma[i] = sd(s[,i])
 cat( sprintf("serie %d:   %f  %f \n", i, media[i], sigma[i] ) )
}

layout(matrix(1:5, 5,1)) # splittiamo lo schermo
mar.orig = par('mar')    # salviamo i vecchi margini per definirne dei nuovi 
par(mar=c(2.5,2,1.5,1))  # b, l, t, r 
for (i in 1:5) plot(s[,i], ty='l')
pausa()
for (i in 6:10)  plot(s[,i], ty='l')
pausa()
# zoom sui primi 200 elementi
for (i in 1:5)  plot(s[1:200,i], ty='l')
pausa()
for (i in 6:10)  plot(s[1:200,i], ty='l')

par(mar=mar.orig) # resettiamo i vecchi margini
layout(1)         # rimettiamo a posto lo schermo