使用rugarch包，在ugarchspec函数里，mean.model = list(include.mean = T, archm = T)，这两个变量都设成T，variance.model = list(model = "sGARCH")，model选择sGARCH，这样参数估计出来的便是GARCH-M model了。