R Linear model predict function for many variables

Question

I have a matrix of variables X and fitted linear model Y ~ X1 + X2 + X3 + X4. Let\'s say I want to predict Y for a vector of means of