\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{fit.gee}\hlkwb{<-}\hlkwd{lm}\hlstd{(Y}\hlopt{~}\hlstd{Z}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X2}\hlopt{+}\hlstd{Z}\hlopt{:}\hlstd{X1}\hlopt{:}\hlstd{X2}\hlopt{+}\hlstd{ps.fit.c}\hlopt{+}\hlstd{ps.fit.c}\hlopt{:}\hlstd{X1}\hlopt{+}\hlstd{ps.fit.c}\hlopt{:}\hlstd{X2}\hlopt{+}\hlstd{ps.fit.c}\hlopt{:}\hlstd{X1}\hlopt{:}\hlstd{X2)}
\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(fit.gee)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+                 Estimate Std. Error    t value Pr(>|t|)
+ (Intercept)     3.477510   0.028840 120.580161 0.000000
+ Z               0.989859   0.202076   4.898455 0.000001
+ ps.fit.c        6.918620   0.237115  29.178343 0.000000
+ Z:X1            1.022784   0.157344   6.500303 0.000000
+ Z:X2            0.982109   0.076041  12.915514 0.000000
+ X1:ps.fit.c    -2.761434   0.175974 -15.692298 0.000000
+ X2:ps.fit.c     1.441217   0.087193  16.528950 0.000000
+ Z:X1:X2         1.027346   0.068479  15.002324 0.000000
+ X1:X2:ps.fit.c -0.107657   0.077801  -1.383748 0.166747
\end{verbatim}
\begin{alltt}
\hlstd{mu0.gee}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit.gee,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{0}\hlstd{,ps.fit.c)))}
\hlstd{mu1.gee}\hlkwb{<-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(fit.gee,}\hlkwc{newdata}\hlstd{=}\hlkwd{data.frame}\hlstd{(X1,X2,}\hlkwc{Z}\hlstd{=}\hlnum{1}\hlstd{,ps.fit.c)))}
\hlkwd{c}\hlstd{(mu0.gee,mu1.gee,mu1.gee}\hlopt{-}\hlstd{mu0.gee)}
\end{alltt}
\begin{verbatim}
+ [1]  4.456889  2.538435 -1.918454
\end{verbatim}
\end{kframe}
\end{knitrout}