作者FengOrz (啦啦啦)
看板Statistics
标题[程式] R，离散cdf底下面积计算与画图
时间Wed Aug 5 20:02:05 2015
[软体程式类别]: R
[程式问题]: 离散分配cdf底下面积计算与绘图
[软体熟悉度]: 低(1~3个月)
[问题叙述]:
我想把计算离散分配的cdf底下面积+画出来，
原始资料有很多个x及其对应的cdf值，以下为简例
例如说原本
0, x<0
F1(x)= 1/2, 0<=x<100 积分後会变成F2(x)= x/2, 0<=x<=100
1, x>=100
G1(x)= 0, x<50 积分後会变成G2(x)= x-50, x>=50
1, x=50
然後想画出F2(x)与G2(x)
因资料较多，没办法逐笔积分出来，也不会画...
其中一个分配的前六笔资料长这样
x prob cdf
1 -0.0436 0.01136364 0.01136364 (=1/88)
2 -0.0406 0.01136364 0.02272727 (=2/88)
3 -0.0322 0.01136364 0.03409091 (=3/88)
4 -0.0284 0.01136364 0.04545455 (=4/88)
5 -0.0108 0.01136364 0.05681818 (=5/88)
6 -0.0066 0.01136364 0.06818182 (=6/88)
另一个分配的前六笔长这样
x prob cdf
1 -0.1274 0.01136364 0.01136364 (=1/88)
2 -0.1251 0.01136364 0.02272727 (=2/88)
3 -0.1106 0.01136364 0.03409091 (=3/88)
4 -0.1054 0.01136364 0.04545455 (=4/88)
5 -0.0776 0.01136364 0.05681818 (=5/88)
6 -0.0755 0.01136364 0.06818182 (=6/88)
最终目的是想比较直到任何一个x为止，这两个分配cdf曲线下的累积面积，是不是一个永远大於另一个
不过这点我觉得知道要怎麽画以後我可以自己弄出来
先感谢各路大神了m(_ _)m
--
┬┬┬╮
⊙ ⊙ ⊙ ⊙ ╬┼┼╡
▼▼▼▼ ▼▼▼▼ ╬┼┼┤
▲▲▲▲ ▲▲▲▲∕ ┼┼╬┤
∕███﹀ ﹀███ ● ┼┼┼╢
︶ ） ︶ ﹨ φuminz ┼┼╬╢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 39.10.169.60
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1438776136.A.FE6.html
1F：→ celestialgod: F2, G2到底是啥= = 08/05 20:14
2F：→ celestialgod: F1, G1的积分吗 08/05 20:15
4F：→ FengOrz: 对!!!!! 请问怎麽弄的~~ 08/05 20:27
6F：→ celestialgod: 我直接积分，然後把F2，G2直接描上去了... 08/05 20:31
7F：→ celestialgod: 你如果cdf没有closed form只能试试看用数值解了 08/05 20:32
8F：→ celestialgod: 更正，CDF的积分没有closed form 08/05 20:32
9F：→ FengOrz: 感谢c大，不过我真正的资料有很多笔x(和对应的cdf值)， 08/05 20:34
10F：→ FengOrz: 不知道有没有更一般化的作法，因为无法一个一个积分出来 08/05 20:35
11F：→ FengOrz: m(_ _)m 08/05 20:35
※ 编辑: FengOrz (27.244.15.242), 08/05/2015 20:35:57
12F：→ celestialgod: 那你一开始问题就不能这样问了= =... 08/05 20:35
13F：→ FengOrz: 对不起问的方法不好，修正一下本文囧rz 08/05 20:37
※ 编辑: FengOrz (27.244.15.242), 08/05/2015 20:44:32
14F：→ celestialgod: 你只能自己用梯形法去算了吧XD 08/05 20:43
15F：→ FengOrz: 我原本是这样弄，不过因为後面还要继续比较F3(x)和G3(x) 08/05 20:45
16F：→ FengOrz: 〒△〒 08/05 20:45
18F：→ celestialgod: 不确定我的梯形法有没有问题Orz 08/05 20:53
19F：→ FengOrz: 好的，我在好好研究，谢啦！ 08/05 21:14