In this paper, we propose a new algorithm for the joint clustering of two sets of statistical units \(\mathcal {N}\) and \(\mathcal {M}\) which are also equipped with an adjacency structure which is represented by a bipartite network. Our model is based on the fuzzy Partition Around Medoids, and it combines it with techniques for community detection in bipartite complex networks based on Barber modularity maximization. The goal is to produce a partition of \(\mathcal {N}\cup \mathcal {M}\) into clusters, each of which is also identified by two medoids, one in \(\mathcal {N}\) and one in \(\mathcal {M}\) , which represent the typical units in the cluster for each set. Such clusters are optimized so that units in the same cluster both have similar values on their attributes and are likely to be adjacent. We test the algorithm on both simulated and real data, to show how it is able to capture a wide range of different interactions between the distribution of the attributes and the network structure.
