Difference between the N-point functions for the high-pass filtered N_side = 64Planck 2015 CMB estimates and the corresponding means estimated from 1000 Monte Carlo simulations. The Stokes parameters Q_r and U_r were locally rotated so that the correlation functions are independent of coordinate frame. The first row shows results for the 2-point function, from left to right, TQ_r, TU_r, Q_rQ_r, and U_rU_r. The second row shows results for the pseudo-collapsed 3-point function, from left to right, TQ_rQ_r, TU_rU_r, Q_rQ_rQ_r, and U_rU_rU_r, and the third row shows results for the equilateral 3-point function, from left to right, TQ_rQ_r, TU_rU_r, Q_rQ_rQ_r and U_rU_rU_r. The red solid, orange triple-dot-dashed, green dashed and blue dot-dashed lines correspond to the Commander, NILC, SEVEM, and SMICA maps, respectively. The shaded dark and light grey regions indicate the 68% and 95% confidence regions, respectively, estimated using SMICA simulations. See Sect. 7.2 for the definition of the separation angle θ.