The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

Testing the principle of equivalence by supernova neutrinos

We study the possible impact of the neutrino oscillation which could be induced by a tiny violation of equivalence principle (VEP) on the observation of neutrinos emitted from supernova driven by gravitational collapse. We show that using supernova neutrinos, one can probe very small values of VEP parameters, delta(tau) less than or similar to O(10(-31)) for massless or degenerated neutrinos and delta(tau) less than or similar to O(10(-16)) x (Deltam(2)/10(-5) eV(2)) for massive neutrinos. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Quantum equivalence between the self-dual and the Maxwell-Chern-Simons models nonlinearly coupled to U(1) scalar fields

The use of master actions to prove duality at quantum level becomes cumbersome if one of the dual fields interacts nonlinearly with other fields. This is the case of the theory considered here consisting of U(1) scalar fields coupled to a self-dual field through a linear and a quadratic term in the self-dual field. Integrating perturbatively over the scalar fields and deriving effective actions for the self-dual and the gauge field we are able to consistently neglect awkward extra terms generated via master action and establish quantum duality up to cubic terms in the coupling constant. The duality holds for the partition function and some correlation functions. The absence of ghosts imposes restrictions on the coupling with the scalar fields.

K-bi-Lipschitz equivalence of real function-germs

In this paper we prove that the set of equivalence classes of germs of real polynomials of degree less than or equal to k, with respect to K-bi-Lipschitz equivalence, is finite.

Statistical methods applied to the analysis of blasting vibrations in a diabase quarry in the city of Campinas, (SP) Brazil

Statistical methods of multiple regression analysis, trend surface analysis and principal components analysis were applied to seismographic data recorded during production blasting at a diabase quarry in the urban area of Campinas (SP), Brazil. The purpose of these analyses was to determine the influence of the following variables: distance (D), charge weight per delay (W), and scaled distance (SD) associated with properties of the rock body (orientation, frequency and angle of geological discontinuities; depth of bedrock and thickness of the soil overburden) in the variation of the peak particle velocity (PPV). This approach yielded variables with larger influences (loads) on the variation of ground vibration, as well as behavior and space tendency of this variation. The results showed a better relationship between PPV and D, with D being the most important factor in the attenuation of the ground vibrations. The geological joints and the depth to bedrock have a larger influence than the explosive charges in the variation of the vibration levels, but frequencies appear to be more influenced by the amount of soil overburden.

CLASSICAL EQUIVALENCE OF LAMBDA-R-PHI-2 THEORIES

It is shown that the action functional S[g, phi] = integral d4 x square-root -g[R/k(1 + klambdaphi2) + partial derivative(mu)phi partial derivative(mu)phi] describes, in general, one and the same classical theory whatever may be the value of the coupling constant lambda.

DIAGNOSIS of REGIONAL CEREBRAL BLOOD FLOW ABNORMALITIES USING SPECT: AGREEMENT BETWEEN INDIVIDUALIZED STATISTICAL PARAMETRIC MAPS and VISUAL INSPECTION BY NUCLEAR MEDICINE PHYSICIANS WITH DIFFERENT LEVELS of EXPERTISE IN NUCLEAR NEUROLOGY

INTRODUCTION: Visual analysis is widely used to interpret regional cerebral blood flow (rCBF) SPECT images in clinical practice despite its limitations. Automated methods are employed to investigate between-group rCBF differences in research Studies but have rarely been explored in individual analyses.OBJECTIVES: To compare visual inspection by nuclear physicians with the automated statistical parametric mapping program using a SPECT dataset of patients with neurological disorders and normal control images.METHODS: Using statistical parametric mapping, 14 SPECT images from patients with various neurological disorders were compared individually with a databank of 32 normal images using a statistical threshold of p<0.05 (corrected for multiple comparisons at the level of individual voxels or clusters). Statistical parametric mapping results were compared with Visual analyses by a nuclear physician highly experienced in neurology (A) as well as a nuclear physician with a general background of experience (B) who independently classified images as normal or altered, and determined the location of changes and the severity.RESULTS: of the 32 images of the normal databank, 4 generated maps showing rCBF abnormalities (p<0.05, corrected). Among the 14 images from patients with neurological disorders, 13 showed rCBF alterations. Statistical parametric mapping and physician A completely agreed on 84.37% and 64.28% of cases from the normal databank and neurological disorders, respectively. The agreement between statistical parametric mapping and ratings of physician B were lower (71.18% and 35.71%, respectively).CONCLUSION: Statistical parametric mapping replicated the findings described by the more experienced nuclear physician. This finding suggests that automated methods for individually analyzing rCBF SPECT images may be a valuable resource to complement visual inspection in clinical practice.

Optimum design of experiments for statistical inference

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Anaerobic threshold estimation by statistical modelling.

Anaerobic threshold (AT) is usually estimated as a change point problem by visual analysis of the cardiorespiratory response to incremental dynamic exercise. In this study, two phase linear (TPL) models of the linear-linear and linear-quadratic type were used for the estimation of AT. The correlation coefficient between the classical and statistical approaches was 0.88, and 0.89 after outlier exclusion. The TPL models provide a simple method for estimating AT that can be easily implemented using a digital computer for the automatic pattern recognition of AT.

On the equivalence of the mechanisms governing switching and hysteresis phenomena

The nucleation and growth model, which is usually applied to switching phenomena, is adapted for explaining surface potential measurements on the P(VDF-TrFE) (polyvinylidene fluoride-trifluoroethylene) copolymer obtained in a constant current corona triode. It is shown that the growth is one-dimensional and that the nucleation rate is unimportant, probably because surface potential measurements take much longer than the switching ones. The surface potential data can therefore be accounted for by a growth model in which the velocity of growth varies exponentially with the electric field. Since hysteresis loops can be obtained from surface potential measurements, it is suggested that similar mechanisms can be used when treating switching and hysteresis phenomena, provided that account is taken of the difference in the time scale of the measurements.

Determination and statistical treatment of short-circuit currents for faults on the lines of networks of 440 kV CESP system

The distribution of short-circuit current is investigated by means of two methods, one direct and the other analytic; both methods consider uniform probability distribution of line faults. In the direct method, the procedure consists of calculating fault currents at equidistant points along the line, starting from one of the end points and considering the other end open. The magnitude of the current is classified according to Brazilian standards (regulation NBR-7118). The analytic method assumes that the distribution of short-circuit currents through the busbar and the distribution of the line length connected to it are known, as well as the independence of values. The method is designed to determine the probability that fault currents through a line will surpass the pre-established magnitude, thus generating frequency distribution curves of short-circuit currents along the lines.

Remarks on fermion-boson equivalence in three dimensions

Starting from a decomposition of the self-dual field in (2 + 1) dimensions, we build up an alternative quantum theory which consists of a self-dual model coupled to a Maxwell-generalized Chern-Simons theory. We discuss the fermion-boson equivalence of this quantum theory by comparing it with the Thirring model. Using these results we were able to compute the mass of the bosonized fermions up to third order in 1/m. Some problems related to the number of poles of the effective propagator are also addressed.

Strangeness content and structure function of the nucleon in a statistical quark model

The strangeness content of the nucleon is determined from a statistical model using confined quark levels, and is shown to have a good agreement with the corresponding values extracted from experimental data. The quark levels are generated in a Dirac equation that uses a linear confining potential (scalar plus vector). With the requirement that the result for the Gottfried sum rule violation, given by the New Muon Collaboration (NMC), is well reproduced, we also obtain the difference between the structure functions of the proton and neutron, and the corresponding sea quark contributions.