\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{mod1} \hlkwb{<-} \hlkwd{lm}\hlstd{(fmodi,}\hlkwc{data}\hlstd{=small.nhanes)}
\hlstd{m1}\hlkwb{<-}\hlkwd{predict}\hlstd{(mod1,nhanes.allsmoke)}
\hlstd{m0}\hlkwb{<-}\hlkwd{predict}\hlstd{(mod1,nhanes.nosmoke)}
\hlstd{ps.trunc}\hlkwb{<-}\hlkwd{pmin}\hlstd{(}\hlkwd{pmax}\hlstd{(}\hlnum{0.001}\hlstd{,ps.lr),}\hlnum{0.999}\hlstd{)}
\hlstd{APO.1}\hlkwb{<-}\hlkwd{mean}\hlstd{(Smoke}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{m1)}\hlopt{/}\hlstd{ps.trunc)}\hlopt{+}\hlkwd{mean}\hlstd{(m1)}
\hlstd{APO.0}\hlkwb{<-}\hlkwd{mean}\hlstd{((}\hlnum{1}\hlopt{-}\hlstd{Smoke)}\hlopt{*}\hlstd{(Y}\hlopt{-}\hlstd{m0)}\hlopt{/}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{ps.trunc))}\hlopt{+}\hlkwd{mean}\hlstd{(m0)}
\hlstd{APO.1}\hlopt{-}\hlstd{APO.0}
\end{alltt}
\begin{verbatim}
+ [1] -1.854776
\end{verbatim}
\begin{alltt}
\hlstd{small.nhanes}\hlopt{$}\hlstd{ps} \hlkwb{<-} \hlstd{ps.lr}
\hlstd{fmodia}\hlkwb{<-}\hlkwd{formula}\hlstd{(BPSysAve}\hlopt{~}\hlstd{SmokeNow}\hlopt{+}\hlstd{Gender}\hlopt{+}\hlstd{Age}\hlopt{+}\hlstd{Race3}\hlopt{+}\hlstd{Education}\hlopt{+}\hlstd{MaritalStatus}\hlopt{+}\hlstd{HHIncome}\hlopt{+}\hlstd{Poverty}\hlopt{+}
                \hlstd{SmokeNow}\hlopt{:}\hlstd{HHIncome}\hlopt{+}\hlstd{SmokeNow}\hlopt{:}\hlstd{Gender}\hlopt{+}\hlstd{SmokeNow}\hlopt{:}\hlstd{Age}\hlopt{+}
                \hlkwd{I}\hlstd{(SmokeNow}\hlopt{/}\hlstd{ps)}\hlopt{+}\hlkwd{I}\hlstd{((}\hlnum{1}\hlopt{-}\hlstd{SmokeNow)}\hlopt{/}\hlstd{(}\hlnum{1}\hlopt{-}\hlstd{ps)))}
\hlstd{mod1.lmX} \hlkwb{<-} \hlkwd{lm}\hlstd{(fmodia,}\hlkwc{data}\hlstd{=small.nhanes)}
\hlstd{nhanes.allsmoke} \hlkwb{<-} \hlstd{nhanes.nosmoke} \hlkwb{<-} \hlstd{small.nhanes;}
\hlstd{nhanes.allsmoke}\hlopt{$}\hlstd{SmokeNow} \hlkwb{<-} \hlnum{1}
\hlstd{nhanes.nosmoke}\hlopt{$}\hlstd{SmokeNow} \hlkwb{<-} \hlnum{0}
\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.lmX,nhanes.allsmoke))}\hlopt{-}\hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.lmX,nhanes.nosmoke))}
\end{alltt}
\begin{verbatim}
+ [1] -1.818062
\end{verbatim}
\end{kframe}
\end{knitrout}