首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
Henrys PA  Brown PE 《Biometrics》2009,65(2):423-430
Summary .  We propose a method to test for significant differences in the levels of clustering between two spatial point processes (cases and controls) while taking into account differences in their first-order intensities. The key advance on earlier methods is that the controls are not assumed to be a Poisson process. Inference and diagnostics are based around the inhomogeneous K -function with confidence envelopes obtained from either resampling events in a nonparametric bootstrap approach, or simulating new events as in a parametric bootstrap. Methods developed are demonstrated using the locations of adult and juvenile trees in a tropical forest. A simulation study briefly examines the accuracy and power of the inferential procedures.  相似文献   

2.
The likelihood ratio test for testing equality of vgE;2 correlated variables is developed. In general, evaluation of the test statistic involves the iterative optimization of a likelihood function with 1 + v(v – 1)/2 parameters. The explicit form of the test statistic is derived in the bivariate case, and an iterative algorithm for determining the maximum likelihood estimates is suggested. A limited Monte Carlo study determines the behavior of the proposed procedure under the null hypothesis and variety of parameter values.  相似文献   

3.
A nationwide health card recording system for dairy cattle was introduced in Norway in 1975 (the Norwegian Cattle Health Services). The data base holds information on mastitis occurrences on an individual cow basis. A reduction in mastitis frequency across the population is desired, and for this purpose risk factors are investigated. In this paper a Bayesian proportional hazards model is used for modelling the time to first veterinary treatment of clinical mastitis, including both genetic and environmental covariates. Sire effects were modelled as shared random components, and veterinary district was included as an environmental effect with prior spatial smoothing. A non-informative smoothing prior was assumed for the baseline hazard, and Markov chain Monte Carlo methods (MCMC) were used for inference. We propose a new measure of quality for sires, in terms of their posterior probability of being among the, say 10% best sires. The probability is an easily interpretable measure that can be directly used to rank sires. Estimating these complex probabilities is straightforward in an MCMC setting. The results indicate considerable differences between sires with regards to their daughters disease resistance. A regional effect was also discovered with the lowest risk of disease in the south-eastern parts of Norway.  相似文献   

4.
5.
    
Swartz MD  Kimmel M  Mueller P  Amos CI 《Biometrics》2006,62(2):495-503
Mapping the genes for a complex disease, such as diabetes or rheumatoid arthritis (RA), involves finding multiple genetic loci that may contribute to the onset of the disease. Pairwise testing of the loci leads to the problem of multiple testing. Looking at haplotypes, or linear sets of loci, avoids multiple tests but results in a contingency table with sparse counts, especially when using marker loci with multiple alleles. We propose a hierarchical Bayesian model for case-parent triad data that uses a conditional logistic regression likelihood to model the probability of transmission to a diseased child. We define hierarchical prior distributions on the allele main effects to model the genetic dependencies present in the human leukocyte antigen (HLA) region of chromosome 6. First, we add a hierarchical level for model selection that accounts for both locus and allele selection. This allows us to cast the problem of identifying genetic loci relevant to the disease into a problem of Bayesian variable selection. Second, we attempt to include linkage disequilibrium as a covariance structure in the prior for model coefficients. We evaluate the performance of the procedure with some simulated examples and then apply our procedure to identifying genetic markers in the HLA region that influence risk for RA. Our software is available on the website http://www.epigenetic.org/Linkage/ssgs-public/.  相似文献   

6.
This paper is concerned with comparing several increasing dose levels (treatments) with a zero dose control when the prior information about the umbrella pattern treatment means is available. The problem of testing whether there is at least one treatment which is better than the control is considered. Multiple test procedures are then proposed for deciding treatments (if any) which are better than the control. Some approximate criticial values of the proposed tests are reported. The results of a Monte Carlo power study are presented.  相似文献   

7.
    
Chen YI 《Biometrics》1999,55(4):1236-1240
We consider identifying the minimum effective dose (MED) in a dose-response study, where the MED is defined to be the lowest dose level producing an effect over that of the zero-dose control. Proposed herein is a nonparametric procedure based on the Mann-Whitney statistic incorporated with the step-down closed testing scheme. A numerical example demonstrates the feasibility of the proposed nonparametric procedure. Finally, the comparative results of a Monte Carlo level and power study for small sample sizes are presented and discussed.  相似文献   

8.
Let categorical data coming from a control group and (r - 1) treated groups be given in an r × c contingency table. A simultaneous test procedure of the (r - 1) hypotheses that the probabilities of all c categories do not differ between the i-th treated group and the control is derived. For small tables and small cell frequencies it is exactly performed by generation of all tables having the given marginal sums. If 2 categories or 2 groups only are given the asymptotic distribution of the test statistic is known; otherwise its distribution may be simulated if the computational expenditure of performing an exact test is too large. By means of a Monte Carlo study it is shown that this method meets its level more reliably and that it has a better power than others.  相似文献   

9.
    
The relative influence of niche vs. neutral processes in ecosystem dynamics is an on‐going debate, but the extent to which they structured the earliest animal communities is unknown. Some of the oldest known metazoan‐dominated paleocommunities occur in Ediacaran age (~ 565 million years old) strata in Newfoundland, Canada and Charnwood Forest, UK. These comprise large and diverse populations of sessile organisms that are amenable to spatial point process analyses, enabling inference of the most likely underlying niche or neutral processes governing community structure. We mapped seven Ediacaran paleocommunities using LiDAR, photogrammetry and a laser line probe. We found that neutral processes dominate these paleocommunities, with niche processes exerting limited influence, in contrast with the niche‐dominated dynamics of modern marine ecosystems. The dominance of neutral processes suggests that early metazoan diversification may not have been driven by systematic adaptations to the local environment, but instead may have resulted from stochastic demographic differences.  相似文献   

10.
Czika W  Weir BS 《Biometrics》2004,60(1):69-74
Disease genes can be mapped on the basis of associations between genetic markers and disease status, with the case-control design having the advantage of not requiring individuals from different generations. When the marker loci have multiple alleles, there has been debate on whether the power of tests for association increases or decreases. We show here that the multiple-allele version of Armitage's trend test has increased power over the two-allele version under the requirement of equifrequent alleles, but not in general. The trend test has the advantage of remaining valid even when the sampled population is not in Hardy-Weinberg equilibrium. A departure from Hardy-Weinberg means that association tests depend on gametic and nongametic linkage disequilibrium between marker and disease loci, and we illustrate the magnitude of these effects with simulated data.  相似文献   

11.
12.
    
Song PX  Gao X  Liu R  Le W 《Biometrics》2006,62(2):545-554
Identifying local extrema of expression profiles is one primary objective in some cDNA microarray experiments. To study the replication dynamics of the yeast genome, for example, local peaks of hybridization intensity profiles correspond to putative replication origins. We propose a nonparametric kernel smoothing (NKS) technique to detect local hybridization intensity extrema across chromosomes. The novelty of our approach is that we base our inference procedures on equilibrium points, namely those locations at which the first derivative of the intensity curve is zero. The proposed smoothing technique provides both point and interval estimation for the location of local extrema. Also, this technique can be used to test for the hypothesis of either one or multiple suspected locations being the true equilibrium points. We illustrate the proposed method on a microarray data set from an experiment designed to study the replication origins in the yeast genome, in that the locations of autonomous replication sequence (ARS) elements are identified through the equilibrium points of the smoothed intensity profile curve. Our method found a few ARS elements that were not detected by the current smoothing methods such as the Fourier convolution smoothing.  相似文献   

13.
    
Strong statistical evidence was found for differences in tolerance to natural infections of Tobacco streak virus (TSV) in sunflower hybrids. Data from 470 plots involving 23 different sunflower hybrids tested in multiple trials over 5 years in Australia were analysed. Using a Bayesian Hierarchical Logistic Regression (BHLR) model for analysis provided: (a) a rigorous method for investigating the relative effects of hybrid, seasonal rainfall and proximity to inoculum source on the incidence of severe TSV disease; (b) a natural method for estimating the probability distributions of disease incidence in different hybrids under historical rainfall conditions; and (c) a method for undertaking all pairwise comparisons of disease incidence between hybrids while controlling the familywise error rate without any drastic reduction in statistical power. The tolerance identified in field trials was effective against the main TSV strain associated with disease outbreaks, TSV‐parthenium. Glasshouse tests indicate this tolerance to also be effective against the other TSV strain found in central Queensland, TSV‐crownbeard. The use of tolerant germplasm is critical to minimise the risk of TSV epidemics in sunflower in this region. We found strong statistical evidence that rainfall during the early growing months of March and April had a negative effect on the incidence of severe infection with greatly reduced disease incidence in years that had high rainfall during this period.  相似文献   

14.
    
Liang LJ  Weiss RE 《Biometrics》2007,63(3):733-741
Phylogenetic modeling is computationally challenging and most phylogeny models fit a single phylogeny to a single set of molecular sequences. Individual phylogenetic analyses are typically performed independently using publicly available software that fits a computationally intensive Bayesian model using Markov chain Monte Carlo (MCMC) simulation. We develop a Bayesian hierarchical semiparametric regression model to combine multiple phylogenetic analyses of HIV-1 nucleotide sequences and estimate parameters of interest within and across analyses. We use a mixture of Dirichlet processes as a prior for the parameters to relax inappropriate parametric assumptions and to ensure the prior distribution for the parameters is continuous. We use several reweighting algorithms for combining completed MCMC analyses to shrink parameter estimates while adjusting for data set-specific covariates. This avoids constructing a large complex model involving all the original data, which would be computationally challenging and would require rewriting the existing stand-alone software.  相似文献   

15.
    
Summary .  We introduce a method of estimating disease prevalence from case–control family study data. Case–control family studies are performed to investigate the familial aggregation of disease; families are sampled via either a case or a control proband, and the resulting data contain information on disease status and covariates for the probands and their relatives. Here, we introduce estimators for overall prevalence and for covariate-stratum-specific (e.g., sex-specific) prevalence. These estimators combine the proportion of affected relatives of control probands with the proportion of affected relatives of case probands and are designed to yield approximately unbiased estimates of their population counterparts under certain commonly made assumptions. We also introduce corresponding confidence intervals designed to have good coverage properties even for small prevalences. Next, we describe simulation experiments where our estimators and intervals were applied to case–control family data sampled from fictional populations with various levels of familial aggregation. At all aggregation levels, the resulting estimates varied closely and symmetrically around their population counterparts, and the resulting intervals had good coverage properties, even for small sample sizes. Finally, we discuss the assumptions required for our estimators to be approximately unbiased, highlighting situations where an alternative estimator based only on relatives of control probands may perform better.  相似文献   

16.
17.
18.
19.
This paper is concerned with testing for umbrella alternatives in a k-sample location problem when the underlying populations have possibly different shapes. Following CHEN and WOLFE (1990b), rank-based modifications of the HETTMANSPERGER-NORTON (1987) tests are considered for both the settings where the peak of the umbrella is known and where it is unknown. The proposed procedures are exactly distribution-free when the continuous populations are identical with any shape. Moreover, the modified test for peak-known umbrella alternatives remains asymptotically distribution-free when the continuous populations are assumed to be symmetric, even if they differ in shapes. Comparative results of a Monte Carlo study are presented.  相似文献   

20.
    
Wang L  Zhou XH 《Biometrics》2007,63(4):1218-1225
Heteroscedastic data arise in many applications. In heteroscedastic regression analysis, the variance is often modeled as a parametric function of the covariates or the regression mean. We propose a kernel-smoothing type nonparametric test for checking the adequacy of a given parametric variance structure. The test does not need to specify a parametric distribution for the random errors. It is shown that the test statistic has an asymptotical normal distribution under the null hypothesis and is powerful against a large class of alternatives. We suggest a simple bootstrap algorithm to approximate the distribution of the test statistic in finite sample size. Numerical simulations demonstrate the satisfactory performance of the proposed test. We also illustrate the application by the analysis of a radioimmunoassay data set.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号