体均数μo+β1。若β1<0,说明服药3个月时的服药人群平均体重低于服药前的平均体重,即 该减肥药有效,反之无疗效;t1=0,t2=1时,p为服药6个月时的体重总体均数uo+B2,而B2<0 和β2>0同样反映该减肥药有效或无效。若P2<β1,说明服药6个月时的服药人群平均体重低 于3个月时的服药人群平均体重到期间。我们同样借助 Stata软件对上述资料用混合模型进 行统计分析,相应的 Stata软件的数据格式如下。 0 1234 000 23451234 1 0000000000 47 48 00000 其中t1和t2为服药时间的自变量,其他与例12-1相同,相应的 Stata操作命令如下: Stata命令为 treg y tI t2,i(no 相应输出结果如下 Random-effects GLs regression Number of obs Group variable (i):no Number of groups R 0.955 Obs per group: min between avg= overall =0. 4602 max Random effects u i gaussian Wald chi2(2 170.18 corr(u i, X) 0(assumed) Prob> chi2 0.0000 Coef. Std. Err z P>z [95% Conf. Interval] -1.2,3829708 3.130.002-1.950609 4493909 -4.8,3829708-12.530.000-5.550609-4.049391 cons 51.61.10453646.720.00049.4351553.76485 sI ga_u|2.394438 60553007 rho|,93989071 (fraction of variance due to u_i 3个月时的体重与6个月时的体重比较的 Stata命令和输出结果如下 test tl=t2(Ho:β2=B1) (1)t1-t2=0.0 chi2(1)=88.36 ob>chi2=0.0000 B1估计值为-1.2(kg),β2估计值为-48(kg),服药前体重总体均数μo的估计值为体均数0+1。若1<0，说明服药 3 个月时的服药人群平均体重低于服药前的平均体重，即 该减肥药有效，反之无疗效；t1=0，t2=1 时，为服药 6 个月时的体重总体均数0+2，而2<0 和2>0 同样反映该减肥药有效或无效。若2<1，说明服药 6 个月时的服药人群平均体重低 于 3 个月时的服药人群平均体重到期间。我们同样借助 Stata 软件对上述资料用混合模型进 行统计分析，相应的 Stata 软件的数据格式如下。 no y t1 t2 1 48 0 0 2 53 0 0 3 52 0 0 4 52 0 0 5 53 0 0 1 46 1 0 2 51 1 0 3 52 1 0 4 51 1 0 5 52 1 0 1 42 0 1 2 47 0 1 3 48 0 1 4 48 0 1 5 49 0 1 其中 t1 和 t2 为服药时间的自变量，其他与例 12-1 相同，相应的 Stata 操作命令如下： Stata 命令为 xtreg y t1 t2 ,i(no) 相应输出结果如下： Random-effects GLS regression Number of obs = 15 Group variable (i) : no Number of groups = 5 R-sq: within = 0.9551 Obs per group: min = 3 between = . avg = 3.0 overall = 0.4602 max = 3 Random effects u_i ~ Gaussian Wald chi2(2) = 170.18 corr(u_i, X) = 0 (assumed) Prob > chi2 = 0.0000 ------------------------------------------------------------------------------ y | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- t1 | -1.2 .3829708 -3.13 0.002 -1.950609 -.4493909 t2 | -4.8 .3829708 -12.53 0.000 -5.550609 -4.049391 _cons | 51.6 1.104536 46.72 0.000 49.43515 53.76485 -------------+---------------------------------------------------------------- sigma_u | 2.394438 sigma_e | .60553007 rho | .93989071 (fraction of variance due to u_i) ------------------------------------------------------------------------------ 3 个月时的体重与 6 个月时的体重比较的 Stata 命令和输出结果如下： test t1=t2 （H0：2＝1） ( 1) t1 - t2 = 0.0 chi2( 1) = 88.36 Prob > chi2 = 0.0000 1 估计值为－1.2(kg)， 2 估计值为－4.8(kg)，服药前体重总体均数0 的估计值为