Approximate and Exact Algorithms for Constrained (Un) Weighted Two-dimensional Two-staged Cutting Stock Problems

In this paper we propose two algorithms for solving both unweighted and weighted constrained two-dimensional two-staged cutting stock problems. The problem is called two-staged cutting problem because each produced (sub)optimal cutting pattern is realized by using two cut-phases. In the first cut-phase, the current stock rectangle is slit down its width (resp. length) into a set of vertical (resp. horizontal) strips and, in the second cut-phase, each of these strips is taken individually and chopped across its length (resp. width).First, we develop an approximate algorithm for the problem. The original problem is reduced to a series of single bounded knapsack problems and solved by applying a dynamic programming procedure. Second, we propose an exact algorithm tailored especially for the constrained two-staged cutting problem. The algorithm starts with an initial (feasible) lower bound computed by applying the proposed approximate algorithm. Then, by exploiting dynamic programming properties, we obtain good lower and upper bounds which lead to significant branching cuts. Extensive computational testing on problem instances from the literature shows the effectiveness of the proposed approximate and exact approaches.     

An Econometric Model of Network Formation With Degree Heterogeneity

I introduce a model of undirected dyadic link formation which allows for assortative matching on observed agent characteristics (homophily) as well as unrestricted agent‐level heterogeneity in link surplus (degree heterogeneity). Like in fixed effects panel data analyses, the joint distribution of observed and unobserved agent‐level characteristics is left unrestricted. Two estimators for the (common) homophily parameter, β₀, are developed and their properties studied under an asymptotic sequence involving a single network growing large. The first, tetrad logit (TL), estimator conditions on a sufficient statistic for the degree heterogeneity. The second, joint maximum likelihood (JML), estimator treats the degree heterogeneity {A_i0}_{i = 1}^N as additional (incidental) parameters to be estimated. The TL estimate is consistent under both sparse and dense graph sequences, whereas consistency of the JML estimate is shown only under dense graph sequences.     

The intergenerational continuity of parental conditional regard and its role in mothers’ and adolescents’ contingent self‐esteem and depressive symptoms

Although research has documented the adverse consequences of parental academic conditional regard in different developmental periods, few studies have examined antecedents and, in particular, the possibility of intergenerational continuity of this parenting dimension. The current study aimed to identify patterns of intergenerational similarity in two types of conditional regard (i.e., positive and negative). Additionally, it examined mothers’ and adolescents’ contingent self‐esteem (CSE) and depressive symptoms as outcomes of this process. In total, 211 mothers and their 10‐ to 16‐year‐old adolescents filled out questionnaires assessing perceived conditional regard in the relationship with their own mother (i.e., grandmothers’, and mothers’ conditional positive (PACPR) and negative regard (PACNR)), contingent self‐esteem (i.e., maternal child‐invested CSE and adolescents’ academic CSE), and maternal and adolescents’ depressive symptoms. Results revealed direct intergenerational similarity for PACPR, but not for PACNR. Within generations, PACPR was related positively to CSE, which, in turn, was related positively to depressive symptoms. Moreover, PACNR was related positively to depressive symptoms. Our findings provide preliminary evidence for the intergenerational continuity of conditional regard, which has important consequences because of its relation to CSE and subsequent depressive symptoms within both generations.     

Generalised quasi‐symmetry models for ordinal contingency tables

This paper extends the ordinary quasi‐symmetry (QS) model for square contingency tables with commensurable classification variables. The proposed generalised QS model is defined in terms of odds ratios that apply to ordinal variables. In particular, we present QS models based on global, cumulative and continuation odds ratios and discuss their properties. Finally, the conditional generalised QS model is introduced for local and global odds ratios. These models are illustrated through the analysis of two data sets.     

A Potts‐mixture spatiotemporal joint model for combined magnetoencephalography and electroencephalography data

We develop a new methodology for determining the location and dynamics of brain activity from combined magnetoencephalography (MEG) and electroencephalography (EEG) data. The resulting inverse problem is ill‐posed and is one of the most difficult problems in neuroimaging data analysis. In our development we propose a solution that combines the data from three different modalities, magnetic resonance imaging (MRI), MEG and EEG, together. We propose a new Bayesian spatial finite mixture model that builds on the mesostate‐space model developed by Daunizeau & Friston [Daunizeau and Friston, NeuroImage 2007; 38, 67–81]. Our new model incorporates two major extensions: (i) We combine EEG and MEG data together and formulate a joint model for dealing with the two modalities simultaneously; (ii) we incorporate the Potts model to represent the spatial dependence in an allocation process that partitions the cortical surface into a small number of latent states termed mesostates. The cortical surface is obtained from MRI. We formulate the new spatiotemporal model and derive an efficient procedure for simultaneous point estimation and model selection based on the iterated conditional modes algorithm combined with local polynomial smoothing. The proposed method results in a novel estimator for the number of mixture components and is able to select active brain regions, which correspond to active variables in a high‐dimensional dynamic linear model. The methodology is investigated using synthetic data and simulation studies and then demonstrated on an application examining the neural response to the perception of scrambled faces. R software implementing the methodology along with several sample datasets are available at the following GitHub repository https://github.com/v2south/PottsMix . The Canadian Journal of Statistics 47: 688–711; 2019 © 2019 Statistical Society of Canada     

基于时间序列模型的海洋溶解氧分析与预测

季节效应是指某序列由于受自然气候等因素的影响,随季节的变化而呈现出周期性的变化规律。季节效应一般以月或季为单位。文中研究的海洋溶解氧含量因为受每日气温周期性变化的影响,以2 h为单位也呈周期性的变化规律,类似于季节效应。采用2 h为采样间隔,尝试模拟乘积季节模型来对溶解氧含量进行分析。利用条件期望预测的差分方程形式对序列进行短期预测,并在实践中比对预测数据和实际数据,证明取得了不错的预测精度。     

The effect of non-normal disturbances and conditional heteroskedasticity on multiple cointegration tests

!t is well-known that Johansen's multiple cointegration tests' results and those of Johansen and Juselius' tests for restricrions on cointegrating vectors and their weights have far-reaching implications for economic modelling and analysis. Therefore, it is important to ensure that the tests have desirable finite sample properties. Although the statistics are derived under Gaussian distribution,the asympotic results are derived under a much wider class of distributions. Using simulation, this paper investigates the effect of non-normal disturbances on these tests in finite samples. Further, ARCH/GARCH type conditional heteroskedasticity is present in many economic and financial time series. This paper examines the finite properties of the tests when the error term follows ARCH/GARCH type processes. From the evidence, it appears that researchers should not be overly concerned by the possibility of small departures from non-normality when using Johansen's suggested techniques even in finite samples. ARCH and GARCH effects may be more problematic, however. In particular it becomes more important ro test whether the restriction implicit in the integrated (or near-integrated) ARCH-type Drocess actually holds in time series for the application of the cointegraiion rank tests and the test for restrictions on cointegrating weights. The tests for restrictions on cointegrating vectors apper to be robust for non-normal errors and for all ARCH and GARCH type processes considered.     

Double threshold autoregressive conditionally heteroscedastic model building by genetic algorithms

Asymmetric behaviour in both mean and variance is often observed in real time series. The approach we adopt is based on double threshold autoregressive conditionally heteroscedastic (DTARCH) model with normal innovations. This model allows threshold nonlinearity in mean and volatility to be modelled as a result of the impact of lagged changes in assets and squared shocks, respectively. A methodology for building DTARCH models is proposed based on genetic algorithms (GAs). The most important structural parameters, that is regimes and thresholds, are searched for by GAs, while the remaining structural parameters, that is the delay parameters and models orders, vary in some pre-specified intervals and are determined using exhaustive search and an Asymptotic Information Criterion (AIC) like criterion. For each structural parameters trial set, a DTARCH model is fitted that maximizes the (penalized) likelihood (AIC criterion). For this purpose the iteratively weighted least squares algorithm is used. Then the best model according to the AIC criterion is chosen. Extension to the double threshold generalized ARCH (DTGARCH) model is also considered. The proposed methodology is checked using both simulated and market index data. Our findings show that our GAs-based procedure yields results that comparable to that reported in the literature and concerned with real time series. As far as artificial time series are considered, the proposed procedure seems to be able to fit the data quite well. In particular, a comparison is performed between the present procedure and the method proposed by Tsay [Tsay, R.S., 1989, Testing and modeling threshold autoregressive processes. Journal of the American Statistical Association, Theory and Methods, 84, 231–240.] for estimating the delay parameter. The former almost always yields better results than the latter. However, adopting Tsay's procedure as a preliminary stage for finding the appropriate delay parameter may save computational time specially if the delay parameter may vary in a large interval.     

Bayesian analysis of spatial data using different variance and neighbourhood structures

In disease mapping, the overall goal is to study the incidence or mortality risk caused by a specific disease in a number of geographical regions. It is common to assume that the response variable follows a Poisson distribution, whose average rate can be explained by a group of covariates and a random effect. For this random effect, it is considered conditional autoregressive (CAR) models, which carry information about the neighbourhood relationship between the regions. The focus of this paper was to explore and compare some CAR models proposed in the literature. An application with epidemiological data was conducted to model the risk of death due to Crohn's Disease and Ulcerative Colitis in the State of São Paulo – Brazil. Finally, a simulation study was done to strengthen the results and assess the performance of the models in the presence of various levels of spatial dependence.