data(cells)
   K <- Kest(cells)
   # K is an object of class "fv"
   plot(K, iso ~ r)                # plots iso against r
   plot(K, sqrt(iso/pi) ~ r)   # plots sqrt(iso/r)  against r
   plot(K, cbind(iso,theo) ~ r)   # plots iso against r  AND theo against r
   plot(K, .  ~ r)            # plots all available estimates of K against r
   plot(K, sqrt(./pi) ~ r)   # plots all estimates of L-function
                             # L(r) = sqrt(K(r)/pi)
   plot(K, cbind(iso,theo) ~ r, col=c(2,3))
                                   # plots iso against r  in colour 2
                                   # and theo against r in colour 3
   plot(K, iso ~ r, subset=quote(r < 0.2))
                                   # plots iso against r for r < 10
   # making a legend by hand
   v <- plot(K, . ~ r, ylab="K(r)", legend=FALSE)
   legend(0.05, 0.15, legend=v$meaning, lty=v$lty, col=v$col)Run the code above in your browser using DataLab