GLM过程的lsmean结果疑问

2024-07-13 13:04| 来源: 网络整理| 查看: 265

有一数据集，包含以下变量，因为男女、城乡、年份等年龄组成不一样，进行TC的均数比较时需要用年龄进行调节，使用GLM过程，用年龄age作为协变量计算lsmean。 [td]#变量类型长度输出格式标签1sex数值81.性别1=男2=女2age数值8 年龄（岁）3tc数值85.2TC4yrs字符4 年份5district字符4$4.城乡6id数值8 ID号

为了进行年份间、城乡间、男女间进行比较，分别运行以下3个glm过程。

libname a "d:\"; proc sort data=a.all;by district sex; run; proc glm data=a.all; title '表1 城乡、男女TC年龄调节均数的年份间比较'; class yrs district sex; model TC=yrs district sex age; by district sex; lsmeans yrs/stderr; run; proc sort data=a.all;by yrs sex; run; proc glm data=a.all; title '表2 不同年份、男女TC年龄调节均数的城乡间比较'; class yrs district sex; model TC=yrs district sex age; by yrs sex; lsmeans district/stderr; run; proc sort data=a.all;by yrs district; run; proc glm data=a.all; title '表3 不同年份、城乡TC年龄调节均数的男女间比较'; class yrs district sex; model TC=yrs district sex age; by yrs district; lsmeans sex/stderr; run;复制代码分别得出结果，按道理，同样是用age年龄进行调节，最终结果应该是一样的，比如，2001年城市男性的调节均数（TC的lsmean）在表1、表2、表3中应该是一样的，但实际不同，为什么？

表1结果：城乡=城市性别1=男2=女=1

[td]yrstc LSMEAN标准误差Pr > |t|20014.704213650.05855886 |t|城市4.635279570.03999276 |t|14.637513810.04075315

【本文地址】

公司简介

联系我们