SAS求多元线性问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人需要利用SAS求出三个特征方程y（i）=c0+c1*x1+c2*x2+c3*x3+c4*x4+c5*5
data discr1;
input x1 x2 x3 x4 x5 g @@;
cards;
461.6569 488.25 0.0288 156.5592 51.0938 1
396.5858 423 0.0338 122.0231 63.103 2
435.6569 462.875 0.0306 126.0641 64.8653 3
493.8284 531.125 0.0274 153.0133 51.065 1
478.1716 510.75 0.0281 143.818 63.3964 2
447.1421 473 0.0297 122.1139 64.3938 3
449.8995 475.875 0.0295 151.4528 50.4853 1
426.4853 451.75 0.0312 131.0329 64.0884 2
530.1421 561.625 0.0251 130.0581 66.7228 3
507.4853 542.125 0.0265 155.2275 49.1891 1
418 445.375 0.032 127.967 66.9045 2
432.5858 461.5 0.031 121.2331 64.9549 3
478.0711 506.125 0.0278 151.6727 50.5903 1
409 436.5 0.0328 129.6025 67.1023 2
442 470.5 0.0303 125.131 65.0238 3
;
proc discrim data=discr1 testdata=discr2
anova manova simple list testout=discr3;
class g;
proc print;
run;

g表示三种不同组别的特征参数，运行后发生错误应该怎么调试

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：性问题 discrim MANOVA estout simple SAS 线性

你前面已经定义过数据库discr2了吗？

另外如果你把错误信息贴出来，应该求助起来更容易。:)

前面没有定义Discr2和Discr3。错误提示为：

ERROR: 文件 WORK.DISCR2.DATA 不存在。
22   class g;

NOTE: SAS 系统由于错误而停止了该步的处理。
WARNING: 数据集 WORK.DISCR3 可能不完整。在停止该步后有 0 个观测和 0 个变量。

我的目的是求三个基于贝叶斯概率的模型函数y(i),求出三个函数以后，用被检验的图像的一组特征分别带入三个模型函数，求出相应的值，其中哪个值最大，就证明这个被检验的图像与这个模型属于一类品种。

[此贴子已经被作者于2009-5-31 15:37:02编辑过]

我觉得应该用logistic analysis,而不是discriminent analysis
简单的建模模式：
p(g=i)=exp(c0+c1*x1+c2*x2+c3*x3+c4*x4+c5*x5)/[1+exp(c0+c1*x1+c2*x2+c3*x3+c4*x4+c5*x5)],i=1,2,3
sas写法：
 proc logistic data=discr1;
 model g=x1--x5;
 run;
有些option可以在帮助里面查询