Modeling threshold conditional heteroscedasticity with regime-dependent skewness and kurtosis