Cortical neurons in vivo generate highly irregular spike sequences. Recently, it was experimentally found that the local variation of interspike intervals, LV, is nearly constant for every spike sequence for the same neurons. On the contrary, the coefficient of variation, CV, varies over different spike sequences. Here, we first show that these characteristic features are also applicable in bursting spike sequences that are obtained from the rat gustatory cortex. Next, we show that the conventional leaky integrate-and-fire model does not fully account for reproducing these statistical features in data of real bursting spike sequences. We resolve this difficulty by proposing an alternative neuron model which is a reduction of the bursting neuron model involving the persistent sodium current. Our study implies that (1) the characteristic features of CV and LV are the results of the endogenous bursting and (2) the bursting behavior in the gustatory cortex is caused mainly by the persistent sodium current.