\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlkwd{library}\hlstd{(splines)}
\hlstd{mod1.PSlm1} \hlkwb{<-} \hlkwd{lm}\hlstd{(BPSysAve}\hlopt{~}\hlstd{SmokeNow}\hlopt{+}\hlstd{ps.lr,}\hlkwc{data}\hlstd{=small.nhanes)}
\hlstd{APO.PSlm1.1} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm1,nhanes.allsmoke))}
\hlstd{APO.PSlm1.0} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm1,nhanes.nosmoke))}
\hlstd{APO.PSlm1.1} \hlopt{-} \hlstd{APO.PSlm1.0}
\end{alltt}
\begin{verbatim}
+ [1] -1.10791
\end{verbatim}
\begin{alltt}
\hlstd{mod1.PSlm2} \hlkwb{<-} \hlkwd{lm}\hlstd{(BPSysAve}\hlopt{~}\hlstd{(}\hlnum{1}\hlopt{+}\hlstd{SmokeNow)}\hlopt{*}\hlstd{(ps.lr}\hlopt{+}\hlkwd{I}\hlstd{(ps.lr}\hlopt{^}\hlnum{2}\hlstd{)),}\hlkwc{data}\hlstd{=small.nhanes)}
\hlstd{APO.PSlm2.1} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm2,nhanes.allsmoke))}
\hlstd{APO.PSlm2.0} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm2,nhanes.nosmoke))}
\hlstd{APO.PSlm2.1} \hlopt{-} \hlstd{APO.PSlm2.0}
\end{alltt}
\begin{verbatim}
+ [1] -1.681372
\end{verbatim}
\begin{alltt}
\hlstd{mod1.PSlm3} \hlkwb{<-} \hlkwd{lm}\hlstd{(BPSysAve}\hlopt{~}\hlstd{(}\hlnum{1}\hlopt{+}\hlstd{SmokeNow)}\hlopt{*}\hlkwd{bs}\hlstd{(ps.lr,}\hlkwc{df}\hlstd{=}\hlnum{4}\hlstd{,}\hlkwc{Boundary.knots} \hlstd{=} \hlkwd{range}\hlstd{(}\hlnum{0}\hlstd{,}\hlnum{1}\hlstd{)),}
                    \hlkwc{data}\hlstd{=small.nhanes)}
\hlstd{APO.PSlm3.1} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm3,nhanes.allsmoke))}
\hlstd{APO.PSlm3.0} \hlkwb{<-} \hlkwd{mean}\hlstd{(}\hlkwd{predict}\hlstd{(mod1.PSlm3,nhanes.nosmoke))}
\hlstd{APO.PSlm3.1} \hlopt{-} \hlstd{APO.PSlm3.0}
\end{alltt}
\begin{verbatim}
+ [1] -1.971314
\end{verbatim}
\end{kframe}
\end{knitrout}