数学建模专题汇总-离散模型

离散模型

§ 1 离散回归模型

一、离散变量

如果我们用0，1，2，3，4，⋯说明企业每年的专利申请数，申请数是一个离散的变量,但是它是间隔尺度变量，该变量类型不在本章的讨论的被解释变量中。但离散变量0和1可以用来说明企业每年是否申请专利的事项，类似表示状态的变量才在本章的讨论中。在专利申请数的问题中，离散变量0，1，2，3 和4 等数字具

有具体的经济含义，不能随意更改；而在是否申请专利的两个选择对象的选择问题中，数字0和1只是用于区别两种不同的选择，是表示一种状态。本专题讨论有序尺度变量和名义尺度变量的被解释变量。

、离散因变量

在讨论家庭是否购房的问题中，可将家庭购买住房的决策用数字1 表示，而将家庭不购买住房的决策用数字0 表示。

1 yes

x

0 no

如果x 作为说明某种具体经济问题的自变量，则应用以前介绍虚拟变量知识就足够了。如果现在考虑某个家庭在一定的条件下是否购买住房问题时，则表示状态的虚拟变量就不再是自变量，而是作为一个被说明对象的因变量出现在经济模型中。因此，需要对以前讨论虚拟变量的分析方法进行扩展，以便使其能够适应分析类似家庭是否购房的问题。因为在家庭是否购房问题中,虚拟因变量的具体取值仅是为了区别不同的状态，所以将通过虚拟因变量讨论备择对象选择的回归模型称为离散选择模型。

三、线性概率模型

现在约定备择对象的0 和1 两项选择模型中，下标i 表示各不同的经济主体，取值

0或l的因变量 y i表示经济主体的具体选择结果，而影响经济主体进行选择的自变量 x i 。如果选择响应YES 的概率为 p(y i 1/ x i ) ，则经济主体选择响应NO 的概率为 1 p(y i 1/ x i)，

则E(y i /x i) 1 p(y i 1/x i) 0 p(y i 0/x i)＝ p(y i 1/x i)。根据经典线性回归，我们知道其总体回归方程是条件期望建立的，这使我们想象可以构造线性概率模型

p(y i 1/ x i) E(y i / x i) x iβ

0 1 x i1 L k x ik u i

描述两个响应水平的线性概率回归模型可推知，根据统计数据得到的回归结果并不一定能够保证回归模型的因变量拟合值界于[0，1]。如果通过回归模型式得到的因变量拟合值完全偏离0或l两个数值，则描述两项选择的回归模型的实际用途就受到很大的限制。为避免出现回归模型的因变量预测值偏离0或1的情形，需要限制因变量的取值范围并对回归模型式进行必要的修正。由于要对其进行修正，那么其模型就会改变，模型改变会导致似然函

数改变，这就是我们下面要讨论的。现在我们讨论的模型与判别分析的目的是一样的，但有区别。

§ 2 二元离散选择模型

一、效用函数

为了使得二元选择问题的有进一步研究可能，首先建立一个效用函数。在讨论家庭是否购房的问题中，可将家庭购买住房的决策用数字1 表示，而将家庭不购买住房的决策用数字0表示。用 U i1表示第i个人选择买房的效用， U i0表示第i个人

选择不买房的效用。其效用均为随机变量，于是有

1 1 1

U i1 1 X iβ1u i1(1)

U i00X iβ0u i0(2)

将(1)-(2),得

U i 1

U i 0

1

2

X i (β1

β0

) (u i 1

u i 0

)

记：

y i *

U i 1

U i 0

12

β*

β1 β0

10

u i *

u i 1

u i 0

则有 **

Y i

*

*

X i

β*

u i *

，格林称该模型为潜回归

这是二元选择模型的切入点。称 Y i *

为过渡变量（潜在的） 测

的。

当效用差 Y i *

大于零，则应该选 “1，”即购房；

，这个变量是不可观

当效用差Y i*小于零，则应该选“0，”即不购房。

故p(Y i 1) p(Y i*0) P(u i*－X iβ*) 1 F －X iβ*

p(Y i 0) p(Y i*0) P(u i*－X iβ*) F－X iβ*

此处已经通过Y i*，将自变量与事件发生的概率联系起来了。为概率提供了一个潜在的结构模型。

现在的问题是F ( )服从何种分布？F ( )既然是分布函数，则必须满足分布函数的条件.

二、两类常用的模型

根据以上的分析，我们的问题已经转化为作为F ( )有什么形状，即密度函数f 具有什么样的函数形式。采用累积标准正态概率分布函数的模型称作Probit 模型，或概率单位模型，用正态分布的累积概率作为Probit 模型的预测概率。另外logistic

函数也能满足这样的要求，采用 logistic 函数的模型称作 logit 模型，或对数单位模 型。注：分布在此时是以 y 轴为对称。

（一） Logit 模型 因为

p(Y i 1) p(Y i

0) P(u i *

X i β*

) P(u i *

X i β* ) F

X i β*

如果我们取 F .）为逻辑函数

LOGIT ),即

(x) F (x)

x

1e xx

1 e 1 e

满足分布函数的条件） ，有

p(Y i 1) F

X i β

(

*

1 X i β*

) 1 e

1* X i β

*

e * X i β* 1 e X i β为了更简化模型

Y i

*

*

X i β*

u i *

，我们令