tag 标签: 最大值经管大学堂:名校名师名课

相关帖子

版块 作者 回复/查看 最后发表
如何每隔5行求取一次最大值和平均数 attachment R语言论坛 xingzhaoh 2013-7-10 19 21311 oliyiyi 2023-9-9 13:28:07
stata一组变量求最大值均值等 Stata专版 gongxuche1991 2013-5-6 10 12304 小点喜欢晴天 2017-6-30 10:37:53
恋爱、婚姻与家庭及其价值动因 哲学与心理学版 choudh 2013-6-28 2 5930 wonderliza 2014-10-29 09:55:39
悬赏 如何实现在stata中挑选出过去若干个最大值 - [悬赏 200 个论坛币] Stata专版 xu_tony77 2013-1-15 8 10760 ttangssong 2014-5-15 09:16:35
【讨论】如何挑出4个最大值和最小值 Stata专版 dxystata 2013-6-21 9 2737 luisluan 2013-12-15 19:49:43
情感强度第二定律 跨学科讨论区 choudh 2013-9-7 1 2690 nannan1107 2013-9-7 17:53:19
概率论与数理统计(12) 学者专栏 fancunhui 2013-8-19 0 7549 fancunhui 2013-8-19 16:17:43
数据库处理,寻找最大值退回变量名问题。 SAS专版 liulubucm 2013-5-23 4 2000 liulubucm 2013-7-16 11:14:31
悬赏 请教高手,求二十元函数的最大值 - [悬赏 2 个论坛币] 计量经济学与统计软件 TimeT 2013-6-17 2 1683 TimeT 2013-6-20 22:27:13
悬赏 200论坛币求请问怎样求和、均数、最大值、最小值 - [!reward_solved!] SAS专版 lnlhckao123 2013-5-4 6 1819 zs236 2013-5-4 18:17:43
求最大值,望指教 attach_img SAS专版 我是守望者 2013-4-9 7 1816 zhengbo8 2013-4-10 13:58:51
如何求最大值 attachment Stata专版 qiaqiao 2013-1-28 2 2524 迅哥儿2012 2013-3-24 20:21:14
一个变量的标准差,方差,最大值最小值等等 SAS专版 jingliang30 2013-2-21 22 11117 webgu 2013-2-22 18:45:28
新手说说SAS和R SAS专版 complicated 2013-1-30 11 12468 ilikeyahoo 2013-2-22 11:57:57
spss分析中缺失值的处理问题!求解 SPSS论坛 595727181 2013-1-14 8 5330 595727181 2013-1-17 21:18:43
[求助]如何求一列变量的最大值? SPSS论坛 思亭 2009-5-26 4 7032 panda56 2011-9-2 12:55:57
[求助]请教如何在同一个变量下一定范围内找最大值 SAS专版 inertialife 2008-8-22 1 2353 inertialife 2008-8-22 19:29:00
如何将横坐标的刻度最大值调整?[R] R语言论坛 hayes 2008-4-2 2 5181 hayes 2008-4-4 11:48:00
请问: mle 找不到最大值 Stata专版 phoebs 2007-9-28 3 3281 phoebs 2007-9-28 22:40:00

相关日志

分享 R语言的一个小范例:数值模拟和绘图
accumulation 2015-5-23 01:48
下午上课,要给经济工程专业的本科生简要介绍一下用于统计分析和绘图的R语言。考虑到他们所在的年级并没有系统的学完统计学和计量经济学,打算不从回归分析入手,先教他们对R的基本操作有个了解。写了如下这段代码,准备课上带同学们做一下。顺便发到网上,供大家参考。 背景: 假定Beta 是 一个由alpha(0~1)决定的量,满足关系 问在alpha的取值空间内,Bata值的走势如何,以及最大值出现在什么地方,并绘图说明。 talk is cheap. show me the code (废话少说,放码过来) 一点一点写,并适当做说明 qiu - function(alpha){(alpha / (1-alpha))^{-alpha}} 这段代码的意思是,定义一个名为qiu的function,输入值alpha之后,会自动求得对应的值。比如 qiu(0.67) 0.622206 现在我们来搞alpha。 x - seq(0.001,0.99, by=0.001) 定义一个数组x,取值从0.01开始,每次增加0.001单位,直到0.999。然后看看x这个数组的情况吧 x length(x) 结果就不贴出来了。现在,把x的值交给qiu来处理,生成的一系列值,我们定义为y。 y - qiu(x) 画个草图来看看分布情况如何 plot(y~x) 显然,随着x从0到1, y并不是单调变化的:先增大,后减小。那么,y的最高值是多少?对应的x值是多大? max(y) x 第一行告诉你y的最大值1.3211 第二行的which.max(y)告诉你当y值最大时,所对应的编号是218。把218带入x中,找到对应的x的值:x =0.218 schx=c(x ) schy=c(max(y)) 分别将y值最大时的x和y值,定义为schx和schy,供下文使用。 重新画张图,前面那张太丑了 dev.off() plot(y~x, ylim=c(0,1.5), xlim=c(0,1), type="l", ylab=expression(beta), xlab=expression(alpha)) 第一行dev.off()告诉R的绘图程序,关闭此前的图。第二行的ylim和xlim定义x和y轴的取值范围。type是告诉plot程序,散点图以line的形式呈现。ylab和xlab是x和y轴的名称:由于我们需要使用希腊字母,因此需要使用expression(alpha)和expression(beta)来打出与。 比原来那张图看着帅一些了是吧。继续调整。 par(new=TRUE) plot(max(y)~x , pch=2, lty=3, col="red",ylim=c(0,1.5), xlim=c(0,1), ylab="", xlab="") 第一行告诉绘图程序,我下面要在原有那张图的基础上,继续添加内容,不要覆盖原图。 第二行中,把y最大值时的点标出来,pch表示这个点的形状(1是空心圆,2是空心三角,等等,你们自己试试看从1到20吧)。col="red",将这个点画成红色。ylab和xlab设为空,ylim和xlim和上文的值相同,使得两张图x轴、y轴重合。 差不多行了。如果你想继续和我一样骚包一点的话,跟我向下继续设置。 text(schx, schy-0.2, substitute(paste("(", xx ,", ", yy, ")"), list(xx=schx,yy=schy))) 这行代码的作用是,在途中加入一段文字。语法text(a,b,c)的意思是,a代表横坐标的位置,b代表纵坐标的位置,c代表所需要加入的语言。a和b加在一起,告诉plot()需要把一段文字c放在(a,b)这个地方。 c需要做进一步的说明。"(", 以及 ")", 表示这其中是有纯文字部分的,plot()直接把它们打出来即可(注意逗号是要保留的)。加在两个小括号中间的是xx和yy,这是两个值。 substitute(paste(),list()): substitute()中包括两部分,一部分是paste(),一部分是list()。paste告诉plot()要加入这段东西了,其中包括两个待赋值的xx和yy。list()为它们分别赋值。 看看效果吧。 看起来还可以。要是能把这个式子也写上去就更完美了。 text(0.4, 0.5, substitute( paste(beta == (frac(alpha, 1-alpha))^{-alpha}, ",", "0", alpha, "1") ) ) 收工。这张图看起来还是有些难看。。。。不过基本信息都已经有了。 (来源:中国统计网)
个人分类: 金融学|0 个评论

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-26 12:35