# Carichiamo dei dati da analizzare da rete. Vediamo prima il file:
readLines("http://macosa.dima.unige.it/om/prg/R/polar1.txt",n=4)
# Ottengo:
# # temperature mese per mese di due anni"
# "41"                                     
# "60"                                     
# "77"
# Devo saltare 1 riga, in cui è contenuto un commento (skip=1)
dati <- scan("http://macosa.dima.unige.it/om/prg/R/polar1.txt",skip=1)
# Ho caricato i dati nella varibile "dati" e ho ottenuto l'informazione
# che essi sono 24:   Read 24 items
max(dati)
# Ottengo: 99. Questo è il valore massimo
# Copio "a scatola nera" questa riga
rad <- function(x) x*pi/180; A <- seq(0,330,30)*pi/180
# Scelgo un intervallo per le x ed uno per le y che si allontani
# dall'origine al più quanto il valore massimo
plot(c(-110,110),c(-110,110),type="n",xlab="",ylab="",axes=FALSE,asp=1)
# Traccio dei cerchi (in questo caso 4, con raggi 25,..., 100
for (i in 1:4) symbols(0,0,circles=i*25, inches=FALSE, add=TRUE, fg="grey60")
# Traccio 12 "raggi"
for(i in 0:6) abline(0,tan(rad(30*i)),lty=3)
# Traccio la "curva" che rappresenta i dati
lines(cos(A)*dati,sin(A)*dati)
# Evidenzio i dati con dei punti, diversamente colorati
points(cos(A)*dati[1:12],sin(A)*dati[1:12],col="red")
points(cos(A)*dati[13:24],sin(A)*dati[13:24],col="blue")
# Metto due scritte
ang <- rad(192); r <- 100; text(cos(ang)*r,sin(ang)*r,"100")
ang <- rad(6); r <- 70; text(cos(ang)*r,sin(ang)*r,"gen")
# Risultato finale: