This paper proposes a new approach for optimizing several mean responses in a general linear model. This approach combines the characteristics of both P-model and chance constraint concepts to formulate an equivalent deterministic problem and then solve the latter using familiar optimization techniques. Parameter estimates of the general linear model are given as a-by-product of the proposed procedure. A numerical example is given to illustrate the proposed algorithm using Mathematica software.