Tengo un conjunto de datos (prueba) como se indica a continuación:

Type    Met1    Met2    Met3    Met4
TypeA   65  43  97  77
TypeA   46  25  76  77
TypeA   44  23  55  46
TypeA   46  44  55  77
TypeA   33  22  55  54
TypeB   66  8   66  47
TypeB   55  76  66  65
TypeB   55  77  88  46
TypeB   36  67  55  44
TypeB   67  55  76  65

He comprobado muchos enlaces en las gráficas de caja, pero aún no he tenido éxito para el tipo de gráfica de caja que quiero. Deseo tener una gráfica de caja con mi eje X que tenga el tipo A (amarillo, naranja) para todos los Mets (Met1, Met2, Met3, Met4). En esencia, quiero algo como lo siguiente (tomado deaquí):

Estoy intentando cosas como

boxplot(formula = len ~ Type , data = test, subset == "TypeA")
boxplot(formula = len ~ Type , data = test, subset == "TypeA", add=TRUE)
Legend(legend = c( "typeA", "typeB" ), fill = c( "yellow", "orange" ) )

Pero no soy capaz de resolverlo con nada de eso. ¿Alguien puede ayudarme a saber cómo puedo hacer tales diagramas de caja en mis datos de prueba de la manera corregida?