Graphical aids to interpretation of Scatchard plots and dose-response curves

General expressions are derived for the limiting slopes and intercepts of graphical representations of experimental binding data in either the Scatchard or the “dose-response” co-ordinate systems. We apply a previously formulated general model that includes heterogeneity and/or cooperativity of receptor affinity. One must establish or assume a physical chemical mechanism in order to fully interpret these limiting slopes and intercepts. However, they do provide satisfactory initial estimates for the binding parameters, for use in a non-linear least squares curve fitting approach using an exact model.     

ESTIMATION OF THE KINETIC PARAMETERS OF UNIDIRECTIONAL TRANSPORT ACROSS THE BLOOD-BRAIN BARRIER

Abstract— The unidirectional transport of metabolic substrates from blood to brain may be defined in terms of Michaelis-Menten saturable ( K m, V _max) and non-saturable ( K _d) components of influx. Various computation procedures have been previously reported to estimate the kinetic parameters when an intracarotid injection technique is used. Transformations of the influx data which allow linear plots to obtain estimates were compared with estimates obtained directly from a best fit on a least means squares criterion for both experimental and simulated data. Large discrepancies were apparent between the various estimates of the kinetic parameters when an equal weight was given to transformed data. For pyruvate (21-day-old rats), K _m, values varied between 1.02 and 6.25 mM and V _max varied between 0.68 and 2.30 μmol g⁻¹ min⁻¹. The estimates were almost equivalent when pyruvate data was re-analysed using a weighting scheme based on the finding that the absolute value of the S.D. of influx increased in proportion to influx. It is recommended that estimates of kinetic parameters be obtained by an iterative, non-linear least squares method to fit appropriately weighted data directly.     

A likelihood approach for functional response models

Functional response is an important determinant of community dynamics, and thus empirical methods for characterizing functional responses are as important in understanding ecological processes. The most commonly used method is based on the sum of squares, and the maximum likelihood method is rarely used. When the likelihood method is used, potentially inappropriate probability distributions such as binomial distributions are typically assumed for the number of prey eaten in experiments. In this study, I present a likelihood approach in which the probability distributions are generated by mechanistic understanding of predation processes using Monte Carlo simulations. An example is given on the Holling type II functional response model, but the method is flexible and allows characterization of a wide variety of functional response models. In the example, the likelihood method consistently resulted in superior estimates than the least squares method.     

Maximum likelihood vs. minimum chi-square--a general comparison with applications to the estimation of recombination fractions in two-point linkage analysis.

Some relationships between the estimates of recombination fraction in two-point linkage analysis obtained by maximum likelihood, minimum chi-square, and general least squares are derived. These theoretical results are based on an approximation for the multinomial distribution. Applications (theoretical and experimental) with RFLP (restriction fragment length polymorphism) markers for a segregating F2 population are given. The minimum chi-square estimate is slightly larger than the maximum likelihood estimate. For applications, however, both estimates must be considered to be approximately equal. The least squares estimates are slightly different (larger or smaller) from these estimates.     

A general framework for constraint approaches to adjusted risk differences

The risk difference is an intelligible measure for comparing disease incidence in two exposure or treatment groups. Despite its convenience in interpretation, it is less prevalent in epidemiological and clinical areas where regression models are required in order to adjust for confounding. One major barrier to its popularity is that standard linear binomial or Poisson regression models can provide estimated probabilities out of the range of (0,1), resulting in possible convergence issues. For estimating adjusted risk differences, we propose a general framework covering various constraint approaches based on binomial and Poisson regression models. The proposed methods span the areas of ordinary least squares, maximum likelihood estimation, and Bayesian inference. Compared to existing approaches, our methods prevent estimates and confidence intervals of predicted probabilities from falling out of the valid range. Through extensive simulation studies, we demonstrate that the proposed methods solve the issue of having estimates or confidence limits of predicted probabilities out of (0,1), while offering performance comparable to its alternative in terms of the bias, variability, and coverage rates in point and interval estimation of the risk difference. An application study is performed using data from the Prospective Registry Evaluating Myocardial Infarction: Event and Recovery (PREMIER) study.     

Correcting for bias in estimation of quantitative trait loci effects

Estimates of quantitative trait loci (QTL) effects derived from complete genome scans are biased, if no assumptions are made about the distribution of QTL effects. Bias should be reduced if estimates are derived by maximum likelihood, with the QTL effects sampled from a known distribution. The parameters of the distributions of QTL effects for nine economic traits in dairy cattle were estimated from a daughter design analysis of the Israeli Holstein population including 490 marker-by-sire contrasts. A separate gamma distribution was derived for each trait. Estimates for both the α and β parameters and their SE decreased as a function of heritability. The maximum likelihood estimates derived for the individual QTL effects using the gamma distributions for each trait were regressed relative to the least squares estimates, but the regression factor decreased as a function of the least squares estimate. On simulated data, the mean of least squares estimates for effects with nominal 1% significance was more than twice the simulated values, while the mean of the maximum likelihood estimates was slightly lower than the mean of the simulated values. The coefficient of determination for the maximum likelihood estimates was five-fold the corresponding value for the least squares estimates.     

Unbiased Estimation of Odds Ratio by Using Negative Binomial Plans

The use of the negative binomial distribution in both the numerator and denominator in prospective studies leads to an unbiased estimate of the odds ratio and an exact expression for its variance. Sample sizes that minimize the variance of odds ratio estimates are specified. The variance of the odds ratio estimate is shown to be close to the Cramér-Rao lower bound.     

Estimation of a common effect parameter from sparse follow-up data

Breslow (1981, Biometrika 68, 73-84) has shown that the Mantel-Haenszel odds ratio is a consistent estimator of a common odds ratio in sparse stratifications. For cohort studies, however, estimation of a common risk ratio or risk difference can be of greater interest. Under a binomial sparse-data model, the Mantel-Haenszel risk ratio and risk difference estimators are consistent in sparse stratifications, while the maximum likelihood and weighted least squares estimators are biased. Under Poisson sparse-data models, the Mantel-Haenszel and maximum likelihood rate ratio estimators have equal asymptotic variances under the null hypothesis and are consistent, while the weighted least squares estimators are again biased; similarly, of the common rate difference estimators the weighted least squares estimators are biased, while the estimator employing "Mantel-Haenszel" weights is consistent in sparse data. Variance estimators that are consistent in both sparse data and large strata can be derived for all the Mantel-Haenszel estimators.     

Log-normal variation belts for growth curves

Prediction (confidence) or tolerance belts compound the uncertainty of sample estimates with the estimated extent of individual variation. The latter is therefore better described by variation belts, in which sample estimates are simply substituted for population parameters. Variation belts can provide valuable graphical indications concerning the goodness of fit of postulated error models. While multiplicative least-squares (MLS) methods appear appropriate in principle for biological growth, they are unsatisfactory in practice when logarithmically transformed data are heteroscedastic. Heteroscedastic multiplicative error models can be fitted by iteratively reweighted multiplicative least squares (IRMLS), but unacceptable negative or infinite residual variance estimates and unreasonably wide variation belts are occasionally obtained. These difficulties can be prevented by constrained iteratively reweighted multiplicative least squares (CIRMLS). Examples are presented concerning the metabolic allometry of white rats, the somatic growth of male elephant seals, and the growth of an experimental population of Paramecium caudatum.     

The Applicability of Ordinary Least Squares to Consistently Short Distances between Taxa in Phylogenetic Tree Construction and the Normal Distribution Test Consequences

Short phylogenetic distances between taxa occur, for example, in studies on ribosomal RNA-genes with slow substitution rates. For consistently short distances, it is proved that in the completely singular limit of the covariance matrix ordinary least squares (OLS) estimates are minimum variance or best linear unbiased (BLU) estimates of phylogenetic tree branch lengths. Although OLS estimates are in this situation equal to generalized least squares (GLS) estimates, the GLS chi-square likelihood ratio test will be inapplicable as it is associated with zero degrees of freedom. Consequently, an OLS normal distribution test or an analogous bootstrap approach will provide optimal branch length tests of significance for consistently short phylogenetic distances. As the asymptotic covariances between branch lengths will be equal to zero, it follows that the product rule can be used in tree evaluation to calculate an approximate simultaneous confidence probability that all interior branches are positive.     

Generalized least squares for the synthesis of correlated information

This paper deals with the synthesis of information from different studies when there is lack of independence in some of the contrasts to be combined. This problem can arise in several different situations in both case-control studies and clinical trials. For efficient estimation we appeal to the method of generalized least squares to estimate the summary effect and its standard error. This method requires estimates of the covariances between those contrasts that are not independent. Although it is not possible to estimate the covariance between effects that have been adjusted for confounding factors we present a method for finding upper and lower bounds for this covariance. In the simplest discussion homogeneity of the relative risks is assumed but the method is then extended to allow for heterogeneity in an overall estimate. We then illustrate the method with several examples from an analysis in which case-control studies of cervical cancer and oral contraceptive use are synthesized.     

One‐inflation and unobserved heterogeneity in population size estimation

We present the one‐inflated zero‐truncated negative binomial (OIZTNB) model, and propose its use as the truncated count distribution in Horvitz–Thompson estimation of an unknown population size. In the presence of unobserved heterogeneity, the zero‐truncated negative binomial (ZTNB) model is a natural choice over the positive Poisson (PP) model; however, when one‐inflation is present the ZTNB model either suffers from a boundary problem, or provides extremely biased population size estimates. Monte Carlo evidence suggests that in the presence of one‐inflation, the Horvitz–Thompson estimator under the ZTNB model can converge in probability to infinity. The OIZTNB model gives markedly different population size estimates compared to some existing truncated count distributions, when applied to several capture–recapture data that exhibit both one‐inflation and unobserved heterogeneity.     

Improved exact confidence intervals for the odds ratio in two independent binomial samples

For two independent binomial samples, the usual exact confidence interval for the odds ratio based on the conditional approach can be very conservative. Recently, Agresti and Min (2002) showed that the unconditional intervals are preferable to conditional intervals with small sample sizes. We use the unconditional approach to obtain a modified interval, which has shorter length, and its coverage probability is closer to and at least the nominal confidence coefficient.     

Calculation of narrower confidence intervals for tree mortality rates when we know nothing but the location of the death/survival events

相似文献   

Point and Interval Estimation of the Population Size Using a Zero‐Truncated Negative Binomial Regression Model

This paper presents the zero‐truncated negative binomial regression model to estimate the population size in the presence of a single registration file. The model is an alternative to the zero‐truncated Poisson regression model and it may be useful if the data are overdispersed due to unobserved heterogeneity. Horvitz–Thompson point and interval estimates for the population size are derived, and the performance of these estimators is evaluated in a simulation study. To illustrate the model, the size of the population of opiate users in the city of Rotterdam is estimated. In comparison to the Poisson model, the zero‐truncated negative binomial regression model fits these data better and yields a substantially higher population size estimate. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The number distribution for involved lymph nodes in cancer

The number of involved lymph nodes exhibits considerable heterogeneity within populations. Here, the implications of population heterogeneity are explored with respect to the kinematics of nodal metastases. Data from the National Cancer Institute's Surveillance, Epidemiology, and End Results program for 224656 breast, 12404 gastric, 18015 rectal, 4117 cervical and 2443 laryngeal cancers as well as 9118 melanomas were used to construct frequency distributions for the number of involved nodes which were then fitted to the negative binomial distribution. The negative binomial distribution described the heterogeneity in nodal involvement well. The patterns of nodal involvement can be explained by either of two models: one where involved nodes could seed further nodal metastases, the other where the number of nodal metastases in any individual was randomly distributed, with the deviations between patients accounted for by population heterogeneity. Since the number of sampled nodes similarly approximated a negative binomial distribution, random involvement with superimposed population heterogeneity would more credibly explain both sets of observations.     

Evaluation of the jackknife technique for fitting multiexponential functions to biochemical data

The jackknife technique was tested by fitting a two-exponential function to the time course of disappearance of radioactivity from the area of a wheat leaf that had been fed ¹⁴CO₂. The function was fitted by both unweighted and weighted least squares, first without and then with the jackknife. Weighting altered the estimates of the function's parameters, but jackknifing did not. Hence jackknifing did not remove any of the bias introduced by incorrect weighting. The confidence limits of the parameters calculated by jackknifing were greater than those estimated from the variance-covariance matrix of the regression, but similar to those derived from replicate experiments. The jackknife also allowed confidence limits for the rate constants and transit time of the underlying two-compartment model to be derived.     

Using the EM algorithm to weight data sets of unknown precision when modelling fish stocks

Stocks of commercial fish are often modelled using sampling data of various types, of unknown precision, and from various sources assumed independent. We want each set to contribute to estimates of the parameters in relation to its precision and goodness of fit with the model. Iterative re-weighting of the sets is proposed for linear models until the weight of each set is found to be proportional to (relative weighting) or equal to (absolute weighting) the set-specific residual invariances resulting from a generalised least squares fit. Formulae for the residual variances are put forward involving fractional allocation of degrees of freedom depending on the numbers of independent observations in each set, the numbers of sets contributing to the estimate of each parameter, and the number of weights estimated. To illustrate the procedure, numbers of the 1984 year-class of North Sea cod (a) landed commercially each year, and (b) caught per unit of trawling time by an annual groundfish survey are modelled as a function of age to estimate total mortality, Z, relative catching power of the two fishing methods, and relative precision of the two sets of observations as indices of stock abundance. It was found that the survey abundance indices displayed residual variance about 29 times higher than that of the annual landings.     

An application of multivariate ratio methods for the analysis of a longitudinal clinical trial with missing data.

This paper presents an analysis of a longitudinal multi-center clinical trial with missing data. It illustrates the application, the appropriateness, and the limitations of a straightforward ratio estimation procedure for dealing with multivariate situations in which missing data occur at random and with small probability. The parameter estimates are computed via matrix operators such as those used for the generalized least squares analysis of catetorical data. Thus, the estimates may be conveniently analyzed by asymptotic regression methods within the same computer program which computes the estimates, provided that the sample size is sufficiently computer program which computes the estimates, provided that the sample size is sufficiently large.     

Non-linear heteroscedastic regression model for determination of methotrexate in human plasma by high-performance liquid chromatography

Generalized least squares regression with variance function estimation was used to derive the calibration function for measurement of methotrexate plasma concentration and its results were compared with weighted least squares regression by usual weight factors and also with that of ordinary least squares method. In the calibration curve range of 0.05 to 100 microM, both heteroscedasticity and non-linearity were present therefore ordinary least squares linear regression methods could result in large errors in the calculation of methotrexate concentration. Generalized least squares regression with variance function estimation worked better than both the weighted regression with the usual weight factors and ordinary least squares regression and gave better estimates for methotrexate concentration.