Skip to main content
U.S. flag

An official website of the United States government

SAS code for analyzing water temperature data

February 9, 2018

This code may be used to fit linear models with multivariate random effects and heterogeneous measurement-level residual variances. The code as written may be used to estimate associations between water temperature ('temp') and continuous year ('yearctr'), study reach (or field station; 'fs'), log-transformed mean July water discharge (in 1000 cms units; 'logmeanJulycms1000'), number of days from a central sampling date (for a given year, days from a standard month and day; 'jdatectr'), time of sampling (in fractions of hours from noon on the given date; 'timectrhr') and interactions thereof.