A new algorithm to obtain the least-squares or MINRES solution in common factor analysis is presented. It is based on the up-and-down Marquardt algorithm developed by the present authors for a general nonlinear least-squares problem. Experiments with some numerical models and some empirical data sets showed that the algorithm worked nicely and that SMC (Squared Multiple Correlation) performed best among four sets of initial values for common variances but that the solution might sometimes be very sensitive to fluctuations in the sample covariance matrix.

We examine the least squares approximation C to a symmetric matrix B, when all diagonal elements get weight w relative to all nondiagonal elements. WhenB has positivity p and C is constrained to be positive semi-definite, our main result states that, when w ≥1/2, then the rank of C is never greater than p, and when w ≤1/2 then the rank of C is at least p. For the problem of approximating a given n × n matrix with a zero diagonal by a squared-distance matrix, it is shown that the sstress criterion leads to a similar weighted least squares solution with w =(n+2)/4; the main result remains true. Other related problems and algorithmic consequences are briefly discussed.

At the ANZIAM conference in Hobart in February 2018, there were several talks on the solution of Laplace problems in multiply connected domains by means of conformal mapping. It appears to be not widely known that such problems can also be solved by the elementary method of series expansions with coefficients determined by least-squares fitting on the boundary. (These are not convergent series; the coefficients depend on the degree of the approximation.) Here we give a tutorial introduction to this method, which converges at an exponential rate if the boundary data are sufficiently well-behaved. The mathematical foundations go back to Runge in 1885 and Walsh in 1929. One of our examples involves an approximate Cantor set with up to 2048 components.

Collocation has become a standard tool for approximation of parameterized systems in the uncertainty quantification (UQ) community. Techniques for least-squares regularization, compressive sampling recovery, and interpolatory reconstruction are becoming standard tools used in a variety of applications. Selection of a collocation mesh is frequently a challenge, but methods that construct geometrically unstructured collocation meshes have shown great potential due to attractive theoretical properties and direct, simple generation and implementation. We investigate properties of these meshes, presenting stability and accuracy results that can be used as guides for generating stochastic collocation grids in multiple dimensions.

We develop a semi-parametric predictor of the IBNR reserve in a macro-model when the claim amount for a certain accident and development year can be expressed in a loglinear form composed of a deterministic part and a random error. We need to make assumptions only on the first two moments of the error, without any specified parametric assumption on its distribution. We give its properties, present its advantages and compare the estimates obtained with various predictors of the IBNR reserve, parametric and non-parametric, using a data set.