what is the r code for simulating PowerGARCH,Threshold GARCH,and GJR GARCH in order to capture heteroscedasticity,volatility clustering,etc,,so that i can have simulation of mean part and simulation on innovation part. thanks [[alternative HTML version deleted]]