dat <- as.data.frame(read.csv("https://raw.githubusercontent.com/tmatis12/datafiles/main/normtemp.csv"))                   
males <- dat[1:65,3]
females=dat[1:65,3]
summary(males)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   58.00   70.00   73.00   73.37   78.00   86.00
sd(males)
## [1] 5.875184
hist(males,main="Histogram of males beats",col="light blue",xlab="males heart beat")

qqnorm(males,main="normal probability plot of males beats",col="green",xlab="theoritical quantity")

summary(females)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   58.00   70.00   73.00   73.37   78.00   86.00
sd(females)
## [1] 5.875184
hist(females,main="history of females beats",col="light yellow",xlab="females heart beat")

qqnorm(females,main="normal probability of females beats",col="red",xlab="theoritical quantity")

boxplot(males,females,main="box plot of males/females heart beats",col="pink",ylab="males/females heart beats",names=c("males","females"))

summary(cars)
##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00

Including Plots

You can also embed plots, for example:

plot(pressure)

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.