The SURVEYMEANS procedure might produce incorrect GEOMEAN results when you use strata as domains


Some of the domain geometric means and their standard errors might be incorrect when the following conditions are met:

The incorrect results occur in the "Geometric Means for Domains" output table that is displayed as well as in data sets saved using the ODS OUTPUT DomainGeoMeans=SAS-data-set statement.

Workaround

You can prevent this issue by creating a new variable that is a copy of the stratification variable and using the new variable instead of the original stratification variable in the DOMAIN statement. 

Here is an example:

data MyData;
  set MyData;
  StrataVarCopy=StrataVar;
run;
proc surveymeans data=MyData geomean;
  var y;
  weight wt;
  strata StrataVar;
  domain DomainVar*StrataVarCopy;
run;