Hello, I am working now with time series and I found your script which was very useful for me: <script> open b-g.gdt garch 1 1 ; Y series e = $uhat series h = $h dataset addobs 10 a0 = $coeff[2] a1 = $coeff[3] b1 = $coeff[4] series hfc = h # set future errors to their expectation e = misszero(e) # forecast the variance hfc = a0 + a1 * e(-1)^2 + b1 * hfc(-1) smpl 1970 ; print e h hfc --byobs </script> I have question if it is possible in gretl to prepare similar prediction with one regressant in variance equation for normal,t and GED distributions. Moreover I would like to ask if the same procedure is possible for GJR model. I will be very gratefull for any suggestion from you. Best regards, Tomasz Bielawski