#import data

dat=read.csv("D:/CAO HOC/HP2/THONG KE MOI TRUONG/climatedat.csv")
str(dat)
## 'data.frame':    126 obs. of  9 variables:
##  $ STATION: int  300114 301449 302343 302369 302789 304155 304642 305076 305773 306733 ...
##  $ ECOSYS : chr  "G-NF" "G-DMG" "G-MG" "G-DMG" ...
##  $ BIOME  : chr  "Grassland" "Grassland" "Grassland" "Grassland" ...
##  $ MAT    : num  2.4 4.5 4.9 5.1 2.8 -0.9 -0.5 2.8 3.3 -3.5 ...
##  $ MWMT   : num  17.5 18.8 18.4 18.3 17.4 17.4 16.3 15.1 17.3 16.5 ...
##  $ MCMT   : num  -24.9 -23.1 -23.3 -24.3 -22.2 -30.4 -29 -21.3 -24.5 -32.3 ...
##  $ MAP    : int  443 415 429 405 431 480 511 550 488 431 ...
##  $ MSP    : int  287 257 258 254 291 292 330 319 309 260 ...
##  $ DRYNESS: int  17 3 12 1 14 45 52 42 22 60 ...
attach(dat)
boxplot(MAT~BIOME)

## R Markdown
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:

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