rubin_rules function

Combine estimates using Rubin's rules