  We study the entanglement entropy of harmonic oscillators in noncommutative
phase space. We propose a new definition of quantum R\'enyi entropy based on
Wigner functions in noncommutative phase space. Using the R\'enyi entropy, we
calculate the entanglement entropy of the ground state of the 2D isotropic
harmonic oscillators. We find that for some values of the noncommutative
parameters, the harmonic oscillators can be entangled in noncommutative phase
space. This is a new entanglement-like effect caused by the noncommutativity of
the phase space.
