Subscript out of bounds error - r

Question

I have an issue when running this script in Rstudio where I get the error
Error in TMX03_site[, 24] : subscript out of bounds
I have 5 scripts in total like this two of which work fine and three come up with the same error.
I have looked through other post related to this issue and have tried the suggestions but have found no solution to my problem yet. I hope someone out there could trouble shoot this for me and find out where I am going wrong.
Many thanks
setwd("G:/All ArcGIS/tmax03/03")
m03_1 <- raster("1320_setnull")
m03_2 <- raster("2320_setnull")
m03_3 <- raster("3320_setnull")
m03_4 <- raster("4320_setnull")
m03_5 <- raster("5320_setnull")
m03_6 <- raster("6320_setnull")
m03_7 <- raster("7320_setnull")
m03_8 <- raster("8320_setnull")
m03_9 <- raster("9320_setnull")
m03_10 <- raster("10320_setnull")
m03_11 <- raster("11320_setnull")
m03_12 <- raster("12320_setnull")
m03_13 <- raster("1380_setnull")
m03_14 <- raster("2380_setnull")
m03_15 <- raster("3380_setnull")
m03_16 <- raster("4380_setnull")
m03_17 <- raster("5380_setnull")
m03_18 <- raster("6380_setnull")
m03_19 <- raster("7380_setnull")
m03_20 <- raster("8380_setnull")
m03_21 <- raster("9380_setnull")
m03_22 <- raster("10380_setnull")
m03_23 <- raster("11380_setnull")
m03_24 <- raster("12380_setnull")
m03.stacked<-stack(m03_1,m03_2,m03_3,m03_4,m03_5,m03_6,m03_7,m03_8,m03_10,
m03_11,m03_12,m03_13,m03_14,m03_15,m03_16,m03_17,m03_18,m03_19,m03_20,
m03_21,m03_22,m03_23,m03_24)
writeRaster(m03.stacked, "m03.stacked.nc",overwrite=TRUE)
write03 <- 'm03.stacked.nc' # link to nc
site_03 <- read.csv("C:/Users/Vaughn/Desktop/suteri/tmax/sites_03.csv")
for (x in 1:nrow(site_03)) {
longlat <- c(site_03[x,1],site_03[x,2])
}
#convert to a matrix and transpose by
site_03 <- t(matrix(longlat))
TMX03 <- brick(write03)
TMX03_site <- extract(TMX03, site_03)
TMX03_site1 <- TMX03_site[, 1]
TMX03_site2 <- TMX03_site[, 2]
TMX03_site3 <- TMX03_site[, 3]
TMX03_site4 <- TMX03_site[, 4]
TMX03_site5 <- TMX03_site[, 5]
TMX03_site6 <- TMX03_site[, 6]
TMX03_site7 <- TMX03_site[, 7]
TMX03_site8 <- TMX03_site[, 8]
TMX03_site9 <- TMX03_site[, 9]
TMX03_site10 <- TMX03_site[, 10]
TMX03_site11 <- TMX03_site[, 11]
TMX03_site12 <- TMX03_site[, 12]
TMX03_site13 <- TMX03_site[, 13]
TMX03_site14 <- TMX03_site[, 14]
TMX03_site15 <- TMX03_site[, 15]
TMX03_site16 <- TMX03_site[, 16]
TMX03_site17 <- TMX03_site[, 17]
TMX03_site18 <- TMX03_site[, 18]
TMX03_site19 <- TMX03_site[, 19]
TMX03_site20 <- TMX03_site[, 20]
TMX03_site21 <- TMX03_site[, 21]
TMX03_site22 <- TMX03_site[, 22]
TMX03_site23 <- TMX03_site[, 23]
TMX03_site24 <- TMX03_site[, 24]
> str(TMX03)
Formal class 'RasterBrick' [package "raster"] with 12 slots
..# file :Formal class '.RasterFile' [package "raster"] with 13 slots
.. .. ..# name : chr "G:\\All ArcGIS\\tmax03\\03\\m03.stacked.nc"
.. .. ..# datanotation: chr "FLT4S"
.. .. ..# byteorder : chr "little"
.. .. ..# nodatavalue : num -3.4e+38
.. .. ..# NAchanged : logi FALSE
.. .. ..# nbands : int 23
.. .. ..# bandorder : chr "BIL"
.. .. ..# offset : int 0
.. .. ..# toptobottom : logi FALSE
.. .. ..# blockrows : int 0
.. .. ..# blockcols : int 0
.. .. ..# driver : chr "netcdf"
.. .. ..# open : logi FALSE
..# data :Formal class '.MultipleRasterData' [package "raster"] with 14
.. .. ..# values : logi[0 , 0 ] slots
.. .. ..# offset : num 0
.. .. ..# gain : num 1
.. .. ..# inmemory : logi FALSE
.. .. ..# fromdisk : logi TRUE
.. .. ..# nlayers : int 23
.. .. ..# dropped : NULL
.. .. ..# isfactor : logi FALSE
.. .. ..# attributes: list()
.. .. ..# haveminmax: logi FALSE
.. .. ..# min : num [1:23] Inf Inf Inf Inf Inf ...
.. .. ..# max : num [1:23] -Inf -Inf -Inf -Inf -Inf ...
.. .. ..# unit : chr ""
.. .. ..# names : chr [1:23] "X1" "X2" "X3" "X4" ...
..# legend :Formal class '.RasterLegend' [package "raster"] with 5 slots
.. .. ..# type : chr(0)
.. .. ..# values : logi(0)
.. .. ..# color : logi(0)
.. .. ..# names : logi(0)
.. .. ..# colortable: logi(0)
..# title : chr "variable"
..# extent :Formal class 'Extent' [package "raster"] with 4 slots
.. .. ..# xmin: num 1903933
.. .. ..# xmax: num 3021933
.. .. ..# ymin: num 5305094
.. .. ..# ymax: num 6757094
..# rotated : logi FALSE
..# rotation:Formal class '.Rotation' [package "raster"] with 2 slots
.. .. ..# geotrans: num(0)
.. .. ..# transfun:function ()
..# ncols : int 1118
..# nrows : int 1452
..# crs :Formal class 'CRS' [package "sp"] with 1 slot
.. .. ..# projargs: chr NA
..# history : list()
..# z :List of 1
.. ..$ unknown: int [1:23(1d)] 1 2 3 4 5 6 7 8 9 10 ...
packages that are loaded are raster and sp