提问人:Diana 提问时间:11/14/2023 最后编辑:Stu SztukowskiDiana 更新时间:11/14/2023 访问量:22
有没有办法打印拟合分布的 Proc 单变量输出?
Is there a way to print Proc Univariate output for Fitted distribution?
问:
我有一个观测数据数据集,我想测试它们是否是 ~ LogNormal (μ;σ)
我使用了以下内容:
Proc Univariate data=myData;
var Var1;
histogram / LogNormal;
run;
作为输出,我有不同的表,我想将“对数正态分布参数”的结果存储在表中。
对数正态分布的参数(示例)
参数 | 象征 | 估计 |
---|---|---|
门槛 | 西塔 | 0 |
规模 | 截塔 | 9.09 |
形状 | 西格马 | 1.04 |
但是,在我看来,我只能存储基本统计数据(例如正态性检验、偏度检验)
有人可以帮忙吗?
提前致谢!
答:
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
上一个:图像无法呈现
评论