\begin{knitrout}
\definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe}
\begin{alltt}
\hlstd{beta.intercept}\hlkwb{<-}\hlnum{110}
\hlstd{beta.smoke}\hlkwb{<-}\hlnum{3}
\hlstd{beta.gender}\hlkwb{<-}\hlnum{3.5}
\hlstd{beta.age}\hlkwb{<-}\hlnum{0.5}
\hlstd{beta.ed.grad}\hlkwb{<-}\hlopt{-}\hlnum{5}
\hlstd{beta.married}\hlkwb{<-}\hlopt{-}\hlnum{2.5}
\hlstd{beta.HH65}\hlkwb{<-}\hlnum{6}
\hlstd{beta.HH75}\hlkwb{<-}\hlnum{7}
\hlstd{beta.HH100}\hlkwb{<-}\hlnum{8}
\hlstd{beta.poverty}\hlkwb{<-}\hlopt{-}\hlnum{2}
\hlstd{beta.agepov}\hlkwb{<-}\hlnum{0.2}

\hlstd{Y}\hlkwb{<-}\hlkwd{rep}\hlstd{(}\hlnum{0}\hlstd{,}\hlkwd{nrow}\hlstd{(small.nhanes))}
\hlstd{Y}\hlkwb{<-}\hlstd{Y}\hlopt{+}\hlstd{beta.intercept}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{SmokeNow}\hlopt{==}\hlnum{1}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{SmokeNow}\hlopt{==}\hlnum{1}\hlstd{]}\hlopt{+}\hlstd{beta.smoke}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{Gender}\hlopt{==}\hlstr{'male'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{Gender}\hlopt{==}\hlstr{'male'}\hlstd{]}\hlopt{+}\hlstd{beta.gender}
\hlstd{Y}\hlkwb{<-}\hlstd{Y}\hlopt{+}\hlstd{beta.age}\hlopt{*}\hlstd{small.nhanes}\hlopt{$}\hlstd{Age}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{Education}\hlopt{==}\hlstr{'College Grad'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{Education}\hlopt{==}\hlstr{'College Grad'}\hlstd{]}\hlopt{+}\hlstd{beta.ed.grad}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{MaritalStatus}\hlopt{==}\hlstr{'Married'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{MaritalStatus}\hlopt{==}\hlstr{'Married'}\hlstd{]}\hlopt{+}\hlstd{beta.married}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'65000-74999'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'65000-74999'}\hlstd{]}\hlopt{+}\hlstd{beta.HH65}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'75000-99999'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'75000-99999'}\hlstd{]}\hlopt{+}\hlstd{beta.HH75}
\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'more 99999'}\hlstd{]}\hlkwb{<-}\hlstd{Y[small.nhanes}\hlopt{$}\hlstd{HHIncome}\hlopt{==}\hlstr{'more 99999'}\hlstd{]}\hlopt{+}\hlstd{beta.HH100}
\hlstd{Y}\hlkwb{<-}\hlstd{Y}\hlopt{+}\hlstd{beta.poverty}\hlopt{*}\hlstd{small.nhanes}\hlopt{$}\hlstd{Poverty}
\hlstd{Y}\hlkwb{<-}\hlstd{Y}\hlopt{+}\hlstd{beta.agepov}\hlopt{*}\hlstd{small.nhanes}\hlopt{$}\hlstd{Age}\hlopt{*}\hlstd{small.nhanes}\hlopt{$}\hlstd{Poverty}    \hlcom{#Not fitted}
\hlstd{Y}\hlkwb{<-}\hlstd{Y}\hlopt{+}\hlkwd{rnorm}\hlstd{(}\hlkwd{nrow}\hlstd{(small.nhanes))}\hlopt{*}\hlnum{3}
\hlstd{small.nhanes}\hlopt{$}\hlstd{BPSysAve}\hlkwb{<-}\hlstd{Y}
\end{alltt}
\end{kframe}
\end{knitrout}