In this paper we apply the directional derivative technique to characterize D-multifunction, quasi D-multifunction and use them to obtain ε-optimality for set valued vector optimization problem with multivalued maps. We introduce the notions of local and partial-ε-minimum (weak) point and study ε-optimality, ε-Lagrangian multiplier theorem and ε-duality results.