{"word_count": 632, "character_count": 6417}
En cours d'actuariat, nous commencerons à parler des tables de mortalité prospectives.Côté théorie, je peux renvoyer vers quelques slides de Pierre Devolder, ici. Je renvoie ici pour un court papier de vulgarisation sur le sujet. Sur l'ajustement de lois de mortalité prospectives, je renvoie au très beau mémoire en ligne ici (ou là pour des choses très proches en anglais). Pour l'approche statique, je renvoie également ici pour une méthode simple d'ajustement de loi de Makeham. Pour les aspects dynamique, il y a des documents techniques ici ou là.Sinon pour les données que nous utiliserons en TD, elles sont ici (pour l'exposition, i.e. le nombre de vivants) et là (pour les décès). Il faut retravailler un peu les données, à cause d'une valeur "110+", mais rien de plus simple,&#62; DEATH=read.table("https://perso.univ-rennes1.fr/arthur.charpentier/R/FR-Deaths_1x1.txt",header=TRUE)&#62; EXPOSURE=read.table("https://perso.univ-rennes1.fr/arthur.charpentier/R/FR-Exposures_1x1.txt",header=TRUE)&#62; DEATH$Age=as.numeric(as.character(DEATH$Age))Message d'avis :NAs introduits lors de la conversion automatique &#62; DEATH$Age[is.na(DEATH$Age)]=110&#62; EXPOSURE$Age=as.numeric(as.character(EXPOSURE$Age))Message d'avis :NAs introduits lors de la conversion automatique &#62; EXPOSURE$Age[is.na(EXPOSURE$Age)]=110Coté code, je renvoie à un précédant billet sur le sujet, ici, et là pour le package (non officiel) de démographie développé par Rob Hyndman.Enfin, à la fin du cours, on parlera un peu de modèles de durée. J'essayerai de taper quelques choses dans les jours à venir, mais côté TD, on utilisera ce que fait Frédéric Planchet (qui connaît fort bien le sujet), et qui se trouve en ligne ici. Les données sont en ligne ici, avec un peu de code R pour retraiter les donnéesici et là, et un peu de code pour des modèles classiques, comme du modèle de Cox (ici) ou du Kaplan Meier (là).Les codes R utilisés en TD sont les suivants (avec en plus les graphiques obtenus),&#62; D=DEATH[DEATH$Year==an,]&#62; E=EXPOSURE[EXPOSURE$Year==an,]&#62; MU = D[,3:5]/E[,3:5]&#62; plot(0:110,log(MU[,1]),type="l",col="red")&#62; lines(0:110,log(MU[,2]),col="blue")Avec ci-dessus les log des taux de décès instantannés pour l'année 1986.&#62; PH=matrix(NA,111,111) # ligne X, colonne H&#62; PF=PH&#62; for(x in 0:110){+ PH[x+1,1:(111-x)]=exp(-cumsum(MU[(x+1):111,2]))+ PF[x+1,1:(111-x)]=exp(-cumsum(MU[(x+1):111,1]))+ }&#62; x=0&#62; plot(1:111,PH[x+1,],ylim=c(0,1),type="l",col="blue")&#62; lines(1:111,PF[x+1,],col="red")pour la fonction de survie à la naissance (ci-dessus) ou à 50 ans (ci-dessous),&#62; somme=function(x){sum(x,na.rm=TRUE)}&#62; EH=EF=rep(NA,111)&#62; EH=apply(PH,1,somme)&#62; EF=apply(PF,1,somme)&#62; plot(0:110,EH,type="l",col="blue")&#62; lines(0:110,EF,col="red")Ce qui correspond aux espérances de vie résiduelles.&#62; for(k in 1:3){+ Z=MU[,k]+ Z[is.nan(Z)==TRUE]=NA+ MU[,k]=Z+ }&#62; AGE=unique(DEATH$Age)&#62; ANNEE=unique(DEATH$Year)&#62; MUF=matrix(MU[,1],length(AGE),length(ANNEE))&#62; MUH=matrix(MU[,2],length(AGE),length(ANNEE))&#62; persp(AGE[1:100],ANNEE,log(MUH[1:100,]),+ theta=-30,col="light green",shade=TRUE)&#62; persp(AGE[1:100],ANNEE,log(MUH[1:100,107:1]),+ theta=-30,col="light green",shade=TRUE)&#62; image(AGE,ANNEE,log(MUH))On a retourné le graph en mettant les années à l'envers, mais on peut aussi regarder la surface pour les femmes, ainsi que l'allure des courbes de niveau&#62; persp(AGE[1:100],ANNEE,log(MUF[1:100,107:1]),+ theta=-30,col="light green",shade=TRUE)&#62; image(AGE,ANNEE,log(MUF))Mais nous n'avons ici qu'une lecture statique des tables, et nous ne suivons pas vraiment un individu: en même temps que la personne vieilli, le temps passe.... Considérons une personne née en 1900, et regardons (ex post) l'évolution de sa probabilité de décès (en regardant non pas la probabilité de décéder à 18 d'une personne de 18 ans en 1900, mais en tenant compte du fait que la personne a fêté ses 18 ans en 1918).&#62; MUFP=MUHP=rep(NA,111)    # MU PROSPECTIF&#62; for(i in 1:111){+     MUHP[i]=MUH[i,colonne+i]+     MUFP[i]=MUF[i,colonne+i]+ }&#62; plot(0:110,log(MUHP),col="purple",type="l")&#62; D=DEATH[DEATH$Year==annee.naissance,]&#62; E=EXPOSURE[EXPOSURE$Year==annee.naissance,]&#62; MU = D[,3:5]/E[,3:5]&#62; lines(0:110,log(MU[,2]),col="blue")On peut aussi l'ajustement par un modèle à la Lee-Carter (via la fonction fournie par LifeMetrics)&#62; source("https://perso.univ-rennes1.fr/arthur.charpentier/fitModels.r")&#62; DEATH=DEATH[DEATH$Age&#60;90,] # virer age&#62;90&#62; EXPOSURE=EXPOSURE[EXPOSURE$Age&#60;90,] # virer age&#62;90&#62; XV=unique(DEATH$Age)&#62; YV=unique(DEATH$Year)&#62; ETF=t(matrix(EXPOSURE[,3],length(XV),length(YV)))&#62; DTF=t(matrix(DEATH[,3],length(XV),length(YV)))&#62; ETH=t(matrix(EXPOSURE[,4],length(XV),length(YV)))&#62; DTH=t(matrix(DEATH[,4],length(XV),length(YV)))&#62; WA=matrix(1,length(YV),length(XV))&#62; LCF = fit701(xv=XV,yv=YV,etx=ETF,dtx=DTF,wa=WA)-1166100 -&#62; -263940.2  -&#62;-204175.2 -188375.9 -&#62; -185822.4  -&#62;-184403.5 -183223.3 -&#62; -183069.8  -&#62;-182350.3 -181701.1 -&#62; -181600.6  -&#62;-181145.5 -180729.6 -&#62; -180671.5  -&#62;-180371.8 -180101.9 -&#62; -180068.6  -&#62;-179871.2 &#62; LCH = fit701(xv=XV,yv=YV,etx=ETH,dtx=DTH,wa=WA)-2564403 -&#62; -1028776  -&#62;-953629.9 -881334.2 -&#62; -868279.3  -&#62;-834413.6 -804946.8 -&#62; -802817.2  -&#62;-786343 -772775 -&#62; -772007.5  -&#62;-764501.7 -758653.8 -&#62; -758399  -&#62;-755212 -752931 -&#62; -752851  -&#62;-751620.2 -750811.1 -&#62; -750786.7  -&#62;-750351.4 &#62; plot(LCF$x,LCF$beta1,type="l",col="red")&#62; lines(LCH$x,LCH$beta1,col="blue")&#62; plot(LCF$x,LCF$beta2,type="l",col="red")&#62; lines(LCH$x,LCH$beta2,col="blue")&#62; plot(LCF$y,LCF$kappa2,type="l",col="red")&#62; lines(LCH$y,LCH$kappa2,col="blue")&#62; MUF      = DTF/ETF&#62; MUH      = DTH/ETH&#62; MUFMODEL = MUF; MUHMODEL = MUH&#62; for(i in 1:nrow(MUF)){+ for(j in 1:ncol(MUF)){+ MUFMODEL[i,j] = exp(LCF$beta1[j]+ +LCF$beta2[j]*LCF$kappa2[i])+ MUHMODEL[i,j] = exp(LCH$beta1[j]+ +LCH$beta2[j]*LCH$kappa2[i])+ }}soit, pour avoir un dessin,&#62; persp(LCH$y,LCH$x,log(MUH),+ theta=45,col="light green",shade=TRUE)&#62; persp(LCH$y,LCH$x,log(MUHMODEL),+ theta=45,col="light blue",shade=TRUE)On peut ainsi comparer les taux bruts (en vert) avec les taux estimés (en bleu)mais aussi regarder où se trouvent les erreurs le plus grandes,