duracionDiasDeCadaFase=15
#Dado el vector fechasCorte de fechas, elegir el valor de fechaCorte más similar a dfLong$Fecha y ponerlo en la nueva columna dfLong$fechaCorte
fecha_mas_cercana <- sapply(dfGPS$Fecha, function(f) {
diffs <- abs(difftime(fechasCorte, f, units = "days"))
fechasCorte[which.min(diffs)]
}) %>% as.Date()
dfGPS$fechaCorteCercano <- fecha_mas_cercana
dfGPS$distCambio=as.numeric(difftime(dfGPS$Fecha,dfGPS$fechaCorteCercano,units="days"))
dfGPS=dfGPS %>%
mutate(distCambio=as.numeric(difftime(dfGPS$Fecha,dfGPS$fechaCorteCercano,units="days"))) %>%
mutate(fase=case_when(distCambio <= -duracionDiasDeCadaFase*3 ~ "1. Mucho antes",
distCambio <= -duracionDiasDeCadaFase ~ "1. Mucho antes",
distCambio <= 0 ~ "2. Antes",
distCambio <= duracionDiasDeCadaFase ~ "3. Después",
TRUE ~ "4. Mucho después"))
variable="Distance"
analisis <- function(df,variable){
modelo=lmer( as.formula(sprintf("%s ~ fase + (1|CODE)",variable)), data = dfGPS)
salida1=modelsummary::modelsummary(modelo,stars = TRUE, statistic = "p.value")
salida2=modelsummary::modelplot(modelo)
salida3=emmeans(modelo,specs = "fase") %>% plot()+coord_flip()+xlab(variable)
return(list(salida3,salida2,salida1))
}
lasVariables=names(dfGPS)[4:10]
for(i in lasVariables){
print(i)
salida=analisis(dfGPS,i)
print(salida[[1]])
print(salida[[2]])
print(salida[[3]])
}
## [1] "Duration"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+-----------+
## | | (1) |
## +======================+===========+
## | (Intercept) | 64.928*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase2. Antes | 2.301* |
## +----------------------+-----------+
## | | (0.037) |
## +----------------------+-----------+
## | fase3. Después | 6.791*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase4. Mucho después | 1.883** |
## +----------------------+-----------+
## | | (0.004) |
## +----------------------+-----------+
## | SD (Intercept CODE) | 1.959 |
## +----------------------+-----------+
## | SD (Observations) | 16.099 |
## +----------------------+-----------+
## | Num.Obs. | 3173 |
## +----------------------+-----------+
## | R2 Marg. | 0.017 |
## +----------------------+-----------+
## | R2 Cond. | 0.031 |
## +----------------------+-----------+
## | AIC | 26666.1 |
## +----------------------+-----------+
## | BIC | 26702.5 |
## +----------------------+-----------+
## | ICC | 0.0 |
## +----------------------+-----------+
## | RMSE | 16.05 |
## +======================+===========+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+===========+
## [1] "Distance"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+-------------+
## | | (1) |
## +======================+=============+
## | (Intercept) | 4598.162*** |
## +----------------------+-------------+
## | | (<0.001) |
## +----------------------+-------------+
## | fase2. Antes | 428.520*** |
## +----------------------+-------------+
## | | (<0.001) |
## +----------------------+-------------+
## | fase3. Después | 280.761*** |
## +----------------------+-------------+
## | | (<0.001) |
## +----------------------+-------------+
## | fase4. Mucho después | 206.007*** |
## +----------------------+-------------+
## | | (<0.001) |
## +----------------------+-------------+
## | SD (Intercept CODE) | 230.671 |
## +----------------------+-------------+
## | SD (Observations) | 1371.927 |
## +----------------------+-------------+
## | Num.Obs. | 3173 |
## +----------------------+-------------+
## | R2 Marg. | 0.009 |
## +----------------------+-------------+
## | R2 Cond. | 0.036 |
## +----------------------+-------------+
## | AIC | 54849.8 |
## +----------------------+-------------+
## | BIC | 54886.2 |
## +----------------------+-------------+
## | ICC | 0.0 |
## +----------------------+-------------+
## | RMSE | 1367.50 |
## +======================+=============+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+=============+
## [1] "Speed_Avg_m_min"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+-----------+
## | | (1) |
## +======================+===========+
## | (Intercept) | 71.342*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase2. Antes | 2.930** |
## +----------------------+-----------+
## | | (0.004) |
## +----------------------+-----------+
## | fase3. Después | -2.837*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase4. Mucho después | 1.088+ |
## +----------------------+-----------+
## | | (0.076) |
## +----------------------+-----------+
## | SD (Intercept CODE) | 2.740 |
## +----------------------+-----------+
## | SD (Observations) | 14.909 |
## +----------------------+-----------+
## | Num.Obs. | 3173 |
## +----------------------+-----------+
## | R2 Marg. | 0.010 |
## +----------------------+-----------+
## | R2 Cond. | 0.042 |
## +----------------------+-----------+
## | AIC | 26192.5 |
## +----------------------+-----------+
## | BIC | 26228.9 |
## +----------------------+-----------+
## | ICC | 0.0 |
## +----------------------+-----------+
## | RMSE | 14.86 |
## +======================+===========+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+===========+
## [1] "Distance_Z4"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+------------+
## | | (1) |
## +======================+============+
## | (Intercept) | 558.484*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase2. Antes | 95.359*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase3. Después | 155.138*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase4. Mucho después | 60.125*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | SD (Intercept CODE) | 77.790 |
## +----------------------+------------+
## | SD (Observations) | 406.935 |
## +----------------------+------------+
## | Num.Obs. | 3173 |
## +----------------------+------------+
## | R2 Marg. | 0.015 |
## +----------------------+------------+
## | R2 Cond. | 0.049 |
## +----------------------+------------+
## | AIC | 47151.5 |
## +----------------------+------------+
## | BIC | 47187.9 |
## +----------------------+------------+
## | ICC | 0.0 |
## +----------------------+------------+
## | RMSE | 405.57 |
## +======================+============+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+============+
## [1] "Distance_Z5_6"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+------------+
## | | (1) |
## +======================+============+
## | (Intercept) | 148.865*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase2. Antes | 31.637*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase3. Después | 16.771* |
## +----------------------+------------+
## | | (0.030) |
## +----------------------+------------+
## | fase4. Mucho después | 0.507 |
## +----------------------+------------+
## | | (0.928) |
## +----------------------+------------+
## | SD (Intercept CODE) | 32.439 |
## +----------------------+------------+
## | SD (Observations) | 135.963 |
## +----------------------+------------+
## | Num.Obs. | 3173 |
## +----------------------+------------+
## | R2 Marg. | 0.005 |
## +----------------------+------------+
## | R2 Cond. | 0.059 |
## +----------------------+------------+
## | AIC | 40211.3 |
## +----------------------+------------+
## | BIC | 40247.6 |
## +----------------------+------------+
## | ICC | 0.1 |
## +----------------------+------------+
## | RMSE | 135.48 |
## +======================+============+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+============+
## [1] "Speed_Max"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+-----------+
## | | (1) |
## +======================+===========+
## | (Intercept) | 26.316*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase2. Antes | 0.680*** |
## +----------------------+-----------+
## | | (<0.001) |
## +----------------------+-----------+
## | fase3. Después | -0.160 |
## +----------------------+-----------+
## | | (0.324) |
## +----------------------+-----------+
## | fase4. Mucho después | 0.151 |
## +----------------------+-----------+
## | | (0.198) |
## +----------------------+-----------+
## | SD (Intercept CODE) | 0.992 |
## +----------------------+-----------+
## | SD (Observations) | 2.852 |
## +----------------------+-----------+
## | Num.Obs. | 3173 |
## +----------------------+-----------+
## | R2 Marg. | 0.005 |
## +----------------------+-----------+
## | R2 Cond. | 0.112 |
## +----------------------+-----------+
## | AIC | 15732.3 |
## +----------------------+-----------+
## | BIC | 15768.7 |
## +----------------------+-----------+
## | ICC | 0.1 |
## +----------------------+-----------+
## | RMSE | 2.84 |
## +======================+===========+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+===========+
## [1] "Ac_dec_2_5"
## Warning: Removed 2 rows containing missing values or values outside the scale range
## (`geom_segment()`).
##
## +----------------------+------------+
## | | (1) |
## +======================+============+
## | (Intercept) | 114.528*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase2. Antes | 13.875*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase3. Después | 8.347*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | fase4. Mucho después | 7.464*** |
## +----------------------+------------+
## | | (<0.001) |
## +----------------------+------------+
## | SD (Intercept CODE) | 18.671 |
## +----------------------+------------+
## | SD (Observations) | 44.248 |
## +----------------------+------------+
## | Num.Obs. | 3173 |
## +----------------------+------------+
## | R2 Marg. | 0.008 |
## +----------------------+------------+
## | R2 Cond. | 0.158 |
## +----------------------+------------+
## | AIC | 33118.6 |
## +----------------------+------------+
## | BIC | 33155.0 |
## +----------------------+------------+
## | ICC | 0.2 |
## +----------------------+------------+
## | RMSE | 44.08 |
## +======================+============+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +======================+============+