www = "http://www.biostatisticaumg.it/dataset/raul.csv"
raul = read.csv(www, header = TRUE)
attach(raul)
tail(raul)

www = "http://www.biostatisticaumg.it/dataset/fresher.csv"
fresher = read.csv( www, header = TRUE )
attach(fresher)
tail(fresher)



relation1 = weight ~ height
model1 = lm(relation1)
summary(model1)



formula1 = Ldhtre ~ Ldhuno
modello1 = lm(formula1)
summary(modello1)


par(mfrow = c(1,2))
plot(relation1, main = "fresher, p < 0.001")
abline(model1)
plot(formula1, main = "raul, p < 0.001")
abline(modello1)