1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
|
lb1<-"[A terme]"
op <- par(mfrow = c(2, 2))# 2 x 2 images
plot(terme$gsp,terme$delfia,col="blue",main=paste(lb1,lb,", N=",
length(terme$agest)))
abline(fitterme,col="red") (1ere case)
lm_tout <- round(coef(fitterme), 3) # coefficient
mtext(bquote(y == .(lm_tout[2])*x + .(lm_tout[1])),adj=0)
x<-format(summary(fitterme)$r.squared,digits=3,width=4)
mtext(bquote(R2 == .(x)),adj=1,padj=-0.2)
rm(x,lm_tout,lb1)
hist(log(terme$delfia), breaks = 15, probability=T,right=F,col="lightblue",
border="blue", main="Delfia (log)")2eme case
points(density(log(terme$delfia), bw=.3), type='l', col="red", lwd=3)
hist(log(terme$gsp), breaks = 15, probability=T,right=F, col="lightblue",
border="blue", main="Gsp (log)")3eme case
points(density(log(terme$gsp), bw=.3), type='l', col='red', lwd=3)
<--- ici !
par(op) |
Partager