I like using gls() to fit regressions b/c it permits use of varFuncs to model heteroscedasticity.
To my dismay, however, I just found out that the predict function used in conjunction with a model object from gls() will not produce se.fit values nor will it generate intervals.
Is there a ready made alternative function, or will I have to write my own to do the job?
Tim
Timothy G. Gregoire
J. P. Weyerhaeuser Professor of Forest Management
School of Forestry & Environmental Studies, Yale University
360 Prospect Street, New Haven, CT 06511-2104 U.S.A.
office: 1.203.432.9398 mobile: 1.203.508.4014, fax: 1.203.432.3809
G&V sampling text: http://crcpress.com/product/isbn/9781584883708