Bonjour,
Mon problème concerne l'utilisation de la fonction Mort2Dsmooth,
cette fonction du package : MortalitySmooth, permet d'effectuer un lissage matriciel par B spline.
voici la fonction telle que décrite dans le package :
avec x, ici mon vecteur d'âge
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6 Mort2Dsmooth(x, y, Z, offset, W, overdispersion=FALSE, ndx = c(floor(length(x)/5), floor(length(y)/5)), deg = c(3, 3), pord = c(2, 2), lambdas = NULL, df = NULL, method = 1, coefstart = NULL, control = list())
y mon vecteur d'années
et Z mon nombre de décès (matrice de mes nombres de décès par âge et par année, source : Human Mortality Databse, France.)
voici mon code :
message d'erreur :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20 DATAdeath=read.table(file.choose(""),header=FALSE,sep="",dec=",",quote="") DATAexpo=read.table(file.choose(""),header=FALSE,sep="",dec=",",quote="") années=1816:2012 ages=18:100 #j'ai bien vérifié la dimension de DATAdeath, DATAexpo, ligne et colonne. colnames(DATAdeath) <- ages rownames(DATAdeath) <- années DATAdeath colnames(DATAexpo) <- ages rownames(DATAexpo) <- années DATAexpo Mort2Dsmooth(années,ages,as.matrix(DATAdeath),as.matrix(DATAexpo), W <- matrix(1,length(années), length(ages)), overdispersion=FALSE, ndx = c(floor(length(x)/5), floor(length(y)/5)), deg = c(3, 3), pord = c(2, 2), df = NULL, method = 1, coefstart = NULL,control=list())
Je ne vois pourquoi mon vecteur x, ici années, n'est pas détecté.Erreur dans Mort2Dsmooth_checker(x = x, y = y, Z = Z, offset = offset, W = W, : objet 'x' introuvable
Merci !
Partager