Rare-earth defect pairs in GaN: LDA+U calculations