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

Is there a way to print Proc Univariate output for Fitted distribution?

提问人:Diana 提问时间:11/14/2023 最后编辑:Stu SztukowskiDiana 更新时间:11/14/2023 访问量:22

问:

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

我使用了以下内容:

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

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

对数正态分布的参数(示例)

参数 象征 估计
门槛 西塔 0
规模 截塔 9.09
形状 西格马 1.04

但是,在我看来,我只能存储基本统计数据(例如正态性检验、偏度检验)

有人可以帮忙吗?

提前致谢!

SAS 分销

评论


答:

1赞 Stu Sztukowski 11/14/2023 #1

启用 ODS 跟踪,您可以找到存储此信息的表,然后使用 输出。例如: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

proc univariate data=sashelp.cars;
    var horsepower;
    histogram / lognormal;
    ods output ParameterEstimates;
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