[程式] R中gamma的survival function
------------------------------------------------------------------------
[軟體程式類別]:
R
[程式問題]:
在R中跑gamma function的survival function
[軟體熟悉度]:
低(1~3個月)
[問題敘述]:
想要把gamma在survival function 圖畫出來
可是因為它的survival function 是1-cdf
然後gamma的cdf沒辦法直接積分出來,所以無法只帶數字進去算
是要做成三組數字帶進去看,曲線是怎麼跑?
可是一直都沒辦法將它們會合在一起
[程式範例]:
myq=seq(0.005,0.995,by=0.005)
q1=qgamma(myq,3,1)
dp1=pgamma(q1,3,1)
S1=1-dp1
q2=qgamma(myq,5,1)
dp2=pgamma(q2,5,1)
S2=1-dp2
Sa=cbind(S1,S2)
plot(S1~q1,type="l",main="Gamma cdf Plot ", ylab="cdf", xlab="x")
windows()
plot(S2~q2,type="l",main="Gamma cdf Plot ", ylab="cdf", xlab="x")
-----------------------------------------------------------------------------
這個是目前同學在幫我修改過的程式碼,會跑出兩張圖
可是沒辦法像課本一樣
雖然只是存活分析最前半段,可是已經寫了一天還有別人的幫忙
但還是畫不出來
希望有這方面的高手可以指點一下
因為還有後面的hazard rate=f(x)/s(x)
可是我已經卡在這邊好久了
拜託大家了!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.41.162.205
※ 編輯: spiritsky 來自: 114.41.162.205 (10/03 13:07)
→
10/03 19:45, , 1F
10/03 19:45, 1F
→
10/03 19:46, , 2F
10/03 19:46, 2F
推
10/03 20:10, , 3F
10/03 20:10, 3F
→
10/03 20:32, , 4F
10/03 20:32, 4F
→
10/04 01:23, , 5F
10/04 01:23, 5F
→
10/04 07:47, , 6F
10/04 07:47, 6F
→
10/04 21:39, , 7F
10/04 21:39, 7F
→
10/04 21:52, , 8F
10/04 21:52, 8F