Hosmer-Lemeshow

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Pierre Julion
Messages : 11
Enregistré le : 22 Oct 2014, 12:51

Hosmer-Lemeshow

Messagepar Pierre Julion » 12 Déc 2014, 15:09

Bonjour à tous,
je cherche à étudier l'adéquation des modèles de régression logistique conditionnelle suivants :

modele1<-clogit(cas~ast+ptb+pest+strata(id),cc)
modele2<-clogit(cas~ecz+ptb+pest+strata(id),cc)
modele3<-clogit(cas~rhin+ptb+pest+strata(id),cc)


Puis pour faire mon Hosmer-Lemeshow j'ai créé la fonction
hosmerlem = function(y, yhat, g=10) {
cutyhat = cut(yhat,
breaks = quantile(yhat, probs=seq(0,
1, 1/g)), include.lowest=TRUE)
obs = xtabs(cbind(1 - y, y) ~ cutyhat)
expect = xtabs(cbind(1 - yhat, yhat) ~ cutyhat)
chisq = sum((obs - expect)^2/expect)
P = 1 - pchisq(chisq, g - 2)
return(list(chisq=chisq,p.value=P))
}


Comme ni fitted.values ni fitted n'existent pour clogit, j'appelle la fonction en faisant

sansna<-na.omit(cc)
hosmerlem(sans$cas,predict(modele1)


Et à chaque fois j'ai une p-value à 1 :

$chisq
[1] -279.1566

$p.value
[1] 1

$table_Obs_Pred
0 1 0 1
[-1.72,-0.618] 52 6 107.549208 -49.549208
(-0.618,-0.427] 42 16 87.616971 -29.616971
(-0.427,-0.282] 45 12 76.879292 -19.879292
(-0.282,-0.122] 46 12 68.895700 -10.895700
(-0.122,0.00484] 38 20 60.845750 -2.845750
(0.00484,0.147] 41 16 52.864398 4.135602
(0.147,0.24] 40 18 46.865963 11.134037
(0.24,0.411] 34 23 38.426574 18.573426
(0.411,0.611] 33 25 28.301339 29.698661
(0.611,1.58] 17 41 8.754806 49.245194


Où est le problème...?

Merci d'avance

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 14 Déc 2014, 11:51

Il nous faudrait les données pour pouvoir reproduire l'exemple. D'où vient clogit ?
Renaud

Pierre Julion
Messages : 11
Enregistré le : 22 Oct 2014, 12:51

Messagepar Pierre Julion » 16 Déc 2014, 12:51

Pas de problèmes mais je ne sais pas comment vous les faire passer.
La fonction clogit vient du package survival

Eric Casellas
Messages : 767
Enregistré le : 06 Jan 2009, 14:59

Messagepar Eric Casellas » 16 Déc 2014, 12:52

Pour mettre des données sur le forum, voir http://forums.cirad.fr/logiciel-R/viewtopic.php?t=3302
Eric

Pierre Julion
Messages : 11
Enregistré le : 22 Oct 2014, 12:51

Messagepar Pierre Julion » 16 Déc 2014, 14:31

OK Merci de l'aide

Je précise que la formule générale pour les régressions logistiques étaient

Code : Tout sélectionner

logast<-clogit(cas~factor(ast)+factor(etud)+factor(fumeur)+factor(alcool)+factor(ptb)+factor(pest)+strata(id),na.omit(log)
logecz<-clogit(cas~factor(ecz)+factor(etud)+factor(fumeur)+factor(alcool)+factor(ptb)+factor(pest)+strata(id),na.omit(log)
lognbreall<-clogit(cas~nbreall+factor(etud)+factor(fumeur)+factor(alcool)+factor(ptb)+factor(pest)+strata(id),na.omit(log)


Voici le jeu de données :

Code : Tout sélectionner

log<-structure(list(id = c(1400006L, 1400008L, 1400013L, 1400014L,
1400015L, 1400022L, 1400023L, 1400027L, 1400028L, 1400029L, 1400032L,
1400034L, 1400036L, 1400038L, 1400039L, 1400040L, 1400041L, 1400042L,
1400046L, 1400051L, 1400053L, 1400055L, 1400056L, 1400058L, 1400059L,
1400064L, 1400069L, 1400073L, 1400098L, 1400100L, 1400109L, 1400111L,
1400112L, 1400115L, 1400122L, 1400128L, 1400131L, 1400132L, 1400133L,
1400138L, 1400143L, 1400147L, 1400152L, 1400153L, 1400154L, 1400155L,
1400156L, 1400159L, 1400160L, 1400164L, 1400166L, 1400006L, 1400008L,
1400013L, 1400014L, 1400015L, 1400022L, 1400023L, 1400027L, 1400028L,
1400029L, 1400032L, 1400034L, 1400036L, 1400038L, 1400039L, 1400040L,
1400041L, 1400042L, 1400046L, 1400051L, 1400053L, 1400055L, 1400056L,
1400058L, 1400059L, 1400064L, 1400069L, 1400073L, 1400098L, 1400100L,
1400109L, 1400111L, 1400112L, 1400115L, 1400122L, 1400128L, 1400131L,
1400132L, 1400133L, 1400138L, 1400143L, 1400147L, 1400152L, 1400153L,
1400154L, 1400155L, 1400156L, 1400159L, 1400160L, 1400164L, 1400166L,
1400006L, 1400008L, 1400013L, 1400014L, 1400015L, 1400022L, 1400023L,
1400027L, 1400028L, 1400029L, 1400032L, 1400034L, 1400036L, 1400038L,
1400039L, 1400040L, 1400041L, 1400042L, 1400046L, 1400051L, 1400053L,
1400055L, 1400056L, 1400058L, 1400059L, 1400064L, 1400069L, 1400073L,
1400098L, 1400100L, 1400109L, 1400111L, 1400112L, 1400115L, 1400122L,
1400128L, 1400131L, 1400132L, 1400133L, 1400138L, 1400143L, 1400147L,
1400152L, 1400153L, 1400154L, 1400155L, 1400156L, 1400159L, 1400160L,
1400164L, 1400166L, 3300090L, 3300098L, 3300099L, 3300100L, 3300108L,
3300111L, 3300128L, 3300132L, 3300133L, 3300141L, 3300155L, 3300950L,
3300951L, 3300960L, 3300964L, 3300968L, 3300971L, 3300972L, 3300977L,
3300978L, 3300981L, 3300982L, 3300983L, 3300984L, 3300985L, 3300990L,
3300999L, 3301004L, 3301005L, 3301018L, 3301019L, 3301024L, 3301027L,
3301038L, 3301056L, 3301057L, 3301059L, 3301061L, 3301063L, 3301067L,
3301070L, 3301072L, 3301073L, 3301075L, 3301079L, 3301080L, 3301092L,
3301096L, 3301097L, 3301099L, 3301104L, 3301108L, 3301111L, 3301113L,
3301115L, 3301117L, 3301125L, 3301126L, 3301128L, 3301135L, 3301142L,
3301149L, 3301151L, 3301152L, 3301153L, 3301157L, 3301159L, 3301169L,
3301171L, 3301175L, 3301184L, 3301185L, 3301187L, 3301216L, 3301269L,
3301324L, 3301326L, 3301330L, 3301342L, 3301354L, 3301365L, 3301370L,
3301375L, 3301385L, 3301394L, 3301430L, 3301432L, 3301437L, 3301440L,
3301445L, 3301449L, 3301453L, 3301457L, 3301458L, 3301459L, 3301467L,
3301474L, 3301475L, 3301479L, 3301496L, 3301526L, 3301546L, 3300090L,
3300098L, 3300099L, 3300100L, 3300108L, 3300111L, 3300128L, 3300132L,
3300133L, 3300141L, 3300155L, 3300950L, 3300951L, 3300960L, 3300964L,
3300968L, 3300971L, 3300972L, 3300977L, 3300978L, 3300981L, 3300982L,
3300983L, 3300984L, 3300985L, 3300990L, 3300999L, 3301004L, 3301005L,
3301018L, 3301019L, 3301024L, 3301027L, 3301038L, 3301056L, 3301057L,
3301059L, 3301061L, 3301063L, 3301067L, 3301070L, 3301072L, 3301073L,
3301075L, 3301079L, 3301080L, 3301092L, 3301096L, 3301097L, 3301099L,
3301104L, 3301108L, 3301111L, 3301113L, 3301115L, 3301117L, 3301125L,
3301126L, 3301128L, 3301135L, 3301142L, 3301149L, 3301151L, 3301152L,
3301153L, 3301157L, 3301159L, 3301169L, 3301171L, 3301175L, 3301184L,
3301185L, 3301187L, 3301216L, 3301269L, 3301324L, 3301326L, 3301330L,
3301342L, 3301354L, 3301365L, 3301370L, 3301375L, 3301385L, 3301394L,
3301430L, 3301432L, 3301437L, 3301440L, 3301445L, 3301449L, 3301453L,
3301457L, 3301458L, 3301459L, 3301467L, 3301474L, 3301475L, 3301479L,
3301496L, 3301526L, 3301546L, 3300090L, 3300098L, 3300099L, 3300100L,
3300108L, 3300111L, 3300128L, 3300132L, 3300133L, 3300141L, 3300155L,
3300950L, 3300951L, 3300960L, 3300964L, 3300968L, 3300971L, 3300972L,
3300977L, 3300978L, 3300981L, 3300982L, 3300983L, 3300984L, 3300985L,
3300990L, 3300999L, 3301004L, 3301005L, 3301018L, 3301019L, 3301024L,
3301027L, 3301038L, 3301056L, 3301057L, 3301059L, 3301061L, 3301063L,
3301067L, 3301070L, 3301072L, 3301073L, 3301075L, 3301079L, 3301080L,
3301092L, 3301096L, 3301097L, 3301099L, 3301104L, 3301108L, 3301111L,
3301113L, 3301115L, 3301117L, 3301125L, 3301126L, 3301128L, 3301135L,
3301142L, 3301149L, 3301151L, 3301152L, 3301153L, 3301157L, 3301159L,
3301169L, 3301171L, 3301175L, 3301184L, 3301185L, 3301187L, 3301216L,
3301269L, 3301324L, 3301326L, 3301330L, 3301342L, 3301354L, 3301365L,
3301370L, 3301375L, 3301385L, 3301394L, 3301430L, 3301432L, 3301437L,
3301440L, 3301445L, 3301449L, 3301453L, 3301457L, 3301458L, 3301459L,
3301467L, 3301474L, 3301475L, 3301479L, 3301496L, 3301526L, 3301546L,
3400002L, 3400007L, 3400015L, 3400016L, 3400018L, 3400030L, 3400034L,
3400039L, 3400041L, 3400055L, 3400060L, 3400076L, 3400083L, 3400086L,
3400090L, 3400092L, 3400093L, 3400096L, 3400098L, 3400100L, 3400101L,
3400107L, 3400108L, 3400121L, 3400123L, 3400125L, 3400150L, 3400151L,
3400152L, 3400154L, 3400155L, 3400160L, 3400164L, 3400166L, 3400002L,
3400007L, 3400015L, 3400016L, 3400018L, 3400030L, 3400034L, 3400039L,
3400041L, 3400055L, 3400060L, 3400076L, 3400083L, 3400086L, 3400090L,
3400092L, 3400093L, 3400096L, 3400098L, 3400100L, 3400101L, 3400107L,
3400108L, 3400121L, 3400123L, 3400125L, 3400150L, 3400151L, 3400152L,
3400154L, 3400155L, 3400160L, 3400164L, 3400166L, 3400002L, 3400007L,
3400015L, 3400016L, 3400018L, 3400030L, 3400034L, 3400039L, 3400041L,
3400055L, 3400060L, 3400076L, 3400083L, 3400086L, 3400090L, 3400092L,
3400093L, 3400096L, 3400098L, 3400100L, 3400101L, 3400107L, 3400108L,
3400121L, 3400123L, 3400125L, 3400150L, 3400151L, 3400152L, 3400154L,
3400155L, 3400160L, 3400164L, 3400166L, 5000003L, 5000004L, 5000011L,
5000019L, 5000021L, 5000028L, 5000034L, 5000038L, 5000040L, 5000041L,
5000042L, 5000044L, 5000046L, 5000052L, 5000057L, 5000065L, 5000070L,
5000075L, 5000080L, 5000081L, 5000083L, 5000092L, 5000096L, 5000097L,
5000104L, 5000113L, 5000115L, 5000116L, 5000117L, 5000119L, 5000127L,
5000003L, 5000004L, 5000011L, 5000019L, 5000021L, 5000028L, 5000034L,
5000038L, 5000040L, 5000041L, 5000042L, 5000044L, 5000046L, 5000052L,
5000057L, 5000065L, 5000070L, 5000075L, 5000080L, 5000081L, 5000083L,
5000092L, 5000096L, 5000097L, 5000104L, 5000113L, 5000115L, 5000116L,
5000117L, 5000119L, 5000127L, 5000003L, 5000004L, 5000011L, 5000019L,
5000021L, 5000028L, 5000034L, 5000038L, 5000040L, 5000041L, 5000042L,
5000044L, 5000046L, 5000052L, 5000057L, 5000065L, 5000070L, 5000075L,
5000080L, 5000081L, 5000083L, 5000092L, 5000096L, 5000097L, 5000104L,
5000113L, 5000115L, 5000116L, 5000117L, 5000119L, 5000127L),
    cas = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
    1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0), maladie = c(1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L), agediag = c(66.302532512, 54.299794661,
    66.997946612, 77.963039014, 68.265571526, 89.831622177, 73.42642026,
    49.368925394, 61.253935661, 40.358658453, 49.207392197, 37.954825462,
    59.186858316, 53.708418891, 38.028747433, 72.676249144, 43.613963039,
    62.976043806, 55.734428474, 53.49486653, 65.059548255, 67.545516769,
    79.126625599, 52.83504449, 70.617385352, 42.258726899, 55.770020534,
    67.945242984, 53.648186174, 58.587268994, 61.336071184, 73.681040383,
    78.288843258, 50.168377823, 43.263518138, 56.396988364, 70.223134839,
    51.362080767, 51.638603696, 50.094455852, 40.566735113, 57.267624914,
    67.052703628, 46.187542779, 65.995893224, 70.499657769, 55.857631759,
    52.10403833, 66.40109514, 62.836413415, 53.872689938, 68.865160849,
    52.62696783, 68.525667351, 79.225188227, 66.499657769, 90.633812457,
    71.572895277, 50.929500342, 60.550308008, 38.694045175, 48.887063655,
    38.047912389, 57.492128679, 55.570157426, 38.976043806, 71.321013005,
    44.689938398, 61.787816564, 55.323750856, 52.585900068, 62.863791923,
    64.722792608, 78.075290897, 52.717316906, 68.533880903, 40.492813142,
    53.242984257, 69.577002053, 53.919233402, 60.646132786, 62.012320329,
    71.852156057, 78.458590007, 49.645448323, 42.778918549, 53.683778234,
    70.99247091, 52.952772074, 53.316906229, 52.536618754, 38.422997947,
    58.283367556, 66.724161533, 44.314852841, 66.970568104, 70.502395619,
    53.785078713, 51.893223819, 65.519507187, 61.620807666, 53.070499658,
    67.121149897, 54.067077344, 67.403148528, 79.211498973, 68.840520192,
    90.206707734, 75.608487337, 49.01026694, 60.815879535, 38.855578371,
    50.020533881, 39.731690623, 56.613278576, 54.548939083, 39.72347707,
    74.855578371, 42.442162902, 62.729637235, 54.715947981, 55.16495551,
    66.937713895, 69.415468857, 78.373716632, 52.654346338, 72.049281314,
    44.533880903, 56.673511294, 66.699520876, 51.375770021, 56.635181383,
    61.577002053, 75.148528405, 75.696098563, 48.947296372, 45.245722108,
    57.199178645, 68.136892539, 50.453114305, 50.028747433, 48.550308008,
    38.324435318, 56.290212183, 68.350444901, 46.255989049, 64.624229979,
    70.39835729, 56.443531828, 50.412046543, 65.021218344, 65.451060917,
    54.305270363, 68.443531828, 67.808350445, 50.214921287, 87.994524298,
    52.377823409, 73.87816564, 51.926078029, 64.388774812, 58.592744695,
    56.084873374, 75.195071869, 64.791238877, 60.062970568, 51.600273785,
    57.63449692, 53.081451061, 59.600273785, 69.859000684, 67.663244353,
    54.045174538, 70.335386721, 61.880903491, 54.579055441, 55.674195756,
    67.857631759, 72.681724846, 54.255989049, 52.824093087, 54.26146475,
    59.561943874, 54.770704997, 68.588637919, 67.942505133, 55.017111567,
    57.952087611, 54.778918549, 65.798767967, 44.692676249, 35.436002738,
    61.021218344, 77.25119781, 54.967830253, 61.785078713, 42.995208761,
    39.794661191, 50.321697467, 49.196440794, 78.970568104, 55.148528405,
    71.077344285, 54.272416153, 75.394934976, 75.989048597, 81.070499658,
    68.648870637, 71.846680356, 38.50513347, 71.542778919, 58.776180698,
    54.346338125, 69.535934292, 61.80698152, 62.524298426, 70.001368925,
    59.271731691, 79.463381246, 46.606433949, 58.064339493, 33.226557153,
    61.541409993, 80.681724846, 75.802874743, 48.629705681, 54.277891855,
    51.775496235, 62.672142368, 57.048596851, 83.227926078, 62.910335387,
    67.813826146, 74.554414784, 57.180013689, 69.587953457, 81.174537988,
    69.467488022, 42.688569473, 46.012320329, 52.328542094, 50.88843258,
    56.980150582, 63.436002738, 68.958247775, 60.955509925, 45.615331964,
    48.153319644, 43.441478439, 63.630390144, 66.924024641, 75.846680356,
    44.391512663, 79.017111567, 80.892539357, 68.41615332, 69.80971937,
    51.989048597, 87.553730322, 50.119096509, 71.331964408, 52.914442163,
    66.984257358, 59.871321013, 53.270362765, 73.585215606, 65.267624914,
    59.033538672, 54.288843258, 60.34770705, 50.962354552, 59.663244353,
    71.397672827, 70.409308693, 56.481861739, 68.739219713, 60.807665982,
    53.418206708, 55.909650924, 70.056125941, 71.419575633, 53.297741273,
    51.868583162, 52.177960301, 58.67761807, 53.771389459, 70.863791923,
    68.117727584, 57.697467488, 55.745379877, 54.614647502, 69.976728268,
    42.811772758, 36.542094456, 59.523613963, 76.024640657, 56.10403833,
    62.231348392, 41.418206708, 37.823408624, 48.881587953, 51.279945243,
    78.338124572, 55.260780287, 68.325804244, 54.088980151, 76.372347707,
    75.682409309, 81.812457221, 67.214236824, 72.821355236, 35.934291581,
    69.297741273, 59.151266256, 51.997262149, 69.708418891, 63.682409309,
    60.848733744, 69.404517454, 61.700205339, 77.489390828, 48.240930869,
    60.095824778, 33.9137577, 61.034907598, 79.509924709, 75.140314853,
    48.032854209, 52.791238877, 51.876796715, 62.266940452, 59.728952772,
    81.744010951, 61.11430527, 65.711156742, 75.496235455, 59.488021903,
    67.531827515, 83.09650924, 69.338809035, 41.80698152, 45.303216975,
    51.277207392, 51.392197125, 57.084188912, 64.386036961, 68.054757016,
    59.797399042, 47.764544832, 47.430527036, 41.667351129, 64.553045859,
    66.006844627, 76.788501027, 45.891854894, 79.214236824, 78.746064339,
    68.780287474, 69.163586585, 50.031485284, 86.639288159, 50.707734428,
    74.455852156, 51.627652293, 65.957563313, 56.665297741, 53.363449692,
    75.61670089, 62.891170431, 61.14715948, 52.205338809, 60.421629021,
    52.807665982, 57.519507187, 70.302532512, 66.283367556, 55.326488706,
    71.263518138, 60.892539357, 52.328542094, 57.678302533, 67.72073922,
    70.327173169, 55.4770705, 50.472279261, 52.769336071, 61.801505818,
    56, 66.069815195, 67.386721424, 56.388774812, 57.516769336,
    53.125256674, 72.320328542, 42.620123203, 34.368240931, 60.082135524,
    79.898699521, 57.004791239, 61.80698152, 44.944558522, 39.121149897,
    49.46201232, 47.206023272, 78.535249829, 52.665297741, 72.950034223,
    51.321013005, 73.94661191, 76.840520192, 81.292265572, 67.386721424,
    71.167693361, 37.475701574, 72.205338809, 57.333333333, 55.40862423,
    68.580424367, 60.180698152, 60.777549624, 67.967145791, 60.501026694,
    79.961670089, 47.561943874, 58.71321013, 30.833675565, 61.700205339,
    78.321697467, 74.587268994, 49.987679671, 56.336755647, 49.990417522,
    62.080766598, 57.930184805, 82.439425051, 63.433264887, 67.841204654,
    73.987679671, 55.811088296, 67.767282683, 80.873374401, 68.580424367,
    44.772073922, 48.112251882, 52.41889117, 49.790554415, 55.693360712,
    65.919233402, 68.027378508, 62.340862423, 43.381245722, 48.186173854,
    39.893223819, 64.328542094, 66.836413415, 73.475701574, 45.516769336,
    77.149897331, 81.273100616, 56.27926078, 67.047227926, 69.875427789,
    69.284052019, 47.649555099, 67.997262149, 68.824093087, 55.660506502,
    80.227241615, 48.733744011, 51.400410678, 66.477754962, 75.778234086,
    72.741957563, 54.412046543, 70.173853525, 74.863791923, 69.568788501,
    54.918548939, 63.496235455, 64.336755647, 63.665982204, 49.503080082,
    52.599589322, 44.706365503, 37.995893224, 84.48733744, 70.113620808,
    50.95687885, 30.29431896, 56.720054757, 61.53045859, 50.472279261,
    40.062970568, 57.804243669, 68.125941136, 69.325119781, 70.592744695,
    46.622861054, 69.360711841, 66.732375086, 53.965776865, 80.372347707,
    46.422997947, 49.492128679, 64.79945243, 76.788501027, 73.650924025,
    53.267624914, 69.303216975, 77.434633812, 67.945242984, 58.913073238,
    62.614647502, 65.508555784, 61.00752909, 57.420944559, 51.436002738,
    43.216974675, 39.997262149, 82.406570842, 69.905544148, 52.48733744,
    30.258726899, 55.819301848, 60.342231348, 53.122518823, 39.82477755,
    54.562628337, 69.07871321, 70.614647502, 71.222450376, 47.069130732,
    67.852156057, 66.496919918, 57.393566051, 80.55578371, 49.817932923,
    53.694729637, 65.155373032, 74.409308693, 70.970568104, 53.185489391,
    67.4770705, 76.358658453, 67.748117728, 56.484599589, 63.942505133,
    61.694729637, 61.160848734, 59.353867214, 50.929500342, 41.916495551,
    38.863791923, 81.691991786, 72.744695414, 52.5229295, 32.881587953,
    55.348391513, 61.667351129, 49.431895962, 41.607118412, 47.107460643,
    55.589322382, 68.38329911, 70.094455852, 61.108829569, 69.719370294,
    69.659137577, 46.373716632, 46.995208761, 50.98973306, 57.229295003,
    73.9137577, 51.030800821, 49.14715948, 50.893908282, 72.832306639,
    67.54825462, 56.613278576, 59.841204654, 58.206707734, 54.182067077,
    53.859000684, 65.623545517, 69.880903491, 57.555099247, 55.170431211,
    69.941136208, 40.336755647, 66.212183436, 45.429158111, 75.882272416,
    45.144421629, 54.151950719, 67.550992471, 71.345653662, 60.553045859,
    67.854893908, 69.900068446, 47.364818617, 47.696098563, 50.598220397,
    59.520876112, 73.80971937, 52.878850103, 50.269678303, 51.159479808,
    70.702258727, 68.156057495, 57.166324435, 57.891854894, 57.894592745,
    52.098562628, 52.8678987, 64.350444901, 68.314852841, 55.701574264,
    54.707734428, 68.492813142, 42.757015743, 64.542094456, 46.729637235,
    74.652977413, 48.07118412, 53.834360027, 66.830937714, 70.447638604,
    61.434633812, 68.20807666, 72.030116359, 48.361396304, 45.188227242,
    50.36550308, 58.338124572, 75.980835044, 50.351813826, 49.869952088,
    50.57357974, 73.615331964, 66.042436687, 57.196440794, 59.381245722,
    58.726899384, 56.848733744, 53.817932923, 67.211498973, 72.413415469,
    56.191649555, 56.342231348, 71.693360712, 40.008213552, 66.948665298,
    43.808350445, 73.336071184), agequest = c(66.7816564, 54.915811088,
    67.359342916, 78.32991102, 68.76386037, 90.39835729, 73.861738535,
    50.658453114, 61.587953457, 40.66255989, 50.001368925, 38.559890486,
    59.455167693, 54.316221766, 38.85284052, 73.002053388, 46.217659138,
    63.671457906, 56.424366872, 53.998631075, 65.45927447, 67.909650924,
    81.60164271, 53.35523614, 70.929500342, 43.227926078, 56.049281314,
    68.410677618, 53.952087611, 58.85284052, 62.140999316, 74.143737166,
    78.606433949, 50.39835729, 43.550992471, 56.629705681, 70.461327858,
    51.553730322, 53.314168378, 50.417522245, 41.034907598, 57.53045859,
    67.244353183, 46.551676934, 66.650239562, 70.99247091, 56.197125257,
    52.793976728, 66.759753593, 63.263518138, 54.409308693, 70.206707734,
    54.12183436, 69.921971253, 82.92128679, 68.046543463, 93.771389459,
    74.718685832, 54.042436687, 61.763175907, 41.065023956, 51.545516769,
    40.396988364, 60.10403833, 57.708418891, 41.620807666, 73.489390828,
    48.038329911, 62.88569473, 56.936344969, 55.789185489, 65.097878166,
    65.815195072, 80.999315537, 54.042436687, 70.360027379, 42.039698836,
    55.148528405, 71.34017796, 55.586584531, 62.874743326, 64.125941136,
    74.12183436, 80.791238877, 51.255304586, 44.73100616, 55.049965777,
    72.325804244, 54.0862423, 55.690622861, 53.831622177, 40.186173854,
    60.273785079, 67.633127995, 45.28678987, 68.41615332, 71.843942505,
    55.411362081, 53.423682409, 66.721423682, 63.230663929, 55.039014374,
    69.782340862, 55.600273785, 69.571526352, 80.574948665, 70.327173169,
    94.269678303, 78.346338125, 51.96440794, 61.820670773, 41.221081451,
    53.210130048, 41.965776865, 57.782340862, 55.739904175, 41.448323066,
    76.561259411, 45.59890486, 65.169062286, 56.021902806, 57.856262834,
    70.288843258, 71.641341547, 81.281314168, 54.844626968, 73.065023956,
    46.633812457, 58.803559206, 68.424366872, 53.464750171, 58.970568104,
    63.318275154, 78.275154004, 78.562628337, 50.729637235, 46.759753593,
    58.595482546, 69.839835729, 52.106776181, 52.167008898, 50.924024641,
    40.219028063, 57.434633812, 69.393566051, 47.8275154, 66.579055441,
    71.542778919, 57.976728268, 52.287474333, 66.143737166, 66.584531143,
    55.427789185, 68.895277207, 68.150581793, 50.685831622, 88.629705681,
    52.856947296, 74.110882957, 52.265571526, 64.854209446, 59.091033539,
    57.177275838, 76.005475702, 65.256673511, 60.410677618, 51.868583162,
    57.921971253, 53.3908282, 60.284736482, 70.721423682, 68.106776181,
    54.288843258, 70.699520876, 62.283367556, 54.787132101, 56.021902806,
    68.235455168, 73.3908282, 55.151266256, 53.106091718, 54.496919918,
    59.926078029, 55.40862423, 69.437371663, 68.577686516, 55.438740589,
    58.590006845, 55.351129363, 66.173853525, 44.958247775, 36.008213552,
    61.314168378, 77.582477755, 56.06844627, 61.932922656, 43.260780287,
    40.260095825, 51.34017796, 49.924709103, 79.238877481, 55.438740589,
    71.819301848, 54.754277892, 75.950718686, 76.52019165, 83.939767283,
    68.985626283, 72.394250513, 40.722792608, 71.8275154, 59.871321013,
    54.568104038, 70.135523614, 62.231348392, 63.345653662, 70.620123203,
    59.79192334, 79.838466804, 48.054757016, 58.458590007, 34.253251198,
    62.001368925, 83.657768652, 76.845995893, 49.541409993, 54.902121834,
    52.202600958, 63.937029432, 57.412731006, 85.075975359, 63.101984942,
    68.153319644, 74.896646133, 57.623545517, 70.436687201, 81.492128679,
    69.798767967, 44.432580424, 46.929500342, 52.824093087, 51.178644764,
    57.535934292, 63.978097194, 69.448323066, 61.505817933, 46.275154004,
    49.316906229, 44.082135524, 64.446269678, 69.119780972, 77.221081451,
    46.36550308, 79.718001369, 82.590006845, 70.130047912, 71.211498973,
    53.122518823, 88.681724846, 51.838466804, 72.167008898, 54.234086242,
    67.759069131, 61.14715948, 54.595482546, 75.041752225, 66.707734428,
    60.323066393, 55.655030801, 62.004106776, 53.308692676, 60.988364134,
    72.640657084, 71.542778919, 57.757700205, 71.567419576, 61.856262834,
    54.551676934, 57.451060917, 71.33744011, 72.988364134, 54.74880219,
    53.25119781, 54.168377823, 61.207392197, 56.060232717, 73.516769336,
    70.258726899, 59.82477755, 56.917180014, 56.988364134, 70.724161533,
    43.997262149, 38.308008214, 60.449007529, 77.689253936, 57.864476386,
    63.20054757, 42.072553046, 40.082135524, 51.526351814, 52.780287474,
    79.26899384, 57.103353867, 70.338124572, 55.498973306, 77.483915127,
    77.127994524, 85.199178645, 68.45174538, 73.492128679, 38.521560575,
    70.149212868, 60.736481862, 52.840520192, 71.816563997, 65.771389459,
    62.455852156, 71.392197125, 63.074606434, 78.554414784, 50.746064339,
    61.368925394, 35.54825462, 62.502395619, 82.784394251, 76.678986995,
    49.919233402, 54.22587269, 53.262149213, 63.958932238, 60.542094456,
    85.149897331, 61.557837098, 66.64202601, 77.185489391, 60.651608487,
    68.802190281, 84.238193018, 70.217659138, 45.431895962, 47.258042437,
    52.150581793, 52.090349076, 57.979466119, 66.579055441, 69.360711841,
    61.204654346, 48.865160849, 49.440109514, 44.180698152, 65.954825462,
    68.446269678, 78.926762491, 48.427104723, 80.761122519, 81.002053388,
    70.15742642, 71.082819986, 52.501026694, 88.281998631, 52.468172485,
    76.257357974, 52.550308008, 67.107460643, 57.790554415, 54.809034908,
    76.876112252, 64.273785079, 62.587268994, 53.730321697, 62.546201232,
    56.005475702, 59.879534565, 72.098562628, 68.911704312, 57.771389459,
    72.503764545, 61.94661191, 53.456536619, 61.32238193, 69.092402464,
    71.33744011, 57.212867899, 52.624229979, 54.633812457, 65.015742642,
    58.655715264, 68.131416838, 68.829568789, 57.330595483, 58.729637235,
    54.179329227, 73.338809035, 44.783025325, 36.043805613, 60.774811773,
    80.977412731, 59.017111567, 62.390143737, 45.67008898, 40.55578371,
    51.104722793, 48.566735113, 79.148528405, 53.032169747, 74.406570842,
    52.002737851, 74.847364819, 77.768651608, 84.38329911, 68.947296372,
    72.936344969, 40.807665982, 73.697467488, 60.344969199, 56.438056126,
    70.617385352, 61.360711841, 62.135523614, 69.32238193, 62.039698836,
    80.736481862, 50.85284052, 59.939767283, 32.394250513, 63.008898015,
    81.768651608, 76.969199179, 51.61670089, 57.746748802, 51.463381246,
    64.336755647, 59.482546201, 85.615331964, 64.13963039, 68.933607118,
    75.613963039, 56.969199179, 69.106091718, 82.072553046, 69.229295003,
    47.279945243, 50.058863792, 53.292265572, 50.707734428, 56.651608487,
    67.345653662, 69.095140315, 64.013689254, 44.73100616, 50.945927447,
    40.906228611, 66.403832991, 69.648186174, 75.290896646, 47.978097194,
    78.379192334, 83.545516769, 56.462696783, 67.386721424, 70.031485284,
    69.683778234, 47.915126626, 70.978781656, 69.185489391, 56.090349076,
    80.66255989, 49.552361396, 51.668720055, 67.665982204, 77.897330595,
    73.388090349, 54.814510609, 71.72073922, 78.099931554, 70.05065024,
    55.247091034, 64.27926078, 64.955509925, 65.349760438, 49.869952088,
    52.993839836, 46.023271732, 39.118412047, 86.269678303, 70.570841889,
    51.315537303, 32.287474333, 57.117043121, 61.84257358, 50.688569473,
    40.446269678, 59.104722793, 70.15742642, 72.394250513, 72.183436003,
    49.226557153, 73.779603012, 67.613963039, 54.778918549, 81.262149213,
    49.040383299, 51.471594798, 67.498973306, 80.271047228, 74.872005476,
    54.406570842, 72.246406571, 81.092402464, 69.136208077, 60.188911704,
    64.804928131, 67.274469541, 63.611225188, 59.359342916, 52.908966461,
    45.681040383, 42.579055441, 85.568788501, 71.030800821, 53.724845996,
    33.401779603, 58.461327858, 61.577002053, 54.105407255, 40.895277207,
    55.983572895, 70.332648871, 71.838466804, 72.73100616, 49.850787132,
    71.143052704, 67.567419576, 58.455852156, 81.927446954, 52.498288843,
    55.123887748, 68.164271047, 77.338809035, 72.364134155, 54.50513347,
    69.77412731, 79.72073922, 68.895277207, 59.066392882, 65.9137577,
    63.701574264, 63.419575633, 60.810403833, 52.41889117, 43.882272416,
    40.791238877, 84.38329911, 74.472279261, 54.450376454, 35.586584531,
    57.374401095, 63.258042437, 50.494182067, 43.238877481, 47.676933607,
    55.819301848, 71.003422313, 70.929500342, 61.681040383, 72.933607118,
    70.198494182, 46.724161533, 47.452429843, 51.40862423, 58.004106776,
    74.683093771, 51.676933607, 50.119096509, 51.520876112, 72.966461328,
    67.879534565, 57.656399726, 60.049281314, 58.54072553, 54.559890486,
    54.280629706, 65.880903491, 70.097193703, 58.105407255, 55.890485969,
    70.29431896, 40.982888433, 66.674880219, 45.683778234, 76.177960301,
    47.455167693, 55.72073922, 70.787132101, 72.925393566, 62.12183436,
    71.613963039, 71.378507871, 49.451060917, 50.203969884, 52.457221081,
    61.724845996, 75.531827515, 54.516084873, 52.030116359, 53.415468857,
    71.301848049, 69.07871321, 60.161533196, 59.843942505, 59.605749487,
    53.144421629, 53.253935661, 66.277891855, 69.77412731, 57.371663244,
    56.336755647, 69.609856263, 44.20807666, 66.154688569, 48.095824778,
    76.194387406, 49.689253936, 55.159479808, 70.004106776, 73.393566051,
    62.88843258, 72.07118412, 75.017111567, 51.041752225, 48.199863107,
    53.01026694, 61.062286105, 77.538672142, 52.087611225, 51.698836413,
    52.813141684, 75.493497604, 66.99247091, 59.79192334, 61.040383299,
    60.941820671, 59.030800821, 54.606433949, 68.854209446, 74.663928816,
    57.87816564, 57.960301164, 72.947296372, 42.028747433, 68.246406571,
    45.793292266, 75.44421629), etud = c(2L, 3L, 1L, 1L, 1L,
    1L, 1L, 3L, 2L, 2L, 4L, 2L, 2L, 4L, 2L, 1L, 1L, 3L, 2L, 3L,
    3L, 3L, 1L, 2L, 1L, 3L, 2L, 3L, 3L, 2L, 3L, 2L, 1L, 2L, 2L,
    4L, 2L, 3L, 2L, 3L, 2L, 2L, 2L, 4L, 3L, 1L, 2L, 4L, 1L, 2L,
    2L, 3L, 4L, 2L, 3L, 4L, 1L, 4L, 2L, 2L, 2L, 2L, 3L, 2L, 1L,
    2L, 1L, 4L, 3L, 4L, 3L, 1L, 1L, 1L, 3L, 2L, 4L, 4L, 4L, 4L,
    3L, 2L, 3L, 1L, 2L, 1L, 4L, 4L, 1L, 4L, 4L, 4L, 1L, 1L, 4L,
    1L, 1L, 3L, 3L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 4L, 1L, 2L,
    3L, 2L, 3L, 3L, 4L, 1L, 4L, 1L, 1L, 1L, 4L, 3L, 1L, 1L, 2L,
    1L, 4L, 4L, 2L, 2L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 4L, 4L,
    1L, 4L, 3L, 1L, 1L, 3L, 2L, 2L, 4L, 1L, 1L, 1L, 1L, 1L, 1L,
    4L, 2L, 4L, 3L, 2L, 2L, 2L, 2L, 1L, 3L, 1L, 4L, 3L, 3L, 3L,
    1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 4L, 4L, 2L, 1L, 1L, 1L,
    3L, 4L, 3L, 1L, 2L, 4L, 2L, 1L, 1L, 2L, 3L, 4L, 4L, 4L, 2L,
    2L, 1L, 3L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 3L, 2L,
    2L, 1L, 1L, 1L, 3L, 2L, 4L, 2L, 2L, 1L, 2L, 2L, 4L, 4L, 4L,
    1L, 3L, 4L, 1L, 4L, 3L, 2L, 3L, 2L, 3L, 2L, 2L, 1L, 4L, 1L,
    2L, 1L, 3L, 2L, 1L, 1L, 1L, 4L, 2L, 2L, 1L, 1L, 1L, 1L, 3L,
    2L, 2L, 1L, 4L, 1L, 2L, 1L, 4L, 3L, 2L, 3L, 1L, 3L, 2L, 3L,
    1L, 2L, 4L, 4L, 3L, 2L, 4L, 2L, 4L, 3L, 3L, 1L, 2L, 3L, 3L,
    4L, 2L, 2L, 2L, 1L, 1L, 3L, 2L, 4L, 2L, 3L, 4L, 1L, 4L, 4L,
    2L, 1L, 1L, 4L, 2L, 1L, 4L, 3L, 4L, 4L, 2L, 1L, 2L, 2L, 2L,
    1L, 3L, 2L, 4L, 2L, 1L, 2L, 2L, 3L, 2L, 2L, 4L, 2L, 4L, 1L,
    4L, 1L, 4L, 2L, 1L, 4L, 3L, 3L, 4L, 3L, 3L, 4L, 4L, 3L, 4L,
    3L, 4L, 4L, 3L, 4L, 3L, 3L, 2L, 1L, 2L, 1L, 3L, 2L, 2L, 3L,
    3L, 4L, 1L, 2L, 1L, 2L, 1L, 3L, 4L, 2L, 3L, 4L, 2L, 1L, 3L,
    2L, 3L, 3L, 4L, 2L, 4L, 4L, 1L, 1L, 3L, 4L, 4L, 4L, 3L, 4L,
    4L, 4L, 1L, 3L, 4L, 3L, 4L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L,
    1L, 4L, 1L, 3L, 2L, 2L, 1L, 2L, 3L, 2L, 1L, 1L, 2L, 2L, 2L,
    3L, 4L, 3L, 2L, 4L, 4L, 2L, 2L, 3L, 2L, 1L, 1L, 2L, 2L, 1L,
    1L, 2L, 2L, 3L, 4L, 4L, 2L, 3L, 1L, 3L, 3L, 2L, 2L, 4L, 4L,
    2L, 2L, 1L, 1L, 4L, 2L, 2L, 2L, 2L, 1L, 4L, 3L, 2L, 2L, 2L,
    3L, 2L, 2L, 1L, 2L, 4L, 2L, 2L, 1L, 2L, 3L, 3L, 4L, 2L, 4L,
    2L, 4L, 3L, 4L, 2L, 1L, 2L, 4L, 4L, 3L, 2L, 2L, 4L, 2L, 2L,
    4L, 2L, 4L, 4L, 3L, 4L, 2L, 4L, 1L, 1L, 2L, 4L, 3L, 4L, 4L,
    4L, 4L, 4L, 2L, 3L, 4L, 4L, 4L, 2L, 2L, 3L, 4L, 2L, 4L, 1L,
    4L, 3L, 2L, 2L, 4L, 3L, 4L, 2L, 3L, 2L, 4L, 4L, 3L, 2L, 2L,
    4L, 2L, 2L, 2L, 4L, 4L, 3L, 2L, 1L, 3L, 2L, 4L, 3L, 2L, 2L,
    4L, 2L, 1L, 1L, 1L, 3L, 2L, 2L, 3L, 2L, 2L, 1L, 1L, 1L, 2L,
    2L, 2L, 2L, 1L, 1L, 4L, 1L, 2L, 3L, 1L, 1L, 4L, 1L, 2L, 4L,
    2L, 4L, 2L, 2L, 2L, 2L, 1L, 1L, 4L, 2L, 4L, 3L, 1L, 1L, 3L,
    2L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 1L, 1L, 1L, 1L, 1L, 3L, 3L,
    1L, 4L, 1L, 2L, 1L, 2L, 3L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 3L,
    2L, 2L, 4L, 1L, 4L, 1L, 1L, 1L, 2L, 1L, 1L, 4L, 2L, 1L, 1L,
    4L, 1L, 3L, 1L), fumeur = c(1L, 0L, 1L, 0L, 0L, 1L, 2L, 0L,
    2L, 2L, 0L, 1L, 0L, 0L, 0L, 0L, 2L, 2L, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L,
    0L, 1L, 2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 2L, 1L, 1L,
    1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 2L, 2L, 0L, 0L, 0L, 0L, 2L,
    0L, 2L, 0L, 0L, 1L, 1L, 0L, 0L, 2L, 2L, 1L, 2L, 1L, 0L, 1L,
    1L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 2L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 2L,
    2L, 1L, 2L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 2L, 0L, 2L,
    0L, 2L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 2L, 0L, 0L, 1L, 1L,
    0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 2L, 1L, 0L,
    0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L,
    0L, 0L, 1L, 0L, 1L, 1L, 2L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L,
    1L, 1L, 0L, 2L, 0L, 0L, 0L, 0L, 2L, 2L, 0L, 0L, 0L, 1L, 0L,
    1L, 1L, 0L, 1L, 1L, 0L, 1L, 2L, 0L, 0L, 1L, 0L, 1L, 0L, 0L,
    1L, 2L, 0L, 0L, 2L, 0L, 1L, 1L, 1L, 1L, 0L, 2L, 0L, 2L, 0L,
    1L, 2L, 0L, 1L, 2L, 1L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 0L, 0L,
    1L, 1L, 1L, 0L, 0L, 2L, 2L, 0L, 0L, 1L, 2L, 0L, 2L, 1L, 2L,
    2L, 1L, 1L, 0L, 0L, 2L, 2L, 2L, 1L, 0L, 2L, 1L, 2L, 0L, 0L,
    0L, 1L, 2L, 2L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 2L, 2L, 1L, 0L,
    0L, 1L, 0L, 0L, 2L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L,
    1L, 0L, 2L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,
    0L, 0L, 0L, 0L, 2L, 0L, 1L, 0L, 1L, 1L, 2L, 0L, 1L, 0L, 0L,
    1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L,
    0L, 0L, 1L, 2L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L,
    1L, 0L, 2L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L,
    2L, 0L, 0L, 2L, 0L, 1L, 1L, 0L, 0L, 1L, 2L, 2L, 0L, 1L, 0L,
    2L, 1L, 1L, 2L, 1L, 0L, 0L, 2L, 0L, 0L, 0L, 0L, 1L, 0L, 2L,
    1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 2L, 2L, 1L, 0L, 0L, 2L, 0L,
    0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 0L, 1L, 1L, 1L, 0L, 1L, 2L, 0L, 2L, 0L, 0L, 0L, 1L, 1L,
    0L, 1L, 0L, 0L, 1L, 2L, 0L, 0L, 0L, 0L, 0L, 1L, 2L, 1L, 2L,
    1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L,
    1L, 2L, 2L, 0L, 2L, 2L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 2L, 1L, 0L, 1L, 2L, 1L, 0L, 2L, 0L, 0L, 1L, 2L, 0L, 0L,
    2L, 2L, 0L, 2L, 2L, 0L, 1L, 0L, 2L, 2L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 2L, 1L, 2L, 1L, 1L, 0L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L,
    0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L,
    0L, 0L, 0L, 2L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L,
    1L), alcool = c(1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, NA, 0L, 1L, 1L, 0L, 0L, NA, 1L, 0L, 0L,
    1L, 1L, 1L, 0L, 1L, 0L, NA, 1L, 0L, 1L, 0L, 2L, NA, 1L, 0L,
    1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 1L, 1L, 1L,
    0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, NA, 1L, 0L, 0L,
    1L, 1L, NA, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, NA, 2L, 2L,
    1L, 0L, 0L, NA, 1L, 2L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L,
    0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L,
    1L, 1L, NA, 0L, 1L, 0L, 1L, 1L, NA, 1L, 1L, 1L, 1L, 1L, 1L,
    1L, 0L, 0L, NA, 1L, 0L, 1L, 1L, 1L, NA, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L,
    1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 2L, 0L, 0L, 0L, 0L, 2L,
    1L, 2L, 2L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L,
    1L, 2L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L,
    NA, 2L, 0L, 1L, 1L, 2L, 0L, 0L, 0L, 2L, 0L, 1L, 1L, 1L, 1L,
    0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 2L, 1L, 0L,
    1L, 0L, 0L, 0L, 2L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 2L, NA,
    NA, 1L, 0L, 0L, 2L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 2L,
    1L, 2L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L,
    0L, 0L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 0L, 0L, 1L,
    1L, 1L, 2L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, NA, 1L, 0L,
    1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 2L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L,
    0L, 2L, 1L, 0L, 1L, 1L, 1L, 0L, 2L, 1L, 0L, NA, NA, 0L, 1L,
    1L, 2L, 1L, 2L, 1L, 0L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 0L, 2L,
    1L, 1L, 1L, 2L, 1L, 1L, 1L, 0L, 1L, 2L, 2L, 0L, 1L, 1L, 1L,
    1L, 2L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 1L,
    1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, NA, 2L, 1L, 0L, 1L, 1L,
    1L, 1L, 1L, 2L, 0L, 1L, 1L, 1L, 2L, 0L, 1L, 1L, 1L, 1L, 1L,
    0L, 0L, 0L, 0L, 1L, 0L, 2L, 0L, 0L, 2L, 0L, 1L, 0L, 1L, 1L,
    1L, 1L, 0L, 1L, 1L, 0L, 0L, 2L, NA, NA, 1L, 1L, 1L, 0L, 0L,
    0L, 0L, 1L, NA, 1L, 1L, 1L, 0L, 1L, 1L, NA, 1L, 0L, 0L, 1L,
    0L, 1L, 1L, 0L, NA, 0L, 2L, 0L, 0L, 0L, 0L, 1L, 0L, 2L, 0L,
    1L, 1L, 0L, 0L, 1L, 2L, 0L, NA, 2L, 2L, 2L, 0L, 1L, 1L, NA,
    1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, NA, 1L, 2L, 2L, 0L, 1L, 1L,
    1L, 2L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, NA, 1L, 2L, 1L,
    1L, 1L, 2L, NA, 0L, 2L, 0L, 0L, 1L, 2L, 0L, 0L, NA, 0L, 1L,
    0L, 1L, 0L, 1L, 1L, 0L, 2L, 0L, 0L, 1L, 1L, 0L, NA, 1L, 1L,
    NA, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 1L, 1L, 1L, 1L, NA, 1L,
    0L, NA, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 2L, 2L, 0L, 0L, 0L,
    1L, 1L, 0L, 0L, 1L, 1L, 2L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, NA,
    0L, 1L, NA, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L,
    0L, 1L, 2L, 1L, 0L, 0L, 1L, 1L, 2L, 0L, 0L, 0L, 1L), ptb = c(0L,
    NA, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA,
    0L, 2L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 1L,
    0L, 0L, 1L, 2L, 1L, 0L, 0L, 1L, 2L, 2L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L,
    2L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L,
    0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, NA, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
    1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L,
    0L, 0L, 0L, 0L, NA, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 2L, 1L,
    2L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L,
    0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, NA, 0L, 0L, 0L, 0L, 0L,
    NA, NA, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 2L, 0L, 0L, 0L, 1L, 0L, 2L, 0L, 0L, 0L, NA,
    0L, 1L, 0L, 1L, 0L, 1L, 0L, NA, 0L, 1L, 1L, 0L, 1L, 1L, 0L,
    1L, 1L, NA, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L,
    0L, 1L, 0L, NA, 0L, 1L, 1L, 1L, NA, 1L, 1L, 0L, 1L, 1L, 1L,
    0L, 1L, 1L, 0L, 0L, 1L, NA, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 2L,
    0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L,
    1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L,
    1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L,
    1L, 0L, NA, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L,
    1L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L,
    0L, 0L, 2L, 0L, 1L, NA, 0L, NA, 0L, 1L, 0L, 0L, 1L, 0L, 1L,
    1L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 2L, 0L, 0L, 0L, 1L,
    1L, 1L, 1L, 1L, 1L, NA, 0L, 0L, 1L, NA, 0L, 0L, 1L, 0L, 0L,
    0L, 2L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L,
    0L, 0L, 1L, 2L, 0L, 1L, 0L, 0L, 1L, 1L, NA, 1L, 1L, 0L, 1L,
    0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, NA, 1L,
    1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    NA, 0L, 1L, 2L, 2L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L,
    1L, 2L, 1L, 1L, 0L, 1L, NA, 2L, 2L, 0L, 1L, 2L, 0L, 1L, 0L,
    1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L,
    0L, 0L, 1L, 1L, NA, NA, 0L, 0L, 0L, 0L, 1L, NA, 1L, 1L, 0L,
    1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    1L, 1L, 0L, 0L, 0L, 1L, 2L, 0L, 0L, 1L, 2L, 1L, 0L, 0L, 2L,
    1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, NA, 0L,
    1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, NA, 0L, 0L, 0L, 1L,
    1L, 0L, 0L, 2L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L,
    2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L), pest = c(1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L), ast = c(0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L), ecz = c(1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, NA, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, NA,
    0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,
    0L, NA, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L,
    0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, NA, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L),
    rhin = c(0L, 0L, 0L, 1L, 0L, 0L, 0L, NA, 1L, 0L, 0L, 0L,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, NA, 0L, 1L, 0L,
    NA, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    1L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, NA, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L,
    0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L,
    0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L,
    1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
    0L, NA, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 0L, 1L,
    0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L,
    0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L,
    1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, NA, 0L, 0L,
    0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, NA, 0L, 0L, 0L,
    1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, NA, 0L, 0L, 1L, 0L, 0L,
    0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L), nbreall = c(1L,
    0L, 1L, 1L, 1L, 0L, 0L, NA, 2L, 0L, 0L, NA, 0L, 0L, 1L, 0L,
    0L, 1L, 0L, 0L, 0L, 2L, 0L, 3L, 0L, 2L, 0L, 0L, 0L, 0L, 1L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 2L, 1L, 0L, 0L, 0L, 0L, 0L,
    3L, 0L, 0L, 0L, 1L, 0L, 2L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    1L, 0L, 0L, 1L, 0L, 0L, 0L, NA, 0L, 2L, 1L, NA, 0L, 0L, 1L,
    3L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 2L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L,
    0L, 1L, 0L, 1L, 1L, 2L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L,
    0L, 0L, 0L, 2L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 2L, 1L, 0L, 2L,
    0L, 2L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 3L, 0L, 1L, 1L, 0L, 0L,
    0L, 1L, 2L, NA, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L,
    2L, 0L, 0L, 1L, 0L, 0L, 2L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L,
    0L, 0L, 0L, 1L, 0L, NA, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 2L, 0L, 1L, 0L, 1L, 2L, 0L, 0L, 1L, 0L, 0L, 1L, 2L, 0L,
    0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L,
    0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 2L,
    0L, 1L, 1L, 1L, 0L, 0L, 0L, 2L, 1L, 0L, 1L, 0L, 0L, 1L, 1L,
    0L, 0L, 1L, 0L, 1L, 1L, 0L, 2L, 1L, 1L, 1L, 1L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, NA, 0L, 0L,
    0L, 0L, 0L, 2L, 0L, 0L, 0L, 3L, 0L, 0L, 0L, 1L, 0L, 0L, 1L,
    1L, 0L, NA, 0L, 0L, 2L, 0L, 0L, NA, 0L, 0L, 0L, 1L, 0L, 0L,
    1L, 0L, 0L, 2L, 2L, 1L, 2L, 0L, 2L, 0L, NA, 0L, 0L, 3L, 0L,
    0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 0L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 2L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
    1L, 2L, 0L, 0L, 0L, 3L, 0L, 2L, 1L, 0L, 0L, 0L, 0L, 0L, 1L,
    1L, 0L, 1L, 1L, 3L, 1L, 1L, 2L, 0L, 0L, 0L, 1L, 1L, 1L, 0L,
    2L, 0L, 1L, 0L, 2L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 1L, 1L,
    0L, NA, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L,
    1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    2L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 3L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 0L, 2L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
    0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,
    0L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
    2L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L,
    0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,
    0L, 1L, 0L, 0L, 1L, 0L, NA, 0L, 1L, 0L, 3L, 2L, 0L, 0L, 0L,
    1L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, NA, 0L, 0L, 1L, 2L, 1L, 0L,
    0L, 0L, 0L, 0L, 2L, 0L, 0L, NA, 0L, 0L, 0L, 1L, 0L, 0L, 0L,
    0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 2L, 0L, 0L, 1L, 0L, 0L, 0L,
    1L, 0L, 0L, 1L, 0L, NA, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L,
    0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L)), .Names = c("id", "cas",
"maladie", "agediag", "agequest", "etud", "fumeur", "alcool",
"ptb", "pest", "ast", "ecz", "rhin", "nbreall"), class = "data.frame", row.names = c(5L,
7L, 11L, 12L, 13L, 17L, 18L, 21L, 22L, 23L, 25L, 26L, 28L, 29L,
30L, 31L, 32L, 33L, 36L, 41L, 43L, 44L, 45L, 46L, 47L, 49L, 52L,
56L, 74L, 76L, 82L, 83L, 84L, 87L, 90L, 94L, 96L, 97L, 98L, 102L,
105L, 108L, 111L, 112L, 113L, 114L, 115L, 117L, 118L, 120L, 122L,
131L, 133L, 137L, 138L, 139L, 143L, 144L, 147L, 148L, 149L, 151L,
152L, 154L, 155L, 156L, 157L, 158L, 159L, 162L, 167L, 169L, 170L,
171L, 172L, 173L, 175L, 178L, 182L, 200L, 202L, 208L, 209L, 210L,
213L, 216L, 220L, 222L, 223L, 224L, 228L, 231L, 234L, 237L, 238L,
239L, 240L, 241L, 243L, 244L, 246L, 248L, 257L, 259L, 263L, 264L,
265L, 269L, 270L, 273L, 274L, 275L, 277L, 278L, 280L, 281L, 282L,
283L, 284L, 285L, 288L, 293L, 295L, 296L, 297L, 298L, 299L, 301L,
304L, 308L, 326L, 328L, 334L, 335L, 336L, 339L, 342L, 346L, 348L,
349L, 350L, 354L, 357L, 360L, 363L, 364L, 365L, 366L, 367L, 369L,
370L, 372L, 374L, 379L, 385L, 386L, 387L, 391L, 392L, 396L, 397L,
398L, 403L, 406L, 409L, 410L, 415L, 419L, 422L, 424L, 425L, 427L,
428L, 430L, 431L, 432L, 433L, 434L, 436L, 440L, 443L, 444L, 451L,
452L, 453L, 455L, 463L, 473L, 474L, 475L, 477L, 479L, 481L, 483L,
485L, 486L, 487L, 488L, 489L, 494L, 497L, 498L, 499L, 501L, 504L,
505L, 506L, 507L, 508L, 514L, 515L, 517L, 518L, 522L, 526L, 527L,
528L, 529L, 531L, 532L, 537L, 539L, 540L, 542L, 543L, 544L, 550L,
553L, 555L, 556L, 557L, 560L, 565L, 566L, 568L, 569L, 575L, 580L,
588L, 590L, 594L, 595L, 599L, 601L, 604L, 607L, 608L, 609L, 614L,
616L, 617L, 618L, 620L, 622L, 624L, 631L, 637L, 638L, 639L, 643L,
644L, 648L, 649L, 650L, 655L, 658L, 661L, 662L, 667L, 671L, 674L,
676L, 677L, 679L, 680L, 682L, 683L, 684L, 685L, 686L, 688L, 692L,
695L, 696L, 703L, 704L, 705L, 707L, 715L, 725L, 726L, 727L, 729L,
731L, 733L, 735L, 737L, 738L, 739L, 740L, 741L, 746L, 749L, 750L,
751L, 753L, 756L, 757L, 758L, 759L, 760L, 766L, 767L, 769L, 770L,
774L, 778L, 779L, 780L, 781L, 783L, 784L, 789L, 791L, 792L, 794L,
795L, 796L, 802L, 805L, 807L, 808L, 809L, 812L, 817L, 818L, 820L,
821L, 827L, 832L, 840L, 842L, 846L, 847L, 851L, 853L, 856L, 859L,
860L, 861L, 866L, 868L, 869L, 870L, 872L, 874L, 876L, 883L, 889L,
890L, 891L, 895L, 896L, 900L, 901L, 902L, 907L, 910L, 913L, 914L,
919L, 923L, 926L, 928L, 929L, 931L, 932L, 934L, 935L, 936L, 937L,
938L, 940L, 944L, 947L, 948L, 955L, 956L, 957L, 959L, 967L, 977L,
978L, 979L, 981L, 983L, 985L, 987L, 989L, 990L, 991L, 992L, 993L,
998L, 1001L, 1002L, 1003L, 1005L, 1008L, 1009L, 1010L, 1011L,
1012L, 1018L, 1019L, 1021L, 1022L, 1026L, 1030L, 1031L, 1032L,
1033L, 1035L, 1036L, 1041L, 1043L, 1044L, 1046L, 1047L, 1048L,
1054L, 1057L, 1059L, 1060L, 1061L, 1064L, 1069L, 1070L, 1072L,
1073L, 1079L, 1084L, 1092L, 1094L, 1098L, 1099L, 1103L, 1105L,
1108L, 1111L, 1112L, 1113L, 1118L, 1120L, 1121L, 1122L, 1124L,
1126L, 1128L, 1136L, 1139L, 1145L, 1146L, 1148L, 1155L, 1158L,
1160L, 1161L, 1171L, 1173L, 1183L, 1186L, 1187L, 1191L, 1193L,
1194L, 1195L, 1197L, 1199L, 1200L, 1204L, 1205L, 1217L, 1218L,
1219L, 1230L, 1231L, 1232L, 1234L, 1235L, 1237L, 1241L, 1242L,
1254L, 1257L, 1263L, 1264L, 1266L, 1273L, 1276L, 1278L, 1279L,
1289L, 1291L, 1301L, 1304L, 1305L, 1309L, 1311L, 1312L, 1313L,
1315L, 1317L, 1318L, 1322L, 1323L, 1335L, 1336L, 1337L, 1348L,
1349L, 1350L, 1352L, 1353L, 1355L, 1359L, 1360L, 1372L, 1375L,
1381L, 1382L, 1384L, 1391L, 1394L, 1396L, 1397L, 1407L, 1409L,
1419L, 1422L, 1423L, 1427L, 1429L, 1430L, 1431L, 1433L, 1435L,
1436L, 1440L, 1441L, 1453L, 1454L, 1455L, 1466L, 1467L, 1468L,
1470L, 1471L, 1473L, 1477L, 1478L, 1491L, 1492L, 1497L, 1500L,
1502L, 1506L, 1510L, 1513L, 1514L, 1515L, 1516L, 1517L, 1519L,
1524L, 1529L, 1535L, 1539L, 1543L, 1547L, 1548L, 1550L, 1557L,
1560L, 1561L, 1568L, 1574L, 1576L, 1577L, 1578L, 1580L, 1585L,
1591L, 1592L, 1597L, 1600L, 1602L, 1606L, 1610L, 1613L, 1614L,
1615L, 1616L, 1617L, 1619L, 1624L, 1629L, 1635L, 1639L, 1643L,
1647L, 1648L, 1650L, 1657L, 1660L, 1661L, 1668L, 1674L, 1676L,
1677L, 1678L, 1680L, 1685L, 1691L, 1692L, 1697L, 1700L, 1702L,
1706L, 1710L, 1713L, 1714L, 1715L, 1716L, 1717L, 1719L, 1724L,
1729L, 1735L, 1739L, 1743L, 1747L, 1748L, 1750L, 1757L, 1760L,
1761L, 1768L, 1774L, 1776L, 1777L, 1778L, 1780L, 1785L))


Merci !

Florent Aubry
Messages : 324
Enregistré le : 25 Juin 2010, 10:21

Messagepar Florent Aubry » 18 Déc 2014, 12:41

As-tu essayé en modifiant le nombre de groupes pour ton test car semble-t-il le choix du nombre de groupes peut largement influencer les résultats du test (cf. "Why I Don’t Trust the Hosmer-Lemeshow Test for Logistic Regression" par P. Allison, http://www.statisticalhorizons.com/hosmer-lemeshow ?

De plus, J. Bartlett fait remarquer que "As far as I have seen, there is little guidance as to how to choose the number of groups g. Hosmer and Lemeshow's conclusions from simulations were based on using g>p+1, suggesting that if we have 10 covariates in the model, we should choose g>11, although this doesn't appear to be mentioned in text books or software packages.
Intuitively, using a small value of g ought to give less opportunity to detect mis-specification. However, if we choose g to large, the numbers in each group may be so small that it will be difficult to determine whether differences between observed and expected are due to chance or indicative or model mis-specification.
A further problem, highlighted by many others (e.g. Paul Allison) is that, for a given dataset, if one changes g, sometimes one obtains a quite different p-value, such that with one choice of g we might conclude our model does not fit well, yet with another we conclude there is no evidence of poor fit. This is indeed a troubling aspect of the test." (http://thestatsgeek.com/2014/02/16/the-hosmer-lemeshow-goodness-of-fit-test-for-logistic-regression/)

Donc, le problème ne vient peur-être pas d'une erreur dans le code de R, mais d'un problème dans le choix de ta stratégie pour détecter les erreurs de modélisation.

Pierre Julion
Messages : 11
Enregistré le : 22 Oct 2014, 12:51

Messagepar Pierre Julion » 04 Jan 2015, 14:34

Oui j'ai essayé en diminuant à 5 groupes, mais le problème reste identique. Je me vois mal descendre plus bas.

Je sais que le test d'H-L n'est pas la panacée mais c'est le seul que je maîtrise. Quelles autres stratégies proposeriez-vous ?


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité