1. Importing data:
nfl2016 <- read.csv("C:/Users/Angel/Desktop/nfl2016.csv", stringsAsFactors=FALSE)
View(nfl2016)

2.Sampling weights of 11 players

my.nfl <- sample(nfl2016$WT, size=11, replace=F)
  1. Making a function
oneSummarizedSample <- function(data, sample.size){
  my.sample <- sample(x=data,size=sample.size)
  my.summary <- c(mean=mean(my.sample),
                 var=var(my.sample),
                 se=sd(my.sample)/sqrt(sample.size)) 
  return(my.summary)
}

oneSummarizedSample(data=nfl2016$WT, sample.size=11)
##       mean        var         se 
##  242.81818 2204.96364   14.15808
  1. Running the function 1000 times:
my.sampling.distribution <- replicate(1000, oneSummarizedSample(data=nfl2016$WT, sample.size=11))
  1. Making a histogram
hist(my.sampling.distribution[1, ], main="Sampling Distribution of Weights of NFL Players", xlab="Average Weight")