\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{fmods}\hlkwb{<-}\hlkwd{formula}\hlstd{(SmokeNow}\hlopt{~}\hlstd{Gender}\hlopt{+}\hlstd{Age}\hlopt{+}\hlstd{Race3}\hlopt{+}\hlstd{Education}\hlopt{+}\hlstd{MaritalStatus}\hlopt{+}\hlstd{HHIncome}\hlopt{+}\hlstd{Poverty)}
\hlstd{ps.mod} \hlkwb{<-} \hlkwd{glm}\hlstd{(fmods,}\hlkwc{data}\hlstd{=small.nhanes,}\hlkwc{family}\hlstd{=}\hlstr{"binomial"}\hlstd{)}
\hlstd{ps.lr} \hlkwb{<-} \hlkwd{predict}\hlstd{(ps.mod,}\hlkwc{type}\hlstd{=}\hlstr{"response"}\hlstd{)}
\hlstd{Y}\hlkwb{<-}\hlstd{small.nhanes}\hlopt{$}\hlstd{BPSysAve}
\hlstd{ps.lr.quints} \hlkwb{<-} \hlkwd{cut}\hlstd{(ps.lr,quints,}\hlkwc{labels}\hlstd{=}\hlnum{1}\hlopt{:}\hlnum{5}\hlstd{)}
\hlstd{p.strat} \hlkwb{<-} \hlkwd{table}\hlstd{(ps.lr.quints)}\hlopt{/}\hlkwd{length}\hlstd{(ps.lr.quints)}
\hlstd{ATE.strat} \hlkwb{<-} \hlkwd{rep}\hlstd{(}\hlnum{NA}\hlstd{,}\hlnum{5}\hlstd{)}
\hlkwa{for}\hlstd{(j} \hlkwa{in} \hlnum{1}\hlopt{:}\hlnum{5}\hlstd{) \{}
    \hlstd{ATE.strat[j]}\hlkwb{<-}\hlkwd{mean}\hlstd{(Y[Smoke}\hlopt{==}\hlnum{1} \hlopt{&} \hlstd{ps.lr.quints}\hlopt{==}\hlstd{j])}\hlopt{-}\hlkwd{mean}\hlstd{(Y[Smoke}\hlopt{==}\hlnum{0} \hlopt{&} \hlstd{ps.lr.quints}\hlopt{==}\hlstd{j])}
\hlstd{\}}
\hlkwd{sum}\hlstd{(ATE.strat}\hlopt{*}\hlstd{p.strat)}
\end{alltt}
\begin{verbatim}
+ [1] -1.816879
\end{verbatim}
\end{kframe}
\end{knitrout}