Sopra copia dei grafici ottenibili con i seguenti file JavaScript: primo, secondo, terzo.
Ecco come si possono ottenere grafici simili con R (vedi):
source("http://macosa.dima.unige.it/r.R") anno = 1900+c(12,14,24,33,34,36,37,41,53,56,57,60,60,60,61,61,61,62,62,63, 71,73,76,76,77,78,80,80,83,83,84,85,85,87,88,89,93) record = 200+c(0, 1, 3, 4, 6, 7, 9,11,12,15,16,17,18,22,23,24,25,26,27,28, 29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45) min(anno); max(anno); min(record); max(record) # 1912 1993 200 245 BF=6; HF=4 Plane(1910,1995, 195,245); AXES(1900,195,"black") gridHC( seq(200,250,5), "red"); gridVC( seq(1915,1995,5), "red") Point(anno,record, "blue") # Plane(1910,1995, 195,245); AXES(1900,195,"black") gridHC( seq(200,250,5), "red"); gridVC( seq(1915,1995,5), "red") polyl(anno,record, "brown"); Point(anno,record, "blue") # Plane(1910,1995, 195,245); AXES(1910,195,"black") gridHC( seq(200,250,5), "red"); gridVC( seq(1915,1995,5), "red") stepl(anno,record, "brown"); Point(anno,record, "blue") # # # ovvero (senza caricare la libreria "macosa.dima.unige.it/r.R" # e recuperando i dati da rete) readLines("http://macosa.dima.unige.it/om/prg/stf/alto1.txt",n=3) # "# Record di salto in alto maschile in funzione della loro data" # "1912,200" # "1914,201" # Ho capito che c'e' una riga di intestazione standard (che R automaticamente salta) # e che i dati sono separati da ",". Metto i dati in una tabella (hm) hm <- read.table("http://macosa.dima.unige.it/om/prg/stf/alto1.txt",sep=",") # Visualizzo gli estremi delle due variabili, che vengono etichettate V1 e V2 min(hm$V1); max(hm$V1); min(hm$V2); max(hm$V2) # 1912 1993 200 245 # Traccio il grafico per punti plot(hm$V1,hm$V2, ylim=c(200,245),col="blue",xlab="",ylab="",pch=20) abline(h=seq(200,245,5),v=seq(1910,2005,5),lty=3,col="grey50") dev.new() # Traccio il grafico congiungendo i punti plot(hm$V1,hm$V2, ylim=c(200,245),col="blue",xlab="",ylab="",pch=20) abline(h=seq(200,250,5),v=seq(1910,2005,5),lty=3,col="grey50") lines(hm$V1,hm$V2) dev.new() # Traccio il grafico e dei segmenti orizzontali che vissualizzano i periodi # in cui erano in vigore i record plot(hm$V1,hm$V2, ylim=c(200,245),col="blue",xlab="",ylab="",pch=20) abline(h=seq(200,250,5),v=seq(1910,2005,5),lty=3,col="grey50") for(i in 1:36) lines(c(hm$V1[i],hm$V1[i+1]),c(hm$V2[i],hm$V2[i]))