dati<-read.table("pesmol.txt",header=T)
names(dati)
dim(dati)
summary(dati)


res<-lm(formula=lunghezza~peso,data=dati)
summary(res)



plot(dati$peso,dati$lunghezza,type="b")
abline(reg=res,col="red")
title("regressione  peso-lunghezza molle")

plot(dati$peso,res$residuals,type="b")
abline(a=0,b=0,col="red")

title("residui peso-lunghezza molle")


confint(res,level=0.99)
pr<-predict(res,level=0.99,
        interval="confidence")
cbind(dati,pr)


matplot(dati$peso,pr,type="l",lty="solid")
title("IC per valori attesi")

#<-rbind(dati,c(7,NaN))
#prnew<-predict(res,datinew,level=0.99,
  #          interval="confidence")
#cbind(datinew,prnew)
dati$peso2<-dati$peso^2
res2<-lm(formula=lunghezza~peso+peso2,data=dati)
summary(res2)

anova(res,res2)