У меня есть следующий набор данных
subject stim TRBV
<chr> <chr> <chr>
1 HAT-19 (2) 1100-2 TRBV20-1
2 DL (HC2) (2) 1100-2 TRBV6-1
3 MB (HC1) (3) BSV18 TRBV20-1
4 HAT-19 (2) 1100-2 TRBV7-6
5 HAT-001 (2) 1100-2 TRBV15
6 HAT-001 (3) BSV18 TRBV6-2
7 HAT-19 (2) 1100-2 TRBV6-4
8 HAT-001 (3) BSV18 TRBV20-1
9 MB (HC1) (2) 1100-2 TRBV20-1
10 HAT-001 (2) 1100-2 TRBV6-4
Я хочу узнать, сколько раз я получаю определенное значение "TRBV" для определенных "Subject" и "Stim"
Например,, если Subject = "HAT-19" и activ = "(2) 1100-2" и TRBV = "TRBV20-1", я сделал следующее:
x <- my_data[which (my_data$subject == "HAT-19" & my_data$stim == "(2) 1100-2" & my_data$TRBV == "TRBV20-1"),]
y<-x$TRBV
z<-length(y)
Это работает.Но становится все более монотонным делать все эти шаги для всех предметов.Как я могу сократить шаги?