  We present a nonparametric approach to reconstruct the interaction between
dark energy and dark matter directly from SNIa Union 2.1 data using Gaussian
processes, which is a fully Bayesian approach for smoothing data. In this
method, once the equation of state ($w$) of dark energy is specified, the
interaction can be reconstructed as a function of redshift. For the decaying
vacuum energy case with $w=-1$, the reconstructed interaction is consistent
with the standard $\Lambda$CDM model, namely, there is no evidence for the
interaction. This also holds for the constant $w$ cases from $-0.9$ to $-1.1$
and for the Chevallier-Polarski-Linder (CPL) parametrization case. If the
equation of state deviates obviously from $-1$, the reconstructed interaction
exists at $95\%$ confidence level. This shows the degeneracy between the
interaction and the equation of state of dark energy when they get constraints
from the observational data.
