Me gustaría crear un diagrama de dispersión simple en R o MATLAB que incluya dos variables $ x $ y $ y $ que tengan errores asociados, $ \ epsilon_x $ y $ \ epsilon_y $.

Sin embargo, en lugar de agregar barras de error, esperaba crear un "cuadro sombreado" alrededor de cada par $ (x, y) $ donde la altura de la caja varía de ($ y - \ epsilon_y $) a ($ y + \ epsilon_y $) y el ancho de la caja varía de ($ x - \ epsilon_y $) a ($ x + \ epsilon_y $).

¿Es esto posible en R o MATLAB? Si es así, qué paquete o código puedo usar para generar estos gráficos. Idealmente, me gustaría que el paquete también admita límites de error asimétricos.