We introduce a Monte Carlo approach to the calculation of more distant renormalized interactions with higher accuracy than is possible with previous methods. We have applied our method to study the effects of multispin interactions, which turn out to be far more important than commonly assumed. Even though the individual multispin interactions usually have smaller coupling constants than two-spin interactions, they can dominate the effects of two-spin interactions because they are so numerous.