Organisation de type tétracoralliaire des rides septales de Palaeacis (Cnidaria, Carbonifére)

Plusquellec, Y., Lafuste, J. & Webb, G. E. 1990 10 15: Organisation de type tétracoralliaire des rides septales de Palaeacis (Cnidaria, Carbonifére). Lethaia, Vol. 23, pp. 385–397. Oslo. ISSN 0024–1164. Precedemment signalée chez Kerforneidictyum (Pleurodictyforme dévonien), la rare disposition des rides septales suivant un mode tétracoralliaire est décrite chez un Palaeack sp. du Viséen d'Algérie et chez P. cuneiformis subsp. A du Viseen d'Australie. Chez Palaeacis sp. sont mises en evidence une ride septale cardinale forte en position externe opposée à une ride antipode déterminant un plan de symétrie bilatérale; de part et d'autre une zone systématiquement confuse marque l'emplacement des rides alaires. La reconnaissance des périantipodes prête à discussion et ce secteur présente une originalité certaine par rapport à celui des Tétracoralliaires. Chez P. cuneiformis subsp. A, une coupe transversale montre la présence de quatre rides proéminentes qui fusionnent dans I'axe du corallite et I'existence d'une symétrie bilatérale. La position systématique du genre Palaeacis est discutée: son attribution aux Tétracoralliaires est exclue, ancun caractére n'est incompatible avec le statut de Tabulé, la création d'une unité taxon-omique nouvelle regroupant Palaeacis, Trachypsammia et les Pleurodictyformes est envisageable mais actuellement prématurée. ? Palaeacis, Tabulata, Rugosa, rides septales, systématique, relations deparenré, Viséen. The tetracoral pattern of septal ridges described in Kerforneidictyum (Devonian Pleurodictyum-like coral) is scar'ce in the Tabulata, but has been found in Palaeacis sp. from the Visean of Algeria and in P. cuneiformis subsp. A from the Visean of Australia. Palaeacis sp. shows a strong cardinal septal ridge in peripheral position, just opposite the counter ridge. thus creating a plane of bilateral symmetry: an area on both sides. which is always blurred, indicates the position of the alar ridges. counter-lateral ridges are not easily identified and yield differing interpretation. Compared with the Rugosa this area presents distinctive features. In P. cuneiformis subsp. A, a transverse section exhibits four well-developed septal ridges that coalesce in the axis of the corallite with obvious bilateral symmetry. The systematic position of the genus Palaeacis is questioned: Rugosa it is not; it may belong to Tabulata: however, erecting a new higher taxon for Palaeacis, Trachypsammia and other Pleurodictyum-like corals is premature at our present state of knowledge. ? Palaeacis, Tabulata, Rugosa, septal ridges, systematic. affinities, Visean.     

Les possibilités de l'application des germes pathogènes dans la lutte contre les chenilles défoliatrices

Résumé L'utilisation des germes pathogènes est une solution du problème de la sélectivité. Des solutions concurrentes existent avec certains pesticides chimiques. L'utilisation des germes pathogènes est donc envisagée d'une manière pragmatique en négligeant les problèmes de recherches fondamentales. Pour cette raison les mycoses et les protozoonoses ne sont pas retenues. Les viroses offrent un exemple de la plus grande spécificité d'action. L'application des virus d'insectes a été surtout développée en sylviculture. Une grande difficulté est la multiplication des virus sur organismes vivants. Les bactérioses peuvent être seules retenues pour l'agriculture et l'arboriculture intensives: elles ont un spectre d'action assez étendu mais suffisamment sélectif; les bactéries qui nous intérresent peuvent être multipliées industriellement. Le statut deBacillus thuringiensis Berliner et les propriétés des préparations commerciales ont fait l'objet de développements récents. Dans une revue critique des applications de ces préparations en arboriculture fruitière, on considère d'abord le cas des Mineuses et des Tordeuses. Presque toutes ces chenilles sont susceptibles àB. thuringiensis. Mais la difficulté d'utilisation des préparations est de nature éthologique. De bons résultats sont obtenus au Canada contreArgyrotaenia et en U.R.S.S. contreHyponomcuta. En U.R.S.S. on utilise des produits mixtes à concentration réduite (0,5%) de préparation bactérienne et à concentration extrêmement faible (0,002%) de D.D.T. Les résultats contre le Carpocapse sont négatifs. Les résultats contre les chenilles défoliatrices sont très satisfaisants: en Hollande et au Canada contreMalacosoma; au Canada contre les Géométrides,Operophtera etAlsophila espèces très susceptibles àB. thuringiensis, à la dose de 1 à 2 pounds de préparation (titrant 30×10⁹ spores viables par gramme) pour 100 gallons, soit environ, 1,5%; en France, ces résultats sont confirmés. En conclusion, l'utilisation en verger des prépararations commerciales deBacillus thuringiensis doit être ?intégrée? dans un programme général de protection comprenant à la fois les traitements chimiques et les méthodes culturales. De tels programmes sont actuellement expérimentés dans certains vergers de Pommiers par l'I.P.O. et le T.N.O. en Hollande.     

Survie de Bemisia tabaci (Homoptera: Aleyrodidae) aux conditions hivernales de la zone nord méditerranée

L’objectif de ce travail est d’étudier l’influence des températures hivernales sur la survie des B. tabaci adultes. Les individus sont soumis à un cycle de températures qui correspond aux températures observées pendant les 24 heures d’une des journées les plus froides de l’hiver du sud de la France (Perpignan). L’évolution de ces températures est reproduite en conditions contrôlées, et les individus sont exposés à un, trois ou quatre de ces cycles. Après un cycle de 24 heures, le taux de survie moyen des individus est de 73,5% alors que celui du lot témoin est de 85%. Après une exposition à 3 cycles, le taux moyen de survie reste encore relativement élevé (44,4%). Les chances de survivre à 4 cycles de 24 heures est de 7,4%. Le Nord méditerranéen semble être une zone tampon pour le maintien en extérieur des populations de B. tabaci car ni les températures minimales régulièrement enregistrées, ni leur durée, ne sont clairement limitantes.     

Une abeille afrotropicale spécialisée dans la récolte du pollen de Graminées (Poaceae): Lipotriches notabilis (Schletterer 1891) (Hymenoptera Apoidea Halictidae)

Dans les zones de savanes de l’Afrique, un genre d’abeille, Lipotriches Gerstaecker 1858, s’est spécialisé dans la collecte du pollen de graminées. Un site de nidification et l’aire de butinage de Lipotriches notabilis ont été suivis pendant trois années dans la région de Ngaoundéré au Cameroun. Le régime alimentaire pollinique est composé presque exclusivement du pollen de graminées, notamment Brachiaria ruziziensis abondant dans cette région. Le maïs est aussi visité et l’abeille contribue indirectement à la pollinisation par la mise en suspension dans l’air du pollen. Comme les graminées n’offrent pas de ressource sucrée, les femelles comme les mâles de L. notabilis butinent de temps en temps les Asteraceae pour la collecte de nectar. La consommation du pollen de graminées par les femelles a été aussi observée. Lespèce niche en bourgades dont la taille varie d’une dizaine à une centaine de nids. Le nid creusé dans un sol horizontal est du type progressif. Il comprend un tumulus, une cheminée verticale, un conduit principal vertical pouvant atteindre 65 cm de profondeur et en moyenne 3 conduits latéraux obliques de 4 à 16 cm aboutissant dans une ou plusieurs cellules successives. Les cellules des conduits latéraux sont approvisionnées simultanément et fermées de manière régressive. On compte un maximum de 10 cellules par nid. Généralement un nid est habité par une seule femelle, mais certains nids sont habités par deux ou trois. Lactivité de butinage de cette espèce est limitée dans la matinée. Dès 7 heures, elle s’envole pour le site de butinage. Après 11 heures, il n’y a généralement plus de nids ouverts. Les mâles n’ont pas été aperçus au niveau du site de nidification et l’accouplement a lieu sur le site de butinage. La période d’activité commence avec la saison des pluies en avril et se termine au début de la saison sèche en décembre, avec la fanaison des graminées.     

Imported and indigenous natural enemies of citrus coccids and aphids in Cyprus,and an assessment of their potential value in integrated control programmes

Résumé L'application d'insecticides de contact sur lesCitrus à Chypre a apporté des perturbations à l'équilibre naturel si bien que certains insectes nuisibles ont pris une importance accrue et de sévères infestations d'insectes normalement sans importance sont survenues. Ceci est particulièrement sensible dans le cas des Coccides. La situation dans les cultures non traitées révèle de fa?on évidente que, à partLepidosaphes beckii newm., les cochenilles nuisibles ont des ennemis naturels capables d'assurer un contr?le économique en l'absence de facteurs adverses susceptibles de freiner leur activité. L'équilibre naturel du Pou rouge (Aonidiella aurantii mask) est très facilement bouleversé. Il devient spécialement dangereux quand les insecticides de contact sont utilisés. Les ennemis les plus fréquents des Cochenilles desCitrus sont répertoriés. Dans quelques cas, l'importance relative des espèces attaquant certains insectes nuisibles dans des situations variées, est précisée. L'importation de 5 espèces d'Aphytis est notée. L'intégration de la lutte chimique qui est nécessaire avec la lutte biologique est discutée. Enfin, les moyens de détermination des programmes de traitements sont esquissés.      

Etude Ultrastructurale du Genre Enteromonas da Fonseca (Zoomastigophorea) et Révision de l'Ordre des Diplomonadida Wenyon

RESUME. Deux espèces d'Enteromonas sont observées, provenant, l'une de l'intestin de Triton, l'autre des crottes du Lapin domestique. La cellule piriforme porte un noyau antérieur et 4 flagelles insérts près du pôle ventral du noyau. Le flagelle récurrent (R) est logé dans une dépression ventrale ou cytostome. Les cinétosomes, disposés en une paire antérieure (#1, #2) et une paire postérieure (#3, R), sont liés entre eux par des microfibrilles. Une fibre microtubulaire située au-dessus du noyau est reliée au cinétosome #1. Une autre fibre microtubulaire sous-nucléaire est homologue de la fibre microtubulaire croisée qui existe chez les cellules de Diplozoaires. Le cytostome est bordé par 2 lèvres: la gauche proéminente et armée par plusieurs rangées de microtubules, la droite contenant seulement une mince fibre microtubulaire associée à des microfibrilles. Le cytostome occupe les 2/3 de la face ventrale. Le flagelle récurrent pénètre dans le cytostome puis dépasse l'extrémite de la cellule. Les Bactéries sont phagocytées au fond du cytostome, entre les 2 lèvres distendues. Elles sont digérées dans les nombreuses vacuoles et les corps résiduels sont évacués par rupture de la membrane cellulaire. L'ergastoplasme est concentré près de la périphérie de la cellule. Il n'y a pas de mitochondrie ni d'appareil de Golgi. Dans les kystes observés la cellule plurinucléée est enfermée dans une enveloppe kystique microfibrillaire, les axonèmes sont libres dans le cytoplasme. Les formes diplomonades sont nombreuses et ressemblent aux cellules d'Hexamita, excepté par le cytostome qui est différent. Dans ces formes, les 2 monades sont souvent disposées selon une symétrie axiale binaire mais quelquefois elles sont associées de façon plus anarchique. La cinétide d'Enteromonas est organisée comme celle d'un zoïde de Diplozoaire. Il est possible que le genre Enteromonas soit à l'origine des Diplomonadida et que l'état diplomonadien transitoire chez Enteromonas se soit stabilisé ensuite chez les Diplomonadida. Enteromonas apparaît plus primitif que les autres genres de Diplomonadida aussi nous proposons de créer 2 sous-ordres: celui des Enteromonadina avec le genre Enteromonas et celui des Diplomonadina avec les genres Trepomonas, Trigonomonas, Hexamita, Spironucleus, Octomitus, Giardia. La disposition des cinétosomes et l'existence du cytostome sont les principaux caractères communs entre Enteromonas et les Retortamonadida, cependant les fibres annexes ne sont pas homologues. Une étude plus complète de la division nucléaire et cellulaire de ces 2 ordres de Zooflagellés est nécessaire pour donner un meilleur schéma évolutif. SYNOPSIS. Fine structure of 2 species of Enteromonas, one from the intestine of the salamander, Triturus vulgaris, and another from the feces of domestic rabbit, Oryctolagus cuniculi, is described. The pyriform cell has an anteriorly located nucleus. The 4 flagella originate from an area near the anterior end of the nucleus. The recurrent flagellum (R) is lodged in a ventral depression or cytostome. The kinetosomes, arranged into 2 pairs, anterior (#1, #2) and posterior (#3, R), are interconnected by microfibrils. One microtubular fiber, connected to kinetosome #1, is situated near the anterior surface of the nucleus. Another, subnuclear, microtubular fiber is homologous to the “crossed'’fiber found in Diplozoa. The cytostome is bordered by 2 lips: the preeminent left lip is equipped with several rows of microtubules, while the right lip contains only a thin microtubular fiber associated with microfibrils. The cytostome occupies 2/3 of the ventral surface. The recurrent flagellum passes over the anterior surface of the cell and then comes to lie in the cytostome. The bacteria are phagocytosed in the bottom part of the cytostome between the 2 distended lips. They are digested in numerous vacuoles. The undigested residual bodies are evacuated by a rupture of the cell membrane. The ergastoplasm is concentrated near the cell periphery. Mitochondria and the Golgi apparatus are absent. In the cyst stage, the multinucleate cell is enclosed in a microfibrillar membrane; the axonemes lie free in the cytoplasm. Diplomonad forms of Enteromonas resembling Hexamita are numerous, except that the cytostome is different in these 2 genera. In such forms, the arrangement of the 2 individuals often has binary axial symmetry, but on occasion they are associated in a more anarchic fashion. The mastigont of Enteromonas is organized like that of a single zooid of a diplozoon. It is possible that the genus Enteromonas is ancestral to Diplomonadida and that the diplomonad state, transitory in Enteromonas, became permanently established in Diplomonadida. Enteromonas appears to be more primitive than the other genera of Diplomonadida. Thus we propose 2 suborders: Enteromonadina, subord. nov. with the genus Enteromonas, and Diplomonadina Wenyon, emend., with the genera Trepomonas, Trigonomonas, Hexamita, Spironucleus, Octomitus, Giardia. The arrangement of the kinetosomes and the existence of a cytostome are the principal characters common to Enteromonas and Retortamonadida, while their “accessory'’fibers are not homologous. A more complete study of division of the 2 zooflagellate orders is necessary for the presentation of a more detailed evolutionary scheme of these groups.     

Adaptation d'Oecophylla longinoda [Formicidae-Formicinae] aux variations spatio-temporelles de la densité de proies

Résumé L'activité prédatrice de 12 sociétés d'Oecophylla longinoda Latreille a été étudiée durant 124 h d'observations pendant lesquelles 2 482 proies ont été capturées. La prédation s'exerce sur les arbres supportant les nids ainsi que sur des territoires ?secondaires? habituellement situés sur le sol. Les 7 situations suivantes sont comparées: 1) En forêt, durant les périodes ?habituelles?; 2) Durant le passage d'une colonne de Magnans; 3) Le territoire de chasse ?secondaire? d'une société installée sur un arbre situé en bordure de rivière correspond à une zone où de nombreux insectes se concentrent pour s'abreuver; 4) L'arbre support est un manguier dont on suit le cycle biologique (inter-saison, floraison, fruits m?rs) et son impact sur la prédation par lesOecophylla; 5) On compare la situation ?habituelle? à ce qui se passe lorsqu'un essaimage de Termites est rabattu sur les territoires d'une société; 6) Le territoire ?secondaire? est situé sous la véranda d'une maison éclairée chaque nuit où de nombreux insectes ailés sont attirés; 7) Les territoires ?secondaires? de 2 sociétés étant contigus, on compare les périodes de combats aux périodes calmes. Les territoires ?secondaires? sont situés sur des zones de concentration en proies (insectes de la litière des forêts, zones où ils viennent s'abreuver, attraction par la lumière, la nuit). LesOecophylla y capturent souvent des insectes aptères, cas relativement rare sur l'arbre support. Ces territoires sont défendus contre les ouvrières d'autres sociétés. Du fait de leur territorialité très marquée avec surveillance permanente du territoire et un recrutement très efficace, lesOecophylla répondent très rapidement à toute prolifération en proies. A cela s'ajoute une grande diversité de proies potentielles de sorte que cette espèce peut être considérée comme très efficace en lutte biologique, les ouvrières réagissant spécifiquement à toute abondance d'un type d'insecte. Pour une société comprenant 12 calies, le nombre de proies capturées par an est de l'ordre de 45 000. Ramené à un nombre équivalent d'individus et de jours d'activité, il est du même ordre que celui desFormica lugubris (Cherix, 1986).      

Biogéographie et systématique des abeilles rubicoles du genre Ceratina Latreille au Turkestan (Hymenoptera,Apoidea, Xylocopinae)

Les affinités fauniques et la délimitation même du Turkestan en tant qu’entité biogéographique sont mal connues ou controversées. Elles sont discutées ici dans le cadre d’une révision systématique des espèces d’abeilles rubicoles du genre Ceratina. Onze espèces, dont une nouvelle: Ceratina (Euceratina) haladai n. sp., y sont répertoriées. Trois espèces sont redécrites et une clé d’identification est proposée. Toutes ces espèces sont soit présentes dans la région ouest-paléarctique, soit de proches parentes d’espèces méditerranéennes. Aucune espèce ne présente d’affinités notables avec les taxons orientaux ou est-paléarctiques. La carte de distribution de chaque espèce est établie et comparée à la distribution des biomes régionaux. Comme elles sont très largement polylectiques, les plantes butinées ne constituent probablement pas un facteur limitant de leur distribution. A l’inverse, leur nidification rubicole les cantonne à des biotopes riches en ronces (Rubus species) ou en autres plantes à tiges creuses ou médullaires (par ex. Verbascum species). Ceci explique leur implantation plutôt synanthrope en montagne et le long des fleuves.     

BIOMETRIC METHODOLOGY

I. Lobach, R. J. Carroll, C. Spinka, M. H. Gail, and N. Chatterjee 673 Haplotype‐Based Regression Analysis and Inference of Case‐‐Control Studies with Unphased Genotypes and Measurement Errors in Environmental Exposures Il est maintenant établi que le risque de développer une maladie dite complexe est sous l’influence conjointe de facteurs de susceptibilité génétiques, de facteurs d’exposition environnementaux et de leurs interactions potentielles. Chatterjee et caroll (2005) ont développé une méthode efficace de maximum de vraisemblance rétrospectif pour l’analyse des études cas‐témoins qui fait l’hypothèse d’indépendance des effets génétiques et environnementaux et qui autorise une distribution totalement non paramétriques des covariables environnementales. Spinka et al (2005) ont étendu cette approche aux études où certains types d’information, comme les phases haplotypiques, sont manquantes chez certains individus. Nous étendons cette approche aux situations où certaines exposition environnementales sont mesurées avec erreur. A l’aide d’un modèle de régression logistique polychotomique, nous permettons au statut maladie d’avoir K+1 niveaux. Nous proposons d’utiliser une pseudo‐vraisemblance and un algorithme EM associé pour l’estimation des paramètres. Nous démontrons leur consistance et dérivons la matrice de covariance asymptotique des estimations des paramètres lorsque la variance de l’erreur de mesure est connue et lorsqu’elle est estimée à l’aide de réplications. Les inférences avec des corrections de l’erreur de mesure sont compliquées par le fait que le test de Wald à des propriétés non satisfaisantes en présence d’un large nombre d’erreurs de mesure. Les techniques du rapport de vraisemblance sont connues pour être une alternative satisfaisante. Cependant, les tests du rapport de vraisemblance ne sont pas techniquement corrects dans cette situation car la fonction de vraisemblance est fondée sur un modèle incorrect, i.e. un modèle prospectif dans un schéma d’échantillonnage rétrospectif. Nous corrigeons les résultats asymptotiques standards pour prendre en compte le fait que le test du rapport de vraisemblance est fondé sur une fonction apparentéà la vraisemblance. Les performances de la méthode proposée sont illustrées avec des études de simulation dans le cas où l’information génétique est de type haplotypique et que les données manquantes viennent d’une ambigüité de phase lors de la reconstruction de ces haplotypes. Enfin, nous appliquons notre méthode à une étude d’association cas‐témoins entre la prise de calcium et le risque d’adénome colorectal. B. Mukherjee and N. Chatterjee 685 Exploiting Gene‐Environment Independence for Analysis of Case‐‐Control Studies: An Empirical Bayes‐Type Shrinkage Estimator to Trade‐Off between Bias and Efficiency L’analyse prospective habituelle de données de type cas‐témoin par la régression logistique conduit souvent à des estimateurs très imprécis des interactions gènes‐environnement, en raison du faible nombre de cas ou de témoins pour des croisements génotype‐exposition. Par contraste, sous l’hypothèse d’indépendance gène‐environnement, par les méthodes actuelles rétrospectives, y compris l’approche «cas seulement », on peut estimer les paramètres d’interaction de manière beaucoup plus précise, mais on peut avoir un biais important lorsque l’hypothèse sous‐jacente d’indépendance gène‐environnement n’est pas respectée. Dans cet article, nous proposons un nouvel estimateur de type Bayésien empirique pour analyser des données cas‐témoin, permettant d’assouplir l’hypothèse d’indépendance gène‐environnement en s’appuyant sur les données. Dans le cas particulier d’un gène binaire et d’une exposition binaire, la méthode conduit à un estimateur du paramètre log odds‐ratio de l’interaction sous forme simple explicite, correspondant à une moyenne pondérée des estimateurs usuels « cas seulement » et cas‐témoin. Nous décrivons également une méthode générale pour l’obtention du nouvel estimateur et de sa variance dans le schéma rétrospectif au maximum de vraisemblance proposé par Chatterjee et Carroll (2005). Des simulations et des exemples sur données réelles montrent que l’estimateur proposé marque bien l’équilibre entre le biais et l’efficacité selon la réelle association gène‐environnement et la taille d’échantillon pour une étude donnée. Z. Cai, M. Kuroki, J. Pearl, and J. Tian 695 Bounds on Direct Effects in the Presence of Confounded Intermediate Variables Cet article traite le problème de l’estimation de l’effet directement attribuable au traitement (average controlled direct effect : ACDE) dans l’analyse d’un critère final, lorsque des facteurs de confusion non mesurés adviennent entre la mesure d’une variable intermédiaire et ce critère final. En effet, même lorsque ces facteurs de confusion ne perturbent pas l’estimation de l’effet total (qui est alors identifiable), ils rendent non identifiable l’effet direct attribuable au traitement. Kaufman et coll. (2005) ont utilisé un logiciel de programmation linéaire pour trouver les minimum et maximum des valeurs possibles de l’ACDE pour des données numériques particulières. En utilisant la méthode de programmation linéaire symbolique de Balke–Pearl (1997), nous dérivons ici des formules analytiques pour les limites supérieure et inférieure de l’ACDE sous diverses hypothèses de monotonie. Ces valeurs‐limites sont toujours valides. Dans le domaine de la recherche clinique, elles permettent d’évaluer, de façon abordable en termes de calcul, l’effet directement attribuable au traitement. Elles permettent ainsi d’éclairer l’expérimentateur sur le sens dans lequel s’exerce l’effet attribuable au traitement, ainsi que sur la précision des estimations. T. J. VanderWeele 702 The Sign of the Bias of Unmeasured Confounding Les variables de confusion non mesurées constituent un problème fréquent lorsqu’on veut faire de l’inférence dans des études observationnelles. Nous proposons un théorème qui permet, dans certaines circonstances, de tirer des conclusions sur le signe du biais généré par une variable de confusion non mesurée. Plus précisément, il est possible de déterminer le signe du biais quand il existe des relations monotones d’une part entre la variable de confusion non mesurée et le traitement et, d’autre part, entre la variable de confusion non mesurée et le critère de jugement. Les conditions d’application du théorème sont discutées ainsi que les avantages et les inconvénients liés à son utilisation pour déterminer le sens du biais de la confusion non mesurée. M. Zhang, A. A. Tsiatis, and M. Davidian 707 Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates L’objectif premier d’un essai clinique randomisé est de comparer deux traitements ou plus. Par exemple, dans un essai à deux bras avec une réponse continue, on peut s’intéresser à la différence des moyennes; avec plus de deux traitements, la comparaison peut être basée sur des différences appariées. Avec des réponses binaires, des odds ratio appariés ou des log odds ratio peuvent être utilisés. En général, les comparaisons peuvent être basées sur des paramètres ayant un sens, dans un modèle statistique pertinent. Des analyses standards pour l’estimation et les tests dans ce contexte, sont typiquement basées sur les données collectées de la réponse et du traitement seulement. Dans beaucoup d’essais, une information sur une covariable auxiliaire peut également être disponible et il est intéressant d’exploiter ces données pour améliorer l’efficacité des inférences. En prenant une perspective de théorie semi‐paramétrique, nous proposons une approche largement applicable à l’ajustement de covariables auxiliaires pour obtenir des estimateurs plus efficaces et des tests pour les paramètres de traitement dans l’analyse d’essais cliniques randomisés. Des simulations et des applications démontrent la performance des ces méthodes. C. J. Lloyd 716 A New Exact and More Powerful Unconditional Test of No Treatment Effect from Binary Matched Pairs Nous considérons le problème du test de la différence des probabilités du succès de paires binaires couplées. Commençant par trois tests standard inexacts, le paramètre de nuisance est d’abord estimé et ensuite la dépendance résiduelle éliminée par maximisation, produisant ce que j’appelle la p‐valeur de E+M. On montre que cette p‐valeur, basée sur la statistique de McNemar, est numériquement supérieure aux suggestions précédentes, y compris les p‐valeurs partiellement maximisées telles qu’elles sont décrites dans Berger et Sidik (2003). La dernière méthode cependant peut avoir des avantages du point de vue des calculs pour de grands échantillons. G. Wei and D. E. Schaubel 724 Estimating Cumulative Treatment Effects in the Presence of Nonproportional Hazards Souvent, dans les études médicales portant sur le délai d’apparition d’un événement, l’effet traitement n’est pas constant au cours du temps. Dans le cas du modèle de Cox, la solution la plus fréquente est d’utiliser un modèle qui suppose que soit l’effet traitement est constant par intervalle, soit qu’il varie lentement au cours du temps (le modèle de Cox à risques non proportionnels). Cette approche a au moins deux limitations majeures. Premièrement, il est généralement difficile d’évaluer si la forme du paramètre choisi pour l’effet traitement est correct. Deuxièmement, en présence de risques non proportionnels, les investigateurs sont habituellement plus intéressés par l’effet traitement cumulé plutôt que par l’effet traitement instantané (c’est‐à‐dire déterminer si et quand les fonctions de survie se croisent). Par conséquent, nous proposons un estimateur pour l’effet global du traitement en présence de risques non proportionnels. Notre estimateur est fondé sur le risque de base cumulé estimé sous un modèle de Cox stratifié. Aucune hypothèse n’est nécessaire quant à la forme de la fonction concernant l’absence de proportionnalité. Les propriétés asymptotiques des estimateurs proposés sont déduites, et les propriétés dans le cas d’échantillons de taille finie sont évaluées par des simulations. L’estimation et les intervalles de confiance de l’estimateur peuvent être calculés simultanément. La méthode proposée est appliquée à des données provenant du registre national des rejets d’organes canadien. B. R. Logan, J. P. Klein, and M.‐J. Zhang 733 Comparing Treatments in the Presence of Crossing Survival Curves: An Application to Bone Marrow Transplantation Dans certaines études cliniques comparant des traitements en terme de courbes de survie, les chercheurs peuvent anticiper que les courbes de survie se croiseront à un certain point, aboutissant à s’intéresser à la comparaison de la survie à long terme. Cependant, la simple comparaison des courbes de survie à un point donné peut être inefficace, et l’utilisation d’un test du log‐rank pondéré peut être trop sensible aux différences initiales dans la survie. Nous formulons le problème comme celui de tester les différences des courbes de survie après un temps pré‐spécifié, et proposons une variété de techniques pour tester cette hypothèse. Nous étudions ces méthodes en utilisant des simulations et les illustrons sur une étude comparant la survie pour des transplantations autologues et allogéniques de moelle. R. Song, M. R. Kosorok, and J. Cai 741 Robust Covariate‐Adjusted Log‐Rank Statistics and Corresponding Sample Size Formula for Recurrent Events Data Des évènements récurrents se rencontrent fréquemment dans les essais cliniques. Cet article développe des statistiques robustes du log‐rank ajustées pour des covariables, appliquées à des évènements récurrents en nombre arbitraire avec censure indépendante et la formule de taille d’échantillon correspondante. Les tests du log‐rank proposés sont robustes pour différents procédés de génération de données et ajustés pour des covariables prédictives. Ils se ramènent au cadre développé par Kong et Slud (1997) dans le cas d’un événement unique. La formule de taille d’échantillon repose sur la normalité asymptotique des statistiques du log‐rank ajustées pour des covariables sous certaines alternatives locales et sous un modèle pour les covariables au baseline dans le contexte d’évènements récurrents. Quand la taille de l’effet est petite et quand les covariables au baseline n’apportent pas une information significative sur les temps de survenue, elle se ramène à la même forme que celle de Schoenfeld (1983) dans les cas d’un événement unique ou de temps de survenue indépendants pour un sujet. Nous réalisons des simulations pour étudier le contrôle de l’erreur de type I et pour comparer la puissance de plusieurs méthodes sur des échantillons finis. La formule de taille d’échantillon proposée est illustrée à partir des données d’un essai de la rhDNase. J. D. Nielsen and C. B. Dean 751 Clustered Mixed Nonhomogeneous Poisson Process Spline Models for the Analysis of Recurrent Event Panel Data Un modèle semi paramétrique flexible pour analyser des données de comptage longitudinales provenant de mélanges est présenté. Les données de comptage en panel se réfèrent ici à des données de comptage sur des évènements récurrents collectés comme le nombre d’évènements qui sont arrivés durant une période de suivi spécifique. Le modèle suppose que les comptages pour chaque sujet sont générés par des mélanges de processus de Poisson non homogènes avec des fonctions d’intensité lissées modélisées par des splines pénalisés. Les effets des covariables dépendant du temps sont aussi incorporés dans le modèle dans l’intensité du processus en utilisant les fonctions splines. Des mélanges discrets de ces modèles splines de processus non homogènes de Poisson extraient l’information fonctionnelle provenant des grappes sous‐jacentes représentant les sous populations cachées. La motivation de ce travail est une expérience pour tester l’efficacité de phéromones en interrompant le processus d’accouplement de papillons tortriciae de l’écorce du cerisier. Les papillons matures proviennent de sous populations cachées mais distinctes et on s’intéressait au contrôle des réponses des sous populations. Des effets aléatoires intra‐grappes sont utilisés pour prendre en compte les structures de corrélation et l’hétérogénéité communes à ce genre de données. Une approche pour faire de l’inférence par équations d’estimation, requérant seulement des hypothèses faibles, est développée et les propriétés à distance finie des fonctions d’estimation proposées sont étudiées empiriquement par simulation. R. M. Elashoff, G. Li, and N. Li 762 A Joint Model for Longitudinal Measurements and Survival Data in the Presence of Multiple Failure Types Dans ce papier, nous étudions un modèle conjoint pour des mesures longitudinales et des données de survie à risques compétitifs. Notre modèle conjoint fournit une approche souple pour prendre en compte les données manquantes possiblement non ignorables dans les mesures longitudinales dues aux sorties d’étude. Il s’agit aussi d’une extension des modèles conjoints précédents avec un seul type d’événements, offrant une voie possible pour modéliser les événements censurés informativement comme un risque compétitif. Notre modèle consiste en un sous‐modèle linéaire à effets mixtes pour le critère longitudinal et un sous‐modèle de fragilitéà risques proportionnels cause spécifique (Prentice et al., 1978, Biometrics 34, 541–554) pour les données de survie à risques compétitifs, liés l’un à l’autre par des effets aléatoires latents. Nous proposons d’obtenir les estimations du maximum de vraisemblance des paramètres par un algorithme EM et d’estimer leurs écarts‐types en utilisant une méthode de vraisemblance profilée. La méthode développée fonctionne bien dans nos études de simulation et est appliquée à un essai clinique sur les atteintes pulmonaires de la sclérodermie. S. Goetgeluk and S. Vansteelandt 772 Conditional Generalized Estimating Equations for the Analysis of Clustered and Longitudinal Data Un problème commun et important dans des plans d’échantillonnage en grappes est que l’effet des expositions intra‐grappes (c’est‐à‐dire des expositions qui varient dans les grappes) sur le résultat peut être confondu avec des facteurs mesurés ou non des facteurs du niveau de la grappe (c’est‐à‐dire des mesures qui ne changent pas à l’intérieur des grappes). Quand certains résultats sont malades/non expliqués, l’estimation de cet effet par l’intermédiaire de modèles de moyenne de la population ou de modèles à effets aléatoires peut introduire un biais. Nous adaptons ce fait en développant une théorie générale de l’analyse de données en grappes qui permet une estimation convergente et asymptotiquement normale des effets (CAN) des expositions intra‐grappes en présence des niveaux de confusion du niveau de grappe. Des estimateurs efficaces semi‐paramétriques sont obtenus en résolvant les équations d’estimation conditionnelles généralisées (CGEE). Nous comparons cette approche à une proposition populaire de Neuhaus et Kalbfleisch (1998) qui séparent l’effet d’exposition en composantes intra‐grappes ou entre grappes à l’intérieur d’un modèle aléatoire avec une constante. Nous constatons que cette dernière approche fournit des estimateurs convergents et efficaces quand le modèle est linéaire, mais qu’elle est moins flexible en termes de spécification de modèles. Avec des modèles non‐linéaires, cette approche peut aboutir à des estimateurs non convergents et inefficaces, bien que le biais reste faible dans la majorité des cas pratiques. Y. Ma, W. Tang, C. Feng, and X. M. Tu 781 Inference for Kappas for Longitudinal Study Data: Applications to Sexual Health Research De nombreuses recherches comportementales, médicales, psycho‐sociales ou de santé utilisent des mesures de la fiabilité des instruments ou de la concordance entre juges pour évaluer les propriétés psychométriques de leurs instruments, le niveau de consensus en matière de diagnostic médical, ou la fiabilité d’une variable de substitution. L’indice de concordance le plus fréquemment utilisé pour des réponses catégorielles est le coefficient kappa de Cohen. Aujourd’hui, dans de nombreuses études, on observe des réponses corrélées, ce qui complique ou invalide l’inférence à partir des méthodes antérieures. En outre, les problèmes posés par les observations incomplètes sont particulièrement gênants dans les études longitudinales, de plus en plus nombreuses, et l’absence de méthodes adaptées entrave le développement des recherches. L’approche que nous présentons repose sur une nouvelle classe d’estimateurs des coefficients kappa qui prend en compte les difficultés que l’on rencontre dans les études longitudinales où interviennent plusieurs évaluateurs, en particulier celles liées aux observations incomplètes. Nous illustrons la méthode à partir de données réelles concernant le sexe et la santé. B. J. Reich and J. S. Hodges 790 Modeling Longitudinal Spatial Periodontal Data: A Spatially Adaptive Model with Tools for Specifying Priors and Checking Fit La perte d’attachement (AL) ou longueur sous la racine dentaire qui n’est plus attachée à l’os par le ligament parodontal est une mesure habituelle de la maladie parodontale. Dans cet article, nous développons un modèle spatio‐temporel de surveillance de la progression de l’AL. Notre modèle est une extension de l’a‐priori conditionnel autorégressif (CAR), qui lisse spatialement les estimations vers leurs voisins. Cependant, puisque l’AL présente souvent des accès de grandes valeurs dans l’espace et dans le temps, nous proposons un modèle spatiotemporel CAR non‐stationnaire qui permet au degré de lissage temporel et spatial de varier dans les différentes régions de la bouche. Pour celà, nous attribuons à chaque site de mesure d’AL son propre ensemble de paramètres de variance, et lissons les variances avec des a‐priori spatiaux. Nous suggérons une heuristique pour mesurer la complexité des variances site‐spécifiques, et nous l’utilisons pour sélectionner les a‐priori qui garantissent que les paramètres du modèle sont bien identifiés. Sur des données d’un essai clinique, le modèle améliore la qualité de l’ajustement par rapport au modèle usuel dynamique CAR pour 90 des 99 mesures d’AL. Y. Guan 800 A KPSS Test for Stationarity for Spatial Point Processes Nous proposons une méthode formelle pour tester la stationnarité d’un processus spatial ponctuel. La statistique de test proposée est basée sur l’intégrale des carrés des écarts entre les nombres d’événements observés et leur moyenne sous l’hypothèse de stationnarité. Nous montrons que la statistique de test obtenue converge en loi vers une fonctionnelle d’un mouvement brownien bidimensionnel. Pour appliquer le test, nous comparons la statistique calculée aux valeurs critiques à l’extrémité supérieure de cette fonctionnelle. Notre méthode suppose seulement une condition de dépendance faible pour le processus, mais ne requiert aucun modèle paramétrique pour cela. De la sorte, elle peut être appliquée à une large classe de modèles de processus spatiaux ponctuels. Nous étudions l’efficacité du test à la fois par des simulations et des applications à deux exemples sur données réelles où la non stationnaritéétait préalablement envisagée en se basant sur des aspects graphiques. Notre test confirme cette non stationnarité suspectée dans les deux cas. W. Zucchini, D. Raubenheimer, and I. L. MacDonald 807 Modeling Time Series of Animal Behavior by Means of a Latent‐State Model with Feedback Nous décrivons une famille de modèles adaptés aux séries temporelles de comportement nutritionnel animal. Les modèles incluent un état non observé pouvant être interprété comme l’état de motivation de l’animal, et un mécanisme pour donner cet état en retour du comportement observé. Nous discutons des méthodes d’évaluation et de maximisation de la vraisemblance d’une série observée de comportements, de l’estimation des paramètres, et de l’inférence pour la plus probable séquence d’états sous‐jacents. Nous indiquons plusieurs extensions de ces modèles, incorporant des effets aléatoires. Nous appliquons ces méthodes à l’analyse du comportement nutritionnel des chenilles Helicoverpa armigera, et nous démontrons ainsi le potentiel de cette famille de modèles dans l’étude des comportements. R. King and S. P. Brooks 816 On the Bayesian Estimation of a Closed Population Size in the Presence of Heterogeneity and Model Uncertainty Nous considérons l’estimation de la taille d’une population fermée, qui présente souvent un intérêt pour des populations animales sauvages, en utilisant une étude de capture‐recapture. L’estimation de la taille totale de population peut être très sensible au choix du modèle utilisé pour l’ajuster aux données. Nous considérons une approche bayésienne, dans laquelle nous considérons chacun des huit modèles plausibles originellement décrits par Otis et al. (1978) dans un contexte simple, incluant des modèles avec une composante individuelle d’hétérogénéité. Nous montrons comment nous pouvons obtenir une estimation moyenne de la population totale, incluant simultanément le paramètre et l’incertitude du modèle. Pour illustrer la méthodologie nous effectuons d’abord une étude de simulation et nous analysons deux ensembles de données où la taille de population est connue, avant de considérer un exemple réel concernant une population de dauphins du nord‐est de l’Ecosse. A. J. Branscum and T. E. Hanson 825 Bayesian Nonparametric Meta‐Analysis Using Polya Tree Mixture Models Un objectif habituel dans une méta‐analyse est l’estimation d’une mesure d’effet unique en utilisant les données de plusieurs études qui sont toutes dessinées pour traiter la même question scientifique. Parce que les études sont typiquement réalisées dans des lieux géographiquement dispersés, les développements récents dans l’analyse statistique des données de méta‐analyses impliquent l’utilisation de modèles avec effets aléatoires qui tiennent compte de la variabilité entre études attribuable aux différences d’environnement, de démographie, de génétique, et d’autres sources d’hétérogeneïté entre les populations. En se fondant sur la théorie asymptotique, les statistiques résumées propres à chaque étude sont modélisées selon des distributions normales représentant les vraies mesures d’effets latents. Dans un deuxième temps une approche paramétrique modélise ces mesures latentes au moyen d’une distribution normale, qui est uniquement une supposition commode dénuée de justification théorique. Pour éliminer l’influence de modèles paramétriques par trop restrictifs sur l’inférence, nous considérons une classe plus large de distributions pour les effets aléatoires. Nous développons un nouveau modèle de mélange d’arbres de Polya Bayesien non paramétrique et hiérarchique. Nous présentons une méthodologie pour tester un modèle de mélange d’arbre de Polya contre un modèle avec effets aléatoires normaux. Ces méthodes fournissent aux chercheurs une approche immédiate pour conduire une analyse de sensibilité de la normalité des effets aléatoires. Une application impliquant une méta‐analyse d’études épidémiologiques pour caractériser l’association entre la consommation d’alcool et le cancer du sein est présentée qui, alliée à des résultats issus de données simulées, met en lumière la performance des mélanges d’arbre de Polya en présence de non‐normalité des mesures d’effet dans la population d’origine. G. Massonnet, P. Janssen, and T. Burzykowski 834 Fitting Conditional Survival Models to Meta‐Analytic Data by Using a Transformation Toward Mixed‐Effects Models Les modèles à fragilité sont communément utilisés pour modéliser des survies sur données en ‘clusters’. Les approches classiques pour ajuster des modèles de fragilité sont basées sur la vraisemblance. Nous proposons une approche alternative dans laquelle le problème initial d'« ajustement d’un modèle de fragilité» est reformulé en le problème d’« ajustement d’un modèle linéaire mixte » grâce à l’utilisation d’un modèle de transformation. Nous montrons que le principe de transformation s’applique aussi aux modèles multivariés à odds proportionnels ainsi qu’aux modèles de risques additifs multivariés. Il permet de plus de faire le lien entre des méthodes séparées en ce sens qu’il fournit un moyen général d’ajuster des modèles conditionnels pour des données de survie multivariées en utilisant la méthodologie des modèles mixtes. Nous nous centrons sur les modèles à fragilité afin d’étudier les caractéristiques spécifiques de la méthode proposée. En nous basant sur une étude de simulation nous montrons que la méthode proposée produit une alternative correcte et simple pour l’ajustement de modèles à fragilité sur des données présentant un nombre suffisamment grand de classes et des tailles d’échantillons modérées à grande à l’intérieur de niveaux de sous‐groupes de covariables. La méthode proposée est appliquée à des données provenant de 27 essais randomisés sur des cancers colorectaux avancés, disponibles par l’intermédiaire du groupe de Meta‐Analyse du Cancer. R. A. Hubbard, L. Y. T. Inoue, and J. R. Fann 843 Modeling Nonhomogeneous Markov Processes via Time Transformation Les études longitudinales constituent un outil puissant pour représenter la progression des maladies chroniques. Ces études sont habituellement menées avec des patients examinés lors de visites périodiques, donnant ainsi des données de type panel. Dans ce schéma d’observation, les temps exacts de transitions entre les états de la maladie, et la séquence suivie des états de la maladie sont inconnus, et les modèles à base de processus de Markov sont souvent utilisés pour décrire la progression. La plupart des applications reposant sur les processus markoviens supposent l’homogénéité par rapport au temps, autrement dit que les taux de transition sont constants au cours du temps. Cette hypothèse n’est pas satisfaite lorsque les taux de transition dépendent du temps à compter du début de la maladie. Cependant, peu d’outils statistiques existent pour prendre en compte cette non‐homogénéité. Nous proposons des modèles dans lesquels l’échelle de temps d’un processus de Markov non‐homogène est transformée en une échelle de temps opérationnelle par rapport à laquelle le processus est homogène. Nous développons une méthode pour l’estimation jointe de la transformation de temps et de la matrice d’intensité de transition pour le processus homogène transformé. Nous établissons une estimation au maximum de vraisemblance utilisant l’algorithme du score de Fisher par des études de simulation, et nous comparons la performance de notre méthode par rapport à des modèles homogènes et des modèles homogènes par morceaux. Nous appliquons notre méthodologie à une étude sur la progression du délire dans une cohorte de receveurs de transplantation de cellules souche, et nous montrons que notre méthode identifie les tendances temporelles dans l’incidence du délire et dans la récupération. C. C. Drovandi and A. N. Pettitt 851 Multivariate Markov Process Models for the Transmission of Methicillin‐Resistant Staphylococcus Aureus in a Hospital Ward L’Aureus Staphylococcus résistant à la méthicilline (SARM) est un pathogène qui continue d’être un problème majeur à l’hôpital. Nous développons des modèles et des algorithmes de calcul basés sur des données observées d’incidences hebdomadaires pour estimer les paramètres de transmission du MRSA. Nous étendons le modèle déterministe de McBryde, Pettitt et McElwain (2007) impliquant une population sous‐jacente de patients colonisés par MRSA et de personnels soignants qui décrit, parmi d’autres processus, la transmission entre patients non colonisés et personnels soignants colonisés et vice‐versa. Nous développons de nouveaux modèles de Markov bivariés et trivariés pour inclure l’incidence de telle façon que les taux de transmission estimés puissent être basés directement sur de nouvelles colonisations plutôt qu’indirectement sur la prévalence. La sensibilité imparfaite de la détection du pathogène est modélisée en utilisant un processus markovien caché. Les avantages de notre approche incluent (i) l’hypothèse d’une valeur discrète pour le nombre de personnels soignants colonisés, (ii) deux paramètres de transmission pouvant être incorporés dans la vraisemblance, (iii) la vraisemblance dépend du nombre de nouveaux cas pour améliorer la précision de l’inférence, (iv) les enregistrements individuels des patients ne sont pas requis, et (v) la possibilité d’une détection imparfaite de la colonisation est incorporée. Nous comparons notre approche avec celle utilisée par McBryde et al (2007), basée sur une approximation qui élimine les personnels soignants du modèle, utilisant une chaîne de Markov par Monte Carlo et des données individuelles des patients. Nous appliquons ces modèles à des données de colonisation de SARM, collectées dans de petites unités de soins intensifs (ICU) à l’hôpital Princesse Alexandra, Brisbane, Australie. A. R. Cook, G. J. Gibson, and C. A. Gilligan 860 Optimal Observation Times in Experimental Epidemic Processes Cet article décrit une méthode de choix de temps d’observation de processus stochastiques destinés à maximiser les espérances des informations de leurs paramètres. On considère deux modèles généralement utilisés pour des processus épidémiologiques: un processus simple de mort et un processus épidémique de susceptibilité d’infection (SI) avec des sources doubles d’infection se développant de l’intérieur ou de l’extérieur de la population. La recherche du plan optimal emploie des méthodes bayésiennes de calcul pour explorer l’espace expérimental commun paramètre‐donnée, combinéà une méthode connue sous le nom de fermeture de moment d’approximation de la vraisemblance pour rendre efficace l’étape d’acceptation. Pour les processus considérés, on montre que choisir un petit nombre d’observations de façon optimale est presqu’aussi bon que collecter des observations de manière plus intensive, pratique générale dans nombre d’expériences épidémiologiques. L’analyse du processus simple de mort permet une comparaison entre l’approche bayésienne intégrale et des dispositifs localement optimaux autour d’une estimation ponctuelle d’a priori basée sur des résultats asymptotiques. On démontre la robustesse de l’approche pour des a priori mal spécifiés du processus épidémique SI, pour lequel le problème de la faisabilité informatique du calcul de la vraisemblance exclut des plans optimaux. Nous montrons que des plans optimaux déduits d’une approche bayésienne sont semblables à des études d’observations d’une épidémie simple et à des études impliquant des épidémies répétées dans des sous‐populations indépendantes. Un optimum différent existe, cependant, quand l’objectif est de maximiser le gain de l’information basée sur des a priori informatifs ou non: ceci a des implications quand une expérience est conçue pour convaincre un observateur naïf ou sceptique plutôt que pour consolider la connaissance d’un observateur mieux informé. On présente quelques prolongements de ces méthodes, comprenant le choix de critères d’information et l’extension à d’autres processus épidémiologiques avec des probabilités de transition. K. F. Lam, J. V. Deshpande, E. H. Y. Lau, U. V. Naik‐Nimbalkar, P. S. F. Yip, and Ying Xu 869 A Test for Constant Fatality Rate of an Emerging Epidemic: With Applications to Severe Acute Respiratory Syndrome in Hong Kong and Beijing L’étiologie, la pathogénie et le pronostic d’une maladie récemment émergente sont généralement inconnus aux cliniciens. Des interventions et des traitements efficaces aussi précoces que possible garantissent la réduction du caractère fatal de la maladie à un minimum, et la suppression de traitements inadéquats. Dans ce contexte, la capacitéà extraire le plus d’informations hors des données disponibles est critique à la possibilité d’importantes prises de décision. L’inefficacité d’un traitement peut être reflétée par une mortalité constante avec le temps tandis qu’un traitement efficace mène normalement à un taux de décès décroissant. Un test statistique pour une mortalité constante avec le temps est proposé dans cet article. La statistique proposée converge vers un mouvement Brownien sous l’hypothèse nulle. Les caractéristiques spécifiques du mouvement Brownien permettent d’analyser également la distribution des premiers temps de passage à partir de tests séquentiels. Ceci permet à l’hypothèse nulle du taux de décès constant d’être rejetée le plus rapidement possible dès que l’évidence statistique devient suffisante. Les études de simulation montrent que les performances du test proposé sont bonnes et qu’il est très sensible pour trouver un taux de mortalité décroissant. Le test proposé est appliqué aux données du syndrome respiratoire aigu sévère (SARS) à Hong Kong et Beijing. J. Peng, C.‐I. C. Lee, K. A. Davis, and W. Wang 877 Stepwise Confidence Intervals for Monotone Dose–Response Studies Dans les études de dose‐réponse un des points les plus importants est l’identification de la dose minimum efficace (minimum effective dose – MED), qui est définie comme la plus petite dose pour laquelle la réponse moyenne est meilleure que celle obtenue pour une dose nulle de contrôle, en terme de différence clinique significative. Les courbes doses‐réponses sont parfois monotones. Pour trouver la MED plusieurs auteurs ont proposé des procédures de test du genre up and down, basées sur des contrastes entre les moyennes. Dans cet article, pour identifier la MED, on propose une amélioration de la méthode de Marcus et Peritz (1976, Journal of the Royal Statistical Society, Series B 38, 157–161) et on développe une application de la méthode DR de Hsu et Berger (1999, Journal of the American Statistical Society 94, 468–482) pour construire, sous des hypothèses de monotonicité, la borne inférieure de l’intervalle de confiance pour la différence entre la réponse à une dose quelconque non nulle et celle du contrôle. La méthode proposée est illustrée par des exemples numériques et on présente des comparaisons de puissance obtenues par simulations. M. Ding, G. L. Rosner, and P. Müller 886 Bayesian Optimal Design for Phase II Screening Trials Dans la littérature, la plupart des plans de screening en phase 2 considèrent un traitement à la fois. Chaque étude est considérée isolément. Nous proposons une approche décisionnelle plus systématique pour le processus de screening de phase 2. Le plan séquentiel permet une plus grande efficacité et un meilleur apprentissage sur les traitements. L’approche incorpore un modèle Bayesien hiérarchique qui permet d’une façon formelle la combinaison d’information entre différentes études reliées entre elles, et améliore l’estimation dans de petits échantillons en exploitant l’information venant des autres traitements. Le plan incorpore une fonction d’utilité qui prend en compte le coût d’échantillonnage et le possible retour sur investissement. Des simulations par ordinateur montrent que cette méthode a une forte probabilité d’écarter des traitements à faible taux de succès et de faire évoluer en phase 3 des traitements à fort taux de succès. A. Vexler, A. Liu, E. Eliseeva, and E. F. Schisterman 895 Maximum Likelihood Ratio Tests for Comparing the Discriminatory Ability of Biomarkers Subject to Limit of Detection Dans cet article, nous étudions la comparaison des aires sous courbes ROC (Receiver Operating Characteristic) corrélées, des biomarqueurs diagnostiques dont les mesures sont soumises à une limite de détection (LOD), source d’erreur de mesure issue de la sensibilité instrumentale en épidémiologie. Nous proposons et nous examinons les tests du rapport des vraisemblances avec des caractéristiques opérationnelles qui sont aisément obtenues sous méthodologie classique du maximum de vraisemblance. R. Li and L. Nie 904 Efficient Statistical Inference Procedures for Partially Nonlinear Models and their Applications Motivés par une analyse d’un jeu de données réelles dans le domaine de l’écologie, nous sommes amenés à considérer une classe de modèles partiellement non linéaires où l’on note à la fois la présence d’une composante non paramétrique et d’une composante paramétrique. Nous développons deux nouvelles procédures d’estimation pour l’estimation des paramètres de la composante paramétrique. La consistance et la normalité asymptotique des estimateurs sont établies. Puis nous proposons ensuite une procédure d’estimation et une procédure de test F généralisé pour la composante non paramétrique dans les modèles partiellement non linéaires. Des propriétés asymptotiques de la nouvelle procédure d’estimation et le test statistique s’en déduisent. Les qualités des échantillons finis des procédures d’inférence proposées sont évaluées par des études de simulations de Monte Carlo. Une application dans le domaine de l’écologie est utilisée pour illustrer la méthode proposée E. Deschepper, O. Thas, and J. P. Ottoy 912 Tests and Diagnostic Plots for Detecting Lack‐of‐Fit for Circular‐Linear Regression Models Les diagnostics pour la régression et les tests de défaut d’ajustement sont principalement focalisés sur les modèles linéaires de régression linéaire. Lorsque les points sont distribués sur la circonférence d’un cercle, des difficultés apparaissent car il n’existe pas de point de départ ou d’origine. La plupart des tests classiques de défauts d’ajustement exigent de choisir arbitrairement une origine, mais des choix différents peuvent conduire à des conclusions différentes. Nous proposons un outil graphique d’ajustement et un test de défaut d’ajustement qui n’exige pas de point naturel initial. La méthode est basée sur des « résidus régionaux » qui sont définis sur des arcs de cercle. La méthode graphique localise formellement et visualise des sous ensembles d’observations dont l’ajustement n’est pas satisfaisant sur le cercle. Un exemple de données issues du domaine alimentaire est utilisé pour exhiber les problèmes mentionnés auparavant avec les tests classiques d’ajustement et pour illustrer la force de la méthode basée sur les « résidus régionaux » pour détecter et localiser les points de départ d’hypothèse de non effet. Une petite étude de simulation montre la bonne performance du test des « résidus régionaux » dans les deux cas d’écart global ou local du modèle. Puis ces idées sont étendues au cas ou il y a plus d’une variable explicative. B. Xie, W. Pan, and X. Shen 921 Variable Selection in Penalized Model‐Based Clustering Via Regularization on Grouped Parameters L’analyse en grappes modèle‐dépendante pénalisée a été proposée dans le cadre de données de haute dimension mais de petite taille d’échantillon comme celles rencontrées dans les études génomiques. En particulier, elle peut être utilisée pour la sélection de variables. Un nouveau schéma de régularisation est proposé pour regrouper ensemble plusieurs paramètres associés à la même variable au sein des grappes. Nous montrons analytiquement et numériquement que ce schéma est plus efficace que la conventionnelle pénalité L₁ pour la sélection de variable. De plus, nous développons une stratégie pour combiner ce schéma de groupement avec le groupement de variables structurées. Des études de simulation et des applications aux données d’expression des micropuces pour la découverte de sous‐types cancéreux illustrent les avantages de notre nouvelle approche par rapport à plusieurs méthodes existantes. S. H. Lee, J. Lim, M. Vannucci, E. Petkova, M. Preter, and D. F. Klein 931 Order‐Preserving Dimension Reduction Procedure for the Dominance of Two Mean Curves with Application to Tidal Volume Curves Le travail présenté ici est issu d’une étude de cas reposant sur des tracés, à dimension et fréquence élevées, de capacité vitale respiratoire effectués pendant des bouffées d’angoisse provoquées. L’objectif est de développer une procédure pour évaluer la significativité de la supériorité d’une courbe moyenne par rapport à une autre. L’idée essentielle de la méthode suggérée est de préserver l’ordre sur les moyennes tout en réduisant la dimensionnalité des données. La matrice des données observées est projetée sur un ensemble de matrices de rang inférieur sous une contrainte de positivité. Une procédure de test multivarié est alors appliquée dans cette dimension inférieure. Nous utilisons des données simulées pour illustrer les propriétés statistiques de la méthode de test proposée. Les résultats obtenus sur l’étude de cas confirment l’hypothèse préliminaire des investigateurs et fournit un support à leur objectif de construction d’un modèle expérimental pour la bouffée d’angoisse clinique chez les sujets normaux.     

Un difficile problème de taxonomie : les Micrasema (Trichoptera : Brachycentridae) des eaux courantes de la Péninsule Ibérique et des Pyrénées

Un des problèmes taxonomiques les plus compliqués concernant les Trichoptères européens est celui de la plupart des espèces ibériques et pyrénéennes du genre Micrasema, habitants caractéristiques du Rhitral, ces espèces restant à ce jour imparfaitement connues et difficiles à déterminer. L´étude d’un grand nombre d´exemplaires (mâles) a permis de trouver une solution pour la plupart des problèmes. M. servatum (Navás 1918), distribuée dans les parties N.O. de la Péninsule; M. cenerentola Schmid 1952, connue du Centre et du N.O. de l´Espagne; et M. salardum Schmid 1952, connue des Pyrénées Orientales espagnoles, sont des espèces distinctes. M. gabusi Schmid 1952 décrite de la Sierra Nevada, est un synonyme de M. longulum McLachlan 1876, une espèce largement distribuée en Espagne. M. difficile Mosely 1933, décrite des Pyrénées françaises, est un synonyme de M. vestitum Navás 1918, espèce endémique des Pyrénées Centrales et Occidentales (France et Espagne). L´identité de M. moestum (Hagen 1868) pose un problème particulier : il est possible que sous ce nom on désigne actuellement une paire d´espèces cryptiques à l´état adulte, qu´il est pratiquement impossible de distinguer par leurs genitalia, mais dont les larves sont bien distinctes.     

BIOMETRIC PRACTICE

Y. K. Cheung 940 Simple Sequential Boundaries for Treatment Selection in Multi‐Armed Randomized Clinical Trials with a Control Dans une situation où, avec des ressources trop limitées pour réaliser des comparaisons systématiques avec un standard, on veut étudier un grand nombre de traitements, il peut être utile de conduire un essai à plusieurs bras pour éliminer précocement les traitements inintéressants. Quand l’étude a des critères de jugement dont l’évaluation est relativement rapide, comme l’évaluation d’un volume de lésion à partir d’une image chez les patients victimes d’une attaque, des analyses intermédiaires fréquentes sont éthiquement et pratiquement souhaitables pour les cliniciens. Dans cet article je propose une classe de frontières pour les essais à plusieurs bras, dans lesquels l’objectif est soit la sélection d’un traitement donnant une amélioration significative par rapport au groupe contrôle, soit la déclaration de la futilité de la recherche si il n’existe pas de tel traitement. Les limites proposées sont faciles àétablir en aveugle et peuvent être appliquées à un calendrier d’examens évolutif. Il est simple de calibrer l’essai par rapport à un niveau de confiance fixéà l’avance et cela peut être fait quand le taux de réponse dans le groupe contrôle n’est que connu que par intervalle. Une des méthodes proposées est appliquée pour redéfinir un essai de sélection avec un critère de jugement fourni par imagerie chez des patients avec une attaque aiguë, on le compare par simulations à un dispositif optimal à deux étapes. En moyenne, la méthode proposée emploie des échantillons de taille plus faible que la méthode en deux étapes, cet avantage devient substantiel quand il existe un traitement supérieur au contrôle. L. Liu, X. Huang, and J. O'Quigley 950 Analysis of Longitudinal Data in the Presence of Informative Observational Times and a Dependent Terminal Event, with Application to Medical Cost Data Dans les études observationnelles longitudinales, les mesures répétées sont souvent recueillies à des temps d’observation informatifs. De plus, il peut exister un événement terminal dépendant tel que le décès qui interrompt le suivi. Par exemple les patients en moins bonne santé chercheront plus à se faire soigner et les coûts pour chaque visite ont tendance àêtre plus élevés. Ils ont aussi des taux de mortalité plus élevés. Dans cet article, nous proposons un modèle à effets aléatoires de mesures répétées en présence à la fois de délais d’observation informatifs et d’un événement terminal dépendant. Trois sous‐modèles sont utilisés respectivement pour (1) l’intensité des temps d’observation récurrents, (2) la quantité d’observations répétées à chaque temps d’observation et (3) le risque instantané de décès. Des effets aléatoires corrélés sont inclus pour lier les trois sous modèles. L’estimation peut s’effectuer facilement par des techniques de quadrature gaussiennes, par exemple avec la procédure NLMIXED de SAS. La méthode proposée est illustrée par une analyse du processus de coûts générés chez des patients souffrant d’insuffisance cardiaque chronique. Les données sont extraites de la base de données cliniques CDR du système de santé de l’Université de Virginie. J. Cao, G. F. Fussmann, and J. O. Ramsay 959 Estimating a Predator‐Prey Dynamical Model with the Parameter Cascades Method Les équations différentielles ordinaires (EDO) sont largement utilisées pour décrire le comportement dynamique des populations en interaction, Cependant les systèmes EDO donnent rarement des solutions quantitatives proches des observations réelles de terrain ou des résultats expérimentaux, parce que les systèmes naturels sont soumis à un bruit environnemental ou démographique et que les écologistes sont dans l’incertitude pour paramétrer correctement leurs données. Dans cet article nous introduisons les « cascades de paramètres » pour améliorer l’estimation des paramètres EDO de manière à bien ajuster les solutions aux données réelles. Cette méthode est basée sur le lissage pénalisé modifié en définissant la pénalité par les EDO et en généralisant l’estimation profilée pour aboutir à une estimation rapide et à une bonne précision des paramètres EDO sur des données affectées de bruit. Cette méthode est appliquée à un ensemble d’EDO prévues à l’origine pour décrire un système prédateur‐proie expérimental soumis à des oscillations. Le nouveau paramétrage améliore considérablement l’ajustement du modèle EDO aux données expérimentales. La méthode révèle en outre que d’importantes hypothèses structurales sous‐jacentes au modèle EDO originel sont essentiellement correctes. La formulation mathématique des deux termes d’interaction non linéaire (réponses fonctionnelles) qui lient les EDO dans le modèle prédateur‐proie est validée par un estimateur non paramétrique des réponses fonctionnelles à partir des données réelles. Nous suggérons deux applications importantes des « cascades de paramètres »à la modélisation écologique. La méthode peut servir à estimer des paramètres quand les données sont soumises à des bruits ou manquantes, ou quand aucune estimée a priori fiable n’est disponible. La méthode peut par ailleurs aider à valider la qualité structurelle du modèle mathématique. D. S. Johnson, D. L. Thomas, J. M. Ver Hoef, and A. Christ 968 A General Framework for the Analysis of Animal Resource Selection from Telemetry Data Nous proposons un cadre général pour analyser des données télémétriques animales en utilisant des distributions pondérées. On montre que plusieurs interprétations du choix des ressources se présentent quand on compare la ressource utilisée à la ressource disponible. Le cadre général proposé permet de montrer que les modèles habituels de choix des ressources sont des cas particuliers du modèle général avec des hypothèses sur les déplacements et le comportement des animaux. On montre que le cadre des distributions pondérées permet de prendre facilement en compte des données télémétriques fortement autocorrélées, comme c’est typiquement le cas avec la technologie GPS assurant le suivi des animaux. Une analyse de données simulées utilisant plusieurs modèles construits dans le cadre proposé est présentée pour illustrer les gains possibles apportés par cette modélisation flexible. Le modèle est appliquéà des données sur l’ours brun du sud‐est de l’Alaska.     

Gestion informatisée de la biologie délocalisée

La prise en charge par le biologiste de certaines analyses délocalisées se trouve grandement facilitée par la mise en place d’un système informatique reliant l’analyseur délocalisé et le laboratorie. Elle apporte un niveau de sécurité nécessaire pour la prise de responsabilité du laboratoire dans l’acte de biologie délocalisé. L’architecture réseau de la plupart de h?pitaux permet maintenant de transmettre très facilement et de manière fiable des données entre plusieurs services voire entre établissements. Les analyseurs de gazométrie sanguine permettent cette connexion, de plus en plus d’industriels développent d’autres paramètres en plus des simples gaz du sang comme dernièrement l’urée et la créatinine. Les analyseurs permettant de doser les marqueurs cardiaques, ceux réalisant l’hémostase sont désormais fréquemment connectables. Chaque industriel du diagnostic proposant des analyseurs de glycémie a ou va prochainement commercialiser un modèle connectable. Il en est de même pour les bandelettes urinaires. Il est donc probable, et les biologistes sont les premiers demandeurs, que la connexion informatique d’analyseurs de biologie installés dans les unités de soins, soit un pré-requis pour la gestion des analyses délocalisées par le biologiste, et sa prise de responsabilité dans les résultats biologiques issus de ces analyseurs.     

Catalogue et affinités géographiques des Odonata des îles voisines de Madagascar (Insecta: Pterygota)

Une analyse bibliographique exhaustive révèle que 71 espèces et sous-espèces d’Odonates sont présentes sur l’ensemble des îles voisines de Madagascar appartenant aux archipels des Mascareignes, des Comores et des Seychelles. Le taux d’endémisme maximal pour l’ensemble des îles voisines de Madagascar est de 34 %. Le taux d’endémisme est sensiblement plus important dans les Mascareignes (34 %) qu’aux Comores (23 %) ou aux Seychelles (19 %). Le peuplement odonatologique y est dominé par les Libellulidae (36 taxons) et les Coenagrionidae (18 taxons). Chaque archipel possède au moins un des 23 taxons endémiques présents sur ces îles. Les Seychelles comptent deux genres monospécifiques endémiques: Allolestes et Leptocnemis. Les Mascareignes abritent deux genres endémiques: Thalassothemis et Coenagriocnemis. Ce dernier a subi une radiation évolutive au sein de l’archipel. Les îles voisines de Madagascar se caractérisent par une radiation évolutive des genres Hemicordulia et Gynacantha à l’échelle de la région. Les affinités géographiques sont principalement afro-tropicales. La colonisation à partir de l’Afrique est le fait de taxons particuliers qui ne sont pas un échantillon aléatoire de l’ensemble du pool continental. Sur l’ensemble des îles voisines de Madagascar, la richesse spécifique des peuplements d’Odonates d’affinité africaine est maximale sur les îles proches du continent.     

Rôle des abeilles sauvages et domestiques (Hymenoptera : Apoidea) dans la pollinisation de la fève (Vicia faba L. var. major) (Fabaceae) en région de Constantine (Algérie)

Résumé. Les insectes pollinisateurs des cultures de Vicia faba L. en région de Constantine (Algérie) ont été étudiés pendant trois périodes de floraisons (2000, 2001 et 2002). Huit espèces d’abeilles ont été inventoriées. L’abeille sauvage Eucera numida Lepeletier et l’abeille domestique Apis mellifera L. en sont les butineurs les plus abondants. Eucera numida est probablement le pollinisateur le plus efficace car toutes ses visites peuvent être pollinisantes. A l’inverse, un tiers des visites de l’abeille domestique sont consacrées au “vol de nectar” à travers des trous faits par les bourdons à la base des corolles. Ces visites sont non pollinisantes. Les visites d’ Eucera numida sont également plus rapides (m = 10 visites par minute) et son aptitude à déclencher les fleurs de fève est plus grande que celle de l’abeille domestique. Les plants accessibles aux pollinisateurs produisent plus de gousses, les gousses contiennent plus de graines et les graines sont plus grosses et mieux formées que sur les plants encagés.     

REGULAR ARTICLES

L. Cowen and C. J. Schwarz 657 Les Radio‐tags, en raison de leur détectabilitéélevée, sont souvent employés dans les études de capture re‐capture. Une qualité essentielle requise est que les Radio tags ne cessent pas de fonctionner pendant l'étude. Le dysfonctionnement d'un Radio tag avant la fin d'une étude peut conduire à sous‐estimer le taux de survie. Nous développons ici un modèle pour incorporer des données issues de radio tags déficients. Ce modèle a été appliqué aux smolts chinook (tshawytscha d'Oncorhynchus) sur la Colombia river WA. Les estimations des taux de survie obtenus à partir de ce modèle sont beaucoup plus grandes que celles observées dans l'analyse standard de Cormack‐Gaie‐Seber. M. A. Suchard , R. E. Weiss , and J. S. Sinsheimer 665 Les facteurs de Bayes pour comparer deux hypothèses ou plus en compétition sont souvent estimés en construisant un échantillon de chaînes de Markov (MCMC) pour explorer l'espace joint des hypothèses. Pour obtenir des estimations efficaces des facteurs de Bayes, Carlin et Chibb (1995) ont suggéré d'ajuster les odds a priori des hypothèses en compétition de telle sorte que les odds a posteriori soient approximativement un et d'estimer alors les facteurs de Bayes par simple division. Un sous‐produit est que l'on produit souvent plusieurs chaînes MCMC indépendantes, une seulement étant utilisée pour l'estimation. Nous étendons cette approche en proposant trois tests statistiques qui incorporent les résultats des différentes chaînes. Le premier suppose des tirages indépendants dans une chaîne et modélise la fonction indicatrice de l'hypothèse en utilisant une régression logistique pour différents choix des odds a priori. Les deux autres modèles plus complexes abandonnent l'hypothèse d'indépendance en permettant des dépendances retard dans les résultats des différentes chaînes. Ces modèles nous permettent d'estimer l'incertitude dans nos calculs des facteurs de Bayes et d'utiliser pleinement les différentes chaînes MCMC, même si les a priori des odds varient de chaîne à chaîne. Nous appliquons ces méthodes de calcul des facteurs de Bayes pour tester la monophyllie dans deux exemples de phylogénie. Le premier exemple explore la relation d'un pathogène inconnu à un ensemble de pathogènes connus. L'identification de la relation monophyllique du pathogène inconnu peut influer sur le choix de l'antibiotique dans une perspective clinique. Le second exemple se concentre sur la détection des recombinaisons du virus VIH. Pour une application clinique potentielle, ces types d'analyse devraient être achevés aussi efficacement que possible. J. Zhu , J. C. Eickhoff , and P. Yan 674 Des observations de réponses spatio‐temporelles multiples se présentent souvent dans les études environnementales et écologiques. Alors que les modèles purement spatiaux de réponses univariées dans les familles exponentielles de distribution ont progressé au cours des dernières années (e.g. Diggle et al., 1998; Zhang, 2002), peu d'outils statistiques sont disponibles pour des réponse multivariées qui ne sont pas nécessairement Gaussiennes. Une exception cependant est le modèle à facteur commun développé pour des données spatiales multivariées par Wang et Wall (2003). Le but de cet article est d'étendre ce modèle multivarié uniquement spatial et de proposer une classe flexible de modèles linéaires généralisés de variables latentes appliqués à des données multivariées spatio‐temporelles. L'inférence statistique repose sur des estimées du maximum de vraisemblance et de leurs écarts‐type obtenues par l'algorithme EM de Monte Carlo. Nous proposons aussi une nouvelle méthode d'ajustement automatique de la taille de l'échantillon Monte Carlo, qui facilite la convergence de l'algorithme EM et présente ainsi quelque intérêt indépendant. La méthodologie est illustrée par une étude écologique du pin rouge en réponse à des attaques de charançon de l'écorce dans une forêt du Wisconsin. Y.‐G. Wang , X. Lin , and M. Zhu 684 Les méthodes robustes sont utiles pour assurer des inférences statistiques fiables en cas de légers écarts aux hypothèses du modèle. La méthode très largement utilisée des équations d'estimation peut être “robustifiée” en remplaçant les résidus standardisés par les M‐résidus. Si les résidus de Pearson sont supposés centrés sans biais, les estimations des paramètres par l'approche robuste sont asymptotiquement biaisés lorsque les distributions des erreurs ne sont pas symétriques. Nous proposons une méthode pour corriger ce biais. Nos études numériques extensives montrent que la méthode proposée peut réduire ce biais substantiellement. On donne des exemples d'illustration. R. J. Cook , W. Wei , and G. Y. Yi 692 Nous proposons des méthodes semi‐paramétriques pour estimer et tester l'effet d'un traitement en présence d'un événement récurrent observé sur des périodes multiples, avec prise en compte de la censure. Ces méthodes reposent sur l'estimation de fonctions avec pour hypothèse de travail un modèle de Poisson mixte sous lequel les effets aléatoires individuels peuvent être éliminés par conditionnement. Des tests statistiques robustes dérivés de pseudo‐scores sont obtenus par estimation “sandwich” de la variance. L' efficacité relative des analyses conditionnelle et marginale est évaluée de façon analytique sous un modèle mixte de Poisson à temps homogènes. La robustesse et la puissance empirique de l'approche semi‐paramétrique sont étudiées au moyen de simulations. Des adaptations pour prendre en compte les événements récurrents survenant dans des essais en cross‐over sont décrites et ces méthodes sont appliquées aux données d'un essai en cross‐over sur deux périodes chez des patients souffrant d'asthme bronchique. X. Song and Y. Huang 702 Plusieurs méthodes de modélisation fonctionnelle ont été proposées, en présence d'erreur de mesure des covariables avec le modèle de taux proportionnels. Sont inclus dans ces méthodes, l'estimateur du score conditionnel (Tsiatis and Davidian, 2001, Biometrika 88, 447–458) l'estimateur avec correction paramétrique (Nakamura, 1992, Biometrics 48, 829–838) et l'estimateur avec correction non paramétrique (Huang and Wang, 2000, Journal of the American Statistical Association 95, 1209–1219) dans l'ordre d'hypothèses de plus en plus faible sur l'erreur. Bien qu'ils soient tous consistants, chacun de ces estimateurs souffre de possibles difficultés en présence de petits échantillons ou d'erreurs de mesure substantielles. Dans cet article, ayant observé que les estimateurs du score conditionnel ou avec corrections paramétriques étaient asymptotiquement équivalents pour une erreur normale, nous avons cherché quelle était leur performance relative pour des échantillons de taille finie et remarqué que le premier était le meilleur. Ceci signifie une approche générale affinée des méthodes de correction paramétriques et non paramétriques. Les corrections affinées des estimateurs sont asymptotiquement équivalentes à leurs correspondantes standards, mais ont des propriétés numériques améliorées et sont plus performantes lorsque les estimations standards n'existent pas ou sont extrêmes. Nous présentons des résultats de simulations et une application de ces méthodes à un essai clinique sur le HIV. A.‐C. Andrei and S. Murray 715 Ce travail étudie l'analyse séquentielle de données de survie appariées en utilisant une nouvelle classe de tests non‐paramétriques, fondés sur un log‐rank pondéré apparié (PWLR), et les tests à partir de Kaplan‐Meier pondérés appariés (PWKM). Pendant le déroulement d'un essai chacun de ces tests peut alternativement être considéré comme la statistique la plus extrême. En suivant l'évolution de PEMAX, le maximum des différences absolues du PWLR standardisé et du PWKM, on a, à la fois, les avantages des tests de rangs et les avantages des autres tests. Par des simulations on montre que l'utilisation de PEMAX contrôle le risque de type I et est quasiment aussi puissant que l'utilisation du meilleur des deux tests. Ceci, que l'hypothèse des hasards proportionnels soit respectée ou non. Désormais, PEMAX préserve la puissance de façon plus robuste que chacun des deux tests pris individuellement, tout en restant relativement simple à mettre en place et à utiliser dans une analyse. Un exemple, à partir de l'étude “Early Treatment Diabetic Retinopathy Study” (ETDRS), est donné. B. Lu 721 Dans les études d'observation avec un traitement dépendant du temps et des covariables dépendant du temps, il est souhaitable d'équilibrer la distribution des covariables à chaque instant. Un score de propension (“propensity score”) dépendant du temps basé sur le modèle de Cox à risques proportionnels est proposé et utilisé pour apparier les sujets exposés au risque. Il est montré que l'appariement sur ce score de propension permet d'obtenir une distribution équilibrée des covariables dans les groupes traité et non traité. Un appariement optimal avec divers plans est mis en place et comparé dans une étude de traitement chirurgical, cytoscopie et hydrodistention, pour une maladie chronique de la vessie, la cystite interstitielle. Des simulations suggèrent aussi que l'analyse statistique après appariement est plus performante que l'analyse sans appariement en terme d'estimation et d'intervalle de confiance. M. Naskar , K. Das , and J. G. Ibrahim 729 Une classe très générale de distributions de survie multivariées est utilisée pour analyser des données groupées de défaillance, sujettes à la censure et à des modes multiples de défaillance. Conditionnellement à des quantités spécifiques au groupe, la distribution jointe du temps de défaillance et de la variable indicatrice de l'évènement peuvent s'exprimer comme un mélange de la distribution des temps de défaillance dûs à un certain type (ou spécifique d'une cause) et de la distribution marginale du type de défaillance. Nous supposons ici que les distributions marginales des différents types de défaillance sont de fonctions logistiques de certaines covariables. Les quantités spécifiques au groupe sont soumises à une certaine distribution inconnue qui cause la fragilité. La distribution de fragilité est modélisée non paramétriquement par un processus de Dirichlet (DP). Dans une telle approche semi paramétrique une méthode hybride d'estimation est proposée basée sur l'algorithme pondéré du restaurant chinois (WCR) qui nous aide à générer des observations à partir de la distribution prédictive de la fragilité. L'algorithme de Monte‐Carlo ECM (MCEMC) joue un rôle vital pour obtenir des estimations des paramètres qui établissent l'importance des effets des facteurs causaux pour les défaillances d'un certain type. Une étude de simulation est conduite pour étudier la consistance de notre méthodologie. La méthodologie proposée est utilisée pour analyser des données réelles de contamination par le virus VIH sur une cohorte de femmes prostituées du Sénégal. F. Bretz , J. C. Pinheiro , and M. Branson 738 Traditionnellement, les stratégies d'analyse de données issues d'études dose‐réponse se partagent en deux grandes catégories : les procédures de comparaisons multiples et la modélisation. L'approche par modélisation suppose une relation fonctionnelle entre la réponse et la dose – celle‐ci est ici considérée comme un facteur quantitatif –, selon un modèle paramétrique spécifiéà l'avance. On utilise alors ce modèle pour estimer la dose permettant d'obtenir un niveau donné de réponse, mais la validité des conclusions dépend fortement de la pertinence du modèle de dose‐réponse, le vrai modèle étant a priori inconnu. A l'inverse, les procédures de comparaisons multiples considèrent la dose comme un facteur qualitatif et font très peu d'hypothèses – voire aucune – sur le modèle de dose‐réponse sous‐jacent. Avec ces procédures, il s'agira souvent d'identifier la dose minimale efficace, statistiquement significative et susceptible de produire un effet biologique considéré comme intéressant. Une des approches possibles est de tester la significativité de contrastes comparant différents niveaux de dose, tout en contrôlant le risque global de Type I. Ce type de procédure est relativement robuste ; cependant, l'inférence ne peut bien sûr pas sélectionner d'autres doses que celles utilisées dans l'étude. Dans cet article, nous décrivons une stratégie unifiée d'analyse de données issues d'études dose‐réponse, stratégie qui combine les techniques de comparaisons multiples et de modélisation. Nous supposons en fait l'existence de plusieurs modèles candidats, tous paramétriques, et utilisons des techniques de comparaisons multiples pour choisir le modèle le plus vraisemblable vis‐à‐vis de la vraie courbe dose‐réponse sous‐jacente, tout en contrôlant le risque global de Type I. On utilise alors le modèle ainsi sélectionné pour cerner, à l'aide d'une approche inférentielle, les doses les mieux appropriées. J. O'Quigley 749 La méthode de ré‐évaluation séquentielle (CRM) est un algorithme pour la recherche de dose qui s'appuie sur une remise à jour dynamique. Comme d'autres méthodes dynamiques, la CRM fournit une estimation de la dose courante qui correspond à un percentile cible. Cette dose, calculée à partir de tous les sujets inclus, sera utilisée pour la prochaine expérience. La mise en oeuvre de cette ré‐évaluation est possible grâce à l'utilisation d'un modèle simple. Dans l'état actuel des choses ni la CRM ni toute autre méthode dynamique permet d'estimer un percentile cible à partir de données rétrospectives, sauf dans le cas exceptionnel où le modèle simplifié génère les données. Dans cet article notre attention se porte sur un problème spécifique, à savoir l'analyse rétrospective par la CRM de données générées par un mécanisme arbitraire. Nous montrons comment procéder de façon cohérente. La méthodologie proposée n'est pas limitée à ce schéma particulier et s'applique à d'autres plans expérimentaux où l'estimation de la dose est faite par une mise à jour dynamique. L'analyse rétrospective pose un certain nombre de difficultés pour ces plans, en particulier la CRM, car ces plans ont la propriété d'exploiter un modèle de travail sous paramétré. On ne peut répondre à la question ‐ quel aurait été le résultat si on avait utilisé un autre plan expérimental. Néanmoins on peut fournir une estimation de ce résultat. Ceci a de nombreuses applications. Z.‐F. Yu and P. J. Catalano 757 Les effets neurotoxiques d'agents chimiques sont souvent analysés dans le cadre d'études contrôlées sur des rongeurs, pour de multiples variables (binaires et continues) collectées en routine. L'un des objectifs est d'évaluer le risque afin de déterminer les doses optimales non toxiques. De telles études sont cependant confrontées à deux problèmes majeurs concernant les variables continues. Premièrement, il est difficile, pour de telles variables, d'évaluer le risque et de définir une dose de référence. Lorsque le risque est associéà la survenue d'un événement, il peut clairement être défini de façon binaire par la survenue ou non de l'événement. Définir le risque de façon similaire pour des variables continues est moins évident. Dans ce cas, l'événement est souvent défini par un ensemble de valeurs inférieures à un certain seuil, la distribution de ces valeurs étant supposée normale ou log normale. Deuxièmement, alors que, pour de telles études, les variables continues sont traditionnellement analysées séparément, des travaux récents conseillent plutôt d'utiliser des modèles multivariés. Nous proposons une méthode pour modéliser et quantifier le risque associéà des variables continues bivariées. Il s'agit d'une extension des méthodes existantes de régression par quantiles. L'approche basée sur une vraisemblance, permet de dissocier les modèles dose‐réponse de chaque variable, tout en tenant compte de la corrélation bivariée et de la caractérisation globale du risque. L'approche utilisée pour estimer une dose de référence est analogue à celle utilisée pour des variables binaires, sans qu'il soit nécessaire de spécifier un seuil arbitraire. Nous illustrons nos méthodes avec des données issues d'une étude de neurotoxicité de l'exposition de rats au triethyltin. L. Huang , M.‐H. Chen , and J. G. Ibrahim 767 Nous proposons des méthodes bayésiennes pour estimer les paramètres de modèles linéaires généralisés (GLM) avec des données de covariables manquantes qu'on ne peut pas négliger. Nous montrons que lorsqu'une distribution uniforme a priori impropre est employée pour les coefficients de régression, ? , du modèle multinomial de sélection du mécanisme d'obtention des données manquantes, la distribution jointe a posteriori sera toujours impropre si (i) toutes les covariables manquantes sont discrètes et si un coefficient constant est inclus dans le modèle de sélection du mécanisme d'obtention de données manquantes, ou (ii) au moins une covariable est continue et non bornée. Ce caractère impropre se manifestera indifféremment que les a priori propres ou non soient spécifiés pour les paramètres de régression, β , du GLM ou des paramètres, α , de la distribution de la covariable. Pour résoudre ce problème, nous proposons une nouvelle classe d'a priori propres pour les coefficients de régression, ? , dans le modèle de sélection du mécanisme de données manquantes. Ces a priori sont robustes et intéressants pour le calcul au sens où les inférences concernant les β ne sont pas sensibles au choix des hyperparamètres des a priori de ? et qu'ils facilitent un schéma d'échantillonnage de Gibbs qui conduit à une convergence accélérée. De surplus, nous étendons le critère d'évaluation du modèle de Chen, Dey et Ibrahim (2004), appelémesure L pondérée, au GLM et aux problèmes des données manquantes de même que nous étendons le Critère d'Information de Déviance (DIC) de Spiegelhalter et al. (2002) pour évaluer si le mécanisme d'obtention des données manquantes peut être ignoré ou non. Un nouvel algorithme de Monte Carlo de chaîne de Markov est aussi développé pour mener à bien les calculs des a posteriori. Plusieurs simulations sont faites pour étudier la performance des critères bayésiens proposés et la sensibilité aux spécifications a priori. Des ensembles de données provenant d'un essai clinique sur une tumeur mélanique et d'une étude sur le cancer du foie sont présentés pour mieux illustrer las méthodes proposées. A. Golightly and D. J. Wilkinson 781 Cet article traite de l'estimation de constantes de taux dans des modèles stochastiques de processus intra‐cellulaire. Le modèle stochastique discret de cinétique sous jacent est remplacé par une approximation (par équation différentielle stochastique) où le bruit blanc modélise le comportement stochastique et on identifie le modèle utilise sur des intervalles de temps réguliers. L'estimation implique l'utilisation de m‐1 observations latentes entre chaque couples de dates d'observation. Les méthodes MCMC sont alors utilisées pour échantillonner la distribution a posteriori du processus latent et les paramètres du modèle. La méthodologie est appliquée à l'estimation des paramètres d'un processus d'auto régulation des gènes chez un procaryote. G. Diao and D. Y. Lin 789 Les méthodes statistiques pour détecter les gènes qui influencent des caractères quantitatifs à l'aide de marqueurs génétiques sont bien développées pour des phénotypes distribués normalement et complètement observés. De nombreuses expériences impliquent des phénotypes liés à un temps d'échec, qui ont des distributions asymétriques, et qui sont habituellement sujets à censure à cause des perdus de vue aléatoire, des échecs pour des causes en compétition ou de la limitation de la durée de l'expérience. Dans cet article nous développons des méthodes statistiques semi‐paramétriques pour cartographier des Locus de Caractères Quantitatifs (LCQ) basés sur des phénotypes à temps d'échec avec censure. Nous formulons les effets du génotype LCQ sur le temps d'échec au moyen du modèle à risques proportionnels de Cox (1972) et dérivons des procédures d'inférence efficaces basées sur la vraisemblance. De plus, nous montrons comment évaluer la signification statistique quand on recherche les LCQ dans plusieurs régions du génome complet. Des études extensives par simulation démontrent que les méthodes proposées ont de bonnes performances en pratique. Les applications à deux études animales sont présentées. W. Briggs and D. Ruppert 799 Doit‐on pratiquer des mammographies de précaution chez des femmes d'âge moyen en bonne santé? Les chirurgiens doivent‐ils utiliser le score TRISS bien connu pour prédire la chance de survie de leurs patients? Ce sont des exemples de questions auxquelles on est confronté lorsqu'on décide d'utiliser une prédiction de type Oui/Non. Pour retenir une prédiction nous devons montrer qu'elle est qu'elle est plus intéressante que serait notre meilleure vision en absence de prédiction. Le calcul d'une mesure permet d'identifier l'erreur lorsque la prédiction est erronée et d'examiner les performances passées de la prédiction compte tenu de cette erreur. On développe un test statistique dans ce but. Les prédictions qui seront acceptées par ce test seront dites avoir de la compétence. Seules les prédictions compétentes devraient être utilisées. On montrera des méthodes graphiques et numériques pour montrer la compétence. L'utilité des mammographies est étudiée par cette voie. O. Thas and J. C. W. Rayner 808 Dans cet article, nous construisons trois tests d'ajustement lissés pour tester une distribution ZIP contre des alternatives générales lisses au sens de Neyman. Nous appliquons nos tests à un ensemble de données auparavant présenté comme ZIP‐distribué, et nous montrons que le modèle ZIP ne décrit pas correctement les données. En cas de rejet de l'hypothèse nulle de modèle ZIP, les composantes individuelles de la statistique de test, qui sont directement liées à des paramètres interprétables dans un modèle lissé, peuvent être utilisées pour l'examen d'une distribution alternative.     

RÔLE DE L'EAU DE BOISSON SUR LA PRISE DE NOURRITURE SÈCHE ET LE DÉVELOPPEMENT OVARIEN DE LOCUSTA MIGRATORIA MIGRATORIOIDES

Pour déterminer la ration minimale d'eau nécessaire à la survie de Locusta migratoria migratorioides (R. et F.) adulte et les repercussions d'un stress hydrique prolongé sur le fonctionnement ovarien, nous avons séparé le facteur eau du facteur aliment en élevant les criquets de l'éclosion à l'âge adulte avec des germinations de blé lyophilisé, selon une technique décrite par Louveaux et al. (1980). Il est montré que la quantité d'aliments ingérée augmente avec la quantité d'eau bue (r = 0,91). La relation linéaire déterminée permet inversement, connaissant les quantités d'aliments ingérés, d'évaluer à 300 μl la ration quotidienne d'un criquet ayant un accès non limité à l'eau de boisson. Les femelles qui n'ont reçu que 50 μl d'eau par jour (j.) sont mortes avant le 9ème jour de traitement. La survie est assurée à partir de 100 μl/j. L'avancement de la vitellogenèse est également en relation avec la quantité d'eau absorbée quotidiennement. Aucun développement ovarien ne se produit jusqu'à 100 μl/j. Une vitellogenèse ralentie commence à partir de 150 μl/j. et elle devient normale à 300 μl/j. avec des pontes comparables à celles des témoins nourris de blé frais.     

Le pré-analytique en hémostase

La phase pré-analytique est celle de tous les dangers, en particulier lorsque le biologiste n’a pas la ma?trise du recueil de l’échantillon et de son acheminement au laboratoire. En hémostase thrombose, elle est essentielle. Les paramètres pré-analytiques sont d’importance inégale selon les différents examens et leur retentissement peut varier d’un réactif à l’autre et selon les patients. De fait, il existe un risque de rendre un résultat erroné, même avec une technique irréprochable, ainsi qu’un problème de responsabilité technico-légale. La littérature sur les conditions pré-analytiques est pauvre et de nouvelles études sont nécessaires, même si elles sont difficiles à mener et peuvent para?tre de peu d’intérêt. En réalité, leur utilité est aujourd’hui indéniable, en raison notamment du regroupement des laboratoires.     

Microcharops anticarsiae (Hym.: Ichneumonidae) parasitoid ofAnticarsia gemmatalis (Lep.: Noctuidae): Host age preference,sex ratio effects and functional response

Résumé L'endoparasito?deMicrocharops anticarsiae (Hym.: Ichneumonidae), est un des agents naturels de régulation d'Anticarsia gemmatalis (Lep.: Noctuidae). La préférence pour un stade de l'h?te, les effets sur la sex ratio et la réponse fonctionnelle ont été étudiés expérimentalement chezM. anticarsiae pour six densités différentes de l'h?te. Dans des expériences de choix impliquant 2 stades de développement consécutifs ou dans les cas où étaient offerts en même temps les quatre premiers stades larvaires, les femelles du parasito?de n'ont présenté aucune préférence pour un stade donné. Ces résultats montrent que les femelles deM. anticarsiae acceptent tout h?te disponible indépendamment de son stade. La sex ratio présente des différences significatives entre les quatre stades de l'h?te (G_3df=18.21, p<.001). La sex ratio du parasito?de est plus favorable aux femelles chez les larves se développant dans le 4^e stade de l'h?te. Une réponse fonctionnelle de type II a été observée chez les parasito?des femelles en réponse à des densités variables de l'h?te. Ces informations fondamentales concernant les effets de la sex ratio ont montré quel est le stade de l'h?te qui doit être parasité pour obtenir la sex ratio la plus favorable aux femelles, tandis que la réponse fonctionnelle a fourni le taux moyen (38.0±7.5) de parasitisme par femelle. Ces informations sont importantes pour la production en masse deM. anticarsiae pour le programme de lutte biologique.      

DEVELOPPEMENT EMBRYONNAIRE IN VITRO,EN MILIEU ARTIFICIEL DEFINI DE DEUX PARASITOIDES OVOLARVIPARES,PHRYXE CAUDATA ET LIXOPHAGA DIATRAEAE (DIPTERA,TACHINIDAE)

Les utérus sont prélevés par dissection des femelles âgées soit de 2 jours (P. caudata), soit de 2 à 4 jours (L. diatraeae); ils contiennent alors des oeufs qui viennent juste d'être fécondés. Ces utérus sont placés aseptiquement sur des milieux de culture pour observer l'embryogénèse éventuelle des oeufs. Un milieu physiologique simple (Y) et un milieu nutritif complet (M) sont testés sous deux formes, liquide ou gélifiée avec 0,5% d'agarose. Chez les deux espèces le développement embryonnaire complet a pu être observé in vitro. Les meilleurs résultats, qui fournissent le même pourcentage de larves qu'in vivo, sont constatés avec le milieu M agarosé. Sur les milieux liquides le rendement en larves est inférieur de moitié environ par rapport aux milieux gélifiés correspondants. Avec les milieux simples Y, mème si un fort pourcentage de larves est obtenu (Y agarosé), celes-ci meurent rapidement. Seul le milieu M permet d'obtenir des larves d'excellente viabilité et aptes à se développer. La composition du milieu, mais aussi les possibilités respiratoires offertes aux oeufs agissent sur l'embryogénése.