Estimating the directed information to infer causal relationships in ensemble neural spike train recordings