Re: [R] non-reproducible eigen() output with MKL

2025-05-04 Thread Eric Berger
t did not help. Has anyone seen this > behavior before? > > > > (When n is 5, the identical() call almost always returns "TRUE". As n > increases, the proportion of FALSE results increases, and it is nearly > always FALSE when n is 50.) > > > > Jesse &

Re: [R] non-reproducible eigen() output with MKL

2025-05-04 Thread smallepsilon via R-help
is > > behavior before? > > > > (When n is 5, the identical() call almost always returns "TRUE". As n > > increases, the proportion of FALSE results increases, and it is nearly > > always FALSE when n is 50.) > > > > Jesse > > > >

[R] Estimating regression with constraints in model coefficients - Follow-up on Constrained Ordinal Model — Optimized via COBYLA

2025-05-04 Thread Gregg Powell via R-help
Hello Christofer, Just writing with a detailed follow-up. Attached is a script I was able to get running with a bit of work. I did not include the script in the ext of this email. It is only attached. Optimization Progress We were initially aiming to solve the dual-slope constrained ordinal mo

Re: [R] non-reproducible eigen() output with MKL

2025-05-04 Thread peter dalgaard
ses, and it is nearly always > FALSE when n is 50.) > > Jesse > > *** > > n <- 50 > set.seed(20250504) > Sigma <- rWishart(1, df = n, Sigma = diag(n))[,,1] > e1 <- eigen(Sigma) > e2 <- eigen(Sigma) > identical(e1, e2) > > __

[R] non-reproducible eigen() output with MKL

2025-05-04 Thread smallepsilon via R-help
ely before each eigen() call, but that did not help. Has anyone seen this behavior before? (When n is 5, the identical() call almost always returns "TRUE". As n increases, the proportion of FALSE results increases, and it is nearly always FALSE when n is 50.) Jesse *** n <- 50 s