Replication Data Management,needs and solutions: an initial evaluation of conceptual approaches for integrating heterogeneous replication study data

Replication Data Management (RDM) aims at enabling the use of data collections from several iterations of an experiment. However, there are several major challenges to RDM from integrating data models and data from empirical study infrastructures that were not designed to cooperate, e.g., data model variation of local data sources. [Objective] In this paper we analyze RDM needs and evaluate conceptual RDM approaches to support replication researchers. [Method] We adapted the ATAM evaluation process to (a) analyze RDM use cases and needs of empirical replication study research groups and (b) compare three conceptual approaches to address these RDM needs: central data repositories with a fixed data model, heterogeneous local repositories, and an empirical ecosystem. [Results] While the central and local approaches have major issues that are hard to resolve in practice, the empirical ecosystem allows bridging current gaps in RDM from heterogeneous data sources. [Conclusions] The empirical ecosystem approach should be explored in diverse empirical environments.

Analysis of the impact of globalization and economic growth on food security in developing countries

A pesar de los importantes avances en la reducción del hambre, la seguridad alimentaria continúa siendo un reto de dimensión internacional. La seguridad alimentaria es un concepto amplio y multidimensional, cuyo análisis abarca distintas escalas y horizontes temporales. Dada su complejidad, la identificación de las causas de la inseguridad alimentaria y la priorización de las medias para abordarlas, son dos cuestiones que suscitan un intenso debate en la actualidad. El objetivo de esta tesis es evaluar el impacto de la globalización y el crecimiento económico en la seguridad alimentaria en los países en desarrollo, desde una perspectiva macro y un horizonte temporal a largo plazo. La influencia de la globalización se aborda de una manera secuencial. En primer lugar, se analiza la relación entre la inversión público-privada en infraestructuras y las exportaciones agrarias. A continuación, se estudia el impacto de las exportaciones agrarias en los indicadores de seguridad alimentaria. El estudio del impacto del crecimiento económico aborda los cambios paralelos en la distribución de la renta, y cómo la inequidad influye en el comportamiento de la seguridad alimentaria nacional. Además, se analiza en qué medida el crecimiento económico contribuye a acelerar el proceso de mejora de la seguridad alimentaria. Con el fin de conseguir los objetivos mencionados, se llevan a cabo varios análisis econométricos basados en datos de panel, en el que se combinan datos de corte transversal de 52 países y datos temporales comprendidos en el periodo 1991-2012. Se analizan tanto variables en niveles como variables en tasas de cambio anual. Se aplican los modelos de estimación de efectos variables y efectos fijos, ambos en niveles y en primeras diferencias. La tesis incluye cuatro tipos de modelos econométricos, cada uno de ellos con sus correspondientes pruebas de robustez y especificaciones. Los resultados matizan la importancia de la globalización y el crecimiento económico como mecanismos de mejora de la seguridad alimentaria en los países en desarrollo. Se obtienen dos conclusiones relativas a la globalización. En primer lugar, los resultados sugieren que la promoción de las inversiones privadas en infraestructuras contribuye a aumentar las exportaciones agrarias. En segundo lugar, se observa que las exportaciones agrarias pueden tener un impacto negativo en los indicadores de seguridad alimentaria. La combinación de estas dos conclusiones sugiere que la apertura comercial y financiera no contribuye por sí misma a la mejora de la seguridad alimentaria en los países en desarrollo. La apertura internacional de los países en desarrollo ha de ir acompañada de políticas e inversiones que desarrollen sectores productivos de alto valor añadido, que fortalezcan la economía nacional y reduzcan su dependencia exterior. En relación al crecimiento económico, a pesar del incuestionable hecho de que el crecimiento económico es una condición necesaria para reducir los niveles de subnutrición, no es una condición suficiente. Se han identificado tres estrategias adicionales que han de acompañar al crecimiento económico con el fin de intensificar su impacto positivo sobre la subnutrición. Primero, es necesario que el crecimiento económico sea acompañado de una distribución más equitativa de los ingresos. Segundo, el crecimiento económico ha de reflejarse en un aumento de inversiones en salud, agua y saneamiento y educación. Se observa que, incluso en ausencia de crecimiento económico, mejoras en el acceso a agua potable contribuyen a reducir los niveles de población subnutrida. Tercero, el crecimiento económico sostenible en el largo plazo parece tener un mayor impacto positivo sobre la seguridad alimentaria que el crecimiento económico más volátil o inestable en el corto plazo. La estabilidad macroeconómica se identifica como una condición necesaria para alcanzar una mayor mejora en la seguridad alimentaria, incluso habiéndose mejorado la equidad en la distribución de los ingresos. Por último, la tesis encuentra que los países en desarrollo analizados han experimentado diferentes trayectorias no lineales en su proceso de mejora de sus niveles de subnutrición. Los resultados sugieren que un mayor nivel inicial de subnutrición y el crecimiento económico son responsables de una respuesta más rápida al reto de la mejora de la seguridad alimentaria. ABSTRACT Despite the significant reductions of hunger, food security still remains a global challenge. Food security is a wide concept that embraces multiple dimensions, and has spatial-temporal scales. Because of its complexity, the identification of the drivers underpinning food insecurity and the prioritization of measures to address them are a subject of intensive debate. This thesis attempts to assess the impact of globalization and economic growth on food security in developing countries with a macro level scale (country) and using a long-term approach. The influence of globalization is addressed in a sequential way. First, the impact of public-private investment in infrastructure on agricultural exports in developing countries is analyzed. Secondly, an assessment is conducted to determine the impact of agricultural exports on food security indicators. The impact of economic growth focuses on the parallel changes in income inequality and how the income distribution influences countries' food security performance. Furthermore, the thesis analyzes to what extent economic growth helps accelerating food security improvements. To address the above mentioned goals, various econometric models are formulated. Models use panel data procedures combining cross-sectional data of 52 countries and time series data from 1991 to 2012. Yearly data are expressed both in levels and in changes. The estimation models applied are random effects estimation and fixed effects estimations, both in levels and in first differences. The thesis includes four families of econometric models, each with its own set of robustness checks and specifications. The results qualify the relevance of globalization and economic growth as enabling mechanisms for improving food security in developing countries. Concerning globalization, two main conclusions can be drawn. First, results showed that enhancing foreign private investment in infrastructures contributes to increase agricultural exports. Second, agricultural exports appear to have a negative impact on national food security indicators. These two conclusions suggest that trade and financial openness per se do not contribute directly to improve food security in development countries. Both measures should be accompanied by investments and policies to support the development of national high value productive sectors, to strengthen the domestic economy and reduce its external dependency. Referring to economic growth, despite the unquestionable fact that income growth is a pre-requisite for reducing undernourishment, results suggest that it is a necessary but not a sufficient condition. Three additional strategies should accompany economic growth to intensifying its impact on food security. Firstly, it is necessary that income growth should be accompanied by a better distribution of income. Secondly, income growth needs to be followed by investments and policies in health, sanitation and education to improve food security. Even if economic growth falters, sustained improvements in the access to drinking water may still give rise to reductions in the percentage of undernourished people. And thirdly, long-term economic growth appears to have a greater impact on reducing hunger than growth regimes that combine periods of growth peaks followed by troughs. Macroeconomic stability is a necessary condition for accelerating food security. Finally, the thesis finds that the developing countries analyzed have experienced different non-linear paths toward improving food security. Results also show that a higher initial level of undernourishment and economic growth result in a faster response for improving food security.

A stochastic approximation algorithm with Markov chain Monte-Carlo method for incomplete data estimation problems

We propose a general procedure for solving incomplete data estimation problems. The procedure can be used to find the maximum likelihood estimate or to solve estimating equations in difficult cases such as estimation with the censored or truncated regression model, the nonlinear structural measurement error model, and the random effects model. The procedure is based on the general principle of stochastic approximation and the Markov chain Monte-Carlo method. Applying the theory on adaptive algorithms, we derive conditions under which the proposed procedure converges. Simulation studies also indicate that the proposed procedure consistently converges to the maximum likelihood estimate for the structural measurement error logistic regression model.

Efeito da suplementação prolongada de grão de soja cru e integral no pré-parto sobre o desempenho produtivo, qualidade oocitária e embrionária, e função imune de vacas leiteiras

Objetivou-se avaliar o efeito da suplementação prolongada de grão de soja cru e integral (GSI) como fonte de ácido graxo Ω6 sobre o desempenho produtivo, perfil metabólico, qualidade oocitária e embrionária e função imune de vacas leiteiras no período de transição e início de lactação. Foram selecionadas 44 vacas da raça Holandesa, multíparas e gestantes, com parto previsto para 90 dias após o início da avaliação e fornecimento das dietas experimentais, porém em razão da ocorrência de enfermidades metabólicas ou infecciosas (3 abortos; 3 deslocamentos de abomaso; 3 enfermidades podais; 4 distocias) 13 animais foram retirados do experimento. As vacas foram distribuídas aleatoriamente em quatro grupos experimentais diferindo entre eles o início do fornecimento de grão de soja cru e integral (GSI) durante o pré-parto. A dieta era baseada na inclusão de 12% de GSI %MS, com aproximadamente 5,1% de extrato etéreo (EE) o início de seu fornecimento foi conforme descrito a seguir: Grupo 0: Animais não receberam dieta contendo GSI no pré-parto; Grupo 30: Início do fornecimento de dieta com GSI nos 30 dias finais da gestação; Grupo 60: Início do fornecimento de dieta com GSI nos 60 dias finais da gestação; Grupo 90: Início do fornecimento de dieta com GSI nos 90 dias finais da gestação. Após o parto, todas as vacas receberam dieta única com 5,1% de EE, baseada na inclusão de 12% de GSI %MS até 90 dias de lactação. Os animais foram arraçoados de acordo com o consumo de matéria seca no dia anterior, de forma a ser mantido porcentual de sobras das dietas, diariamente, entre 5 e 10%. As amostras dos alimentos e sobras foram coletadas diariamente e armazenadas a -20ºC. Semanalmente as amostras coletadas diariamente foram misturadas e foi retirada uma amostra composta referente a um período de uma semana, a fim de mensurar o consumo de matéria seca e nutrientes. Amostras de fezes foram coletadas nos dias -56, -21, 21, 56 e 84 dias em relação ao parto, com o propósito de mensurar a digestibilidade da matéria seca e nutrientes. A produção de leite foi mensurada diariamente e para a composição dos teores de gordura, proteína, lactose e perfil de ácidos graxos amostras foram coletadas semanalmente. As amostras de sangue para análise dos metabólitos sanguíneos foram coletadas semanalmente. Amostras de sangue para mensurar a atividade do sistema imune foram coletadas na semanas -8, -4, -2, -1 em relação ao parto, parto, +1, +2, +4 e +8 semanas no período pós-parto. Nos dias 21, 42, 63 e 84 do período pós-parto foram realizadas aspirações foliculares, com posterior fertilização in vitro dos oócitos. Todas as variáveis mensuradas foram analisadas pelo procedimento PROC MIXED do SAS 9.4 através de regressão polinomial, utilizando efeito fixo de tratamento, semana, interação tratamento*semana e efeito de animal dentro de tratamento como aleatório. Utilizou nível de 5% de significância. Foi observado efeito (P<0,05) linear crescente para CEE no pré-parto. Não foi observado diferenças no CMS e nutrientes no pós-parto. Não houve alteração da digestibilidade nos períodos pré e pós-parto. Não houve alteração no balanço de energia e nitrogênio nos periodos pré e pós-parto. Não foi observado diferença na produção, composição e teor dos componentes totais do leite. No perfil de ácidos graxos do leite houve efeito (P<0,05) linear descrescente para as concentrações de C16:1cis, C18:1 cis, total de C:18 insaturado, total de AG monoinsaturados, insaturados e a relação do total de AGS:AGI. Foi observado efeito linear (P<0,05) crescente para o total de AG aturado e efeito (P<0,05) quadrático para C18:2, CLAcis9-trans11, e total de AGPI. Foi observado efeito linear crescente (P<0,05) para colesterol total, LDL no préparto e linear decrescente (P<0,05) para GGT nos períodos pré e pós-parto. Foi observado efeito quadrático (P<0,05) para HDL no pré-parto e AST no pós-parto. Em relação a atividade do sistema imune foi observado efeito linear (P<0,05) crescente para o percentual de CD3+ ativos no pós-parto, para o percentual de monócitos que produziram espécie reativa de oxigênio (ERO) no pós-parto quando foram estimulados por S.aureus e E.coli e para a intensidade de imunofluorescência de ERO para ganulócitos no pós-parto quando estimulados por S.aureus. Foi observado efeito (P<0,05) quadrático para o percentual de granulócitos, mononucleares, CD8+ ativos no pós-parto e para o percentual de granulócitos que produziram ERO no pós-parto quando estimulados por E.coli. A suplementação prolongada com GSI no pré-parto melhora a atividade do sistema imune, não melhora a qualidade oocitária e embrionária bem como não influencia negativamente os parametros produtivos de vacas leiteiras no período de transição e início de lactação

Confirming the Association Between Affective Disturbance and Perceived Cognitive Dysfunction Among Breast Cancer Patients: A Meta-Analysis of Available Literature

Purpose: Breast cancer is the most frequently diagnosed cancer among women worldwide. While undergoing chemotherapy treatment for breast cancer, patients often report experiencing "chemobrain." Previous literature reports correlations between psychological distress and these perceived cognitive problems. The aim of the present study was to examine the strength of the association between affective disturbance and subjective cognitive dysfunction.Methods: This study included a meta-analysis of the literature reporting a correlation between mood and subjective cognitive dysfunction. Eight studies with 1344 breast cancer patients treated with chemotherapy were selected based on stringent study inclusion criteria. Studies reporting a correlation coefficient between mood and subjective cognitive dysfunction were included.Results: In these data, there was no significant correlation between affective disturbance and subjective cognitive dysfunction. A random effects model yielded an overall weighted mean effect size of 0.12.Conclusion: Although this meta-analysis did not confirm the correlation between mood and subjective cognitive dysfunction, there was a clear association between these factors in the original disaggregated analyses, and they are clearly impactful from the time of diagnosis through long-term after care. The clinical implications of the present study and future directions for research are discussed.

Do agricultural subsidies crowd out or stimulate rural credit institutions? The Case of CAP Payments. Factor Markets Working Paper No. 4, September 2011

In this paper we estimate the impact of subsidies from the EU’s common agricultural policy on farm bank loans. According to the theoretical results, if subsidies are paid at the beginning of the growing season they may reduce bank loans, whereas if they are paid at the end of the season they increase bank loans, but these results are conditional on whether farms are credit constrained and on the relative cost of internal and external financing. In the empirical analysis, we use farm-level panel data from the Farm Accountancy Data Network to test the theoretical predictions for the period 1995–2007. We employ fixed-effects and generalised method of moment models to estimate the impact of subsidies on farm loans. The results suggest that subsidies influence farm loans and the effects tend to be non-linear and indirect. The results also indicate that both coupled and decoupled subsidies stimulate long-term loans, but the long-term loans of large farms increase more than those of small farms, owing to decoupled subsidies. Furthermore, the results imply that short-term loans are affected only by decoupled subsidies, and they are altered by decoupled subsidies more for small farms than for large farms; however, when controlling for endogeneity, only the decoupled payments affect loans and the relationship is non-linear.

Governança corporativa e custo de agência: uma visão do mercado brasileiro

Esse trabalho investiga empiricamente a relação entre custo de agência e as medidas de monitoramento interno disponíveis aos investidores brasileiros nas empresas nacionais, utilizando amostras de companhias abertas entre os anos de 2010 e 2014, totalizando 134 empresas analisadas e 536 observações. Para medir tal relação, foram utilizadas, como variáveis de monitoramento interno, informações sobre a remuneração variável dos executivos, entre elas o uso de outorgas de opções de compra de ações, a composição do conselho de administração, dando ênfase à representatividade dos conselheiros independentes e à dualidade entre Chairman e CEO, e o percentual do capital social das companhias que está sob propriedade dos executivos. Como proxy para custo de agência, foram utilizados os indicadores Asset Turnover Ratio e General & Administrative Expenses (G&A) como percentual da Receita Líquida. Neste contexto, foram estabelecidas duas hipóteses de pesquisa e estimados modelos de regressão em painel controlados por efeitos fixos de tempo e empresa, empregando como variável dependente as variáveis proxy do custo de agência e utilizando as variáveis endividamento e tamanho das empresas como variáveis de controle. Os resultados dos modelos demonstram que, na amostra selecionada, há uma relação positiva e significativa entre o percentual da remuneração variável e as proxies de custo de agência, comportamento este contrário ao esperado originalmente. Conclui-se assim que as empresas que apresentam uma maior composição variável no total remunerado ao executivo, incorrem em um maior custo de agência, o que leva à conclusão de que tais ferramentas não são boas estratégias de alinhamento de interesses entre executivos e acionistas. As demais variáveis de monitoramento interno não apresentaram significância.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Characteristics and benefits of professional work - Assessment of their importance over a 30-year career

The Professions in Australia Study is the first longitudinal investigation of the professions in Australia; it spans 33 years. Self-administered questionnaires were distributed on at least eight occasions between 1965 and 1998 to cohorts of students and later practitioners from the professions of engineering, law and medicine. The longitudinal design of this study has allowed for an investigation of individual change over time of three archetypal characteristics of the professions, service, knowledge and autonomy and two of the benefits of professional work, financial rewards and prestige. A cumulative logit random effects model was used to statistically assess changes in the ordinal response scores for measuring importance of the characteristics and benefits through stages of the career path. Individuals were also classified by average trends in response scores over time and hence professions are described through their members' tendency to follow a particular path in attitudes either of change or constancy, in relation to the importance of the five elements (characteristics and benefits). Comparisons in trends are also made between the three professions.

Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

Genome-wide linkage scan for loci influencing plasma triglycerides

Background: Plasma triglyceride concentration is known to be a significant risk factor for cardiovascular disease (CVD). Previous studies have found that the level of triglycerides is strongly influenced by genetic factors. Methods: To identify quantitative trait loci influencing triglycerides, we conducted a genome-wide linkage scan on data from 485 Australian adult dizygotic twin pairs. Prior to linkage analysis, triglyceride values were adjusted for the effects of covariates including age, sex, time since last meal, time of blood collection (CT) and time to plasma separation. Results: The heritability estimate for ln(triglyceride) adjusted for all above fixed effects was 0.49. The highest multipoint LOD score observed was 2.94 (genome-wide p=0.049) on chromosome 7 (at 65cM). This 7p region contains several candidate genes. Two other regions with suggestive multipoint LOD scores were also identified on chromosome 4 (LOD score=2.26 at 62cM) and chromosome X (LOD score=2.01 at 81cM). Conclusions: The linkage peaks found represent newly identified regions for more detailed study, in particular the significant linkage observed on chromosome 7p13. \ (c) 2006 Elsevier B.V. All rights reserved.

The importance of forest area and configuration relative to local habitat factors for conserving forest mammals: A case study of koalas in Queensland, Australia

The loss and fragmentation of forest habitats by human land use are recognised as important factors influencing the decline of forest-dependent fauna. Mammal species that are dependent upon forest habitats are particularly sensitive to habitat loss and fragmentation because they have highly specific habitat requirements, and in many cases have limited ability to move through and utilise the land use matrix. We addressed this problem using a case study of the koala (Phascolarctos cinereus) surveyed in a fragmented rural-urban landscape in southeast Queensland, Australia. We applied a logistic modelling and hierarchical partitioning analysis to determine the importance of forest area and its configuration relative to site (local) and patch-level habitat variables. After taking into account spatial auto-correlation and the year of survey, we found koala occurrence increased with the area of all forest habitats, habitat patch size and the proportion of primary Eucalyptus tree species; and decreased with mean nearest neighbour distance between forest patches, the density of forest patches, and the density of sealed roads. The difference between the effect of habitat area and configuration was not as strong as theory predicts, with the configuration of remnant forest becoming increasingly important as the area of forest habitat declines. We conclude that the area of forest, its configuration across the landscape, as well as the land use matrix, are important determinants of koala occurrence, and that habitat configuration should not be overlooked in the conservation of forest-dependent mammals, such as the koala. We highlight the implications of these findings for koala conservation. (c) 2006 Elsevier Ltd. All rights reserved.

Methods for Categorical Longitudinal Survey Data: Understanding Employment Status of Australian Women

Many variables that are of interest in social science research are nominal variables with two or more categories, such as employment status, occupation, political preference, or self-reported health status. With longitudinal survey data it is possible to analyse the transitions of individuals between different employment states or occupations (for example). In the statistical literature, models for analysing categorical dependent variables with repeated observations belong to the family of models known as generalized linear mixed models (GLMMs). The specific GLMM for a dependent variable with three or more categories is the multinomial logit random effects model. For these models, the marginal distribution of the response does not have a closed form solution and hence numerical integration must be used to obtain maximum likelihood estimates for the model parameters. Techniques for implementing the numerical integration are available but are computationally intensive requiring a large amount of computer processing time that increases with the number of clusters (or individuals) in the data and are not always readily accessible to the practitioner in standard software. For the purposes of analysing categorical response data from a longitudinal social survey, there is clearly a need to evaluate the existing procedures for estimating multinomial logit random effects model in terms of accuracy, efficiency and computing time. The computational time will have significant implications as to the preferred approach by researchers. In this paper we evaluate statistical software procedures that utilise adaptive Gaussian quadrature and MCMC methods, with specific application to modeling employment status of women using a GLMM, over three waves of the HILDA survey.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.

Macro determinants of FDI inflows to Japan:an analysis of source country characteristics

This paper examines the source country determinants of FDI into Japan. The paper highlights certain methodological and theoretical weaknesses in the previous literature and offers some explanations for hitherto ambiguous results. Specifically, the paper highlights the importance of panel data analysis, and the identification of fixed effects in the analysis rather than simply pooling the data. Indeed, we argue that many of the results reported elsewhere are a feature of this mis-specification. To this end, pooled, fixed effects and random effects estimates are compared. The results suggest that FDI into Japan is inversely related to trade flows, such that trade and FDI are substitutes. Moreover, the results also suggest that FDI increases with home country political and economic stability. The paper also shows that previously reported results, regarding the importance of exchange rates, relative borrowing costs and labour costs in explaining FDI flows, are sensitive to the econometric specification and estimation approach. The paper also discusses the importance of these results within a policy context. In recent years Japan has sought to attract FDI, though many firms still complain of barriers to inward investment penetration in Japan. The results show that cultural and geographic distance are only of marginal importance in explaining FDI, and that the results are consistent with the market-seeking explanation of FDI. As such, the attitude to risk in the source country is strongly related to the size of FDI flows to Japan. © 2007 The Authors Journal compilation © 2007 Blackwell Publishing Ltd.