Kraemer and Jacklin (1979) proposed a method of analysis of univariate dyadic social interactions or relational data, and Mendoza and Graziano (1982) extended this method to multivariate relations. Their approach is based on an analysis-of- variance-type model that contains parameters characterizing the behavior of actors and partners and their interactions on each relation. The techniques presented in this article offer an alternative approach to the multivariate analysis of social interactions by realizing that many relations yield discrete-valued data and thus are better modeled by using methods designed for categorical data. This alternative approach is also more general because it allows more types of models to be fit. We illustrate, using the same data analyzed by the earlier methods.