May be function uses Welch-Satterthwaite correction?