# Somewhere else, part 135          et un peu de lecture en français

• “Aux U.S., la vie est un actif financier comme les autres” http://lemonde.fr/economie … ht @econoclaste (une ‘découverte’ de produits qui datent)
• “Les universitaires ont intérêt à contribuer à Wikipédia — qui a intérêt à ce qu’ils contribuent…” http://framablog.org/…
• RT @coulmont Montre tes Saints ! le jour de la Saint-N il naît 1,35 fois plus de bébés prénommés N qu’habituellement Did I miss something?

# Modèles de Prévision, suite

La seconde partie de transparents pour le cours ACT6420 est en ligne. Il s’agit d’un document qui permettra à tout le monde garder des traces, mais je n’entends pas forcément les suivre… Par exemple, ce soir, je propose de jouer un peu avec les données d’incendies à Chicago

> chicago=read.table("http://freakonometrics.free.fr/chicago.txt",
> Y=chicago$Fire > X1=chicago$X_1
> X2=chicago$X_2 > X3=chicago$X_3
> base=data.frame(Y,X1,X2,X3)

ou une base sur les états américains (cf un ancien billet)

> US=read.table("http://freakonometrics.free.fr/US.txt", + header=TRUE,sep=";")

Sinon, la semaine passée, on avait évoqué la matrice chapeau (dont on reparlera davantage lorsque l’on évoquera la phase de diagnostique). L’idée était de noter que

avec

On peut visualiser cette matrice avec le code suivant

> X=cbind(1,cars\$speed)
> M=X %*% solve( t(X) %*% X) %*% t(X)
> Q=matrix(as.numeric(cut(M,c(-10,quantile(M,(1:5)/6),10),labels=1:6)),50,50)
> plot(1:50,1:50,col="white",axes=FALSE,xlab="",ylab="",ylim=c(50,0))
> library(RColorBrewer)
> couleur=rev(brewer.pal(6,"RdBu"))
> for(i in 1:50){for(j in 1:50){ points(i,j,col=couleur[Q[i,j]],pch=15,cex=1)}}
> axis(1)
> axis(2)

avec des cases rouges pour les valeurs importantes (et positive) et des case bleues pour des valeurs importantes (mais négatives) 