作者spiritsky (精灵)
看板Statistics
标题[程式] R中gamma的survival function
时间Sun Oct 3 13:07:08 2010
------------------------------------------------------------------------
[软体程式类别]:
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)
1F：→ arbiteron:pgamma就是CDF 如果三条曲线画在同一张图上的话 10/03 19:45
2F：→ arbiteron:後面的就要用 lines或points的指令 10/03 19:46
3F：推 wasken:你可以用 par(new=T) 也可以 10/03 20:10
4F：→ lin15:add=TRUE 10/03 20:32
5F：→ arbiteron:add=TRUE 似乎只适用在curve(...) 10/04 01:23
6F：→ lin15:plot也可以 10/04 07:47
7F：→ arbiteron:thx 10/04 21:39
8F：→ lin15:後来用好像不行...印象中有在plot下用过才对=.=||| 10/04 21:52