In this paper we deal with the problem of the joint optimization of the precoders, equalizers and relay beamformer of a multiple-input multiple-output interfering relay channel. This network can be regarded az a generalized model for both one-way and two-way relay channels with/without direct interfering links. Unlike the conventional design procedures, we assume that the Channel State Information (CSI) is not known perfectly. The imperfect CSI is described using the norm bounded error framework. We use a system-wide Sum Mean Square Error (SMSE) based problem formulation which is constrained using the transmit power of the terminals and the relay node. The problem at hand, from a worst-case design perspective, is a multilinear, and hence, a nonconvex problem which is also semiinfinite in its constraints. We use a generalized version of the Peterson’s lemma to handle the semi-infiniteness and reduce the original problem to a single Linear Matrix Inequality (LMI). However, this LMI is not convex, and to resolve this issue we propose an iterative algorithm based on the alternating convex search methodology to solve the aforementioned problem. Finally simulation results, i.e., the convergence of the proposed algorithm and the SMSE properties, are included to asses the performance of the proposed algorithm.