liamml
mydata <-
mutate(
data,
t = 1:length(date),
dow = wday(date),
PM2.5lag1 = lag(PM2.5, 1),
PM2.5lag2 = lag(PM2.5, 2),
PM2.5lag3 = lag(PM2.5, 3),
PM2.5lag4 = lag(PM2.5, 4),
PM2.5lag5 = lag(PM2.5, 5),
PM2.5lag6 = lag(PM2.5, 6),
PM2.5lag7 = lag(PM2.5, 7),
PM2.5lag8 = lag(PM2.5, 8),
PM2.5lag9 = lag(PM2.5, 9),
PM2.5lag10 = lag(PM2.5, 11),
PM2.5lag12 = lag(PM2.5, 12),
PM2.5lag13 = lag(PM2.5, 13),
PM2.5lag14 = lag(PM2.5, 14),
PM2.5lag15 = lag(PM2.5, 15),
PM2.5lag16 = lag(PM2.5, 16),
PM2.5lag17 = lag(PM2.5, 17),
PM2.5lag18 = lag(PM2.5, 18),
PM2.5lag19 = lag(PM2.5, 19),
PM2.5lag20 = lag(PM2.5, 20),
PM2.5lag21 = lag(PM2.5, 21),
PM2.5lag22 = lag(PM2.5, 22),
PM2.5lag23 = lag(PM2.5, 23),
PM2.5lag24 = lag(PM2.5, 24),
PM2.5lag25 = lag(PM2.5, 25),
PM2.5lag26 = lag(PM2.5, 26),
PM2.5lag27 = lag(PM2.5, 27),
PM2.5lag28 = lag(PM2.5, 28),
PM2.5ave03 = runMean(PM2.5, 0:3),
PM2.5lag12 = runMean(PM2.5, 1:2),
PM2.5ave13 = runMean(PM2.5, 1:3),
PM2.5ave02 = runMean(PM2.5, 0:2),
PM2.5ave01 = runMean(PM2.5, 0:1),
PM2.5ave23 = runMean(PM2.5, 2:3),
PM10lag1 = lag(PM10, 1),
PM10lag2 = lag(PM10, 2),
PM10lag3 = lag(PM10, 3),
PM10lag4 = lag(PM10, 4),
PM10lag5 = lag(PM10, 5),
PM10lag6 = lag(pM10, 6),
PM10lag7 = lag(PM10, 7),
PM10lag8 = lag(PM10, 8),
PM10lag9 = lag(PM10, 9),
PM10lag10 = lag(PM10, 11),
PM10lag12 = lag(PM10, 12),
PM10lag13 = lag(PM10, 13),
PM10lag14 = lag(PM10, 14),
PM10lag15 = lag(PM10, 15),
PM10lag16 = lag(PM10, 16),
PM10lag17 = lag(PM10, 17),
PM10lag18 = lag(PM10, 18),
PM10lag19 = lag(PM10, 19),
PM10lag20 = lag(PM10, 20),
PM10lag21 = lag(PM10, 21),
PM10lag22 = lag(PM10, 22),
PM10lag23 = lag(PM10, 23),
PM10lag24 = lag(PM10, 24),
PM10lag25 = lag(PM10, 25),
PM10lag26 = lag(PM10, 26),
PM10lag27 = lag(PM10, 27),
PM10lag28 = lag(PM10, 28),
PM10ave03 = runMean(PM10, 0:3),
PM10ave12 = runMean(PM10, 1:2),
PM10ave13 = runMean(PM10, 1:3),
PM10ave02 = runMean(PM10, 0:2),
PM10ave01 = runMean(PM10, 0:1),
PM10ave23 = runMean(PM10, 2:3),
NO2lag1 = lag(NO2, 1),
NO2lag2 = lag(NO2, 2),
NO2lag3 = lag(NO2, 3),
NO2lag4 = lag(NO2, 4),
NO2lag5 = lag(NO2, 5),
NO2lag6 = lag(NO2, 6),
NO2lag7 = lag(NO2, 7),
NO2lag8 = lag(NO2, 8),
NO2lag9 = lag(NO2, 9),
NO2lag10 = lag(NO2, 11),
NO2lag12 = lag(NO2, 12),
NO2lag13 = lag(NO2, 13),
NO2lag14 = lag(NO2, 14),
NO2lag15 = lag(NO2, 15),
NO2lag16 = lag(NO2, 16),
NO2lag17 = lag(NO2, 17),
NO2lag18 = lag(NO2, 18),
NO2lag19 = lag(NO2, 19),
NO2lag20 = lag(NO2, 20),
NO2lag21 = lag(NO2, 21),
NO2lag22 = lag(NO2, 22),
NO2lag23 = lag(NO2, 23),
NO2lag24 = lag(NO2, 24),
NO2lag25 = lag(NO2, 25),
NO2lag26 = lag(NO2, 26),
NO2lag27 = lag(NO2, 27),
NO2lag28 = lag(NO2, 28),
NO2ave03 = runMean(NO2, 0:3),
NO2ave12 = runMean(NO2, 1:2),
NO2ave13 = runMean(NO2, 1:3),
NO2ave02 = runMean(NO2, 0:2),
NO2ave01 = runMean(NO2, 0:1),
NO2ave23 = runMean(PM10, 2:3),
SO2lag1 = lag(SO2, 1),
SO2lag2 = lag(SO2, 2),
SO2lag3 = lag(SO2, 3),
SO2lag4 = lag(SO2, 4),
SO2lag5 = lag(SO2, 5),
SO2lag6 = lag(SO2, 6),
SO2lag7 = lag(SO2, 7),
SO2lag8 = lag(SO2, 8),
SO2lag9 = lag(SO2, 9),
SO2lag10 = lag(SO2, 11),
SO2lag12 = lag(SO2, 12),
SO2lag13 = lag(SO2, 13),
SO2lag14 = lag(SO2, 14),
SO2lag15 = lag(SO2, 15),
SO2lag16 = lag(SO2, 16),
SO2lag17 = lag(SO2, 17),
SO2lag18 = lag(SO2, 18),
SO2lag19 = lag(SO2, 19),
SO2lag20 = lag(SO2, 20),
SO2lag21 = lag(SO2, 21),
SO2lag22 = lag(SO2, 22),
SO2lag23 = lag(SO2, 23),
SO2lag24 = lag(SO2, 24),
SO2lag25 = lag(SO2, 25),
SO2lag26 = lag(SO2, 26),
SO2lag27 = lag(SO2, 27),
SO2lag28 = lag(SO2, 28),
SO2ave03 = runMean(SO2, 0:3),
SO2ave12 = runMean(SO2, 1:2),
SO2ave13 = runMean(SO2, 1:3),
SO2ave02 = runMean(SO2, 0:2),
SO2ave01 = runMean(SO2, 0:1),
SO2ave23 = runMean(SO2, 2:3),
COlag1 = lag(CO, 1),
COlag2 = lag(CO, 2),
COlag3 = lag(CO, 3),
COlag4 = lag(CO, 4),
COlag5 = lag(CO, 5),
COlag6 = lag(CO, 6),
COlag7 = lag(CO, 7),
COlag8 = lag(CO, 8),
COlag9 = lag(CO, 9),
COlag10 = lag(CO, 11),
COlag12 = lag(CO, 12),
COlag13 = lag(CO, 13),
COlag14 = lag(CO, 14),
COlag15 = lag(CO, 15),
COlag16 = lag(CO, 16),
COlag17 = lag(CO, 17),
COlag18 = lag(CO, 18),
COlag19 = lag(CO, 19),
COlag20 = lag(CO, 20),
COlag21 = lag(CO, 21),
COlag22 = lag(CO, 22),
COlag23 = lag(CO, 23),
COlag24 = lag(CO, 24),
COlag25 = lag(CO, 25),
COlag26 = lag(CO, 26),
COlag27 = lag(CO, 27),
COlag28 = lag(CO, 28),
COave03 = runMean(CO, 0:3),
COave12 = runMean(CO, 1:2),
COave13 = runMean(CO, 1:3),
COave02 = runMean(CO, 0:2),
COave01 = runMean(CO, 0:1),
COave23 = runMean(CO, 2:3),
O3lag1 = lag(O3, 1),
O3lag2 = lag(O3, 2),
O3lag3 = lag(O3, 3),
O3lag4 = lag(O3, 4),
O3lag5 = lag(O3, 5),
O3lag6 = lag(O3, 6),
O3lag7 = lag(O3, 7),
O3lag8 = lag(O3, 8),
O3lag9 = lag(O3, 9),
O3lag10 = lag(O3, 11),
O3lag12 = lag(O3, 12),
O3lag13 = lag(O3, 13),
O3lag14 = lag(O3, 14),
O3lag15 = lag(O3, 15),
O3lag16 = lag(O3, 16),
O3lag17 = lag(O3, 17),
O3lag18 = lag(O3, 18),
O3lag19 = lag(O3, 19),
O3lag20 = lag(O3, 20),
O3lag21 = lag(O3, 21),
O3lag22 = lag(O3, 22),
O3lag23 = lag(O3, 23),
O3lag24 = lag(O3, 24),
O3lag25 = lag(O3, 25),
O3lag26 = lag(O3, 26),
O3lag27 = lag(O3, 27),
O3lag28 = lag(O3, 28),
O3ave03 = runMean(O3, 0:3),
O3ave12 = runMean(O3, 1:2),
O3ave13 = runMean(O3, 1:3),
O3ave02 = runMean(O3, 0:2),
O3ave01 = runMean(O3, 0:1),
O3ave23 = runMean(O3, 2:3)
))
这就是源代码