13_29

(1) Come rappresentesti graficamente i seguenti dati:
Il comune XX è suddiviso in 4 circoscrizioni in cui vi sono, rispettivamente 23253, 41647 (scusate: avevo messo per sbaglio 416471 - ma l'esercizio si svolge allo stesso modo), 35829, 56936 abitanti.

(2) Carica i seguenti file via via tra un # e l'altro. Commenta, brevemente, i comandi dati, spiegando quanto "fanno". Completa la fase (F) e illustra, inoltre, le motivazioni delle fasi (A)-(F)

# (A)
readLines("http://macosa.dima.unige.it/R/pop_ita_2011.txt",n=4)
dati <- read.table("http://macosa.dima.unige.it/R/pop_ita_2011.txt",sep=",",header=TRUE)
dim(dati)
dati[c(1,15),1:6]
# (B)
Italia <- NULL; for (i in 1:15) Italia <- c(Italia, sum(dati[i,2:6]) )
dati2 <- cbind(dati,Italia)
dati2
# (C)
par( mai = c(0.5,0.5,0.1,0.1) )
plot(dati2$anno,dati2$Italia,ylim=c(0,60000) )
abline(h=axTicks(2),v=axTicks(1),lty=3)
lines(dati2$anno,dati2$Italia,col="blue")
lines(dati2$anno,dati2$I.N.Occ.,col="red")
lines(dati2$anno,dati2$I.N.Ori.,col="green")
lines(dati2$anno,dati2$I.Cen.,col="brown")
lines(dati2$anno,dati2$I.Mer.,col="magenta")
lines(dati2$anno,dati2$I.Ins.,col="orange")
# (D)
windows()
par( mai = c(0.5,0.5,0.1,0.1) )
plot(dati2$anno,dati2$Italia/dati2[1,7]*100)
abline(h=axTicks(2),v=axTicks(1),lty=3)
lines(dati2$anno,dati2$Italia/dati2[1,7]*100,col="blue")
lines(dati2$anno,dati2$I.N.Occ./dati2[1,2]*100,col="red")
lines(dati2$anno,dati2$I.N.Ori./dati2[1,3]*100,col="green")
lines(dati2$anno,dati2$I.Cen./dati2[1,4]*100,col="brown")
lines(dati2$anno,dati2$I.Mer./dati2[1,5]*100,col="magenta")
lines(dati2$anno,dati2$I.Ins./dati2[1,6]*100,col="orange")
# nel 1860 probabilmente i dati non tengono conto dell'intera
# popolazione del Veneto (che è stato in parte aggregato dopo)
# (E)
windows()
par( mai = c(0.5,0.5,0.1,0.1) )
plot(dati2$anno,dati2$Italia/dati2[9,7]*100,ylim=c(20,135) )
abline(h=axTicks(2),v=axTicks(1),lty=3)
lines(dati2$anno,dati2$Italia/dati2[9,7]*100,lty=2)
lines(dati2$anno,dati2$I.N.Occ./dati2[9,2]*100,col="red")
points(dati2$anno,dati2$I.N.Occ./dati2[9,2]*100,col="red",pch=20)
lines(dati2$anno,dati2$I.N.Ori./dati2[9,3]*100,col="green")
points(dati2$anno,dati2$I.N.Ori./dati2[9,3]*100,col="green",pch=20)
lines(dati2$anno,dati2$I.Cen./dati2[9,4]*100,col="brown")
points(dati2$anno,dati2$I.Cen./dati2[9,4]*100,col="brown",pch=20)
lines(dati2$anno,dati2$I.Mer./dati2[9,5]*100,col="magenta")
points(dati2$anno,dati2$I.Mer./dati2[9,5]*100,col="magenta",pch=20)
lines(dati2$anno,dati2$I.Ins./dati2[9,6]*100,col="blue")
points(dati2$anno,dati2$I.Ins./dati2[9,6]*100,col="blue",pch=20)
text(1975,70,"INOcc",col="red")
text(1925,50,"INOri",col="green")
text(1975,50,"ICen",col="brown")
text(1925,30,"IMer",col="magenta")
text(1975,30,"IIns",col="blue")
# (F)
readLines("http://macosa.dima.unige.it/R/pop_lig_2011.txt",n=3)
# Scrivi i comandi per aggiungere la visualizzazione dell'andamento
# dei dati relativi alla Liguria

(3) Calcola il valore (in forma frazionaria) della somma dei 4 valori che Cbin(n,n-2)/n assume per n pari a 25, 50, 100 e 200. Esplicita come hai svolto il calcolo.

[soluzioni nel forum]