tập con() của một vector trong R

Question

Tôi đã viết các chức năng sau đây dựa trên subset(), mà tôi tìm thấy tiện dụng:

ss <- function (x, subset, ...) 
{ 
    r <- eval(substitute(subset), data.frame(.=x), parent.frame()) 
    if (!is.logical(r)) 
     stop("'subset' must be logical") 
    x[r & !is.na(r)] 
} 

Vì vậy, tôi có thể viết:

ss(myDataFrame$MyVariableName, 500 < . & . < 1500) 

thay vì

myDataFrame$MyVariableName[ 500 < myDataFrame$MyVariableName 
           & myDataFrame$MyVariableName < 1500] 

Điều này có vẻ như điều gì đó mà những người khác có thể đã phát triển các giải pháp, mặc dù - bao gồm nội dung nào đó trong RI lõi có thể đã bỏ lỡ. Có gì ngoài kia không?

Answer 1

Cảm ơn bạn đã chia sẻ Ken.

Bạn có thể sử dụng:

x <- myDataFrame$MyVariableName; x[x > 100 & x < 180] 

Yours có thể yêu cầu ít gõ nhưng mã là ít khái quát hóa cho người khác nếu bạn đang chia sẻ mã. Tôi có một vài chức năng tiết kiệm thời gian giống như vậy nhưng sử dụng chúng một cách tiết kiệm vì chúng có thể làm chậm mã của bạn (các bước bổ sung) và yêu cầu bạn cũng bao gồm mã đó cho hàm đó khi bạn chia sẻ tệp với người khác.

So sánh độ dài của văn bản. Gần như giống nhau chiều dài:

ss(mtcars$hp, 100 < . & . < 180) 
x <- mtcars$hp; x[x > 100 & x < 180] 

Hãy so sánh thời gian trên 1000 lần nhắc lại.

library(rbenchmark) 
benchmark(
     tyler = x[x > 100 & x < 180], 
     ken = ss(mtcars$hp, 100 <. & . < 180), 
replications=1000) 

    test replications elapsed relative user.self sys.self user.child sys.child 
2 ken   1000 0.56 18.66667  0.36  0.03   NA  NA 
1 tyler   1000 0.03 1.00000  0.03  0.00   NA  NA 

Vì vậy, tôi đoán nó phụ thuộc vào việc bạn cần tốc độ và/hoặc khả năng chia sẻ so với tiện lợi. Nếu nó chỉ dành cho bạn trên một tập dữ liệu nhỏ tôi sẽ nói nó có giá trị.

EDIT: ĐIỂM CHUẨN MỚI

> benchmark(
+  tyler = {x <- mtcars$hp; x[x > 100 & x < 180]}, 
+  ken = ss(mtcars$hp, 100 <. & . < 180), 
+  ken2 = ss2(mtcars$hp, 100 <. & . < 180), 
+  joran = with(mtcars,hp[hp>100 & hp< 180 ]), 
+ replications=10000) 

    test replications elapsed relative user.self sys.self user.child sys.child 
4 joran  10000 0.83 2.677419  0.69  0.00   NA  NA 
2 ken  10000 3.79 12.225806  3.45  0.02   NA  NA 
3 ken2  10000 0.67 2.161290  0.35  0.00   NA  NA 
1 tyler  10000 0.31 1.000000  0.20  0.00   NA  NA 

Answer 2

Tôi nhận ra rằng giải pháp Ken cung cấp là tổng quát hơn hơn là chỉ chọn các mục bên trong dãy (vì nó sẽ làm việc trên bất kỳ biểu thức logic) nhưng điều này đã nhắc nhở tôi rằng Greg tuyết có nhà khai thác ghi vào so sánh trong gói Dạy Demo của mình: (. = x)

library(TeachingDemos) 
x0 <- rnorm(100) 
x0[ 0 %<% x0 %<% 1.5 ]