Вы можете использовать% в% следующим образом
year <- 2000:2019
n <- round(rnorm(length(year), 0))
name <- c(rep("alex", 5), rep("ben", 10), rep("mary", 5))
babyname <- data.frame(year, n, name)
myname <- "alex"
bn <- babyname[babyname$name %in% myname, ]
bnl
year n name
1 2000 0 alex
2 2001 1 alex
3 2002 -1 alex
4 2003 2 alex
5 2004 -1 alex
Внутри вашей функции и с несколькими значениями в качестве myname это будет
nametrend <- function(myname) {
bn <- babyname[babyname$name %in% myname, ]
bn$name <- droplevels(bn$name)
plot(x = bn$year, y = bn$n, main = myname, col = bn$name)
legend("topleft", legend= levels(bn$name), pch=16, col=unique(bn$name))
}
myname <- c("alex", "mary")
nametrend(myname)