[R] roc and lattice
rhk
mailing-lists at rhkoning.com
Wed Jan 10 08:31:54 CET 2007
Hello, I am afraid I do not fully understand all intricacies of
programming in lattice plots. In the code below I try to plot an ROC
curve, following R-news 4(1). When I condition on the variable 'group' I
get the error message below, when I plot the curve for all data (i.e., y
~ pred.prob), I get the plot I want. Can someone point out why
conditioning gives that message? Thanks, Ruud
> plot.a <- xyplot(y ~ pred.prob|group, data=x.df,
+ xlim=c(0,1),xlab="1-specificiteit",
+ ylab="sensitiviteit",
+ panel=function(x,y,subscripts,...){
+ DD <- table(-x,y)
+ sens <- cumsum(DD[,2])/sum(DD[,2])
+ mspec <- cumsum(DD[,1])/sum(DD[,1])
+ panel.xyplot(mspec,sens,type="l",...)
+ panel.abline(0,1)
+ })
> print(plot.a)
Error in panel(x = c(0.000265710002003069, 0.000345712857778025,
0.000265710002003069, :
subscript out of bounds
More information about the R-help
mailing list