- 博客(18)
- 收藏
- 关注
原创 计量经济学笔记3-Eviews操作-多元线性回归
案例:中国国债发行额国债发行总量(DEBT,亿元)应该与经济总规模,财政赤字的多少,每年的还本付息能力有关系。选择3个解释变量,国内生产总值(GDP:百亿元),财政赤字额(DEF:亿元),年还本付息额(REPAY:亿元),数据来源:国家统计局。这是一个时间序列数据。由于国债数据是1980到2018,所以range到2018。但是REPAY官方数据只到2005年,因此sample区间1980-2005.对各自变量与因变量画出散点图拟合不同的模型以上操作之前的文章里讲过,就不赘述了关于F检验和t检
2022-04-14 19:32:55
5386
原创 计量经济学笔记2-Eviews操作
案例:研究电影票房(income)与制作成本(cost)以及片方分账(producer)之间的关系将这三个变量都选中,然后画图可以手动把图例拖到上面用最小二乘法拟合模型从estimate处可以查看方程其中income是被解释变量,其他都是解释变量如果只用cost来回归,结果如下在实际应用中,要根据研究问题选择合适的变量over...
2022-04-14 19:03:16
1181
原创 运筹学决策分析实例
决策分析可以看成人与自然的博弈。决策者需要在几个选择中做一个决策,自然环境也会有一个状态。每一个决策和自然环境状态的组合都会产生收益(可正可负),列在收益矩阵里面。最终的目标是根据一定的评判标准做出最优决策。接下来会通过英文教材上的一个例题介绍几种评判标准。目录最小收益最大化准则极大似然准则贝叶斯准则敏感性分析借助试验的决策决策树最小收益最大化准则某公司有块地,分析师认为地里有0.25的概率有石油(oil),0.25的概率没有(dry)。另一家公司想买这块地,但这家公司在犹豫是卖掉(sell)还是自
2021-06-08 16:22:16
4133
原创 R语言-运筹学非线性规划实例
若目标函数或约束条件中含有自变量的非线性函数,则这样的规划问题就属于非线性规划。这里简单介绍4种情况下非线性规划的R语言实现。目录无约束规划只含不等式约束的规划系数可变的规划约束条件含等式的规划无约束规划无约束规划中没有约束条件,只需要对目标函数取极值。R语言的nloptr包默认求解的是极小值,因此若要求最大值只需加个负号就行。以这个目标函数为例,我们求它的极大值library(nloptr) eval_f = function(x){ return(list('objective'=-(
2021-06-07 19:45:55
4156
原创 运筹学-图论实例
目录最短路问题Dijkstra算法R语言实现最大流与最小割割集的定义割集的容量最大流-最小割问题求最大流的标号算法最小费用流问题最短路问题求连通图任意两点间总权最小的路Dijkstra算法对每一个顶点给出一组标号,其中第一位代表从始点到该点最短路权的上界,第二位代表该点的上一个点,第三位是T或P,T表示试探性标号,是一种临时标 号,P表示永久性标号,即这个标号不会改变。标有P的点,其第一位标号即为最短路权。算法的每一步都把某一点的T标号改为P标号,当终点被标上P时,全部计算结束。以这个图为例,我们
2021-06-06 10:35:28
4255
2
原创 R语言-运筹学线性规划实例
通过一个简单的小例子来复习线性规划的上机操作。标题例子影子价格灵敏度分析例子首先写出该线性规划问题的数学模型用R语言求解library(lpSolve)direction = 'max'objective.vec = c(2,1) #目标函数系数a1 = c(0,5) #每行约束条件中各变量的系数a2 = c(6,2)a3 = c(1,1)a = rbind(a1,a2,a3) #系数矩阵a.dir = rep('<=',3) #约束方向a.rhs = c(15
2021-06-05 10:20:35
4868
9
原创 R语言-Metropolis Hasting抽样算法
这种算法可以用抽样的方法模拟任何一种分布,计算其均值、方差等特征,前提是知道该分布的密度函数,但不必知道其分布函数。目录引例一引例二Metropolis Hasting算法马尔科夫链(Markov chain)随机游走(Random Walk)应用检验引例一先来看一个简单的小例子,即给出一些均匀分布的随机数,我们可以模拟指数分布。模拟过程如下:首先生成一些[0,1]均匀分布的随机数,可以看作是随机生成的指数分布的概率密度,因此代入指数分布函数(cdf)的反函数,即可求出服从指数分布的随机数。可以再
2021-05-13 21:46:27
3555
1
原创 运筹学-动态规划实例
动态规划是解决多阶段决策问题最优化的一种方法,目标是达到整个过程的整体最优。其基本思想是从最后一阶段开始,用逆序递推方法求解。三个小例子分配医疗队雇佣工人工厂生产分配医疗队世界卫生组织要求将五支医疗队分配到三个国家,分配的医疗队的数目对各国家的效益如表11.1所示。(注:效益指延长该国家人的寿命)设:sis_isi为状态变量,即能够分配给第i个国家(以及剩下的国家)的医疗队数目xix_ixi为决策变量,即实际分配给第i个国家的医疗队数目xi⋆x_i^\starxi⋆为最优的xifn(
2021-05-10 19:02:34
4803
原创 R语言-回归系数的极大似然估计
老师要求我们对回归方程中的回归系数进行极大似然估计,回归方程如下:正常的步骤应该是这样的:步骤一:写出极大似然函数log(β),其中的β为(β0,β1,β2)t(β_0, β_1, β_2)^t(β0,β1,β2)t步骤二:解出log’(β)=0的根但是在写代码的时候不用将极大似然函数表示出来,只需要对log’(β)=0求根即可,这里用的是牛顿法求根。首先来简单介绍一下牛顿算法:这是一种通过不断迭代来求根的算法,先随便给一个初始的x值,再给出迭代原则,最终就能找到方程的根。例如图中是求g
2021-05-05 14:39:07
5428
3
原创 R语言-运筹学指派和运输问题实例
这是我自己的运筹学笔记,由于老师讲的实在一言难尽,导致自己不得不对着英文教材自学,以下的例子就出自英文教材。目录问题提出第一种方法:运输问题第二种方法:指派问题问题提出一家公司打算用三个基站生产四种产品,各基站的生产的成本如表9.27所示。除第2个基站不能生产第3种产品外,其他基站的生产都没有限制。表的最后一行是每种产品要求的生产量,最后一列是基站生产能力的限制。现在问题是如何进行生产分配。我们有两种方法:1.允许一种产品在两个基站分开生产;2.不允许分开生产接下来将对这两种方法依次进行分析
2021-04-29 15:16:45
2366
1
原创 R语言-运筹学虚拟变量运输问题实例
这是我自己的运筹学笔记,由于老师讲的实在一言难尽,导致自己不得不对着英文教材自学,以下的例子就出自英文教材。目录问题介绍R语言解决问题介绍有3条河流向4个城市供水,表格9.10中显示了运费,没有运费则表示无法供水。城市有其最小需求(刚需,第一行),还有奢侈性需求(刚需+额外需求,第二行),首先满足城市的刚需,在此基础上尽量满足其奢侈性需求,并使得费用最小化。难点在于每个城市的需求是不确定的,有一个上界和下界,上界为其他城市的刚需满足后剩下的,下界即为刚需。因此Hollyglass的上界就是(50+
2021-04-27 22:23:51
2042
7
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人