A finite-population revenue management model and a risk-ratio procedure for the joint estimation of population size and parameters

Many dynamic revenue management models divide the sale period into a finite number of periods T and assume, invoking a fine-enough grid of time, that each period sees at most one booking request. These Poisson-type assumptions restrict the variability of the demand in the model, but researchers and practitioners were willing to overlook this for the benefit of tractability of the models. In this paper, we criticize this model from another angle. Estimating the discrete finite-period model poses problems of indeterminacy and non-robustness: Arbitrarily fixing T leads to arbitrary control values and on the other hand estimating T from data adds an additional layer of indeterminacy. To counter this, we first propose an alternate finite-population model that avoids this problem of fixing T and allows a wider range of demand distributions, while retaining the useful marginal-value properties of the finite-period model. The finite-population model still requires jointly estimating market size and the parameters of the customer purchase model without observing no-purchases. Estimation of market-size when no-purchases are unobservable has rarely been attempted in the marketing or revenue management literature. Indeed, we point out that it is akin to the classical statistical problem of estimating the parameters of a binomial distribution with unknown population size and success probability, and hence likely to be challenging. However, when the purchase probabilities are given by a functional form such as a multinomial-logit model, we propose an estimation heuristic that exploits the specification of the functional form, the variety of the offer sets in a typical RM setting, and qualitative knowledge of arrival rates. Finally we perform simulations to show that the estimator is very promising in obtaining unbiased estimates of population size and the model parameters.

Familiality and SNP heritability of age at onset and episodicity in major depressive disorder.

BACKGROUND: Strategies to dissect phenotypic and genetic heterogeneity of major depressive disorder (MDD) have mainly relied on subphenotypes, such as age at onset (AAO) and recurrence/episodicity. Yet, evidence on whether these subphenotypes are familial or heritable is scarce. The aims of this study are to investigate the familiality of AAO and episode frequency in MDD and to assess the proportion of their variance explained by common single nucleotide polymorphisms (SNP heritability). METHOD: For investigating familiality, we used 691 families with 2-5 full siblings with recurrent MDD from the DeNt study. We fitted (square root) AAO and episode count in a linear and a negative binomial mixed model, respectively, with family as random effect and adjusting for sex, age and center. The strength of familiality was assessed with intraclass correlation coefficients (ICC). For estimating SNP heritabilities, we used 3468 unrelated MDD cases from the RADIANT and GSK Munich studies. After similarly adjusting for covariates, derived residuals were used with the GREML method in GCTA (genome-wide complex trait analysis) software. RESULTS: Significant familial clustering was found for both AAO (ICC = 0.28) and episodicity (ICC = 0.07). We calculated from respective ICC estimates the maximal additive heritability of AAO (0.56) and episodicity (0.15). SNP heritability of AAO was 0.17 (p = 0.04); analysis was underpowered for calculating SNP heritability of episodicity. CONCLUSIONS: AAO and episodicity aggregate in families to a moderate and small degree, respectively. AAO is under stronger additive genetic control than episodicity. Larger samples are needed to calculate the SNP heritability of episodicity. The described statistical framework could be useful in future analyses.

Mossbauer studies of amorphous FeSi compositionally modulated thin films

Conversion electron Mossbauer spectra of composition modulated FeSi thin films have been analysed within the framework of a quasi shape independent model in which the distribution function for the hyperfine fields is assumed to be given by a binomial distribution. Both the hyperfine field and the hyperfine field distribution depend on the modulation characteristic length.

A Lei de Lotka na bibliometria brasileira

Usando os dados reportados em artigos publicados em revistas brasileiras e trabalhos apresentados em congressos nacionais, replicaram-se as aplicações da Lei de Lotka à literatura brasileira em 10 campos diferentes. Utilizou-se o modelo do poder inverso pelos métodos do mínimo quadrado e probabilidade máxima. Das 10 literaturas nacionais analisadas, somente a literatura de medicina, siderurgia, jaca e biblioteconomia ajustaram-se ao modelo do poder inverso generalizado pelo método dos mínimos quadrados. No entanto, só duas literaturas (veterinária e cartas do Arquivo Privado de Getúlio Vargas) não se ajustaram ao modelo quando se usou o método da máxima probabilidade. Para ambas literaturas, tentaram-se modelos diferentes. A literatura de veterinária ajustou-se à distribuição binomial negativa, e as cartas do Arquivo Privado de Getúlio Vargas ajustaram-se melhor à distribuição Gauss-Poisson Inversa Generalizada.

Migration Flows and Quality of Life in a Metropolitan Area: the Case of Barcelona-Spain

The phenomenon of human migration is certainly not new and it has been studied from a variety of perspectives. Yet, the attention on human migration and its determinant has not been fading over time as confirmed by recent contributions (see for instance Cushing and Poot 2004 and Rebhun and Raveh 2006). In this paper we combine the recent theoretical contributions by Douglas (1997) and Wall (2001) with the methodological advancements of Guimarães et al. (2000, 2003) to model inter-municipal migration flows in the Barcelona area. In order to do that, we employ two different types of count models, i.e. the Poisson and negative binomial and compare the estimations obtained. Our results show that, even after controlling for the traditional migration factors, QoL (measured with a Composite Index which includes numerous aspects and also using a list of individual variables) is an important determinant of short distance migration movements in the Barcelona area.

Which change talk dimensions are associated with subsequent drinking among young men receiving a brief motivational intervention in a non-clinical setting?

Background: Patient change talk (CT) during brief motivational interventions (BMI) has been linked with subsequent changes in drinking in clinical settings but this link has not been clearly established among young people in non-clinical populations. Objective: To determine which of several CT dimensions assessed during an effective BMI delivered in a non-clinical setting to 20-year old men are associated with drinking 6 months later. Methods: Of 125 individuals receiving a face-to-face BMI session (15.8 ± 5.4 minutes), we recorded and coded a subsample of 42 sessions using the Motivational Interviewing Skill Code 2.1. Each patient change talk utterance was categorized as `Reason´, `Ability´, `Desire´, `Need´, `Commitment´, `Taking steps´, or `Other´. Each utterance was graded according to its strength (absolute value from 1 to 3) and direction (i.e. towards (positive sign) or away (negative sign) from change/in favor of status quo). `Ability´, `Desire´, and `Need´ to change (`ADN´) were grouped together since these codes were too scarce to conduct analyses. Mean strength scores over the entire session were computed for each dimension and later dichotomized in towards change (i.e. mean core > 0) and away from change/in favor of status quo. Negative binomial regression models were used to assess the relationship between CT dimensions and drinking 6 months later, adjusting for drinking at baseline. Results: Compared to subjects with a `Taking steps´ score away from change/in favor of status quo, subjects with a positive `Taking steps´ score reported significantly less drinking 6 months later (Incidence Rate Ration [IRR] for drinks per week: 0.56, 95% Confidence Interval [CI] 0.31, 1.00). IRR (95%CI) for subjects with a positive `ADN´ score was 0.58, (0.32, 1.03). For subjects with a positive `Reason´, `Commitment´, and `Other´ scores, IRR (95%CI) were 1.28 (0.77; 2.12) 1.63 (0.85; 3.14) and 1.03 (0.61; 1.72), respectively. Conclusion: A change talk dimension reflecting steps taken towards change (`Taking steps´) is associated with less drinking 6 months later among young men receiving a BMI in a non-clinical setting. Encouraging patients to take steps towa change may be a worthy objective for clinicians and may explain BMI efficacy.

Efficacy Of Canakinumab (Acz885), A Fully Human Anti-Interleukin (Il)-1beta Monoclonal Antibody, In The Prevention Of Flares In Gout Patients Initiating Allopurinol Therapy

Background: Gout patients initiating urate lowering therapy have an increased risk of flares. Inflammation in gouty arthritis is induced by IL-1b. Canakinumab targets and inhibits IL-1b effectively in clinical studies. This study compared different doses of canakinumab vs colchicine in preventing flares in gout patients initiating allopurinol therapy.Methods: In this 24 week double blind study, gout patients (20-79 years) initiating allopurinol were randomized (1:1:1:1:1:1:2) to canakinumab s.c. single doses of 25, 50, 100, 200, 300 mg, or 150 mg divided in doses every 4 weeks (50+50+25+25 mg [q4wk]) or colchicine 0.5 mg p.o. daily for 16 weeks. Primary outcome was to determine the canakinumab dose giving comparable efficacy to colchicine with respect to the number of gout flares occurring during first 16 weeks. Secondary outcomes included number of patients with gout flares and C-reactive protein (CRP) levels during the first 16 weeks.Results: 432 patients were randomized and 391 (91%) completed the study. All canakinumab doses were better than colchicine in preventing flares and therefore, a canakinumab dose comparable to colchicine could not be determined. Based on a negative binomial model, all canakinumab groups, except 25 mg, reduced the flare rate ratio per patient significantly compared to colchicine group (rate ratio estimates 25 mg 0.60, 50 mg 0.34, 100 mg 0.28, 200 mg 0.37, 300 mg 0.29, q4wk 0.38; p<=0.05). The percentage of patients with flares was lower for all canakinumab groups (25 mg 27.3%, 50 mg 16.7%, 100 mg 14.8%, 200 mg 18.5%, 300 mg 15.1%, q4wk 16.7%) compared to colchicine group (44.4%). All patients taking canakinumab were significantly less likely to experience at least one gout flare than patients taking colchicine (odds ratio range [0.22 - 0.47]; p<=0.05 for all). The median baseline CRP levels were 2.86 mg/L for 25 mg, 3.42 mg/L for 50 mg, 1.76 mg/L for 100 mg, 3.66 mg/L for 200 mg, 3.21 mg/L for 300 mg, 3.23 mg/L for q4wk canakinumab groups and 2.69 mg/L for colchicine group. In all canakinumab groups with median CRP levels above the normal range at baseline, median levels declined within 15 days of treatment and were maintained at normal levels (ULN=3 mg/L) throughout the 16 week period. Adverse events (AEs) occurred in 52.7% (25 mg), 55.6% (50 mg), 51.9% (100 mg), 51.9% (200 mg), 54.7% (300 mg), and 58.5% (q4wk) of patients on canakinumab vs 53.7% of patients on colchicine. Serious AEs (SAE) were reported in 2 (3.6%; 25 mg), 2 (3.7%, 50 mg), 3 (5.6%, 100 mg), 3 (5.6%, 200 mg), 3 (5.7%, 300 mg) and 1 (1.9%, q4wk) patients on canakinumab and in 5 (4.6%) patients on colchicine. One fatal SAE (myocardial infarction, not related to study drug) occurred in colchicine group.Conclusion: In this large randomized, double-blind active controlled study of flare prevention in gout patients initiating allopurinol therapy, treatment with canakinumab led to a statistically significant reduction in flares compared with colchicine (standard of care), and was well tolerated.

Are niche-based species distribution models transferable in space?

Aim To assess the geographical transferability of niche-based species distribution models fitted with two modelling techniques. Location Two distinct geographical study areas in Switzerland and Austria, in the subalpine and alpine belts. Methods Generalized linear and generalized additive models (GLM and GAM) with a binomial probability distribution and a logit link were fitted for 54 plant species, based on topoclimatic predictor variables. These models were then evaluated quantitatively and used for spatially explicit predictions within (internal evaluation and prediction) and between (external evaluation and prediction) the two regions. Comparisons of evaluations and spatial predictions between regions and models were conducted in order to test if species and methods meet the criteria of full transferability. By full transferability, we mean that: (1) the internal evaluation of models fitted in region A and B must be similar; (2) a model fitted in region A must at least retain a comparable external evaluation when projected into region B, and vice-versa; and (3) internal and external spatial predictions have to match within both regions. Results The measures of model fit are, on average, 24% higher for GAMs than for GLMs in both regions. However, the differences between internal and external evaluations (AUC coefficient) are also higher for GAMs than for GLMs (a difference of 30% for models fitted in Switzerland and 54% for models fitted in Austria). Transferability, as measured with the AUC evaluation, fails for 68% of the species in Switzerland and 55% in Austria for GLMs (respectively for 67% and 53% of the species for GAMs). For both GAMs and GLMs, the agreement between internal and external predictions is rather weak on average (Kulczynski's coefficient in the range 0.3-0.4), but varies widely among individual species. The dominant pattern is an asymmetrical transferability between the two study regions (a mean decrease of 20% for the AUC coefficient when the models are transferred from Switzerland and 13% when they are transferred from Austria). Main conclusions The large inter-specific variability observed among the 54 study species underlines the need to consider more than a few species to test properly the transferability of species distribution models. The pronounced asymmetry in transferability between the two study regions may be due to peculiarities of these regions, such as differences in the ranges of environmental predictors or the varied impact of land-use history, or to species-specific reasons like differential phenotypic plasticity, existence of ecotypes or varied dependence on biotic interactions that are not properly incorporated into niche-based models. The lower variation between internal and external evaluation of GLMs compared to GAMs further suggests that overfitting may reduce transferability. Overall, a limited geographical transferability calls for caution when projecting niche-based models for assessing the fate of species in future environments.

Canakinumab reduces the risk of acute gouty arthritis flares during initiation of allopurinol treatment: results of a double-blind, randomised study.

Objective This study assessed the efficacy and safety of canakinumab, a fully human anti-interleukin 1 beta monoclonal antibody, for prophylaxis against acute gouty arthritis flares in patients initiating urate-lowering treatment.Methods In this double-blind, double-dummy, dose-ranging study, 432 patients with gouty arthritis initiating allopurinol treatment were randomised 1:1:1:1:1:1:2 to receive: a single dose of canakinumab, 25, 50, 100, 200, or 300 mg subcutaneously; 4 x 4-weekly doses of canakinumab (50 + 50 + 25 + 25 mg subcutaneously); or daily colchicine 0.5 mg orally for 16 weeks. Patients recorded details of flares in diaries. The study aimed to determine the canakinumab dose having equivalent efficacy to colchicine 0.5 mg at 16 weeks.Results A dose-response for canakinumab was not apparent with any of the four predefined dose-response models. The estimated canakinumab dose with equivalent efficacy to colchicine was below the range of doses tested. At 16 weeks, there was a 62% to 72% reduction in the mean number of flares per patient for canakinumab doses >= 50 mg versus colchicine based on a negative binomial model (rate ratio: 0.28-0.38, p <= 0.0083), and the percentage of patients experiencing >= 1 flare was significantly lower for all canakinumab doses (15% to 27%) versus colchicine (44%, p<0.05). There was a 64% to 72% reduction in the risk of experiencing >= 1 flare for canakinumab doses >= 50 mg versus colchicine at 16 weeks (hazard ratio (HR): 0.28-0.36, p <= 0.05). The incidence of adverse events was similar across treatment groups.Conclusions Single canakinumab doses >= 50 mg or four 4-weekly doses provided superior prophylaxis against flares compared with daily colchicine 0.5 mg.

Efficacy of Canakinumab (ACZ885), a Fully Human Anti-Interleukin (IL)-1beta Monoclonal Antibody, in the Prevention of Flares in Gout Patients Initiating Allopurinol Therapy

Background: Gout patients initiating urate lowering therapy have an increased risk of flares. Inflammation in gouty arthritis is induced by IL-1b. Canakinumab targets and inhibits IL-1b effectively in clinical studies. This study compared different doses of canakinumab vs colchicine in preventing flares in gout patients initiating allopurinol therapy.Methods: In this 24 week double blind study, gout patients (20-79 years) initiating allopurinol were randomized (1:1:1:1:1:1:2) to canakinumab s.c. single doses of 25, 50, 100, 200, 300 mg, or 150 mg divided in doses every 4 weeks (50+50+25+25 mg [q4wk]) or colchicine 0.5 mg p.o. daily for 16 weeks. Primary outcome was to determine the canakinumab dose giving comparable efficacy to colchicine with respect to the number of gout flares occurring during first 16 weeks. Secondary outcomes included number of patients with gout flares and C-reactive protein (CRP) levels during the first 16 weeks.Results: 432 patients were randomized and 391 (91%) completed the study. All canakinumab doses were better than colchicine in preventing flares and therefore, a canakinumab dose comparable to colchicine could not be determined. Based on a negative binomial model, all canakinumab groups, except 25 mg, reduced the flare rate ratio per patient significantly compared to colchicine group (rate ratio estimates 25 mg 0.60, 50 mg 0.34, 100 mg 0.28, 200 mg 0.37, 300 mg 0.29, q4wk 0.38; p<=0.05). The percentage of patients with flares was lower for all canakinumab groups (25 mg 27.3%, 50 mg 16.7%, 100 mg 14.8%, 200 mg 18.5%, 300 mg 15.1%, q4wk 16.7%) compared to colchicine group (44.4%). All patients taking canakinumab were significantly less likely to experience at least one gout flare than patients taking colchicine (odds ratio range [0.22 - 0.47]; p<=0.05 for all). The median baseline CRP levels were 2.86 mg/L for 25 mg, 3.42 mg/L for 50 mg, 1.76 mg/L for 100 mg, 3.66 mg/L for 200 mg, 3.21 mg/L for 300 mg, 3.23 mg/L for q4wk canakinumab groups and 2.69 mg/L for colchicine group. In all canakinumab groups with median CRP levels above the normal range at baseline, median levels declined within 15 days of treatment and were maintained at normal levels (ULN=3 mg/L) throughout the 16 week period. Adverse events (AEs) occurred in 52.7% (25 mg), 55.6% (50 mg), 51.9% (100 mg), 51.9% (200 mg), 54.7% (300 mg), and 58.5% (q4wk) of patients on canakinumab vs 53.7% of patients on colchicine. Serious AEs (SAE) were reported in 2 (3.6%; 25 mg), 2 (3.7%, 50 mg), 3 (5.6%, 100 mg), 3 (5.6%, 200 mg), 3 (5.7%, 300 mg) and 1 (1.9%, q4wk) patients on canakinumab and in 5 (4.6%) patients on colchicine. One fatal SAE (myocardial infarction, not related to study drug) occurred in colchicine group.Conclusion: In this large randomized, double-blind active controlled study of flare prevention in gout patients initiating allopurinol therapy, treatment with canakinumab led to a statistically significant reduction in flares compared with colchicine (standard of care), and was well tolerated.

Social factors related to the clinical severity of influenza cases in Spain during the A(H1N1)2009 virus pandemic

Background During the 2009 influenza pandemic, a change in the type of patients most often affected by influenza was observed. The objective of this study was to assess the role of individual and social determinants in hospitalizations due to influenza A (H1N1) 2009 infection. Methods We studied hospitalized patients (cases) and outpatients (controls) with confirmed influenza A (H1N1) 2009 infection. A standardized questionnaire was used to collect data. Variables that might be related to the hospitalization of influenza cases were compared by estimation of the odds ratio (OR) and 95% confidence intervals (CI) and the variables entered into binomial logistic regression models. Results Hospitalization due to pandemic A (H1N1) 2009 influenza virus infections was associated with non-Caucasian ethnicity (OR: 2.18, 95% CI 1.17 − 4.08), overcrowding (OR: 2.84, 95% CI 1.20 − 6.72), comorbidity and the lack of previous preventive information (OR: 2.69, 95% CI: 1.50 − 4.83). Secondary or higher education was associated with a lower risk of hospitalization (OR 0.56, 95% CI: 0.36 − 0.87) Conclusions In addition to individual factors such as comorbidity, other factors such as educational level, ethnicity or overcrowding were associated with hospitalization due to A (H1N1) 2009 influenza virus infections.

Factors associated with 56-day non-return rate in dairy cattle

The objective of this work was to identify factors associated with the 56-day non-return rate (56-NRR) in dairy herds in the Galician region, Spain, and to estimate it for individual Holstein bulls. The experiment was carried out in herds originated from North-West Spain, from September 2008 to August 2009. Data of the 76,440 first inseminations performed during this period were gathered. Candidate factors were tested for their association with the 56-NRR by using a logistic model (binomial). Afterwards, 37 sires with a minimum of 150 first performed inseminations were individually evaluated. Logistic models were also estimated for each bull, and predicted individual 56-NRR rate values were calculated as a solution for the model parameters. Logistic regression found four major factors associated with 56-NRR in lactating cows: age at insemination, days from calving to insemination, milk production level at the time of insemination, and herd size. First-service conception rate, when a particular sire was used, was higher for heifers (0.71) than for lactating cows (0.52). Non-return rates were highly variable among bulls. Asignificant part of the herd-level variation of 56-NRR of Holstein cattle seems attributable to the service sire. High correlation level between observed and predicted 56-NRR was found.

Predictors for hospitalization and outpatient visits in patients with inflammatory bowel disease: results from the Swiss Inflammatory Bowel Disease Cohort Study.

OBJECTIVES: Patients with inflammatory bowel disease (IBD) have a high resource consumption, with considerable costs for the healthcare system. In a system with sparse resources, treatment is influenced not only by clinical judgement but also by resource consumption. We aimed to determine the resource consumption of IBD patients and to identify its significant predictors. MATERIALS AND METHODS: Data from the prospective Swiss Inflammatory Bowel Disease Cohort Study were analysed for the resource consumption endpoints hospitalization and outpatient consultations at enrolment [1187 patients; 41.1% ulcerative colitis (UC), 58.9% Crohn's disease (CD)] and at 1-year follow-up (794 patients). Predictors of interest were chosen through an expert panel and a review of the relevant literature. Logistic regressions were used for binary endpoints, and negative binomial regressions and zero-inflated Poisson regressions were used for count data. RESULTS: For CD, fistula, use of biologics and disease activity were significant predictors for hospitalization days (all P-values <0.001); age, sex, steroid therapy and biologics were significant predictors for the number of outpatient visits (P=0.0368, 0.023, 0.0002, 0.0003, respectively). For UC, biologics, C-reactive protein, smoke quitters, age and sex were significantly predictive for hospitalization days (P=0.0167, 0.0003, 0.0003, 0.0076 and 0.0175 respectively); disease activity and immunosuppressive therapy predicted the number of outpatient visits (P=0.0009 and 0.0017, respectively). The results of multivariate regressions are shown in detail. CONCLUSION: Several highly significant clinical predictors for resource consumption in IBD were identified that might be considered in medical decision-making. In terms of resource consumption and its predictors, CD and UC show a different behaviour.

Propuesta de muestreo para Neohydatothrips signifer (Thysanoptera: Thripidae) en el cultivo de maracuyá

El objetivo de este trabajo fue desarrollar una propuesta de muestreo de Neohydatothrips signifer en el cultivo de maracuyá. Se correlacionó la infestación promedio del trips en el cultivo con las variables meteorológicas de humedad relativa, temperatura y precipitación, con uso del coeficiente de correlación de Pearson. Así mismo, se determinó el índice de agregación de la población en el cultivo, con la prueba chi cuadrado (X²) para la distribución de Poisson, la prueba de la razón V²/µ, y la distribución de la binomial negativa. La única variable climática que se correlacionó significativamente con la densidad promedio del trips fue la temperatura. Se presentó mayor incidencia y densidad poblacional del insecto en los terminales vegetativos en comparación con los botones florales. El patrón de disposición espacial que se presentó en todos los estados de desarrollo fenológico del cultivo fue de tipo agregado. Con el índice de agregación y el tamaño promedio de la población, se calculó que se deben muestrear 22 terminales nuevos por hectárea, sin importar su posición en la planta. El muestreo del trips en el cultivo de maracuyá debe hacerse sobre los terminales vegetativos de la planta, sin importar la posición del terminal en ella.

Spatial predictions of land-use transitions and associated threats to biodiversity: the case of forest regrowth in mountain grasslands

Question Can we predict where forest regrowth caused by abandonment of agricultural activities is likely to occur? Can we assess how it may conflict with grassland diversity hotspots? Location Western Swiss Alps (4003210m a.s.l.). Methods We used statistical models to predict the location of land abandonment by farmers that is followed by forest regrowth in semi-natural grasslands of the Western Swiss Alps. Six modelling methods (GAM, GBM, GLM, RF, MDA, MARS) allowing binomial distribution were tested on two successive transitions occurring between three time periods. Models were calibrated using data on land-use change occurring between 1979 and 1992 as response, and environmental, accessibility and socio-economic variables as predictors, and these were validated for their capacity to predict the changes observed from 1992 to 2004. Projected probabilities of land-use change from an ensemble forecast of the six models were combined with a model of plant species richness based on a field inventory, allowing identification of critical grassland areas for the preservation of biodiversity. Results Models calibrated over the first land-use transition period predicted the second transition with reasonable accuracy. Forest regrowth occurs where cultivation costs are high and yield potential is low, i.e. on steeper slopes and at higher elevations. Overlaying species richness with land-use change predictions, we identified priority areas for the management and conservation of biodiversity at intermediate elevations. Conclusions Combining land-use change and biodiversity projections, we propose applied management measures for targeted/identified locations to limit the loss of biodiversity that could otherwise occur through loss of open habitats. The same approach could be applied to other types of land-use changes occurring in other ecosystems.