Repeatability across platforms

2 replies [Last post]
Steve's picture
Offline
Joined: 07/30/2009

I have run into an interesting problem. A model (which is admittedly of low stability) that I used in class today ran on the Mac (with warnings, but reasonably close estimates) but not on the PC (non-invertible). Both are Intel processors, both running R 2.9.2 and OpenMx 0.2.2. The model is of simulated data with boundary constraints on residuals because variance explained is very high.

I attach the script and data as a test case for the limits of similarities between architectures.

AttachmentSize
LDEUnivariateExample091116.R2.89 KB
LDEUnivariateExample091116.txt2.88 KB
carey's picture
Offline
Joined: 10/19/2009
from my (humble)

from my (humble) understanding of R, cross platform uniformity is not guaranteed. there can be subtle differences in C, fortran, and other compilers that can generate differences so that the compiled modules do not always give exact answers across platforms.
numerical optimization with constraints is one of those situations where slight differences in compiler settings can generate discrepant results within a platform, let alone across platforms. (most numerical guru's that i've talked to always recommend unconstrained optimization for likelihood problems). suggest trying a model without constraints and then comparing the results across platforms.
if you cannot get a model without constraints to converge, then there is probably something the matter with the model, given the data.
greg

pdeboeck's picture
Offline
Joined: 08/04/2009
This code requires

This code requires "GLLAfunctions.R" to run.