有没有办法打印拟合分布的 Proc 单变量输出？

Question

我有一个观察数据的数据集，我想测试它们是否 ~ LogNormal (μ;σ)

我使用过以下内容：

Proc Univariate data=myData;
var Var1;
histogram / LogNormal;
run;

我有不同的表作为输出，我想将“对数正态分布参数”的结果存储在表中。

对数正态分布的参数（示例）

参数	符号	估计
阈值	西塔	0
规模	泽塔	9.09
形状	西格玛	1.04

但是，在我看来，我只能存储基础统计数据（例如正态性测试，偏度测试）

有人可以帮忙吗？

提前致谢！

Answer 1

启用ODS Trace，您可以找到存储此信息的表，然后使用

ods output

输出。例如：

ods trace on;
proc univariate data=sashelp.cars;
    var horsepower;
    histogram / lognormal;
run;
ods trace off;

检查日志，你会看到：

Output Added:
-------------
Name:       ParameterEstimates
Label:      Parameter Estimates
Template:   base.univariate.FitParms
Path:       Univariate.Horsepower.Histogram.Lognormal.ParameterEstimates
-------------

使用

ods output

获取参数估计值并将其保存到名为

lognormal_est

的表中：

proc univariate data=sashelp.cars;
    var horsepower;
    histogram / lognormal;
    ods output ParameterEstimates=lognormal_est;
run;

VarName     Histogram   Distribution    Parameter   Symbol  Estimate
Horsepower  1           Lognormal       Threshold   Theta   0
Horsepower  1           Lognormal       Scale       Zeta    5.321345
Horsepower  1           Lognormal       Shape       Sigma   0.328634
Horsepower  1           Lognormal       Mean                216.0145
Horsepower  1           Lognormal       Std Dev             72.95021