A novel method to estimate social effect coefficients in the popular so-called linear-in-means regression model in the Social Sciences is presented here that utilizes non-experimental multidimensional network data. The procedure can accommodate social interactions that correlate with the error in the model by making use of a different set of networks links among the same observations that are exogenous in the traditional sense. In particular the full observability of a two-layered mph{multiplex} network data structure is assumed here to propose a new Generalized 3-Stage Least Squares (G3SLS) estimator that is consistent, asymptotically normally distributed, and also easy to implement using widely-used existing statistical software because of its closed-form definition. The underlying assumptions are general enough to accommodate common problems with observational data such as measurement error, simultaneity, and unobserved heterogeneity. Monte Carlo exercises confirm the good small sample performance of the proposed G3SLS estimator in these scenarios. An empirical application finds a positive and significant peer effects in citations among research articles published in top general-interest journals in economics.