q-analogs of multiple zeta functions and their special values