\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{mydata}\hlkwb{<-}\hlkwd{data.frame}\hlstd{(X1,X2,Z,Y)}
\hlkwd{library}\hlstd{(drgee)}

\hlcom{#Correct specification}
\hlstd{dr.est} \hlkwb{<-} \hlkwd{drgee}\hlstd{(}\hlkwc{oformula} \hlstd{=} \hlkwd{formula}\hlstd{(Y}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),} \hlkwc{eformula} \hlstd{=} \hlkwd{formula}\hlstd{(Z}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),}
\hlkwc{iaformula} \hlstd{=} \hlkwd{formula}\hlstd{(}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),} \hlkwc{olink} \hlstd{=} \hlstr{"identity"}\hlstd{,} \hlkwc{elink} \hlstd{=} \hlstr{"logit"}\hlstd{,}
\hlkwc{data} \hlstd{= mydata,} \hlkwc{estimation.method} \hlstd{=} \hlstr{"dr"}\hlstd{)}
\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(dr.est)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+         Estimate Std. Error   z value Pr(>|z|)
+ Z       0.960527   0.186412  5.152697        0
+ Z:X1    1.042272   0.142903  7.293565        0
+ Z:X2    0.971475   0.070635 13.753405        0
+ Z:X1:X2 1.031380   0.062547 16.489619        0
\end{verbatim}
\begin{alltt}
\hlcom{#Mis-specification of PS model, correct specification of the OR nuisance model}
\hlstd{dr.est} \hlkwb{<-} \hlkwd{drgee}\hlstd{(}\hlkwc{oformula} \hlstd{=} \hlkwd{formula}\hlstd{(Y}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),} \hlkwc{eformula} \hlstd{=} \hlkwd{formula}\hlstd{(Z}\hlopt{~}\hlstd{X1),}
\hlkwc{iaformula} \hlstd{=} \hlkwd{formula}\hlstd{(}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),} \hlkwc{olink} \hlstd{=} \hlstr{"identity"}\hlstd{,} \hlkwc{elink} \hlstd{=} \hlstr{"logit"}\hlstd{,}
\hlkwc{data} \hlstd{= mydata,} \hlkwc{estimation.method} \hlstd{=} \hlstr{"dr"}\hlstd{)}
\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(dr.est)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+         Estimate Std. Error   z value Pr(>|z|)
+ Z       0.976057   0.184511  5.289964        0
+ Z:X1    1.033960   0.141403  7.312173        0
+ Z:X2    0.979746   0.068582 14.285742        0
+ Z:X1:X2 1.026245   0.059186 17.339432        0
\end{verbatim}
\begin{alltt}
\hlcom{#Correct specification of PS model, mis-specification of the OR nuisance model}
\hlstd{dr.est} \hlkwb{<-} \hlkwd{drgee}\hlstd{(}\hlkwc{oformula} \hlstd{=} \hlkwd{formula}\hlstd{(Y}\hlopt{~}\hlstd{X1),} \hlkwc{eformula} \hlstd{=} \hlkwd{formula}\hlstd{(Z}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),}
\hlkwc{iaformula} \hlstd{=} \hlkwd{formula}\hlstd{(}\hlopt{~}\hlstd{X1}\hlopt{*}\hlstd{X2),} \hlkwc{olink} \hlstd{=} \hlstr{"identity"}\hlstd{,} \hlkwc{elink} \hlstd{=} \hlstr{"logit"}\hlstd{,}
\hlkwc{data} \hlstd{= mydata,} \hlkwc{estimation.method} \hlstd{=} \hlstr{"dr"}\hlstd{)}
\hlkwd{round}\hlstd{(}\hlkwd{coef}\hlstd{(}\hlkwd{summary}\hlstd{(dr.est)),}\hlnum{6}\hlstd{)}
\end{alltt}
\begin{verbatim}
+         Estimate Std. Error  z value Pr(>|z|)
+ Z       1.137266   0.876159 1.298013 0.194283
+ Z:X1    1.058655   0.728725 1.452750 0.146293
+ Z:X2    0.960641   0.333634 2.879326 0.003985
+ Z:X1:X2 1.141582   0.337615 3.381314 0.000721
\end{verbatim}
\end{kframe}
\end{knitrout}