Influence diagnostics in nonlinear mixed-effects elliptical models

In this work we propose and analyze nonlinear elliptical models for longitudinal data, which represent an alternative to gaussian models in the cases of heavy tails, for instance. The elliptical distributions may help to control the influence of the observations in the parameter estimates by naturally attributing different weights for each case. We consider random effects to introduce the within-group correlation and work with the marginal model without requiring numerical integration. An iterative algorithm to obtain maximum likelihood estimates for the parameters is presented, as well as diagnostic results based on residual distances and local influence [Cook, D., 1986. Assessment of local influence. journal of the Royal Statistical Society - Series B 48 (2), 133-169; Cook D., 1987. Influence assessment. journal of Applied Statistics 14 (2),117-131; Escobar, L.A., Meeker, W.Q., 1992, Assessing influence in regression analysis with censored data, Biometrics 48, 507-528]. As numerical illustration, we apply the obtained results to a kinetics longitudinal data set presented in [Vonesh, E.F., Carter, R.L., 1992. Mixed-effects nonlinear regression for unbalanced repeated measures. Biometrics 48, 1-17], which was analyzed under the assumption of normality. (C) 2009 Elsevier B.V. All rights reserved.

A Lagrangian relaxation approach to a coupled lot-sizing and cutting stock problem

Industrial production processes involving both lot-sizing and cutting stock problems are common in many industrial settings. However, they are usually treated in a separate way, which could lead to costly production plans. In this paper, a coupled mathematical model is formulated and a heuristic method based on Lagrangian relaxation is proposed. Computational results prove its effectiveness. (C) 2009 Elsevier B.V. All rights reserved.

QCD sum rules for the production of the X(3872) as a mixed molecule-charmonium state in B meson decay

We use QCD sum rules to calculate the branching ratio for the production of the meson X(3872) in the decay B -> X(3872)K, assumed to be a mixture between charmonium and exotic molecular vertical bar c (q) over bar vertical bar vertical bar q (c) over bar vertical bar states with J(PC) = 1(++). We find that in a small range for the values of the mixing angle, 5 degrees <= theta <= 13 degrees, we get the branching ratio B(B -> XK) = (1.00 +/- 0.68) x 10(-5), which is in agreement with the experimental upper limit. This result is compatible with the analysis of the mass and decay width of the mode J/psi(n pi) and the radiative decay mode J/psi gamma performed in the same approach. (C) 2011 Elsevier B.V. All rights reserved.

Raman spectroscopy analysis of structural photoinduced changes in GeS(2) + Ga(2)O(3) thin films

Photoexpansion and photobleaching effects have been observed in amorphous GeS(2) + Ga(2)O(3) (GGSO) thin films, when their surfaces were exposed to UV light. The photoinduced changes on the surface of the samples are indications that the structure has been changed as a result of photoexcitation. In this paper, micro-Raman, energy dispersive X-ray analysis (EDX) and backscattering electrons (BSE) microscopy were the techniques used to identify the origin of these effects. Raman spectra revealed that these phenomena are a consequence of the Ge-S bonds` breakdown and the formation of new Ge-O bonds, with an increase of the modes associated with Ge-O-Ge bonds and mixed oxysulphide tetrahedral units (S-Ge-O). The chemical composition measured by EDX and BSE microscopy images indicated that the irradiated area is oxygen rich. So, the present paper provides fundamental insights into the influence of the oxygen within the glass matrix on the considered photoinduced effects. (C) 2010 Elsevier B.V. All rights reserved.

Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

A note on influence diagnostics in nonlinear mixed-effects elliptical models

This paper provides general matrix formulas for computing the score function, the (expected and observed) Fisher information and the A matrices (required for the assessment of local influence) for a quite general model which includes the one proposed by Russo et al. (2009). Additionally, we also present an expression for the generalized leverage on fixed and random effects. The matrix formulation has notational advantages, since despite the complexity of the postulated model, all general formulas are compact, clear and have nice forms. (C) 2010 Elsevier B.V. All rights reserved.

Designing Novel Antitrypanosomal Agents from a Mixed Graph-Theoretical Substructural Approach

Chagas disease is nowadays the most serious parasitic health problem. This disease is caused by Trypanosoma cruzi. The great number of deaths and the insufficient effectiveness of drugs against this parasite have alarmed the scientific community worldwide. In an attempt to overcome this problem, a model for the design and prediction of new antitrypanosomal agents was obtained. This used a mixed approach, containing simple descriptors based on fragments and topological substructural molecular design descriptors. A data set was made up of 188 compounds, 99 of them characterized an antitrypanosomal activity and 88 compounds that belong to other pharmaceutical categories. The model showed sensitivity, specificity and accuracy values above 85%. Quantitative fragmental contributions were also calculated. Then, and to confirm the quality of the model, 15 structures of molecules tested as antitrypanosomal compounds (that we did not include in this study) were predicted, taking into account the information on the abovementioned calculated fragmental contributions. The model showed an accuracy of 100% which means that the ""in silico"" methodology developed by our team is promising for the rational design of new antitrypanosomal drugs. (C) 2009 Wiley Periodicals, Inc. J Comput Chem 31: 882-894. 2010

Probing the dependence of the properties of cellulose acetates and their films on the degree of biopolymer substitution: use of solvatochromic indicators and thermal analysis

Although cellulose acetates, CAs, are extensively employed there is scant information about the systematic dependence of their properties on their degree of substitution, DS; this is the subject of the present work. Nine CAs samples, DS from 0.83 to 3.0 were synthesized; their films were prepared. The following solvatochromic probes have been employed in order to determine the empirical polarity, E (T)(33); ""acidity, alpha""; ""basicity, beta"", and ""dipolarity/polarizability, pi*"" of the casted films: 2,6-dichloro-4-(2,4,6-triphenyl-pyridinium-1-yl) phenolate, WB; 4-nitroaniline; 4-nitroanisole; 4-nitro-N,N-dimethylaniline; 2,6-diphenyl-4-(2,4,6-triphenyl-pyridinium-1-yl)phenolate, RB. Additionally, two systems, ethanol plus ethyl acetate (EtOH-EtAc), and cellulose plus cellulose triacetate, CTA, were employed as models for CAs of different DS. Regarding the model systems, the following was observed: (i) For EtOH-EtAc, the dependence of all solvatochromic parameters on the ""equivalent-DS"" of the binary mixture was non-linear because of preferential solvation; (ii) The dependence of E (T)(33) on equivalent DS of the cellulose-CTA films is linear, but the slope is smaller than that of the corresponding plot for CAs. This is attributed to the more efficient hydrogen bonding in the model system, a conclusion corroborated by IR measurements. The dependence of solvatochromic parameters of CAs on their DS is described by the simple equations; a consequence of the substitution of the OH by the ester group. The thermal properties of bulk CAs samples were investigated by DSC and TGA; their dependence on DS is described by simple equations. The relevance of these data to the processing and applications of CAs is briefly discussed.

Financial time series analysis : Chaos and neurodynamics approach

This work aims at combining the Chaos theory postulates and Artificial Neural Networks classification and predictive capability, in the field of financial time series prediction. Chaos theory, provides valuable qualitative and quantitative tools to decide on the predictability of a chaotic system. Quantitative measurements based on Chaos theory, are used, to decide a-priori whether a time series, or a portion of a time series is predictable, while Chaos theory based qualitative tools are used to provide further observations and analysis on the predictability, in cases where measurements provide negative answers. Phase space reconstruction is achieved by time delay embedding resulting in multiple embedded vectors. The cognitive approach suggested, is inspired by the capability of some chartists to predict the direction of an index by looking at the price time series. Thus, in this work, the calculation of the embedding dimension and the separation, in Takens‘ embedding theorem for phase space reconstruction, is not limited to False Nearest Neighbor, Differential Entropy or other specific method, rather, this work is interested in all embedding dimensions and separations that are regarded as different ways of looking at a time series by different chartists, based on their expectations. Prior to the prediction, the embedded vectors of the phase space are classified with Fuzzy-ART, then, for each class a back propagation Neural Network is trained to predict the last element of each vector, whereas all previous elements of a vector are used as features.

Energy Efficiency through Thermal Energy Storage : Possibilities for the Swedish Building Stock

The need for heating and cooling in buildings constitutes a considerable part of the total energy use in a country and reducing this need is of outmost importance in order to reach national and international goals for reducing energy use and emissions. One important way of reaching these goals is to increase the proportion of renewable energy used for heating and cooling of buildings. Perhaps the largest obstacle with this is the often occurring mismatch between the availability of renewable energy and the need for heating or cooling, hindering this energy to be used directly. This is one of the problems that can be solved by using thermal energy storage (TES) in order to save the heat or cold from when it is available to when it is needed. This thesis is focusing on the combination of TES techniques and buildings to achieve increased energy efficiency for heating and cooling. Various techniques used for TES as well as the combination of TES in buildings have been investigated and summarized through an extensive literature review. A survey of the Swedish building stock was also performed in order to define building types common in Sweden. Within the scope of this thesis, the survey resulted in the selection of three building types, two single family houses and one office building, out of which the two residential buildings were used in a simulation case study of passive TES with increased thermal mass (both sensible and latent). The second case study presented in the thesis is an evaluation of an existing seasonal borehole storage of solar heat for a residential community. In this case, real measurement data was used in the evaluation and in comparisons with earlier evaluations. The literature reviews showed that using TES opens up potential for reduced energy demand and reduced peak heating and cooling loads as well as possibilities for an increased share of renewable energy to cover the energy demand. By using passive storage through increased thermal mass of a building it is also possible to reduce variations in the indoor temperature and especially reduce excess temperatures during warm periods, which could result in avoiding active cooling in a building that would otherwise need it. The analysis of the combination of TES and building types confirmed that TES has a significant potential for increased energy efficiency in buildings but also highlighted the fact that there is still much research required before some of the technologies can become commercially available. In the simulation case study it was concluded that only a small reduction in heating demand is possible with increased thermal mass, but that the time with indoor temperatures above 24 °C can be reduced by up to 20%. The case study of the borehole storage system showed that although the storage system worked as planned, heat losses in the rest of the system as well as some problems with the system operation resulted in a lower solar fraction than projected. The work presented within this thesis has shown that TES is already used successfully for many building applications (e.g. domestic hot water stores and water tanks for storing solar heat) but that there still is much potential in further use of TES. There are, however, barriers such as a need for more research for some storage technologies as well as storage materials, especially phase change material storage and thermochemical storage.

Digital spiral analysis for objective assessment of fine motor timing variability in Parkinson's disease

OBJECTIVES: To develop a method for objective assessment of fine motor timing variability in Parkinson’s disease (PD) patients, using digital spiral data gathered by a touch screen device. BACKGROUND: A retrospective analysis was conducted on data from 105 subjects including65 patients with advanced PD (group A), 15 intermediate patients experiencing motor fluctuations (group I), 15 early stage patients (group S), and 10 healthy elderly subjects (HE) were examined. The subjects were asked to perform repeated upper limb motor tasks by tracing a pre-drawn Archimedes spiral as shown on the screen of the device. The spiral tracing test was performed using an ergonomic pen stylus, using dominant hand. The test was repeated three times per test occasion and the subjects were instructed to complete it within 10 seconds. Digital spiral data including stylus position (x-ycoordinates) and timestamps (milliseconds) were collected and used in subsequent analysis. The total number of observations with the test battery were as follows: Swedish group (n=10079), Italian I group (n=822), Italian S group (n = 811), and HE (n=299). METHODS: The raw spiral data were processed with three data processing methods. To quantify motor timing variability during spiral drawing tasks Approximate Entropy (APEN) method was applied on digitized spiral data. APEN is designed to capture the amount of irregularity or complexity in time series. APEN requires determination of two parameters, namely, the window size and similarity measure. In our work and after experimentation, window size was set to 4 and similarity measure to 0.2 (20% of the standard deviation of the time series). The final score obtained by APEN was normalized by total drawing completion time and used in subsequent analysis. The score generated by this method is hence on denoted APEN. In addition, two more methods were applied on digital spiral data and their scores were used in subsequent analysis. The first method was based on Digital Wavelet Transform and Principal Component Analysis and generated a score representing spiral drawing impairment. The score generated by this method is hence on denoted WAV. The second method was based on standard deviation of frequency filtered drawing velocity. The score generated by this method is hence on denoted SDDV. Linear mixed-effects (LME) models were used to evaluate mean differences of the spiral scores of the three methods across the four subject groups. Test-retest reliability of the three scores was assessed after taking mean of the three possible correlations (Spearman’s rank coefficients) between the three test trials. Internal consistency of the methods was assessed by calculating correlations between their scores. RESULTS: When comparing mean spiral scores between the four subject groups, the APEN scores were different between HE subjects and three patient groups (P=0.626 for S group with 9.9% mean value difference, P=0.089 for I group with 30.2%, and P=0.0019 for A group with 44.1%). However, there were no significant differences in mean scores of the other two methods, except for the WAV between the HE and A groups (P<0.001). WAV and SDDV were highly and significantly correlated to each other with a coefficient of 0.69. However, APEN was not correlated to neither WAV nor SDDV with coefficients of 0.11 and 0.12, respectively. Test-retest reliability coefficients of the three scores were as follows: APEN (0.9), WAV(0.83) and SD-DV (0.55). CONCLUSIONS: The results show that the digital spiral analysis-based objective APEN measure is able to significantly differentiate the healthy subjects from patients at advanced level. In contrast to the other two methods (WAV and SDDV) that are designed to quantify dyskinesias (over-medications), this method can be useful for characterizing Off symptoms in PD. The APEN was not correlated to none of the other two methods indicating that it measures a different construct of upper limb motor function in PD patients than WAV and SDDV. The APEN also had a better test-retest reliability indicating that it is more stable and consistent over time than WAV and SDDV.

Parameter estimation of Poisson generalized linear mixed models based on three different statistical principles : a simulation study

Generalized linear mixed models are flexible tools for modeling non-normal data and are useful for accommodating overdispersion in Poisson regression models with random effects. Their main difficulty resides in the parameter estimation because there is no analytic solution for the maximization of the marginal likelihood. Many methods have been proposed for this purpose and many of them are implemented in software packages. The purpose of this study is to compare the performance of three different statistical principles - marginal likelihood, extended likelihood, Bayesian analysis-via simulation studies. Real data on contact wrestling are used for illustration.

Women's acceptability of misoprostol treatment for incomplete abortion by midwives and physicians : secondary outcome analysis from a randomized controlled equivalence trial at district level in Uganda

OBJECTIVE: This study aimed to assess women´s acceptability of diagnosis and treatment of incomplete abortion with misoprostol by midwives, compared with physicians. METHODS: This was an analysis of secondary outcomes from a multi-centre randomized controlled equivalence trial at district level in Uganda. Women with first trimester incomplete abortion were randomly allocated to clinical assessment and treatment with misoprostol by a physician or a midwife. The randomisation (1:1) was done in blocks of 12 and stratified for health care facility. Acceptability was measured in expectations and satisfaction at a follow up visit 14-28 days following treatment. Analysis of women's overall acceptability was done using a generalized linear mixed-effects model with an equivalence range of -4% to 4%. The study was not masked. The trial is registered at ClinicalTrials.org, NCT 01844024. RESULTS: From April 2013 to June 2014, 1108 women were assessed for eligibility of which 1010 were randomized (506 to midwife and 504 to physician). 953 women were successfully followed up and included in the acceptability analysis. 95% (904) of the participants found the treatment satisfactory and overall acceptability was found to be equivalent between the two study groups. Treatment failure, not feeling calm and safe following treatment, experiencing severe abdominal pain or heavy bleeding following treatment, were significantly associated with non-satisfaction. No serious adverse events were recorded. CONCLUSIONS: Treatment of incomplete abortion with misoprostol by midwives and physician was highly, and equally, acceptable to women. TRIAL REGISTRATION: ClinicalTrials.gov NCT01844024.

The effects of Ikea and external shopping centers on regional trade and retail growth : A comparative analysis of Ikea’s establishment in Kalmar and Karlstad

This paper examines the effects of Ikea store establishment in Kalmar and Karlstad on the trade and retail inside the two cities, and as well on the trade and retail in the close neighboring municipalities and in further peripheral municipalities in both regions. After the establishment of Ikea store, Kalmar and Karlstad have experienced significant growth in trade and retail. The question, however, is how big this growth is in both cities? And how different locations on different distances from Ikea have been affected? What impact there was on different segments of the retail? How different business branches have been affected? How large the catchment area for the emerging new large-scale retail locations is? These questions, in addition to few others, are investigated in this paper. The thesis starts with an introduction chapter containing a background of the topic, problem description, the investigated questions, the purpose, and the outline of the paper. The next chapter includes the frame of reference which consists of literature review and theoretical framework about the external shopping centers and their impact on retail and regional trade development. It includes also information gathered from previous studies technical reports and other available sources about the subject. The third chapter includes description for the methods used to collect the primary and secondary data needed for the purpose of this study. Then the empirical framework which demonstrates the results of the conducted research followed by analysis and concluded in discussion and conclusion. Mixed methods are used as research strategy in this thesis, and the method to conduct the research is based on telephone interviews for the primary (qualitative) data, and documents and desk research for the secondary (quantitative) data. The gathered data is analyzed and designed in a way that allows the usage of comparative analysis technique to present the findings and draw conclusions. The results showed that new established Ikea retail store outside the city boundaries results with many effects on the city center and on the neighboring municipalities as well. The city center seems not to be affected negatively, but on the contrary positive effects were witnessed in both regions, these positive effects are linked to the increase inflow of customers from the external retail area which is known as spillover effect. III On the other hand, the neighboring towns and municipalities are more negatively affected especially with the trade of con-convenience goods as the consumers in these towns and municipalities start to go to the area of Ikea and the large external retail center to do their purchasing, the substitution effect is then said to be occurred. Moreover, the further far municipalities do not seem to be significantly affected by the establishment of Ikea. These effects whether positive or negative could be monitored by looking to few trade parameters such as the turnover, the sales index, and the consumers’ expenditure, these parameters can be very useful to measure the developments and changes in the trade and retail in a given place. 

The tax treatment of employee stock options : generous to a fault

In the 2000 budgets, both the federal and Ontario governments introduced changes to the tax treatment of employee stock options for the explicit purpose of making their tax treatment in Canada similar to or more favourable than that in the United States. The federal budget added a deferral, similar to that currently applicable to options granted by Canadian-controlled private corporations, for up to $100,000 per year of public company stock options. The Ontario budget introduced an exemption from tax for employees involved in research and development on the first $100,000 per year of employee benefits arising on the exercise of qualified stock options or on eligible capital gains arising from the sale of shares acquired by the exercise of eligible stock options. These proposals reflect the apparent acceptance by the two governments that there is a “brain drain” from Canada to the United States of knowledge workers in the “new” economy and that reductions in Canadian taxes should stem this drain. In the author’s view, the tax treatment of employee stock options, even without these changes, is overly generous. Both the federal and provincial proposals ignore the fact that most employee stock options are taxed more favourably in Canada than in the United States in any event. In particular, most employee stock option benefits in Canada are taxed at capital gains tax rates, whereas in the United States most are taxed at full rates. While the US Internal Revenue Code does provide capital gains tax treatment for certain employee stock option benefits, a number of preconditions must be met. Most important, the shares acquired pursuant to the options must be held for a minimum of one year after the option is exercised. In addition, there are monetary limits on the amount of options that qualify for capital gains treatment. In Canada, there are generally no holding period requirements or monetary limits that apply in order for the option holder to benefit from capital gains tax rates. Empirical evidence indicates that the vast majority of employees in the United States exercise their options and immediately sell the shares acquired. These “cashless exercises” do not benefit from capital gains treatment in the United States, whereas similar cashless exercises in Canada generally do. This empirical evidence suggests not only that the 2000 budget proposals are unwarranted, but also that the existing treatment of employee stock options in Canada is already more generous than that in the United States. This article begins with a theoretical “benchmark” for the taxation of employee stock options. The author suggests that employee stock options should be treated in the same manner as other income from employment. In theory, the value of the benefit should be included in income when the option is granted or vests. However, owing to the practical difficulty of valuing employee stock options, the theoretical benchmark proposed is that the value of the benefit (the difference between the fair market value of the shares acquired and the strike price under the option) be taxed when the shares are acquired, and the employer be entitled to a corresponding deduction. The employee stock option rules in Canada and the United States are then compared and contrasted with each other and the benchmark treatment. The article then examines the arguments that have been made for favourable treatment of employee stock options. Included in this critique is a review of the recent empirical work on the Canadian brain drain. Empirical studies suggest that the brain drain—if it exists at all—is small and that, despite what many newspapers and right-wing think-tanks would have us believe, lower taxes in the United States are not the cause. One study, concluding that taxes do have an effect on migration, suggests that even if Canada adopted a tax system identical to that in the United States, the brain drain would be reduced by a mere 10 percent. Indeed, even if Canada eliminated income tax altogether, it would not stop the brain drain. If governments here want to spend money in order to stem the brain drain, they should focus on other areas. For example, Canada produces fewer university graduates in the fields of mathematics, sciences, and engineering than any other G7 country except Italy. The short supply of university graduates in these fields, the apparent loss of top-calibre academics to US
universities, and the consequent lower levels of university research in these areas (an important spawning ground for new ideas in the “new” knowledge-based economy) suggest that Canada may be better served by devoting more resources to its university institutions, particularly in post-graduate programs, rather than continuing the current trend of budget cuts that universities have endured and may further endure if taxes are reduced.
As far as employee stock options are concerned, if Canada does want to look to the United States for guidance on tax reform (which it seems to do with increasing frequency of late), it should adopt the US rules applicable to nonstatutory options, which are close to the proposed benchmark treatment. In the absence of preferential tax treatment, employee stock options would still be included in compensation packages provided that there were sound business reasons for their use. No persuasive evidence has been put forward that the use of stock options, in the absence of tax incentives, is suboptimal. Indeed, the US experience suggests quite the opposite.