我正在尝试使用SAS来估计变量的sigma。
未加权我使用以下代码:
proc univariate data=test1;
var grams; **unweighted but gives both estimates;
histogram / lognormal(zeta=EST sigma=EST);
run;
但由于这是调查数据,我需要包括权重。 我不知道如何在log surveymeans语句中包含一个带有日志正态分布的变量来估计sigma。或者,如果我应该使用proc调查手段。
proc surveymeans data = test plots = histogram;
weight WTDRD1;
cluster Masked_variance_pseudo_PSU;
strata Masked_variance_pseudo_stratum;
var grams;
run;
ods graphics off;
代码仅使用调查权重估算平均值。不是使用对数正态分布的sigma。