[程式] R 離散cdf底下面積計算與畫圖

看板Statistics作者 (啦啦啦)時間10年前 (2015/08/05 20:02), 10年前編輯推噓0(0019)
留言19則, 2人參與, 最新討論串1/1
[軟體程式類別]: 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://www.ptt.cc/bbs/Statistics/M.1438776136.A.FE6.html

08/05 20:14, , 1F
F2, G2到底是啥= =
08/05 20:14, 1F

08/05 20:15, , 2F
F1, G1的積分嗎
08/05 20:15, 2F

08/05 20:22, , 3F
08/05 20:22, 3F

08/05 20:27, , 4F
對!!!!! 請問怎麼弄的~~
08/05 20:27, 4F

08/05 20:28, , 5F

08/05 20:31, , 6F
我直接積分,然後把F2,G2直接描上去了...
08/05 20:31, 6F

08/05 20:32, , 7F
你如果cdf沒有closed form只能試試看用數值解了
08/05 20:32, 7F

08/05 20:32, , 8F
更正,CDF的積分沒有closed form
08/05 20:32, 8F

08/05 20:34, , 9F
感謝c大,不過我真正的資料有很多筆x(和對應的cdf值),
08/05 20:34, 9F

08/05 20:35, , 10F
不知道有沒有更一般化的作法,因為無法一個一個積分出來
08/05 20:35, 10F

08/05 20:35, , 11F
m(_ _)m
08/05 20:35, 11F
※ 編輯: FengOrz (27.244.15.242), 08/05/2015 20:35:57

08/05 20:35, , 12F
那你一開始問題就不能這樣問了= =...
08/05 20:35, 12F

08/05 20:37, , 13F
對不起問的方法不好,修正一下本文囧rz
08/05 20:37, 13F
※ 編輯: FengOrz (27.244.15.242), 08/05/2015 20:44:32

08/05 20:43, , 14F
你只能自己用梯形法去算了吧XD
08/05 20:43, 14F

08/05 20:45, , 15F
我原本是這樣弄,不過因為後面還要繼續比較F3(x)和G3(x)
08/05 20:45, 15F

08/05 20:45, , 16F
〒△〒
08/05 20:45, 16F

08/05 20:53, , 17F
只能這樣畫給你 http://pastebin.com/TAW0jyQq
08/05 20:53, 17F

08/05 20:53, , 18F
不確定我的梯形法有沒有問題Orz
08/05 20:53, 18F

08/05 21:14, , 19F
好的,我在好好研究,謝啦!
08/05 21:14, 19F
文章代碼(AID): #1LmVj8_c (Statistics)