This study focuses on relational data obtained through object relations. Traditional analysis of relational data often ignores attribute information. Therefore, Mikawa et al. proposed a method to estimate the latent structure of continuous relational data using a generative model and parameter estimation. However, real-world relational data can be discrete, and therefore, we propose a new model for binary relational data using a generative model based on the Bernoulli distribution and the Monte Carlo Expectation-Maximization (EM) algorithm for parameter estimation. We also clarify the effectiveness of the proposed model through simulation experiments using artificial data and real data.