在R中是否有一个用于制作灵敏度/特异性与截止点图的方案? 以glm对象作为参数的优选方法?
这样的东西就是我要找的东西:
PS:这个剧情有名字吗?
答案 0 :(得分:3)
相反,您可以使用ROC曲线绘制针对1-特异性的灵敏度。
library(Epi)
x <- rnorm( 100 )
z <- rnorm( 100 )
w <- rnorm( 100 )
tigol <- function( x ) 1 - ( 1 + exp( x ) )^(-1)
y <- rbinom( 100, 1, tigol( 0.3 + 3*x + 5*z + 7*w ) )
ROC( form = y ~ x + z , plot="sp" )
ROCR
和epicalc
包产生类似的情节。 latter应该接受glm对象。