Sopra copia dei grafici ottenibili con i seguenti file JavaScript primosecondoterzo.

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]))