Using team data for the season of 2016, I tested the correlaion b/t _Hit- and Run Scored.

library(Lahman)
dat<-subset(Teams,yearID==2016,select=c(R,H))
cor(dat$H,dat$R)
## [1] 0.6948651

I found the a positive relationship b/t the two variavles. I plotted the relationship

plot(dat$H,dat$R)