Okay, I checked with the nice person from R that wrote the function and it appears that he missed squared (for SD), so the manual calculation is correct. I can close this long thread
Thanks Dason and Greta for your help
This is the best statistics forum.