作者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/m.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