# 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