# carica le seguenti righe: otterrai l'Italia dati1 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/penisola.txt",header=FALSE,sep =",") dati2 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/sicilia.txt",header=FALSE,sep =",") dati3 <- read.table("http://macosa.dima.unige.it/om/prg/R/rmacosa/sardegna.txt",header=FALSE,sep =",") mo <- min(dati1$V1,dati2$V1,dati3$V1);Mo <- max(dati1$V1,dati2$V1,dati3$V1) mv <- min(dati1$V2,dati2$V2,dati3$V2);Mv <- max(dati1$V2,dati2$V2,dati3$V2) plot(c(mo,Mo),c(mv,Mv),type="n",xlab="", ylab="") abline(h=seq(36,48,1),v=seq(6,19,1),lty=3,col="grey") lines(dati1,col="black") lines(dati2,col="brown") lines(dati3,col="blue") # # se vuoi evidenziare i punti puoi aggiungere: points(dati1,col="orange",pch=".", cex=3) points(dati2,col="orange",pch=".", cex=3) points(dati3,col="orange",pch=".", cex=3) # # se vuoi una scala che sia monometrica nel centro d'Italia usa il seguente plot (man mano # che ci si allontana dal centro la scala non è più corretta poiché la superficie della # terra è sferica e non è rappresentabile senza deformazioni su una superficie piana) plot(c(mo,Mo),c(mv,Mv),type="n",xlab="", ylab="", asp=1.31) abline(h=seq(36,48,1),v=seq(6,19,1),lty=3,col="grey") lines(dati1,col="black") lines(dati2,col="brown") lines(dati3,col="blue") #