Foro de discusión y soporte para usuarios de R

Buenas, tengo una tabla de datos con este formato:

function x1 x2 x3
2.232 1 1 1.00
2.242 1 1 1.01
2.732 1 1 1.02
2.770 1 2 1.00
1.932 1 2 1.01
2.132 1 2 1.02
3.222 1.2 1 1
..... ... .. ..


Que representa los valores de function(x1, x2, x3) para cada combinación de valores de x1, x2, x3.
Me gustaría generar un gráfico donde cada punto tenga las coordenadas
x=x1, y=x2, z=(media de function(x1, x2) para todo x3).
¿Cómo puedo hacerlo?
Por ejemplo, para los datos de arriba el primer punto sería
x=x1=1, y=x2=1, z=(2.232+2.242+2.732)/3

En realidad mi tabla tiene más columnas y quiero promediar para aquellas variables que no van a ir en los ejes del gráfico.por ejemplo graficar function(x1, x2) cogiendo las medias para las diferentes x3, x4, x5.

Otra opción sería no promediar sino coger el valor maximo de
function(x1, x2) de entre todos los x1, x2, x3
¿Cómo se hace eso?

Y por último, ¿Cómo grafico function(x1, x2, x3) con x=x1, y=x2, z=x3, y los difrerentes valores de function los represento con diferentes colores?

gracias