The Scalpel and the Butterfly: The Conflict between Animal Research and Animal Protection Deborah Rudacille

ABSTRACT

Although there is a substantial body of research on inter-human empathy and inter-animal empathy, there is a dearth of research comparing humans' empathic reactions to humans and animals. To address this issue, three experiments were conducted in which participants read a scenario about a human or animal abuse victim in need of medical attention, and indicated the degree of empathy they had on an emotional response scale. In Experiment 1, women had significantly more empathy for animals than for humans, whereas men tended to express more empathy for humans than for animals. In Experiment 2, adult women expressed the same degree of empathy for a child as for a puppy. Similarly, in Experiment 3, adult men and women expressed the same degree of empathy for a baby as for a puppy. Overall, results indicated that people feel at least as much empathy for animals as for humans. We suggest that an animal target elicits a great deal of empathy partly because it is perceived as not being responsible for having caused the need situation. Future research will show whether empathy for animals translates to prosocial behavior toward them as well.     

Cloning Vectors for Rice

We developed various binary vectors that can be used for expressing a foreign gene in rice. Vectors pGA3426, pGA3436, and pGA3626 are intended for overexpression of a gene using the maize Ubiquitin promoter, whereas pGA3780 is for rather mild expression of a gene using the rice Actin1 promoter. Vector pGA3777 is for expressing two genes simultaneously. We also developed binary vectors for expressing a fusion protein with a tag. Four vectors (pGA3427, pGA3428, pGA3429, and pGA3438) are for protein tags with sGFP, HA, His, and Myc, respectively. Vector pGA3383 is for analyzing promoter activity using the GUS reporter. In this vector, multiple cloning sites in front of GUS can be utilized for accepting a promoter fragment. We also generated transient expression vectors for studying the subcellular localization of a protein. Vectors pGA3452, pGA3651, and pGA3652 are for GFP fusion; pGA3574 for RFP fusion; pGA3697 for Myc tag; and pGA3698 for HA tag. In addition, we generated pGA3506, pGA3516, pGA3592, and pGA3593 for facilitating the subcloning of full-length cDNA clones into our binary vectors.     

Models and estimators linking individual-based and sample-based rarefaction,extrapolation and comparison of assemblages

Aims In ecology and conservation biology, the number of species counted in a biodiversity study is a key metric but is usually a biased underestimate of total species richness because many rare species are not detected. Moreover, comparing species richness among sites or samples is a statistical challenge because the observed number of species is sensitive to the number of individuals counted or the area sampled. For individual-based data, we treat a single, empirical sample of species abundances from an investigator-defined species assemblage or community as a reference point for two estimation objectives under two sampling models: estimating the expected number of species (and its unconditional variance) in a random sample of (i) a smaller number of individuals (multinomial model) or a smaller area sampled (Poisson model) and (ii) a larger number of individuals or a larger area sampled. For sample-based incidence (presence–absence) data, under a Bernoulli product model, we treat a single set of species incidence frequencies as the reference point to estimate richness for smaller and larger numbers of sampling units.Methods The first objective is a problem in interpolation that we address with classical rarefaction (multinomial model) and Coleman rarefaction (Poisson model) for individual-based data and with sample-based rarefaction (Bernoulli product model) for incidence frequencies. The second is a problem in extrapolation that we address with sampling-theoretic predictors for the number of species in a larger sample (multinomial model), a larger area (Poisson model) or a larger number of sampling units (Bernoulli product model), based on an estimate of asymptotic species richness. Although published methods exist for many of these objectives, we bring them together here with some new estimators under a unified statistical and notational framework. This novel integration of mathematically distinct approaches allowed us to link interpolated (rarefaction) curves and extrapolated curves to plot a unified species accumulation curve for empirical examples. We provide new, unconditional variance estimators for classical, individual-based rarefaction and for Coleman rarefaction, long missing from the toolkit of biodiversity measurement. We illustrate these methods with datasets for tropical beetles, tropical trees and tropical ants.Important findings Surprisingly, for all datasets we examined, the interpolation (rarefaction) curve and the extrapolation curve meet smoothly at the reference sample, yielding a single curve. Moreover, curves representing 95% confidence intervals for interpolated and extrapolated richness estimates also meet smoothly, allowing rigorous statistical comparison of samples not only for rarefaction but also for extrapolated richness values. The confidence intervals widen as the extrapolation moves further beyond the reference sample, but the method gives reasonable results for extrapolations up to about double or triple the original abundance or area of the reference sample. We found that the multinomial and Poisson models produced indistinguishable results, in units of estimated species, for all estimators and datasets. For sample-based abundance data, which allows the comparison of all three models, the Bernoulli product model generally yields lower richness estimates for rarefied data than either the multinomial or the Poisson models because of the ubiquity of non-random spatial distributions in nature.     

Mapping quantitative trait Loci underlying fitness-related traits in a free-living sheep population

We searched for quantitative trait loci (QTL) underlying fitness-related traits in a free-living pedigree of 588 Soay sheep in which a genetic map using 251 markers with an average spacing of 15 cM had been established previously. Traits examined included birth date and weight, considered both as maternal and offspring traits, foreleg length, hindleg length, and body weight measured on animals in August and jaw length and metacarpal length measured on cleaned skeletal material. In some cases the data were split to consider different age classes separately, yielding a total of 15 traits studied. Genetic and environmental components of phenotypic variance were estimated for each trait and, for those traits showing nonzero heritability (N= 12), a QTL search was conducted by comparing a polygenic model with a model including a putative QTL. Support for a QTL at genome-wide significance was found on chromosome 11 for jaw length; suggestive QTL were found on chromosomes 2 and 5 (for birth date as a trait of the lamb), 8 (birth weight as a trait of the lamb), and 15 (adult hindleg length). We discuss the prospects for refining estimates of QTL position and effect size in the study population, and for QTL searches in free-living pedigrees in general.     

A Major Gene for Resting Metabolic Rate Unassociated with Body Composition: Results from the Québec Family Study

A major gene hypothesis for resting metabolic rate (RMR) was investigated using segregation analysis (POINTER) of data on families participating in Phase 2 of the Québec Family Study. Complete analyses were conducted on RMR adjusted for age, and also on RMR adjusted for age and other covariates, primarily fat mass (FM) and fat-free mass (FFM). Prior to adjustment for covariates, support for a major gene hypothesis was equivocal — i.e., there was evidence for either a major gene or a multifactorial component (i.e., polygenic and/or familial environment). The multifactorial model was preferred over the major gene model, although the latter did segregate according to Mendelian expectations. However, after the effects of FM and FFM were accounted for, a major gene effect was unambiguous and compelling. The putative locus accounted for 57% of the variance, affected 7% of the sample, and led to high values of RMR. The lack of a significant multifactorial effect suggested that the familial etiology of RMR adjusted for FM and FFM was likely to be entirely a function of the major locus. Comparing the RMR results from pre- and post-adjustment for FM and FFM suggests a plausible hypothesis. We know from earlier studies in this sample that there is a putative major gene for FM and a major non-Mendelian effect for FFM. The current study leads us to speculate that: (1) the gene(s) affecting body size and body composition also may have an effect on RMR, and further (2) removal of the effect of the major gene(s) for body size and composition allowed for detection of an additional major gene affecting only the RMR. Thus, RMR appears to be an oligogenic trait.     

Repetitive mating and female fitness in Dysdercus cardinalis (Hemiptera: Pyrrhocoridae)

Using three mating-treatment groups–pairing with a male for life, pairing with a male from adult emergence to first oviposition, and no pairing (virgin)–of Dysdercus cardinalis females in a laboratory study, the following results were obtained. Most of the virgin females did not lay eggs; those that did oviposited for the first time at a considerably older age than females in the other two groups. In a lifetime, females pairing with a male for life and those pairing with a male up to first oviposition produced essentially the same number of eggs, and this was more than 3.5 times the number of eggs laid by a virgin female. Virgin females produced an average of 1.3 clutches in a lifetime, compared with about 4.5 clutches produced by females in the other two groups. Reproductive span was significantly shorter, and fecundity per day of reproductive span significantly greater, for females pairing with males for life than for those pairing with a male up to first oviposition. Age at death was significantly different amongst females in different mating-treatment groups. Virgin females had the longest life, followed by females pairing with a male up to first oviposition and females pairing with a male for life, in that order. There was a positive correlation between lifetime gross reproduction and age at death for females pairing with a male for life. There was no relationship between these two characters for females pairing with a male up to first oviposition. Both females pairing with a male for life and those pairing with a male to first oviposition exhibited a significant inverse relationship between fecundity per day and age at death. The results obtained indicate that (1) mating is a prerequisite for normal gonadal activity in Dysdercus cardinalis , and (2) repetitive mating increases the rate of reproduction. This would reduce the mean age of parenthood which is inversely related to the intrinsic rate of increase.     

Screening for functional expression and overexpression of a family of diiron-containing interfacial membrane proteins using the univector recombination system

The large number of uncharacterized genes emerging from genome sequencing projects has resulted in a need for quick and reliable screening methods for protein expression parameters. We have utilized the univector plasmid recombination system (as previously reported) to develop a series of vectors for rapid screening for expression in Escherichia coli. A high level of recombinant protein expression is a requirement for purification of protein for structural determination and other purposes. In other applications, successful complementation of a missing enzyme activity in E. coli, as well as directed evolution studies and metabolic engineering, often require a much lower level of protein expression. In this report we describe the construction of a number of new pHOST vectors that can be screened for both low- and high-level expression. We isolated a mutant vector for MBP fusions that exhibited a more optimal level of expression for complementation of aerobic respiration in hemA(-) E. coli, our functional assay for the alternative oxidase. We then demonstrated the use of our system to rapidly screen for both optimal functional expression and optimal overexpression of the alternative oxidase as well as two other members of a family of membrane-bound diiron carboxylate proteins, the plastid terminal oxidase and 5-demethoxyquinone hydroxylase.     

Standardizing cancer biomarkers criteria: data elements as a foundation for a database. Inflammatory mediator/M-CSF as model marker

The purpose of this position article was to design a set of criteria (data elements) for a wide range of cancer biomarkers (CBs) in an attempt to standardize biomarkers features through a common language as a foundation for a database. Data elements are described as a set of generic criteria, which should characterize nearly all biomarkers introduced in the literature. Data elements were extracted from the review of prominent features that biomarkers represent within various categories. The extracted characteristics of biomarkers produced a short list of shared and unique generic features such as biological nature and history; stage/phase of study; sensitivity and specificity; modes of action; risk assessment; validation status; technology, and recommendation status for diversified biomarkers. To tailor data elements on specific markers, a cytokine, such as macrophage-colony stimulating factor (M-CSF), which has been proposed as a ‹potentially suitable biomarker’ for diagnosis of ovarian, lung, breast, pancreatic, and colorectal cancers, was selected as a Model biomarker. Small scale clinical studies suggested the superior usefulness of M-CSF compared with traditional markers for cancer detection. A key criterion for selecting Model marker and tailoring data elements for detection of cancer was the comparison of data on its specificity and sensitivity with traditional markers. The design of data elements for standardizing CBs criteria is considered a Research Tool and a foundation for developing a comprehensive CBs database useful for oncology researchers for a wide range of biomarkers. Validation, integration and proper packaging, data visualization and recommendation of suitability of CBs, by a panel of experts, for technology development are important challenging next steps toward developing a reliable database, which would allow professionals to effectively retrieve and study integrated information on potentially useful markers; identify important knowledge gaps and limitations of data; and assess state of technologies and commercialization of markers at a point of need. Appropriate use of integrated information on biomarkers in clinical practices would eventually account for more cost-effective characteristics of an individual’s state of health.     

A reduced muscle model and planar musculoskeletal model fit for the simulation of whole-body movements

Musculoskeletal models are made to reflect the capacities of the human body in general, and often a specific subject in particular. It remains challenging to both model the musculoskeletal system and then fit the modelled muscles to a specific human subject. We present a reduced muscle model, a planar musculoskeletal model, and a fitting method that can be used to find a feasible set of active and passive muscle parameters for a specific subject. At a minimum, the fitting method requires inverse dynamics data of the subject, a scalar estimate of the peak activation reached during the movement, and a plausible initial estimate for the strength and flexibility of that subject. While additional data can be used to result in a more accurate fit, this data is not required for the method solve for a feasible fit. The minimal input requirements of the proposed fitting method make it well suited for subjects who cannot undergo a maximum voluntary contraction trial, or for whom recording electromyographic data is not possible. To evaluate the model and fitting method we adjust the musculoskeletal model so that it can perform an experimentally recorded stoop-lift of a 15 kg box.     

Site-dependent Life-Cycle Impact Assessment in Sweden (5 pp)

Goal, Scope and Background Although LCA is traditionally a site-independent tool, there is currently a trend towards making LCA more site-dependent if not site-specific. For Europe, site-dependent impact factors have been calculated on a country basis for acidification, terrestrial eutrophication and toxicological impacts. It is, however, an open question whether this is the optimum level for site-dependent factors. The aim of this paper is to develop site-dependent characterisation factors for different parts of Sweden for air emissions of NOx, SOx and particulates regarding ecosystem and human health impacts. Based on experiences from a case-study, the usability of the site-dependent factors for LCA are discussed, as well as the appropriate level of site-dependency for ecosystem and human health impacts. Method logy. The Ecosense model is used for calculating site-dependent factors for some atmospheric pollutants. Characterisation factors are calculated for four different places in Sweden with two different stack heights. Results and Conclusions The characterisation factors for ecosystem impacts show fairly small differences between different parts of Sweden (within a factor of two). For health impacts, the differences between different parts of the country were larger and more significant (up to one order of magnitude). Also the difference between low and high stack heights may be relevant, especially in densely populated areas. These results suggest that for ecosystems, site-dependent characterisation factors for the considered atmospheric pollutants on a country level may be sufficient for most applications. However, for health impacts, site-dependent factors on a country level may be inappropriate. Beside LCA, the calculated factors and the methodology used should also be useful for other environmental system analysis tools, such as Strategic Environmental Assessment, Cost-Benefit Analysis and Environmental Management Systems.     

Insulin production from hiPSC-derived pancreatic cells in a novel wicking matrix bioreactor

Clinical use of pancreatic β islets for regenerative medicine applications requires mass production of functional cells. Current technologies are insufficient for large-scale production in a cost-efficient manner. Here, we evaluate advantages of a porous cellulose scaffold and demonstrate scale-up to a wicking matrix bioreactor as a platform for culture of human endocrine cells. Scaffold modifications were evaluated in a multiwell platform to find the optimum surface condition for pancreatic cell expansion followed by bioreactor culture to confirm suitability. Preceding scale-up, cell morphology, viability, and proliferation of primary pancreatic cells were evaluated. Two optimal surface modifications were chosen and evaluated further for insulin secretion, cell morphology, and viable cell density for human-induced pluripotent stem cell-derived pancreatic cells at different stages of differentiation. Scale-up was accomplished with uncoated, amine-modified cellulose in a miniature bioreactor, and insulin secretion and cell metabolic profiles were determined for 13 days. We achieved 10-fold cell expansion in the bioreactor along with a significant increase in insulin secretion compared with cultures on tissue culture plastic. Our findings define a new method for expansion of pancreatic cells a on wicking matrix cellulose platform to advance cell therapy biomanufacturing for diabetes.     

A generalized heterozygote deficiency assessed with microsatellites in French common ash populations

Common ash is a temperate forest tree with a colonizing behaviour, a discontinuous spatial distribution and a peculiar and poorly known mating system. Microsatellite markers were used to study the genetic structure in natural populations of common ash. Twelve populations located in northeastern France were analysed at five loci. Levels of genetic variability within and among stands were estimated for the seedling and adult stages. As expected for a forest tree, our results reveal high levels of intrapopulation diversity and a low genetic differentiation between stands. However, a general and significant heterozygote deficiency was found, with a mean F(IS) of 0.163 for the seedlings and of 0.292 for the adult trees. The different explanations for such an excess homozygosity are discussed: a nonMendelian inheritance of alleles, the presence of null alleles, a Wahlund effect and assortative mating.     

Quantitative T2 mapping for detection and quantification of thrombophlebitis in a rabbit model

Short peripheral catheter thrombophlebitis (SPCT), a sterile inflammation of the vein wall, is the most common complication associated with short peripheral catheters (SPCs) and affects up to 80% of hospitalized patients receiving IV therapy. Extensive research efforts have been devoted for improvement and optimization of the catheter material, but means for examination of any novel design are limited, inaccurate and require costly comprehensive pre-clinical and clinical trials. Therefore, there is a conclusive need for a reliable quantitative method for evaluation of SPCT, in particular for research purposes examining the thrombophlebitis-related symptoms of any novel catheter design. In this study, we developed for the first time a quantitative MRI based tool for evaluation of SPCT. The extent and severity of SPCT caused by two different commercially available SPCs with known predisposition for thrombophlebitis, were studied in a rabbit model. MRI analysis was consistent with the standardized pathology evaluation and showed remarkable difference in the percent of edema between the experimental groups. These differences were in line with previous studies and provide evidence that this type of analysis may be useful for future assessment of SPCT in vivo. As a non-invasive method, it may constitute a cost effective solution for examination of new catheters and other medical devices, thereby reducing the need for animal sacrifice.     

Loss of gametophytic self-incompatibility with evolution of inbreeding depression

Gametophytic self-incompatibility (SI) in plants is a widespread mechanism preventing self-fertilization and the ensuing inbreeding depression, but it often evolves to self-compatibility. We analyze genetic mechanisms for the breakdown of gametophytic SI, incorporating a dynamic model for the evolution of inbreeding depression allowing for partial purging of nearly recessive lethal mutations by selfing, and accounting for pollen limitation and sheltered load linked to the S-locus. We consider two mechanisms for the breakdown of gametophytic SI: a nonfunctional S-allele and an unlinked modifier locus that inactivates the S-locus. We show that, under a wide range of conditions, self-compatible alleles can invade a self-incompatible population. Conditions for invasion are always less stringent for a nonfunctional S-allele than for a modifier locus. The spread of self-compatible genotypes is favored by extremely high or low selfing rates, a small number of S-alleles, and pollen limitation. Observed parameter values suggest that the maintenance of gametophytic SI is caused by a combination of high inbreeding depression in self-incompatible populations coupled with intermediate selfing rates of the self-compatible genotypes and sheltered load linked to the S-locus.     

一个可用于复杂大群体的可视化系谱绘制和分析软件

对于在遗传研究和家系研究中大的系谱结构图还很难分析。系谱的绘制通常是遗传性状的分析研究的第一步。系图可以反映整个群体的结构、每个个体之间的相互关系以及基因流的走向,便于理解遗传性状的本质。因为所用家系数目的增大和复杂性的增加,绘制1个清晰的系谱有时变得十分困难。因此开发了1种名为Pedigraph软件,可以解决这个问题。Pedigraph能够完成对于大的复杂的群体的系谱绘制工作,并能进行相应的系谱分析。初步的测试表明这个软件在研究动植物的遗传育种中是1个有用的工具,同时它也可以用于人类的群体和历史等方面的研究。     

A Systems Model Approach to Determining Resilience Surrogates for Case Studies

Resilience theory offers a framework for understanding the dynamics of complex systems. However, operationalizing resilience theory to develop and test empirical hypotheses can be difficult. We present a method in which simple systems models are used as a framework to identify resilience surrogates for case studies. The process of constructing a systems model for a particular case offers a path for identifying important variables related to system resilience, including the slowly-changing variables and thresholds that often are keys to understanding the resilience of a system. We develop a four-step process for identifying resilience surrogates through development of systems models. Because systems model development is often a difficult step, we summarize four basic existing systems models and give examples of how each may be used to identify resilience surrogates. The construction and analysis of simple systems models provides a useful basis for guiding and directing the selection of surrogate variables that will offer appropriate empirical measures of resilience.     

Consistent Sets of Spectrophotometric Chlorophyll Equations for Acetone, Methanol and Ethanol Solvents

A set of equations for determining chlorophyll a (Chl a) and accessory chlorophylls b, c ₂ , c ₁ + c ₂ and the special case of Acaryochloris marina, which uses Chl d as its primary photosynthetic pigment and also has Chl a, have been developed for 90% acetone, methanol and ethanol solvents. These equations for different solvents give chlorophyll assays that are consistent with each other. No algorithms for Chl c compounds (c ₂ , c ₁ + c ₂) in the presence of Chl a have previously been published for methanol or ethanol. The limits of detection (and inherent error, ± 95% confidence limit), for chlorophylls in all organisms tested, was generally less than 0.1 μg/ml. The Chl a and b algorithms for green algae and land plants have very small inherent errors (< 0.01 μg/ml). Chl a and d algorithms for Acaryochloris marina are consistent with each other, giving estimates of Chl d/a ratios which are consistent with previously published estimates using HPLC and a rarely used algorithm originally published for diethyl ether in 1955. The statistical error structure of chlorophyll algorithms is discussed. The relative error of measurements of chlorophylls increases hyperbolically in diluted chlorophyll extracts because the inherent errors of the chlorophyll algorithms are constants independent of the magnitude of absorbance readings. For safety reasons, efficient extraction of chlorophylls and the convenience of being able to use polystyrene cuvettes, the algorithms for ethanol are recommended for routine assays of chlorophylls. The methanol algorithms would be convenient for assays associated with HPLC work.     

An optimal shape problem related to the realistic design of river fishways

A river fishway is a hydraulic structure that facilitates fish in overcoming obstacles (dams, waterfalls, etc.) to their spawning and other migrations in rivers. In this work we present a mathematical formulation of an optimal design problem for a vertical slot fishway, where the state system is given by the 2D shallow water equations fixing the height and velocity of water, the design variables are the geometry of the slots, and the objective function is determined by the existence of rest areas for fish and of a water velocity suitable for fish swimming capability. We also derive an expression for the gradient of the objective function via the adjoint system. From the numerical point of view, we present a characteristic-Galerkin method for solving the shallow water equations, and a direct search algorithm for the computation of the optimal design variables. Finally, we give numerical results obtained for a standard ten pools channel.     

中国西南季风常绿阔叶林不同恢复阶段凋落物动态分析

为探索季风常绿阔叶林不同恢复阶段群落凋落物产量及其动态变化规律,于云南普洱地区通过设置凋落物承接网并定期收集网内的凋落物,对中国西南季风常绿阔叶林区不同恢复阶段群落凋落物产量及其动态变化进行研究。结果表明:(1)不同恢复阶段季风常绿阔叶林年凋落物总量在8 133.1~8 798.3kg/hm2之间,年凋落物总量大小关系为恢复30年群落老龄林群落恢复40年群落。其中叶凋落量最高,其次为枝凋落量,两者贡献量超过总凋落量的3/4。(2)季风常绿阔叶林不同恢复阶段群落中凋落物随时间的动态变化趋势大致相同,2月份达到高峰值,随后逐渐减少,在9月份降至最低,随后又有所升高,为单峰或多峰曲线;在不同凋落物组分凋落量时间动态上,不同恢复阶段群落叶凋落量随月份变化均为单峰曲线;枝凋落量在恢复群落中为单峰曲线,而在老龄林中则为多峰曲线;皮凋落量随月份的变化在恢复30年及老龄林群落间均为单峰曲线,但在恢复40年群落中为平缓曲线,月份间变化不大;繁殖体凋落物在恢复30年及老龄林群落间均为多峰曲线,但在恢复40年群落中为单峰曲线;半分解物凋落量在恢复30年及老龄林群落中随月份呈单峰曲线,在恢复40年群落中则为多峰曲线。(3)在短刺栲、刺栲和红木荷3种优势物种中,短刺栲叶片年凋落量在所有群落中均最大(分别占恢复30年群落的53.93%、恢复40年群落的47.83%、老龄林的28.32%),红木荷次之(分别占恢复30年群落的8.45%、恢复40年群落的10.71%、老龄林的31.69%),刺栲最少(分别占恢复30年群落的6.1%、恢复40年群落的7.53%、老龄林的6.36%)。短刺栲叶凋落量随月份的变化在恢复群落中呈单峰曲线,而在老龄林中则呈现多峰曲线;红木荷在3种群落中则均为单峰曲线;刺栲则是在恢复30年及老龄林中呈单峰曲线,而在恢复40年群落中呈多峰曲线。     

Observations on the Staining of Corynebacterium Diphtheriae

Albert's method, of staining diphtheria cultures consists of staining a fixed smear for one minute (some laboratories stain for five minutes) with a solution containing toluidine blue and malachite (or methyl) green, washing with water, and applying Albert's iodine for one minute. This procedure is discussed and criticized, and in addition the mechanism of the stain is elucidated. Also, the procedure which involves staining a fixed smear for one minute with Loeffler's alkaline methylene blue solution is discussed and criticized.

To overcome the objections to the above staining methods, a different method is proposed. This consists of staining a fixed smear with an acid solution of toluidine blue, washing with water, applying Albert's iodine for one minute, washing with water, and finally applying a safranin solution for 15-20 seconds. The theoretical basis for this method is presented.