statsmodels.sandbox.stats.multicomp.varcorrection_pairs_unequal(var_all, nobs_all, df_all) [source]

return joint variance from samples with unequal variances and unequal sample sizes for all pairs

something is wrong

  • var_all (array_like) – The variance for each sample
  • nobs_all (array_like) – The number of observations for each sample
  • df_all (array_like) – degrees of freedom for each sample
  • varjoint (array) – joint variance.
  • dfjoint (array) – joint Satterthwait’s degrees of freedom


(copy, paste not correct) variance is

1/k * sum_i 1/n_i

where k is the number of samples and summation is over i=0,…,k-1. If all n_i are the same, then the correction factor is 1.

This needs to be multiplies by the joint variance estimate, means square error, MSE. To obtain the correction factor for the standard deviation, square root needs to be taken.

TODO: something looks wrong with dfjoint, is formula from SPSS

© 2009–2012 Statsmodels Developers
© 2006–2008 Scipy Developers
© 2006 Jonathan E. Taylor
Licensed under the 3-clause BSD License.