hellenthal-group-UCL / fastGLOBETROTTER

faster haplotype-based admixture inference and dating software
Other
8 stars 1 forks source link

one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly #3

Closed Shaileshdesai76 closed 1 month ago

Shaileshdesai76 commented 1 month ago

i am testing fastGlobetrotter in my machine, using provided data and file, without any change, however running final stage, i encounter error.

Step: R < fastGLOBETROTTER.R tutorial/paramfile.txt tutorial/samplefile.txt tutorial/recomfile.txt 3 --no-save > tutorial/AllFrenchYoruba30gen50prop.fastGT.out

~Downloads/fastGLOBETROTTER/shailesh/tutorial$ R < fastGLOBETROTTER.R tutorial/paramfile.txt tutorial/samplefile.txt tutorial/recomfile.txt 3 --no-save > tutorial/AllFrenchYoruba30gen50prop.fastGT.out Warning message: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly Warning messages: 1: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly 2: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly 3: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly 4: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly 5: In optim(init, ourlhood, method = method) : one-dimensional optimization by Nelder-Mead is unreliable: use "Brent" or optimize() directly

Please can you guide me what to do in this case

Shaileshdesai76 commented 1 month ago

as i am trying to solve this error, my self, i tried with Brent, but didnt work, as i was re-installing again, i found that why compiling, there is some warming like, (attached file) fastglobetrotter_error.txt, but all command is working fine except last one. Step: R < fastGLOBETROTTER.R tutorial/paramfile.txt tutorial/samplefile.txt tutorial/recomfile.txt 3 --no-save > tutorial/AllFrenchYoruba30gen50prop.fastGT.out

Thank you.

Shaileshdesai76 commented 1 month ago

i am closing this issue, as i found if we simply ignore the warning, though the result will be useful and reliable.