10.08
코드 오류 수정, 코드 결과값 수정, 결론 수정 완료
is.na
:변수 내에서 결측치 논리문 판단 (T,F).
sum(is.na())
:변수 내 결측치 개수를 알려줌.
is.na가 TRUE인 값(결측치)에 대해 평균값을 지정했음. 결측치를 완전히 제거할 경우 오차 발생 가능성이 크므로, 결측치를 평균값으로 대체하고자 함.
<요약>
결측치가 적은 경우 데이터 결과에 큰 문제가 없었지만 결측치가 많을 수록 raw data와의 결과가 달라짐.
결측치가 많은 자료는 median값이 크게 달라지는 경향을 보임. 그 외 min, max값에는 영향을 주지 못하였고, 1분위수나 3분위수에 조금 영향을 미치는 결과. 예상 외로 mean 값은 눈에 띄게 변동하지 않았음.
결측치가 적은 경우는 mean imputation을 사용해도 될 것으로 판단.
R 코드
w12<-read.csv("C:/Bigdata/water_data/data_2012.csv",header = T)
w13<-read.csv("C:/Bigdata/water_data/data_2013.csv",header = T)
w14<-read.csv("C:/Bigdata/water_data/data_2014.csv",header = T)
w15<-read.csv("C:/Bigdata/water_data/data_2015.csv",header = T)
w16<-read.csv("C:/Bigdata/water_data/data_2016.csv",header = T)
w17<-read.csv("C:/Bigdata/water_data/data_2017.csv",header = T)
w18<-read.csv("C:/Bigdata/water_data/data_2018.csv",header = T)
w19<-read.csv("C:/Bigdata/water_data/data_2019.csv",header = T)
w20<-read.csv("C:/Bigdata/water_data/data_2020.csv",header = T)
w21<-read.csv("C:/Bigdata/water_data/data_2021.csv",header = T)
w22<-read.csv("C:/Bigdata/water_data/data_2022.csv",header = T)
r12<-read.csv("C:/Bigdata/rf_data/rf_2012.csv",header = T)
r13<-read.csv("C:/Bigdata/rf_data/rf_2013.csv",header = T)
r14<-read.csv("C:/Bigdata/rf_data/rf_2014.csv",header = T)
r15<-read.csv("C:/Bigdata/rf_data/rf_2015.csv",header = T)
r16<-read.csv("C:/Bigdata/rf_data/rf_2016.csv",header = T)
r17<-read.csv("C:/Bigdata/rf_data/rf_2017.csv",header = T)
r18<-read.csv("C:/Bigdata/rf_data/rf_2018.csv",header = T)
r19<-read.csv("C:/Bigdata/rf_data/rf_2019.csv",header = T)
r20<-read.csv("C:/Bigdata/rf_data/rf_2020.csv",header = T)
r21<-read.csv("C:/Bigdata/rf_data/rf_2021.csv",header = T)
r22<-read.csv("C:/Bigdata/rf_data/rf_2022.csv",header = T)
#결측치 평균값으로 대체
#2012
w12$ swl12<-w12$ swl
mean(w12$ swl)
mean(w12$ swl,na.rm = TRUE)
w12$ swl12[is.na(w12$ swl12)]=25.04245
summary(w12$ swl)
summary(w12$ swl12)
w12$ inf12<-w12$ inf
mean(w12$ inf)
mean(w12$ inf,na.rm = TRUE)
w12$ inf12[is.na(w12$ inf12)]=674.1987
summary(w12$ inf)
summary(w12$ inf12)
w12$ sfw12<-w12$ sfw
mean(w12$ sfw)
mean(w12$ sfw,na.rm = TRUE)
w12$ sfw12[is.na(w12$ sfw12)]=229.2945
summary(w12$ sfw)
summary(w12$ sfw12)
w12$ ecpc12<-w12$ ecpc
mean(w12$ ecpc)
mean(w12$ ecpc,na.rm = TRUE)
w12$ ecpc12[is.na(w12$ ecpc12)]=14.7066
summary(w12$ ecpc)
summary(w12$ ecpc12)
w12$ tototf12<-w12$ tototf
mean(w12$ tototf)
mean(w12$ tototf,na.rm = T)
w12$ tototf12[is.na(w12$ tototf12)]=675.9133
summary(w12$ tototf)
summary(w12$ tototf12)
w12$ tide_level12<-w12$ tide_level
mean(w12$ tide_level)
mean(w12$ tide_level,na.rm = T)
w12$ tide_level12[is.na(w12$ tide_level12)]=354.2489
summary(w12$ tide_level)
summary(w12$ tide_level12)
str(w12$ tide_level)
str(w12$ tide_level12)
#2013
w13$ tide_level13<-w13$ tide_level
mean(w13$ tide_level)
mean(w13$ tide_level,na.rm = T)
w13$ tide_level13[is.na(w13$ tide_level13)]=352.2339
summary(w13$ tide_level)
summary(w13$ tide_level13)
#2014
w14$ swl14<-w14$ swl
mean(w14$ swl)
mean(w14$ swl,na.rm = TRUE)
w14$ swl14[is.na(w14$ swl14)]=25.15173
summary(w14$ swl)
summary(w14$ swl14)
w14$ inf14<-w14$ inf
mean(w14$ inf)
mean(w14$ inf,na.rm = TRUE)
w14$ inf14[is.na(w14$ inf14)]=230.5956
summary(w14$ inf)
summary(w14$ inf14)
w14$ sfw14<-w14$ sfw
mean(w14$ sfw)
mean(w14$ sfw,na.rm = TRUE)
w14$ sfw14[is.na(w14$ sfw14)]=231.489
summary(w14$ sfw)
summary(w14$ sfw14)
w14$ ecpc14<-w14$ ecpc
mean(w14$ ecpc)
mean(w14$ ecpc,na.rm = TRUE)
w14$ ecpc14[is.na(w14$ ecpc14)]=12.51196
summary(w14$ ecpc)
summary(w14$ ecpc14)
w14$ tototf14<-w14$ tototf
mean(w14$ tototf)
mean(w14$ tototf,na.rm = T)
w14$ tototf14[is.na(w14$ tototf14)]=675.9133
summary(w14$ tototf)
summary(w14$ tototf14)
w14$ ecpc14<-w14$ ecpc
mean(w14$ ecpc)
mean(w14$ ecpc,na.rm = TRUE)
w14$ ecpc14[is.na(w14$ ecpc14)]=12.51196
summary(w14$ ecpc)
summary(w14$ ecpc14)
w14$ tototf14<-w14$ tototf
mean(w14$ tototf)
mean(w14$ tototf,na.rm = T)
w14$ tototf14[is.na(w14$ tototf14)]=228.3617
summary(w14$ tototf)
summary(w14$ tototf14)
w14$ tide_level14<-w14$ tide_level
mean(w14$ tide_level)
mean(w14$ tide_level,na.rm = T)
w14$ tide_level14[is.na(w14$ tide_level14)]=342.6716
summary(w14$ tide_level)
summary(w14$ tide_level14)
w14$ fw_1018680_14<-w14$ fw_1018680
mean(w14$ fw_1018680)
mean(w14$ fw_1018680,na.rm = T)
##null data##
#2015
w15$ swl15<-w15$ swl
mean(w15$ swl)
mean(w15$ swl,na.rm = TRUE)
w15$ swl15[is.na(w15$ swl15)]=25.19713
summary(w15$ swl)
summary(w15$ swl15)
w15$ inf15<-w15$ inf
mean(w15$ inf)
mean(w15$ inf,na.rm = TRUE)
w15$ inf15[is.na(w15$ inf15)]=139.738
summary(w15$ inf)
summary(w15$ inf15)
w15$ sfw15<-w15$ sfw
mean(w15$ sfw)
mean(w15$ sfw,na.rm = TRUE)
w15$ sfw15[is.na(w15$ sfw15)]=228.3777
summary(w15$ sfw)
summary(w15$ sfw15)
w15$ ecpc15<-w15$ ecpc
mean(w15$ ecpc)
mean(w15$ ecpc,na.rm = TRUE)
w15$ ecpc15[is.na(w15$ ecpc15)]=15.62358
summary(w15$ ecpc)
summary(w15$ ecpc15)
w15$ tototf15<-w15$ tototf
mean(w15$ tototf)
mean(w15$ tototf,na.rm = T)
w15$ tototf15[is.na(w15$ tototf15)]=136.775
summary(w15$ tototf)
summary(w15$ tototf15)
w15$ tide_level15<-w15$ tide_level
mean(w15$ tide_level)
mean(w15$ tide_level,na.rm = T)
w15$ tide_level15[is.na(w15$ tide_level15)]=341.752
summary(w15$ tide_level)
summary(w15$ tide_level15)
w15$ fw_1018662_15<-w15$ fw_1018662
mean(w15$ fw_1018662)
mean(w15$ fw_1018662,na.rm = T)
w15$ fw_1018662_15[is.na(w15$ fw_1018662_15)]=477.4466
summary(w15$ fw_1018662)
summary(w15$ fw_1018662_15)
w15$ fw_1018680_15<-w15$ fw_1018680
mean(w15$ fw_1018680)
mean(w15$ fw_1018680,na.rm = T)
##fw680 is null data
w15$ fw_1018683_15<-w15$ fw_1018683
mean(w15$ fw_1018683)
mean(w15$ fw_1018683,na.rm = T)
w15$ fw_1018683_15[is.na(w15$ fw_1018683_15)]=169.327
summary(w15$ fw_1018683)
summary(w15$ fw_1018683_15)
#2016
w16$ swl16<-w16$ swl
mean(w16$ swl)
mean(w16$ swl,na.rm = TRUE)
w16$ swl16[is.na(w16$ swl16)]=25.1527
summary(w16$ swl)
summary(w16$ swl16)
w16$ inf16<-w16$ inf
mean(w16$ inf)
mean(w16$ inf,na.rm = TRUE)
w16$ inf16[is.na(w16$ inf16)]=398.9414
summary(w16$ inf)
summary(w16$ inf16)
w16$ sfw16<-w16$ sfw
mean(w16$ sfw)
mean(w16$ sfw,na.rm = TRUE)
w16$ sfw16[is.na(w16$ sfw16)]=216.9966
summary(w16$ sfw)
summary(w16$ sfw16)
w16$ ecpc16<-w16$ ecpc
mean(w16$ ecpc)
mean(w16$ ecpc,na.rm = TRUE)
w16$ ecpc16[is.na(w16$ ecpc16)]=27.0047
summary(w16$ ecpc)
summary(w16$ ecpc16)
w16$ tototf16<-w16$ tototf
mean(w16$ tototf)
mean(w16$ tototf,na.rm = T)
w16$ tototf16[is.na(w16$ tototf16)]=389.2531
summary(w16$ tototf)
summary(w16$ tototf16)
w16$ tide_level16<-w16$ tide_level
mean(w16$ tide_level)
mean(w16$ tide_level,na.rm = T)
w16$ tide_level16[is.na(w16$ tide_level16)]=344.6495
summary(w16$ tide_level)
summary(w16$ tide_level16)
w16$ fw_1018662_16<-w16$ fw_1018662
mean(w16$ fw_1018662)
mean(w16$ fw_1018662,na.rm = T)
w16$ fw_1018662_16[is.na(w16$ fw_1018662_16)]=675.8494
summary(w16$ fw_1018662)
summary(w16$ fw_1018662_16)
w16$ fw_1018680_16<-w16$ fw_1018680
mean(w16$ fw_1018680)
w16$ fw_1018680
mean(w16$ fw_1018680,na.rm = T)
#fw680 is null data
w16$ fw_1018683_16<-w16$ fw_1018683
mean(w16$ fw_1018683)
mean(w16$ fw_1018683,na.rm = T)
w16$ fw_1018683_16[is.na(w16$ fw_1018683_16)]=429.0785
summary(w16$ fw_1018683)
summary(w16$ fw_1018683_16)
#2017
w17$ swl17<-w17$ swl
mean(w17$ swl)
mean(w17$ swl,na.rm = TRUE)
w17$ swl17[is.na(w17$ swl17)]=25.12655
summary(w17$ swl)
summary(w17$ swl17)
w17$ inf17<-w17$ inf
mean(w17$ inf)
mean(w17$ inf,na.rm = TRUE)
w17$ inf17[is.na(w17$ inf17)]=637.1397
summary(w17$ inf)
summary(w17$ inf17)
w17$ sfw17<-w17$ sfw
mean(w17$ sfw)
mean(w17$ sfw,na.rm = TRUE)
w17$ sfw17[is.na(w17$ sfw17)]=216.0446
summary(w17$ sfw)
summary(w17$ sfw17)
w17$ ecpc17<-w17$ ecpc
mean(w17$ ecpc)
mean(w17$ ecpc,na.rm = TRUE)
w17$ ecpc17[is.na(w17$ ecpc17)]=27.9559
summary(w17$ ecpc)
summary(w17$ ecpc17)
w17$ tototf17<-w17$ tototf
mean(w17$ tototf)
mean(w17$ tototf,na.rm = T)
w17$ tototf17[is.na(w17$ tototf17)]=636.8505
summary(w17$ tototf)
summary(w17$ tototf17)
w17$ tide_level17<-w17$ tide_level
mean(w17$ tide_level)
mean(w17$ tide_level,na.rm = T)
w17$ tide_level17[is.na(w17$ tide_level17)]=339.8704
summary(w17$ tide_level)
summary(w17$ tide_level17)
w17$ wl_1018662_17<-w17$ wl_1018662
mean(w17$ wl_1018662)
mean(w17$ wl_1018662,na.rm = T)
w17$ wl_1018662_17[is.na(w17$ wl_1018662_17)]=336.2544
summary(w17$ wl_1018662)
summary(w17$ wl_1018662_17)
w17$ fw_1018662_17<-w17$ fw_1018662
mean(w17$ fw_1018662)
mean(w17$ fw_1018662,na.rm = T)
w17$ fw_1018662_17[is.na(w17$ fw_1018662_17)]=878.699
summary(w17$ fw_1018662)
summary(w17$ fw_1018662_17)
w17$ wl_1018680_17<-w17$ wl_1018680
mean(w17$ wl_1018680)
mean(w17$ wl_1018680,na.rm = T)
w17$ wl_1018680_17[is.na(w17$ wl_1018680_17)]=315.5513
summary(w17$ wl_1018680)
summary(w17$ wl_1018680_17)
w17$ fw_1018680_17<-w17$ fw_1018680
mean(w17$ fw_1018680)
mean(w17$ fw_1018680,na.rm = T)
#fw680 is null data
w17$ wl_1018683_17<-w17$ wl_1018683
mean(w17$ wl_1018683)
mean(w17$ wl_1018683,na.rm = T)
w17$ wl_1018683_17[is.na(w17$ wl_1018683_17)]=317.8243
summary(w17$ wl_1018683)
summary(w17$ wl_1018683_17)
w17$ fw_1018683_17<-w17$ fw_1018683
mean(w17$ fw_1018683)
mean(w17$ fw_1018683,na.rm = T)
w17$ fw_1018683_17[is.na(w17$ fw_1018683_17)]=641.0775
summary(w17$ fw_1018683)
summary(w17$ fw_1018683_17)
w17$ wl_1019630_17<-w17$ wl_1019630
mean(w17$ wl_1019630)
mean(w17$ wl_1019630,na.rm = T)
w17$ wl_1019630_17[is.na(w17$ wl_1019630_17)]=306.4256
summary(w17$ wl_1019630)
summary(w17$ wl_1019630_17)
w17$ fw_1019630_17<-w17$ fw_1019630
mean(w17$ fw_1019630)
mean(w17$ fw_1019630,na.rm = T)
w17$ fw_1019630_17[is.na(w17$ fw_1019630_17)]=1171.493
summary(w17$ fw_1019630)
summary(w17$ fw_1019630_17)
#2018
w18$ swl18<-w18$ swl
mean(w18$ swl)
mean(w18$ swl,na.rm = TRUE)
w18$ swl18[is.na(w18$ swl18)]=25.04118
summary(w18$ swl)
summary(w18$ swl18)
w18$ inf18<-w18$ inf
mean(w18$ inf)
mean(w18$ inf,na.rm = TRUE)
w18$ inf18[is.na(w18$ inf18)]=631.382
summary(w18$ inf)
summary(w18$ inf18)
w18$ sfw18<-w18$ sfw
mean(w18$ sfw)
mean(w18$ sfw,na.rm = TRUE)
w18$ sfw18[is.na(w18$ sfw18)]=212.9452
summary(w18$ sfw)
summary(w18$ sfw18)
w18$ ecpc18<-w18$ ecpc
mean(w18$ ecpc)
mean(w18$ ecpc,na.rm = TRUE)
w18$ ecpc18[is.na(w18$ ecpc18)]=31.05483
summary(w18$ ecpc)
summary(w18$ ecpc18)
w18$ tototf18<-w18$ tototf
mean(w18$ tototf)
mean(w18$ tototf,na.rm = T)
w18$ tototf18[is.na(w18$ tototf18)]=629.3792
summary(w18$ tototf)
summary(w18$ tototf18)
w18$ tide_level18<-w18$ tide_level
mean(w18$ tide_level)
mean(w18$ tide_level,na.rm = T)
w18$ tide_level18[is.na(w18$ tide_level18)]=347.7946
summary(w18$ tide_level)
summary(w18$ tide_level18)
w18$ fw_1018662_18<-w18$ fw_1018662
mean(w18$ fw_1018662)
mean(w18$ fw_1018662,na.rm = T)
w18$ fw_1018662_18[is.na(w18$ fw_1018662_18)]=838.3288
summary(w18$ fw_1018662)
summary(w18$ fw_1018662_18)
##fw680 is null data
w18$ fw_1018683_18<-w18$ fw_1018683
mean(w18$ fw_1018683)
mean(w18$ fw_1018683,na.rm = T)
w18$ fw_1018683_18[is.na(w18$ fw_1018683_18)]=653.9321
summary(w18$ fw_1018683)
summary(w18$ fw_1018683_18)
#2019
w19$ swl19<-w19$ swl
mean(w19$ swl)
mean(w19$ swl,na.rm = TRUE)
w19$ swl19[is.na(w19$ swl19)]=25.02327
summary(w19$ swl)
summary(w19$ swl19)
w19$ inf19<-w19$ inf
mean(w19$ inf)
mean(w19$ inf,na.rm = TRUE)
w19$ inf19[is.na(w19$ inf19)]=268.5078
summary(w19$ inf)
summary(w19$ inf19)
w19$ sfw19<-w19$ sfw
mean(w19$ sfw)
mean(w19$ sfw,na.rm = TRUE)
w19$ sfw19[is.na(w19$ sfw19)]=212.2494
summary(w19$ sfw)
summary(w19$ sfw19)
w19$ ecpc19<-w19$ ecpc
mean(w19$ ecpc)
mean(w19$ ecpc,na.rm = TRUE)
w19$ ecpc19[is.na(w19$ ecpc19)]=31.75061
summary(w19$ ecpc)
summary(w19$ ecpc19)
w19$ tototf19<-w19$ tototf
mean(w19$ tototf)
mean(w19$ tototf,na.rm = T)
w19$ tototf19[is.na(w19$ tototf19)]=250.5232
summary(w19$ tototf)
summary(w19$ tototf19)
w19$ tide_level19<-w19$ tide_level
mean(w19$ tide_level)
mean(w19$ tide_level,na.rm = T)
w19$ tide_level19[is.na(w19$ tide_level19)]=344.6186
summary(w19$ tide_level)
summary(w19$ tide_level19)
w19$ fw_1018662_19<-w19$ fw_1018662
mean(w19$ fw_1018662)
mean(w19$ fw_1018662,na.rm = T)
w19$ fw_1018662_19[is.na(w19$ fw_1018662_19)]=535.2696
summary(w19$ fw_1018662)
summary(w19$ fw_1018662_19)
##fw680 is null data
#2020
w20$ swl20<-w20$ swl
mean(w20$ swl)
mean(w20$ swl,na.rm = TRUE)
w20$ swl20[is.na(w20$ swl20)]=24.92906
summary(w20$ swl)
summary(w20$ swl20)
w20$ inf20<-w20$ inf
mean(w20$ inf)
mean(w20$ inf,na.rm = TRUE)
w20$ inf20[is.na(w20$ inf20)]=1233.032
summary(w20$ inf)
summary(w20$ inf20)
w20$ sfw20<-w20$ sfw
mean(w20$ sfw)
mean(w20$ sfw,na.rm = TRUE)
w20$ sfw20[is.na(w20$ sfw20)]=208.768
summary(w20$ sfw)
summary(w20$ sfw20)
w20$ ecpc20<-w20$ ecpc
mean(w20$ ecpc)
mean(w20$ ecpc,na.rm = TRUE)
w20$ ecpc20[is.na(w20$ ecpc20)]=35.232
summary(w20$ ecpc)
summary(w20$ ecpc20)
w20$ tototf20<-w20$ tototf
mean(w20$ tototf)
mean(w20$ tototf,na.rm = T)
w20$ tototf20[is.na(w20$ tototf20)]=1223.577
summary(w20$ tototf)
summary(w20$ tototf20)
w20$ tide_level20<-w20$ tide_level
mean(w20$ tide_level)
mean(w20$ tide_level,na.rm = T)
w20$ tide_level20[is.na(w20$ tide_level20)]=353.0349
summary(w20$ tide_level)
summary(w20$ tide_level20)
w20$ wl_1018662_20<-w20$ wl_1018662
mean(w20$ wl_1018662)
mean(w20$ wl_1018662,na.rm = T)
w20$ wl_1018662_20[is.na(w20$ wl_1018662_20)]=366.8673
summary(w20$ wl_1018662)
summary(w20$ wl_1018662_20)
w20$ fw_1018662_20<-w20$ fw_1018662
mean(w20$ fw_1018662)
mean(w20$ fw_1018662,na.rm = T)
w20$ fw_1018662_20[is.na(w20$ fw_1018662_20)]=1013.296
summary(w20$ fw_1018662)
summary(w20$ fw_1018662_20)
w20$ wl_1018680_20<-w20$ wl_1018680
mean(w20$ wl_1018680)
mean(w20$ wl_1018680,na.rm = T)
w20$ wl_1018680_20[is.na(w20$ wl_1018680_20)]=366.5208
summary(w20$ wl_1018680)
summary(w20$ wl_1018680_20)
#fw680 is null data
w20$ wl_1018683_20<-w20$ wl_1018683
mean(w20$ wl_1018683)
mean(w20$ wl_1018683,na.rm = T)
w20$ wl_1018683_20[is.na(w20$ wl_1018683_20)]=348.7886
summary(w20$ wl_1018683)
summary(w20$ wl_1018683_20)
w20$ fw_1018683_20<-w20$ fw_1018683
mean(w20$ fw_1018683)
mean(w20$ fw_1018683,na.rm = T)
w20$ fw_1018683_20[is.na(w20$ fw_1018683_20)]=1398.137
summary(w20$ fw_1018683)
summary(w20$ fw_1018683_20)
w20$ wl_1019630_20<-w20$ wl_1019630
mean(w20$ wl_1019630)
mean(w20$ wl_1019630,na.rm = T)
w20$ wl_1019630_20[is.na(w20$ wl_1019630_20)]=323.1091
summary(w20$ wl_1019630)
summary(w20$ wl_1019630_20)
w20$ fw_1019630_20<-w20$ fw_1019630
mean(w20$ fw_1019630)
mean(w20$ fw_1019630,na.rm = T)
w20$ fw_1019630_20[is.na(w20$ fw_1019630_20)]=1829.523
summary(w20$ fw_1019630)
summary(w20$ fw_1019630_20)
#2021
w21$ swl21<-w21$ swl
mean(w21$ swl)
mean(w21$ swl,na.rm = TRUE)
w21$ swl21[is.na(w21$ swl21)]=24.94015
summary(w21$ swl)
summary(w21$ swl21)
w21$ inf21<-w21$ inf
mean(w21$ inf)
mean(w21$ inf,na.rm = TRUE)
w21$ inf21[is.na(w21$ inf21)]=374.3588
summary(w21$ inf)
summary(w21$ inf21)
w21$ sfw21<-w21$ sfw
mean(w21$ sfw)
mean(w21$ sfw,na.rm = TRUE)
w21$ sfw21[is.na(w21$ sfw21)]=209.2311
summary(w21$ sfw)
summary(w21$ sfw21)
w21$ ecpc21<-w21$ ecpc
mean(w21$ ecpc)
mean(w21$ ecpc,na.rm = TRUE)
w21$ ecpc21[is.na(w21$ ecpc21)]=34.76887
summary(w21$ ecpc)
summary(w21$ ecpc21)
w21$ tototf21<-w21$ tototf
mean(w21$ tototf)
mean(w21$ tototf,na.rm = T)
w21$ tototf21[is.na(w21$ tototf21)]=368.5098
summary(w21$ tototf)
summary(w21$ tototf21)
w21$ tide_level21<-w21 tide_level)
mean(w21$ tide_level,na.rm = T)
w21$ tide_level21[is.na(w21$ tide_level21)]=352.2062
summary(w21$ tide_level)
summary(w21$ tide_level21)
w21$ fw_1018662_21<-w21$ fw_1018662
mean(w21$ fw_1018662)
mean(w21$ fw_1018662,na.rm = T)
w21$ fw_1018662_21[is.na(w21$ fw_1018662_21)]=504.1115
summary(w21$ fw_1018662)
summary(w21$ fw_1018662_21)
##fw680 is null data
w21$ fw_1018683_21<-w21$ fw_1018683
mean(w21$ fw_1018683)
mean(w21$ fw_1018683,na.rm = T)
w21$ fw_1018683_21[is.na(w21$ fw_1018683_21)]=411.2394
summary(w21$ fw_1018683)
summary(w21$ fw_1018683_21)
#2022
w22$ swl22<-w22$ swl
mean(w22$ swl)
mean(w22$ swl,na.rm = TRUE)
w22$ swl22[is.na(w22$ swl22)]=25.05185
summary(w22$ swl)
summary(w22$ swl22)
w22$ inf22<-w22$ inf
mean(w22$ inf)
mean(w22$ inf,na.rm = TRUE)
w22$ inf22[is.na(w22$ inf22)]=543.6389
summary(w22$ inf)
summary(w22$ inf22)
w22$ sfw22<-w22$ sfw
mean(w22$ sfw)
mean(w22$ sfw,na.rm = TRUE)
w22$ sfw22[is.na(w22$ sfw22)]=213.3036
summary(w22$ sfw)
summary(w22$ sfw22)
w22$ ecpc22<-w22$ ecpc
mean(w22$ ecpc)
mean(w22$ ecpc,na.rm = TRUE)
w22$ ecpc22[is.na(w22$ ecpc22)]=30.69637
summary(w22$ ecpc)
summary(w22$ ecpc22)
w22$ tototf22<-w22$ tototf
mean(w22$ tototf)
mean(w22$ tototf,na.rm = T)
w22$ tototf22[is.na(w22$ tototf22)]=492.8109
summary(w22$ tototf)
summary(w22$ tototf22)
w22$ tide_level22<-w22$ tide_level
mean(w22$ tide_level)
mean(w22$ tide_level,na.rm = T)
w22$ tide_level22[is.na(w22$ tide_level22)]=346.2235
summary(w22$ tide_level)
summary(w22$ tide_level22)
w22$ fw_1018662_22<-w22$ fw_1018662
mean(w22$ fw_1018662)
mean(w22$ fw_1018662,na.rm = T)
w22$ fw_1018662_22[is.na(w22$ fw_1018662_22)]=688.7929
summary(w22$ fw_1018662)
summary(w22$ fw_1018662_22)
##fw680 is null data.
