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