<\/p><\/div>
\nWe have already introduced binary variables as a special type of discrete variable that can be used to indicate whether or not a subject has a characteristic of interest, such as gender for a person or ownership of a captive insurance company for a firm. Binary variables also describe whether or not an event of interest has occurred, such as an accident. A model with a binary dependent variable allows one to predict whether an event has occurred or a subject has a characteristic of interest.<\/p>\n
Example: MEPS Expenditures.<\/strong> Section 11.4 will describe an extensive database from the Medical Expenditure Panel Survey (MEPS) on hospitalization utilization and expenditures. For these data, we will consider Like the linear regression techniques introduced in prior chapters, we are interested in using characteristics of a person, such as their age, sex, education, income and prior health status, to help explain the dependent variable y<\/em>. Unlike the prior chapters, now the dependent variable is discrete and not even approximately normally distributed. In limited circumstances, linear regression can be used with binary dependent variables – this application is known as a linear probability model<\/em>.<\/p>\n
\n\\begin{eqnarray*}\\small
\ny_i = \\left\\{
\n\\begin{array}{ll}
\n1 & i\\text{th person was hospitalized during the period} \\\\
\n0 & \\text{otherwise}
\n\\end{array}
\n\\right. .
\n\\end{eqnarray*} There are n<\/em>=2,000 persons in this sample, distributed as:
\n\\begin{matrix}{\\text{Table 11.1 Hospitalization by Gender}} \\\\ \\small
\n\\begin{array}{ll|ll}\\hline
\n& & \\text{Male} & \\text{Female} \\\\ \\hline
\n\\text{Not hospitalized} & y=0 & 902 (95.3\\%) & ~~~~941 (89.3\\%) \\\\
\n\\text{Hospitalized} & y=1 & ~~44 ( 4.7\\%) & ~~~~113 (10.7\\%) \\\\
\n\\text{Total} & & 946 & 1,054 \\\\ \\hline
\n\\end{array}
\n\\end{matrix} Table 11.1 suggests that gender has an important influence on whether someone becomes hospitalized.<\/p>\nR Code and Output for Table 11.1<\/strong><\/i><\/a><\/h2>\n
R-Code<\/strong>\r\nHexpend <- read.table(file=\"http:\/\/instruction.bus.wisc.edu\/jfrees\/jfreesbooks\/Regression%20Modeling\/BookWebDec2010\/CSVData\/HealthExpend.csv\", header=TRUE, sep=\",\")\r\nattach(Hexpend)\r\nPOSEXP <- 1*(EXPENDIP>0)\r\nmytable <- table(POSEXP,GENDER)\r\nmytable \r\nprop.table(mytable, 2)<\/pre>\n
R-Code Output<\/strong>\r\n GENDER\r\nPOSEXP 0 1\r\n 0 902 941\r\n 1 44 113\r\n\r\n GENDER\r\nPOSEXP 0 1\r\n 0 0.9535 0.8928\r\n 1 0.0465 0.1072\r\n<\/pre>\n<\/div>\n