source("http://macosa.dima.unige.it/r.R")
#
# sistema di var.casuali A (colonna 1) e B (colonna 2)
# dist. da GE di al,at,bg,bs,co,cr,cu,mn,mi,no,pv,so,to,va,vc
# (distanza tra i capoluoghi di provincia)
# A: dist. in linea d'aria in km (arrotondata alle cinquine)
# B: dist. stradale in km (arrotondata alle cinquine)
AriaStrada=c(
55, 85,
70, 115,
155, 205,
160, 230,
155, 185,
115, 185,
105, 145,
165, 240,
110, 140,
115, 155,
85, 125,
165, 290,
110, 170,
155, 195,
105, 135)
# Conto i dati e poi li metto in due variabili distinte
n=length(AriaStrada); n
## [1] 30
Aria=NULL; Strada=NULL
# li metto in Aria e Strada, lunghe la metà (che prima ho inizializzato)
for(i in 1:(n/2)) {Aria[i]=AriaStrada[i*2-1]; Strada[i]=AriaStrada[i*2]}
Plane(0,max(Aria),0,max(Strada)); Point(0,0, "blue")
for( i in 1:(n/2) ) POINT(Aria[i],Strada[i],"brown")
p = mean(Aria*Strada)/mean(Aria^2); p
# 1.42046
f=function(x) p*x
graph1(f, 0,200, "blue")
#
# Potresti usare anche (vedi)
regression(Aria,Strada, 0,0)
## 1.42 * x
f1=function(x) 1.42*x
graph1(f1, 0,200, "blue")