Marginal probability density and cumulative distribution functions are presented for multidimensional variables defined by non-singular affine transformations of vectors of independent two-piece normal variables, the most important subclass of Ferreira and Steel’s general multivariate skewed di… stributions. The marginal functions are obtained by first expressing the joint density as a mixture of Arellano-Valle and Azzalini’s unified skew-normal densities and then using the property of closure under marginalization of the latter class.