SAS运行 3-way anova ，显示“满足收敛准则，但最终的 Hessian 不是正定的”，怎么办

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

dm 'log;clear;output;clear';

options nodate nocenter pageno=1;

title1 'a';

data one;

input co $ mo $ tre $ per;

cards;

C1  30% 30%-Unfilled  72

C1  30% 30%-Unfilled  77

C1  30% 30%-Unfilled  72

C1  30% 30%-Unfilled  78

C1  30% 30%-Unfilled  75

C1  30% 30%-Unfilled  71

C1  30% 30%-Soil   33

C1  30% 30%-Soil   29

C1  30% 30%-Soil   0

C1  30% 30%-Soil   47

C1  30% 30%-Soil   50

C1  30% 30%-Soil   62

C1  30% 30%-clay   83

C1  30% 30%-clay   76

C1  30% 30%-clay   78

C1  30% 30%-clay   77

C1  30% 30%-clay   80

C1  30% 30%-clay   84

C1  60% 60%-Unfilled  95

C1  60% 60%-Unfilled  96

C1  60% 60%-Unfilled  91

C1  60% 60%-Unfilled  90

C1  60% 60%-Unfilled  94

C1  60% 60%-Unfilled  91

C1  60% 60%-Soil   100

C1  60% 60%-Soil   79

C1  60% 60%-Soil   93

C1  60% 60%-Soil   80

C1  60% 60%-Soil   95

C1  60% 60%-Soil   85

C1  60% 60%-clay   92

C1  60% 60%-clay   93

C1  60% 60%-clay   90

C1  60% 60%-clay   86

C1  60% 60%-clay   95

C1  60% 60%-clay   82

C2  30% 30%-Unfilled  0

C2  30% 30%-Unfilled  0

C2  30% 30%-Unfilled  0

C2  30% 30%-Unfilled  0

C2  30% 30%-Unfilled  0

C2  30% 30%-Unfilled  0

C2  30% 30%-Soil   0

C2  30% 30%-Soil   0

C2  30% 30%-Soil   0

C2  30% 30%-Soil   0

C2  30% 30%-Soil   0

C2  30% 30%-Soil   0

C2  30% 30%-clay   79

C2  30% 30%-clay   0

C2  30% 30%-clay   77

C2  30% 30%-clay   0

C2  30% 30%-clay   75

C2  30% 30%-clay   69

C2  60% 60%-Unfilled  86

C2  60% 60%-Unfilled  76

C2  60% 60%-Unfilled  76

C2  60% 60%-Unfilled  81

C2  60% 60%-Unfilled  75

C2  60% 60%-Unfilled  72

C2  60% 60%-Soil   85

C2  60% 60%-Soil   72

C2  60% 60%-Soil   80

C2  60% 60%-Soil   64

C2  60% 60%-Soil   65

C2  60% 60%-Soil   51

C2  60% 60%-clay   72

C2  60% 60%-clay   82

C2  60% 60%-clay   82

C2  60% 60%-clay   89

C2  60% 60%-clay   93

C2  60% 60%-clay   87

C3  30% 30%-Unfilled  18

C3  30% 30%-Unfilled  1

C3  30% 30%-Unfilled  44

C3  30% 30%-Unfilled  0

C3  30% 30%-Unfilled  0

C3  30% 30%-Unfilled  0

C3  30% 30%-Soil   0

C3  30% 30%-Soil   0

C3  30% 30%-Soil   0

C3  30% 30%-Soil   0

C3  30% 30%-Soil   0

C3  30% 30%-Soil   0

C3  30% 30%-clay   76

C3  30% 30%-clay   73

C3  30% 30%-clay   79

C3  30% 30%-clay   72

C3  30% 30%-clay   60

C3  30% 30%-clay   61

C3  60% 60%-Unfilled  76

C3  60% 60%-Unfilled  74

C3  60% 60%-Unfilled  84

C3  60% 60%-Unfilled  84

C3  60% 60%-Unfilled  90

C3  60% 60%-Unfilled  82

C3  60% 60%-Soil   7

C3  60% 60%-Soil   7

C3  60% 60%-Soil   5

C3  60% 60%-Soil   20

C3  60% 60%-Soil   13

C3  60% 60%-Soil   4

C3  60% 60%-clay   81

C3  60% 60%-clay   82

C3  60% 60%-clay   92

C3  60% 60%-clay   87

C3  60% 60%-clay   73

C3  60% 60%-clay   75

;

run;

proc mixed data=one;

class co mo tre;

model per=mo tre co mo*tre mo*co tre*co mo*tre*co/htype=3;

random co;

lsmeans mo*tre/pdiff adjust=tukey;

run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Hessian ANOVA HESS anov Way

你用了3-way full model, 但是并不是所有的组合都有观测值，所以数据不支持这个模型，看看d.f. 值或者 用 proc freq 就能发现了. 事实上 mo 和tre nested, 你应该从模型中去掉 mo 和 random statement, 用2-way ANOVA.

gerryli 发表于 2018-2-22 02:31
你用了3-way full model, 但是并不是所有的组合都有观测值，所以数据不支持这个模型，看看d.f. 值或者 用 p ...

非常感谢你的回答