rm(list = ls())
#Load in libraries
library(Matrix)
library(readr)
library(lme4)
library(car)
## Loading required package: carData
## Registered S3 methods overwritten by 'car':
## method from
## influence.merMod lme4
## cooks.distance.influence.merMod lme4
## dfbeta.influence.merMod lme4
## dfbetas.influence.merMod lme4
library(emmeans)
library(sjstats)
library(MuMIn)
library(lmerTest)
##
## Attaching package: 'lmerTest'
## The following object is masked from 'package:lme4':
##
## lmer
## The following object is masked from 'package:stats':
##
## step
setwd("c:/users/Paul/Documents/Rwork")
scatted_soilsdata<- read.csv(file="scatted_soilsdata.csv")
scatted_soilsdata$ScatF<-as.factor(scatted_soilsdata$Scat)
scatted_soilsdata$TimeF<-as.factor(scatted_soilsdata$Time)
scatted_soilsdata$HoneySF<-as.factor(scatted_soilsdata$Honeysuckle)
str(scatted_soilsdata)
## 'data.frame': 96 obs. of 13 variables:
## $ Subject : int 1 1 2 2 1 1 2 2 1 1 ...
## $ Time : int 1 1 1 1 2 2 2 2 3 3 ...
## $ Scat : chr "P" "P" "A" "A" ...
## $ Honeysuckle : chr "H" "NH" "H" "NH" ...
## $ TotalN : num 2.19 1.8 3.26 3.36 2.11 1.58 4.23 4.71 3.99 3.03 ...
## $ TotalNO3 : num 1.05 0.66 2.1 2.24 1.66 1.12 3.75 4.29 3.49 2.5 ...
## $ Ncummulative : num 2.19 3.99 3.26 6.61 6.09 ...
## $ Nitratecummulative: num 1.05 1.71 2.1 4.34 3.37 ...
## $ Nmin : num 2.19 3.99 3.26 6.61 3.05 3.84 5.43 7.78 3.89 4.9 ...
## $ NNitr : num 1.05 1.71 2.1 4.34 1.69 2.25 4.05 6.19 2.66 3.49 ...
## $ ScatF : Factor w/ 2 levels "A","P": 2 2 1 1 2 2 1 1 2 2 ...
## $ TimeF : Factor w/ 24 levels "1","2","3","4",..: 1 1 1 1 2 2 2 2 3 3 ...
## $ HoneySF : Factor w/ 2 levels "H","NH": 1 2 1 2 1 2 1 2 1 2 ...
summary(scatted_soilsdata)
## Subject Time Scat Honeysuckle
## Min. :1.0 Min. : 1.00 Length:96 Length:96
## 1st Qu.:1.0 1st Qu.: 6.75 Class :character Class :character
## Median :1.5 Median :12.50 Mode :character Mode :character
## Mean :1.5 Mean :12.50
## 3rd Qu.:2.0 3rd Qu.:18.25
## Max. :2.0 Max. :24.00
##
## TotalN TotalNO3 Ncummulative Nitratecummulative
## Min. : 1.580 Min. : 0.660 Min. : 2.19 Min. : 1.05
## 1st Qu.: 9.783 1st Qu.: 9.225 1st Qu.: 79.67 1st Qu.: 71.54
## Median :17.565 Median :16.895 Median : 273.51 Median : 259.44
## Mean :19.142 Mean :18.459 Mean : 349.98 Mean : 334.52
## 3rd Qu.:30.142 3rd Qu.:29.595 3rd Qu.: 564.89 3rd Qu.: 540.06
## Max. :38.640 Max. :37.720 Max. :1129.13 Max. :1096.02
##
## Nmin NNitr ScatF TimeF HoneySF
## Min. : 2.19 Min. : 1.05 A:48 1 : 4 H :48
## 1st Qu.:12.98 1st Qu.:11.76 P:48 2 : 4 NH:48
## Median :21.70 Median :20.57 3 : 4
## Mean :21.91 Mean :20.66 4 : 4
## 3rd Qu.:28.47 3rd Qu.:27.14 5 : 4
## Max. :47.05 Max. :45.67 6 : 4
## (Other):72
scatmodel1<-lmer(Nmin~ScatF*TimeF+(1|Subject),data=scatted_soilsdata)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## unable to evaluate scaled gradient
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge: degenerate Hessian with 1 negative eigenvalues
## Warning in as_lmerModLT(model, devfun): Model may not have converged with 1
## eigenvalue close to zero: 5.3e-09
anova(scatmodel1)
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## ScatF 41.6 41.57 1 48 25.014 8.023e-06 ***
## TimeF 11715.9 509.39 23 48 306.497 < 2.2e-16 ***
## ScatF:TimeF 1568.0 68.18 23 48 41.021 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eta_sq(scatmodel1,partial=TRUE)
## Warning in as_lmerModLT(model, devfun, tol = tol): Model may not have converged
## with 1 eigenvalue close to zero: 5.3e-09
## term partial.etasq
## 1 ScatF 0.343
## 2 TimeF 0.993
## 3 ScatF:TimeF 0.952
emmeans(scatmodel1, list(pairwise~ScatF), adjust="bonferroni")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of ScatF`
## ScatF emmean SE df lower.CL upper.CL
## A 19.4 0.72 10775 18 20.8
## P 24.5 0.72 10775 23 25.9
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF`
## 1 estimate SE df t.ratio p.value
## A - P -5.09 1.02 10775 -5.001 <.0001
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
emmeans(scatmodel1, list(pairwise~ScatF), adjust="tukey")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of ScatF`
## ScatF emmean SE df lower.CL upper.CL
## A 19.4 0.72 10775 18 20.8
## P 24.5 0.72 10775 23 25.9
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF`
## 1 estimate SE df t.ratio p.value
## A - P -5.09 1.02 10775 -5.001 <.0001
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
emmeans(scatmodel1, list(pairwise~TimeF), adjust="bonferroni")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of TimeF`
## TimeF emmean SE df lower.CL upper.CL
## 1 4.01 0.811 120 2.41 5.62
## 2 5.03 0.811 120 3.42 6.63
## 3 6.60 0.811 120 4.99 8.21
## 4 8.52 0.811 120 6.91 10.13
## 5 10.08 0.811 120 8.47 11.68
## 6 11.37 0.811 120 9.76 12.97
## 7 12.49 0.811 120 10.89 14.10
## 8 13.53 0.811 120 11.92 15.14
## 9 14.97 0.811 120 13.37 16.58
## 10 16.82 0.811 120 15.22 18.43
## 11 18.97 0.811 120 17.37 20.58
## 12 21.34 0.811 120 19.74 22.95
## 13 23.57 0.811 120 21.97 25.18
## 14 25.60 0.811 120 24.00 27.21
## 15 27.48 0.811 120 25.87 29.08
## 16 29.22 0.811 120 27.61 30.82
## 17 30.75 0.811 120 29.15 32.36
## 18 32.09 0.811 120 30.49 33.70
## 19 33.26 0.811 120 31.65 34.87
## 20 34.28 0.811 120 32.67 35.89
## 21 35.22 0.811 120 33.61 36.83
## 22 36.10 0.811 120 34.49 37.70
## 23 36.91 0.811 120 35.31 38.52
## 24 37.68 0.811 120 36.07 39.29
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of TimeF`
## 1 estimate SE df t.ratio p.value
## 1 - 2 -1.012 0.912 48 -1.111 1.0000
## 1 - 3 -2.587 0.912 48 -2.838 1.0000
## 1 - 4 -4.508 0.912 48 -4.945 0.0027
## 1 - 5 -6.065 0.912 48 -6.653 <.0001
## 1 - 6 -7.353 0.912 48 -8.066 <.0001
## 1 - 7 -8.482 0.912 48 -9.305 <.0001
## 1 - 8 -9.518 0.912 48 -10.441 <.0001
## 1 - 9 -10.960 0.912 48 -12.023 <.0001
## 1 - 10 -12.812 0.912 48 -14.055 <.0001
## 1 - 11 -14.960 0.912 48 -16.411 <.0001
## 1 - 12 -17.330 0.912 48 -19.011 <.0001
## 1 - 13 -19.560 0.912 48 -21.457 <.0001
## 1 - 14 -21.590 0.912 48 -23.684 <.0001
## 1 - 15 -23.462 0.912 48 -25.738 <.0001
## 1 - 16 -25.205 0.912 48 -27.650 <.0001
## 1 - 17 -26.742 0.912 48 -29.336 <.0001
## 1 - 18 -28.082 0.912 48 -30.806 <.0001
## 1 - 19 -29.247 0.912 48 -32.084 <.0001
## 1 - 20 -30.267 0.912 48 -33.203 <.0001
## 1 - 21 -31.207 0.912 48 -34.234 <.0001
## 1 - 22 -32.085 0.912 48 -35.197 <.0001
## 1 - 23 -32.903 0.912 48 -36.094 <.0001
## 1 - 24 -33.667 0.912 48 -36.933 <.0001
## 2 - 3 -1.575 0.912 48 -1.728 1.0000
## 2 - 4 -3.495 0.912 48 -3.834 0.1013
## 2 - 5 -5.053 0.912 48 -5.543 0.0003
## 2 - 6 -6.340 0.912 48 -6.955 <.0001
## 2 - 7 -7.470 0.912 48 -8.195 <.0001
## 2 - 8 -8.505 0.912 48 -9.330 <.0001
## 2 - 9 -9.947 0.912 48 -10.912 <.0001
## 2 - 10 -11.800 0.912 48 -12.945 <.0001
## 2 - 11 -13.947 0.912 48 -15.300 <.0001
## 2 - 12 -16.317 0.912 48 -17.900 <.0001
## 2 - 13 -18.547 0.912 48 -20.346 <.0001
## 2 - 14 -20.578 0.912 48 -22.573 <.0001
## 2 - 15 -22.450 0.912 48 -24.628 <.0001
## 2 - 16 -24.192 0.912 48 -26.539 <.0001
## 2 - 17 -25.730 0.912 48 -28.226 <.0001
## 2 - 18 -27.070 0.912 48 -29.696 <.0001
## 2 - 19 -28.235 0.912 48 -30.974 <.0001
## 2 - 20 -29.255 0.912 48 -32.093 <.0001
## 2 - 21 -30.195 0.912 48 -33.124 <.0001
## 2 - 22 -31.073 0.912 48 -34.086 <.0001
## 2 - 23 -31.890 0.912 48 -34.983 <.0001
## 2 - 24 -32.655 0.912 48 -35.822 <.0001
## 3 - 4 -1.920 0.912 48 -2.106 1.0000
## 3 - 5 -3.478 0.912 48 -3.815 0.1075
## 3 - 6 -4.765 0.912 48 -5.227 0.0010
## 3 - 7 -5.895 0.912 48 -6.467 <.0001
## 3 - 8 -6.930 0.912 48 -7.602 <.0001
## 3 - 9 -8.373 0.912 48 -9.185 <.0001
## 3 - 10 -10.225 0.912 48 -11.217 <.0001
## 3 - 11 -12.373 0.912 48 -13.573 <.0001
## 3 - 12 -14.742 0.912 48 -16.172 <.0001
## 3 - 13 -16.973 0.912 48 -18.619 <.0001
## 3 - 14 -19.003 0.912 48 -20.846 <.0001
## 3 - 15 -20.875 0.912 48 -22.900 <.0001
## 3 - 16 -22.617 0.912 48 -24.811 <.0001
## 3 - 17 -24.155 0.912 48 -26.498 <.0001
## 3 - 18 -25.495 0.912 48 -27.968 <.0001
## 3 - 19 -26.660 0.912 48 -29.246 <.0001
## 3 - 20 -27.680 0.912 48 -30.365 <.0001
## 3 - 21 -28.620 0.912 48 -31.396 <.0001
## 3 - 22 -29.497 0.912 48 -32.359 <.0001
## 3 - 23 -30.315 0.912 48 -33.255 <.0001
## 3 - 24 -31.080 0.912 48 -34.095 <.0001
## 4 - 5 -1.558 0.912 48 -1.709 1.0000
## 4 - 6 -2.845 0.912 48 -3.121 0.8416
## 4 - 7 -3.975 0.912 48 -4.361 0.0189
## 4 - 8 -5.010 0.912 48 -5.496 0.0004
## 4 - 9 -6.452 0.912 48 -7.078 <.0001
## 4 - 10 -8.305 0.912 48 -9.111 <.0001
## 4 - 11 -10.453 0.912 48 -11.466 <.0001
## 4 - 12 -12.822 0.912 48 -14.066 <.0001
## 4 - 13 -15.053 0.912 48 -16.512 <.0001
## 4 - 14 -17.082 0.912 48 -18.739 <.0001
## 4 - 15 -18.955 0.912 48 -20.794 <.0001
## 4 - 16 -20.698 0.912 48 -22.705 <.0001
## 4 - 17 -22.235 0.912 48 -24.392 <.0001
## 4 - 18 -23.575 0.912 48 -25.862 <.0001
## 4 - 19 -24.740 0.912 48 -27.140 <.0001
## 4 - 20 -25.760 0.912 48 -28.259 <.0001
## 4 - 21 -26.700 0.912 48 -29.290 <.0001
## 4 - 22 -27.578 0.912 48 -30.252 <.0001
## 4 - 23 -28.395 0.912 48 -31.149 <.0001
## 4 - 24 -29.160 0.912 48 -31.988 <.0001
## 5 - 6 -1.288 0.912 48 -1.412 1.0000
## 5 - 7 -2.417 0.912 48 -2.652 1.0000
## 5 - 8 -3.453 0.912 48 -3.787 0.1170
## 5 - 9 -4.895 0.912 48 -5.370 0.0006
## 5 - 10 -6.747 0.912 48 -7.402 <.0001
## 5 - 11 -8.895 0.912 48 -9.758 <.0001
## 5 - 12 -11.265 0.912 48 -12.358 <.0001
## 5 - 13 -13.495 0.912 48 -14.804 <.0001
## 5 - 14 -15.525 0.912 48 -17.031 <.0001
## 5 - 15 -17.398 0.912 48 -19.085 <.0001
## 5 - 16 -19.140 0.912 48 -20.996 <.0001
## 5 - 17 -20.677 0.912 48 -22.683 <.0001
## 5 - 18 -22.017 0.912 48 -24.153 <.0001
## 5 - 19 -23.183 0.912 48 -25.431 <.0001
## 5 - 20 -24.203 0.912 48 -26.550 <.0001
## 5 - 21 -25.142 0.912 48 -27.581 <.0001
## 5 - 22 -26.020 0.912 48 -28.544 <.0001
## 5 - 23 -26.837 0.912 48 -29.441 <.0001
## 5 - 24 -27.602 0.912 48 -30.280 <.0001
## 6 - 7 -1.130 0.912 48 -1.240 1.0000
## 6 - 8 -2.165 0.912 48 -2.375 1.0000
## 6 - 9 -3.607 0.912 48 -3.957 0.0688
## 6 - 10 -5.460 0.912 48 -5.990 0.0001
## 6 - 11 -7.607 0.912 48 -8.345 <.0001
## 6 - 12 -9.977 0.912 48 -10.945 <.0001
## 6 - 13 -12.207 0.912 48 -13.392 <.0001
## 6 - 14 -14.238 0.912 48 -15.618 <.0001
## 6 - 15 -16.110 0.912 48 -17.673 <.0001
## 6 - 16 -17.852 0.912 48 -19.584 <.0001
## 6 - 17 -19.390 0.912 48 -21.271 <.0001
## 6 - 18 -20.730 0.912 48 -22.741 <.0001
## 6 - 19 -21.895 0.912 48 -24.019 <.0001
## 6 - 20 -22.915 0.912 48 -25.138 <.0001
## 6 - 21 -23.855 0.912 48 -26.169 <.0001
## 6 - 22 -24.733 0.912 48 -27.131 <.0001
## 6 - 23 -25.550 0.912 48 -28.028 <.0001
## 6 - 24 -26.315 0.912 48 -28.867 <.0001
## 7 - 8 -1.035 0.912 48 -1.135 1.0000
## 7 - 9 -2.478 0.912 48 -2.718 1.0000
## 7 - 10 -4.330 0.912 48 -4.750 0.0052
## 7 - 11 -6.478 0.912 48 -7.106 <.0001
## 7 - 12 -8.848 0.912 48 -9.706 <.0001
## 7 - 13 -11.078 0.912 48 -12.152 <.0001
## 7 - 14 -13.107 0.912 48 -14.379 <.0001
## 7 - 15 -14.980 0.912 48 -16.433 <.0001
## 7 - 16 -16.723 0.912 48 -18.344 <.0001
## 7 - 17 -18.260 0.912 48 -20.031 <.0001
## 7 - 18 -19.600 0.912 48 -21.501 <.0001
## 7 - 19 -20.765 0.912 48 -22.779 <.0001
## 7 - 20 -21.785 0.912 48 -23.898 <.0001
## 7 - 21 -22.725 0.912 48 -24.929 <.0001
## 7 - 22 -23.602 0.912 48 -25.892 <.0001
## 7 - 23 -24.420 0.912 48 -26.789 <.0001
## 7 - 24 -25.185 0.912 48 -27.628 <.0001
## 8 - 9 -1.442 0.912 48 -1.582 1.0000
## 8 - 10 -3.295 0.912 48 -3.615 0.1986
## 8 - 11 -5.442 0.912 48 -5.970 0.0001
## 8 - 12 -7.812 0.912 48 -8.570 <.0001
## 8 - 13 -10.043 0.912 48 -11.017 <.0001
## 8 - 14 -12.072 0.912 48 -13.243 <.0001
## 8 - 15 -13.945 0.912 48 -15.298 <.0001
## 8 - 16 -15.688 0.912 48 -17.209 <.0001
## 8 - 17 -17.225 0.912 48 -18.896 <.0001
## 8 - 18 -18.565 0.912 48 -20.366 <.0001
## 8 - 19 -19.730 0.912 48 -21.644 <.0001
## 8 - 20 -20.750 0.912 48 -22.763 <.0001
## 8 - 21 -21.690 0.912 48 -23.794 <.0001
## 8 - 22 -22.567 0.912 48 -24.756 <.0001
## 8 - 23 -23.385 0.912 48 -25.653 <.0001
## 8 - 24 -24.150 0.912 48 -26.492 <.0001
## 9 - 10 -1.853 0.912 48 -2.032 1.0000
## 9 - 11 -4.000 0.912 48 -4.388 0.0172
## 9 - 12 -6.370 0.912 48 -6.988 <.0001
## 9 - 13 -8.600 0.912 48 -9.434 <.0001
## 9 - 14 -10.630 0.912 48 -11.661 <.0001
## 9 - 15 -12.502 0.912 48 -13.715 <.0001
## 9 - 16 -14.245 0.912 48 -15.627 <.0001
## 9 - 17 -15.783 0.912 48 -17.313 <.0001
## 9 - 18 -17.122 0.912 48 -18.783 <.0001
## 9 - 19 -18.288 0.912 48 -20.061 <.0001
## 9 - 20 -19.308 0.912 48 -21.180 <.0001
## 9 - 21 -20.247 0.912 48 -22.211 <.0001
## 9 - 22 -21.125 0.912 48 -23.174 <.0001
## 9 - 23 -21.942 0.912 48 -24.071 <.0001
## 9 - 24 -22.707 0.912 48 -24.910 <.0001
## 10 - 11 -2.147 0.912 48 -2.356 1.0000
## 10 - 12 -4.518 0.912 48 -4.956 0.0026
## 10 - 13 -6.747 0.912 48 -7.402 <.0001
## 10 - 14 -8.777 0.912 48 -9.629 <.0001
## 10 - 15 -10.650 0.912 48 -11.683 <.0001
## 10 - 16 -12.393 0.912 48 -13.594 <.0001
## 10 - 17 -13.930 0.912 48 -15.281 <.0001
## 10 - 18 -15.270 0.912 48 -16.751 <.0001
## 10 - 19 -16.435 0.912 48 -18.029 <.0001
## 10 - 20 -17.455 0.912 48 -19.148 <.0001
## 10 - 21 -18.395 0.912 48 -20.179 <.0001
## 10 - 22 -19.273 0.912 48 -21.142 <.0001
## 10 - 23 -20.090 0.912 48 -22.039 <.0001
## 10 - 24 -20.855 0.912 48 -22.878 <.0001
## 11 - 12 -2.370 0.912 48 -2.600 1.0000
## 11 - 13 -4.600 0.912 48 -5.046 0.0019
## 11 - 14 -6.630 0.912 48 -7.273 <.0001
## 11 - 15 -8.502 0.912 48 -9.327 <.0001
## 11 - 16 -10.245 0.912 48 -11.239 <.0001
## 11 - 17 -11.783 0.912 48 -12.925 <.0001
## 11 - 18 -13.123 0.912 48 -14.395 <.0001
## 11 - 19 -14.287 0.912 48 -15.673 <.0001
## 11 - 20 -15.307 0.912 48 -16.792 <.0001
## 11 - 21 -16.247 0.912 48 -17.823 <.0001
## 11 - 22 -17.125 0.912 48 -18.786 <.0001
## 11 - 23 -17.942 0.912 48 -19.683 <.0001
## 11 - 24 -18.707 0.912 48 -20.522 <.0001
## 12 - 13 -2.230 0.912 48 -2.446 1.0000
## 12 - 14 -4.260 0.912 48 -4.673 0.0067
## 12 - 15 -6.133 0.912 48 -6.727 <.0001
## 12 - 16 -7.875 0.912 48 -8.639 <.0001
## 12 - 17 -9.412 0.912 48 -10.325 <.0001
## 12 - 18 -10.752 0.912 48 -11.795 <.0001
## 12 - 19 -11.918 0.912 48 -13.073 <.0001
## 12 - 20 -12.938 0.912 48 -14.192 <.0001
## 12 - 21 -13.877 0.912 48 -15.224 <.0001
## 12 - 22 -14.755 0.912 48 -16.186 <.0001
## 12 - 23 -15.572 0.912 48 -17.083 <.0001
## 12 - 24 -16.337 0.912 48 -17.922 <.0001
## 13 - 14 -2.030 0.912 48 -2.227 1.0000
## 13 - 15 -3.902 0.912 48 -4.281 0.0244
## 13 - 16 -5.645 0.912 48 -6.193 <.0001
## 13 - 17 -7.183 0.912 48 -7.879 <.0001
## 13 - 18 -8.523 0.912 48 -9.349 <.0001
## 13 - 19 -9.688 0.912 48 -10.627 <.0001
## 13 - 20 -10.707 0.912 48 -11.746 <.0001
## 13 - 21 -11.648 0.912 48 -12.777 <.0001
## 13 - 22 -12.525 0.912 48 -13.740 <.0001
## 13 - 23 -13.342 0.912 48 -14.637 <.0001
## 13 - 24 -14.107 0.912 48 -15.476 <.0001
## 14 - 15 -1.873 0.912 48 -2.054 1.0000
## 14 - 16 -3.615 0.912 48 -3.966 0.0671
## 14 - 17 -5.152 0.912 48 -5.652 0.0002
## 14 - 18 -6.492 0.912 48 -7.122 <.0001
## 14 - 19 -7.657 0.912 48 -8.400 <.0001
## 14 - 20 -8.678 0.912 48 -9.519 <.0001
## 14 - 21 -9.617 0.912 48 -10.550 <.0001
## 14 - 22 -10.495 0.912 48 -11.513 <.0001
## 14 - 23 -11.312 0.912 48 -12.410 <.0001
## 14 - 24 -12.078 0.912 48 -13.249 <.0001
## 15 - 16 -1.742 0.912 48 -1.912 1.0000
## 15 - 17 -3.280 0.912 48 -3.598 0.2087
## 15 - 18 -4.620 0.912 48 -5.068 0.0018
## 15 - 19 -5.785 0.912 48 -6.346 <.0001
## 15 - 20 -6.805 0.912 48 -7.465 <.0001
## 15 - 21 -7.745 0.912 48 -8.496 <.0001
## 15 - 22 -8.623 0.912 48 -9.459 <.0001
## 15 - 23 -9.440 0.912 48 -10.356 <.0001
## 15 - 24 -10.205 0.912 48 -11.195 <.0001
## 16 - 17 -1.538 0.912 48 -1.687 1.0000
## 16 - 18 -2.877 0.912 48 -3.157 0.7610
## 16 - 19 -4.043 0.912 48 -4.435 0.0148
## 16 - 20 -5.062 0.912 48 -5.554 0.0003
## 16 - 21 -6.003 0.912 48 -6.585 <.0001
## 16 - 22 -6.880 0.912 48 -7.547 <.0001
## 16 - 23 -7.697 0.912 48 -8.444 <.0001
## 16 - 24 -8.463 0.912 48 -9.283 <.0001
## 17 - 18 -1.340 0.912 48 -1.470 1.0000
## 17 - 19 -2.505 0.912 48 -2.748 1.0000
## 17 - 20 -3.525 0.912 48 -3.867 0.0914
## 17 - 21 -4.465 0.912 48 -4.898 0.0031
## 17 - 22 -5.343 0.912 48 -5.861 0.0001
## 17 - 23 -6.160 0.912 48 -6.757 <.0001
## 17 - 24 -6.925 0.912 48 -7.597 <.0001
## 18 - 19 -1.165 0.912 48 -1.278 1.0000
## 18 - 20 -2.185 0.912 48 -2.397 1.0000
## 18 - 21 -3.125 0.912 48 -3.428 0.3468
## 18 - 22 -4.003 0.912 48 -4.391 0.0171
## 18 - 23 -4.820 0.912 48 -5.288 0.0008
## 18 - 24 -5.585 0.912 48 -6.127 <.0001
## 19 - 20 -1.020 0.912 48 -1.119 1.0000
## 19 - 21 -1.960 0.912 48 -2.150 1.0000
## 19 - 22 -2.837 0.912 48 -3.113 0.8613
## 19 - 23 -3.655 0.912 48 -4.010 0.0584
## 19 - 24 -4.420 0.912 48 -4.849 0.0037
## 20 - 21 -0.940 0.912 48 -1.031 1.0000
## 20 - 22 -1.817 0.912 48 -1.994 1.0000
## 20 - 23 -2.635 0.912 48 -2.891 1.0000
## 20 - 24 -3.400 0.912 48 -3.730 0.1398
## 21 - 22 -0.877 0.912 48 -0.963 1.0000
## 21 - 23 -1.695 0.912 48 -1.859 1.0000
## 21 - 24 -2.460 0.912 48 -2.699 1.0000
## 22 - 23 -0.818 0.912 48 -0.897 1.0000
## 22 - 24 -1.583 0.912 48 -1.736 1.0000
## 23 - 24 -0.765 0.912 48 -0.839 1.0000
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## P value adjustment: bonferroni method for 276 tests
emmeans(scatmodel1, list(pairwise~TimeF), adjust="tukey")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of TimeF`
## TimeF emmean SE df lower.CL upper.CL
## 1 4.01 0.811 120 2.41 5.62
## 2 5.03 0.811 120 3.42 6.63
## 3 6.60 0.811 120 4.99 8.21
## 4 8.52 0.811 120 6.91 10.13
## 5 10.08 0.811 120 8.47 11.68
## 6 11.37 0.811 120 9.76 12.97
## 7 12.49 0.811 120 10.89 14.10
## 8 13.53 0.811 120 11.92 15.14
## 9 14.97 0.811 120 13.37 16.58
## 10 16.82 0.811 120 15.22 18.43
## 11 18.97 0.811 120 17.37 20.58
## 12 21.34 0.811 120 19.74 22.95
## 13 23.57 0.811 120 21.97 25.18
## 14 25.60 0.811 120 24.00 27.21
## 15 27.48 0.811 120 25.87 29.08
## 16 29.22 0.811 120 27.61 30.82
## 17 30.75 0.811 120 29.15 32.36
## 18 32.09 0.811 120 30.49 33.70
## 19 33.26 0.811 120 31.65 34.87
## 20 34.28 0.811 120 32.67 35.89
## 21 35.22 0.811 120 33.61 36.83
## 22 36.10 0.811 120 34.49 37.70
## 23 36.91 0.811 120 35.31 38.52
## 24 37.68 0.811 120 36.07 39.29
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of TimeF`
## 1 estimate SE df t.ratio p.value
## 1 - 2 -1.012 0.912 48 -1.111 1.0000
## 1 - 3 -2.587 0.912 48 -2.838 0.4357
## 1 - 4 -4.508 0.912 48 -4.945 0.0020
## 1 - 5 -6.065 0.912 48 -6.653 <.0001
## 1 - 6 -7.353 0.912 48 -8.066 <.0001
## 1 - 7 -8.482 0.912 48 -9.305 <.0001
## 1 - 8 -9.518 0.912 48 -10.441 <.0001
## 1 - 9 -10.960 0.912 48 -12.023 <.0001
## 1 - 10 -12.812 0.912 48 -14.055 <.0001
## 1 - 11 -14.960 0.912 48 -16.411 <.0001
## 1 - 12 -17.330 0.912 48 -19.011 <.0001
## 1 - 13 -19.560 0.912 48 -21.457 <.0001
## 1 - 14 -21.590 0.912 48 -23.684 <.0001
## 1 - 15 -23.462 0.912 48 -25.738 <.0001
## 1 - 16 -25.205 0.912 48 -27.650 <.0001
## 1 - 17 -26.742 0.912 48 -29.336 <.0001
## 1 - 18 -28.082 0.912 48 -30.806 <.0001
## 1 - 19 -29.247 0.912 48 -32.084 <.0001
## 1 - 20 -30.267 0.912 48 -33.203 <.0001
## 1 - 21 -31.207 0.912 48 -34.234 <.0001
## 1 - 22 -32.085 0.912 48 -35.197 <.0001
## 1 - 23 -32.903 0.912 48 -36.094 <.0001
## 1 - 24 -33.667 0.912 48 -36.933 <.0001
## 2 - 3 -1.575 0.912 48 -1.728 0.9838
## 2 - 4 -3.495 0.912 48 -3.834 0.0527
## 2 - 5 -5.053 0.912 48 -5.543 0.0003
## 2 - 6 -6.340 0.912 48 -6.955 <.0001
## 2 - 7 -7.470 0.912 48 -8.195 <.0001
## 2 - 8 -8.505 0.912 48 -9.330 <.0001
## 2 - 9 -9.947 0.912 48 -10.912 <.0001
## 2 - 10 -11.800 0.912 48 -12.945 <.0001
## 2 - 11 -13.947 0.912 48 -15.300 <.0001
## 2 - 12 -16.317 0.912 48 -17.900 <.0001
## 2 - 13 -18.547 0.912 48 -20.346 <.0001
## 2 - 14 -20.578 0.912 48 -22.573 <.0001
## 2 - 15 -22.450 0.912 48 -24.628 <.0001
## 2 - 16 -24.192 0.912 48 -26.539 <.0001
## 2 - 17 -25.730 0.912 48 -28.226 <.0001
## 2 - 18 -27.070 0.912 48 -29.696 <.0001
## 2 - 19 -28.235 0.912 48 -30.974 <.0001
## 2 - 20 -29.255 0.912 48 -32.093 <.0001
## 2 - 21 -30.195 0.912 48 -33.124 <.0001
## 2 - 22 -31.073 0.912 48 -34.086 <.0001
## 2 - 23 -31.890 0.912 48 -34.983 <.0001
## 2 - 24 -32.655 0.912 48 -35.822 <.0001
## 3 - 4 -1.920 0.912 48 -2.106 0.8928
## 3 - 5 -3.478 0.912 48 -3.815 0.0555
## 3 - 6 -4.765 0.912 48 -5.227 0.0008
## 3 - 7 -5.895 0.912 48 -6.467 <.0001
## 3 - 8 -6.930 0.912 48 -7.602 <.0001
## 3 - 9 -8.373 0.912 48 -9.185 <.0001
## 3 - 10 -10.225 0.912 48 -11.217 <.0001
## 3 - 11 -12.373 0.912 48 -13.573 <.0001
## 3 - 12 -14.742 0.912 48 -16.172 <.0001
## 3 - 13 -16.973 0.912 48 -18.619 <.0001
## 3 - 14 -19.003 0.912 48 -20.846 <.0001
## 3 - 15 -20.875 0.912 48 -22.900 <.0001
## 3 - 16 -22.617 0.912 48 -24.811 <.0001
## 3 - 17 -24.155 0.912 48 -26.498 <.0001
## 3 - 18 -25.495 0.912 48 -27.968 <.0001
## 3 - 19 -26.660 0.912 48 -29.246 <.0001
## 3 - 20 -27.680 0.912 48 -30.365 <.0001
## 3 - 21 -28.620 0.912 48 -31.396 <.0001
## 3 - 22 -29.497 0.912 48 -32.359 <.0001
## 3 - 23 -30.315 0.912 48 -33.255 <.0001
## 3 - 24 -31.080 0.912 48 -34.095 <.0001
## 4 - 5 -1.558 0.912 48 -1.709 0.9857
## 4 - 6 -2.845 0.912 48 -3.121 0.2684
## 4 - 7 -3.975 0.912 48 -4.361 0.0121
## 4 - 8 -5.010 0.912 48 -5.496 0.0003
## 4 - 9 -6.452 0.912 48 -7.078 <.0001
## 4 - 10 -8.305 0.912 48 -9.111 <.0001
## 4 - 11 -10.453 0.912 48 -11.466 <.0001
## 4 - 12 -12.822 0.912 48 -14.066 <.0001
## 4 - 13 -15.053 0.912 48 -16.512 <.0001
## 4 - 14 -17.082 0.912 48 -18.739 <.0001
## 4 - 15 -18.955 0.912 48 -20.794 <.0001
## 4 - 16 -20.698 0.912 48 -22.705 <.0001
## 4 - 17 -22.235 0.912 48 -24.392 <.0001
## 4 - 18 -23.575 0.912 48 -25.862 <.0001
## 4 - 19 -24.740 0.912 48 -27.140 <.0001
## 4 - 20 -25.760 0.912 48 -28.259 <.0001
## 4 - 21 -26.700 0.912 48 -29.290 <.0001
## 4 - 22 -27.578 0.912 48 -30.252 <.0001
## 4 - 23 -28.395 0.912 48 -31.149 <.0001
## 4 - 24 -29.160 0.912 48 -31.988 <.0001
## 5 - 6 -1.288 0.912 48 -1.412 0.9987
## 5 - 7 -2.417 0.912 48 -2.652 0.5643
## 5 - 8 -3.453 0.912 48 -3.787 0.0595
## 5 - 9 -4.895 0.912 48 -5.370 0.0005
## 5 - 10 -6.747 0.912 48 -7.402 <.0001
## 5 - 11 -8.895 0.912 48 -9.758 <.0001
## 5 - 12 -11.265 0.912 48 -12.358 <.0001
## 5 - 13 -13.495 0.912 48 -14.804 <.0001
## 5 - 14 -15.525 0.912 48 -17.031 <.0001
## 5 - 15 -17.398 0.912 48 -19.085 <.0001
## 5 - 16 -19.140 0.912 48 -20.996 <.0001
## 5 - 17 -20.677 0.912 48 -22.683 <.0001
## 5 - 18 -22.017 0.912 48 -24.153 <.0001
## 5 - 19 -23.183 0.912 48 -25.431 <.0001
## 5 - 20 -24.203 0.912 48 -26.550 <.0001
## 5 - 21 -25.142 0.912 48 -27.581 <.0001
## 5 - 22 -26.020 0.912 48 -28.544 <.0001
## 5 - 23 -26.837 0.912 48 -29.441 <.0001
## 5 - 24 -27.602 0.912 48 -30.280 <.0001
## 6 - 7 -1.130 0.912 48 -1.240 0.9998
## 6 - 8 -2.165 0.912 48 -2.375 0.7526
## 6 - 9 -3.607 0.912 48 -3.957 0.0380
## 6 - 10 -5.460 0.912 48 -5.990 0.0001
## 6 - 11 -7.607 0.912 48 -8.345 <.0001
## 6 - 12 -9.977 0.912 48 -10.945 <.0001
## 6 - 13 -12.207 0.912 48 -13.392 <.0001
## 6 - 14 -14.238 0.912 48 -15.618 <.0001
## 6 - 15 -16.110 0.912 48 -17.673 <.0001
## 6 - 16 -17.852 0.912 48 -19.584 <.0001
## 6 - 17 -19.390 0.912 48 -21.271 <.0001
## 6 - 18 -20.730 0.912 48 -22.741 <.0001
## 6 - 19 -21.895 0.912 48 -24.019 <.0001
## 6 - 20 -22.915 0.912 48 -25.138 <.0001
## 6 - 21 -23.855 0.912 48 -26.169 <.0001
## 6 - 22 -24.733 0.912 48 -27.131 <.0001
## 6 - 23 -25.550 0.912 48 -28.028 <.0001
## 6 - 24 -26.315 0.912 48 -28.867 <.0001
## 7 - 8 -1.035 0.912 48 -1.135 1.0000
## 7 - 9 -2.478 0.912 48 -2.718 0.5181
## 7 - 10 -4.330 0.912 48 -4.750 0.0037
## 7 - 11 -6.478 0.912 48 -7.106 <.0001
## 7 - 12 -8.848 0.912 48 -9.706 <.0001
## 7 - 13 -11.078 0.912 48 -12.152 <.0001
## 7 - 14 -13.107 0.912 48 -14.379 <.0001
## 7 - 15 -14.980 0.912 48 -16.433 <.0001
## 7 - 16 -16.723 0.912 48 -18.344 <.0001
## 7 - 17 -18.260 0.912 48 -20.031 <.0001
## 7 - 18 -19.600 0.912 48 -21.501 <.0001
## 7 - 19 -20.765 0.912 48 -22.779 <.0001
## 7 - 20 -21.785 0.912 48 -23.898 <.0001
## 7 - 21 -22.725 0.912 48 -24.929 <.0001
## 7 - 22 -23.602 0.912 48 -25.892 <.0001
## 7 - 23 -24.420 0.912 48 -26.789 <.0001
## 7 - 24 -25.185 0.912 48 -27.628 <.0001
## 8 - 9 -1.442 0.912 48 -1.582 0.9942
## 8 - 10 -3.295 0.912 48 -3.615 0.0917
## 8 - 11 -5.442 0.912 48 -5.970 0.0001
## 8 - 12 -7.812 0.912 48 -8.570 <.0001
## 8 - 13 -10.043 0.912 48 -11.017 <.0001
## 8 - 14 -12.072 0.912 48 -13.243 <.0001
## 8 - 15 -13.945 0.912 48 -15.298 <.0001
## 8 - 16 -15.688 0.912 48 -17.209 <.0001
## 8 - 17 -17.225 0.912 48 -18.896 <.0001
## 8 - 18 -18.565 0.912 48 -20.366 <.0001
## 8 - 19 -19.730 0.912 48 -21.644 <.0001
## 8 - 20 -20.750 0.912 48 -22.763 <.0001
## 8 - 21 -21.690 0.912 48 -23.794 <.0001
## 8 - 22 -22.567 0.912 48 -24.756 <.0001
## 8 - 23 -23.385 0.912 48 -25.653 <.0001
## 8 - 24 -24.150 0.912 48 -26.492 <.0001
## 9 - 10 -1.853 0.912 48 -2.032 0.9201
## 9 - 11 -4.000 0.912 48 -4.388 0.0112
## 9 - 12 -6.370 0.912 48 -6.988 <.0001
## 9 - 13 -8.600 0.912 48 -9.434 <.0001
## 9 - 14 -10.630 0.912 48 -11.661 <.0001
## 9 - 15 -12.502 0.912 48 -13.715 <.0001
## 9 - 16 -14.245 0.912 48 -15.627 <.0001
## 9 - 17 -15.783 0.912 48 -17.313 <.0001
## 9 - 18 -17.122 0.912 48 -18.783 <.0001
## 9 - 19 -18.288 0.912 48 -20.061 <.0001
## 9 - 20 -19.308 0.912 48 -21.180 <.0001
## 9 - 21 -20.247 0.912 48 -22.211 <.0001
## 9 - 22 -21.125 0.912 48 -23.174 <.0001
## 9 - 23 -21.942 0.912 48 -24.071 <.0001
## 9 - 24 -22.707 0.912 48 -24.910 <.0001
## 10 - 11 -2.147 0.912 48 -2.356 0.7645
## 10 - 12 -4.518 0.912 48 -4.956 0.0019
## 10 - 13 -6.747 0.912 48 -7.402 <.0001
## 10 - 14 -8.777 0.912 48 -9.629 <.0001
## 10 - 15 -10.650 0.912 48 -11.683 <.0001
## 10 - 16 -12.393 0.912 48 -13.594 <.0001
## 10 - 17 -13.930 0.912 48 -15.281 <.0001
## 10 - 18 -15.270 0.912 48 -16.751 <.0001
## 10 - 19 -16.435 0.912 48 -18.029 <.0001
## 10 - 20 -17.455 0.912 48 -19.148 <.0001
## 10 - 21 -18.395 0.912 48 -20.179 <.0001
## 10 - 22 -19.273 0.912 48 -21.142 <.0001
## 10 - 23 -20.090 0.912 48 -22.039 <.0001
## 10 - 24 -20.855 0.912 48 -22.878 <.0001
## 11 - 12 -2.370 0.912 48 -2.600 0.6011
## 11 - 13 -4.600 0.912 48 -5.046 0.0014
## 11 - 14 -6.630 0.912 48 -7.273 <.0001
## 11 - 15 -8.502 0.912 48 -9.327 <.0001
## 11 - 16 -10.245 0.912 48 -11.239 <.0001
## 11 - 17 -11.783 0.912 48 -12.925 <.0001
## 11 - 18 -13.123 0.912 48 -14.395 <.0001
## 11 - 19 -14.287 0.912 48 -15.673 <.0001
## 11 - 20 -15.307 0.912 48 -16.792 <.0001
## 11 - 21 -16.247 0.912 48 -17.823 <.0001
## 11 - 22 -17.125 0.912 48 -18.786 <.0001
## 11 - 23 -17.942 0.912 48 -19.683 <.0001
## 11 - 24 -18.707 0.912 48 -20.522 <.0001
## 12 - 13 -2.230 0.912 48 -2.446 0.7067
## 12 - 14 -4.260 0.912 48 -4.673 0.0047
## 12 - 15 -6.133 0.912 48 -6.727 <.0001
## 12 - 16 -7.875 0.912 48 -8.639 <.0001
## 12 - 17 -9.412 0.912 48 -10.325 <.0001
## 12 - 18 -10.752 0.912 48 -11.795 <.0001
## 12 - 19 -11.918 0.912 48 -13.073 <.0001
## 12 - 20 -12.938 0.912 48 -14.192 <.0001
## 12 - 21 -13.877 0.912 48 -15.224 <.0001
## 12 - 22 -14.755 0.912 48 -16.186 <.0001
## 12 - 23 -15.572 0.912 48 -17.083 <.0001
## 12 - 24 -16.337 0.912 48 -17.922 <.0001
## 13 - 14 -2.030 0.912 48 -2.227 0.8373
## 13 - 15 -3.902 0.912 48 -4.281 0.0153
## 13 - 16 -5.645 0.912 48 -6.193 <.0001
## 13 - 17 -7.183 0.912 48 -7.879 <.0001
## 13 - 18 -8.523 0.912 48 -9.349 <.0001
## 13 - 19 -9.688 0.912 48 -10.627 <.0001
## 13 - 20 -10.707 0.912 48 -11.746 <.0001
## 13 - 21 -11.648 0.912 48 -12.777 <.0001
## 13 - 22 -12.525 0.912 48 -13.740 <.0001
## 13 - 23 -13.342 0.912 48 -14.637 <.0001
## 13 - 24 -14.107 0.912 48 -15.476 <.0001
## 14 - 15 -1.873 0.912 48 -2.054 0.9125
## 14 - 16 -3.615 0.912 48 -3.966 0.0372
## 14 - 17 -5.152 0.912 48 -5.652 0.0002
## 14 - 18 -6.492 0.912 48 -7.122 <.0001
## 14 - 19 -7.657 0.912 48 -8.400 <.0001
## 14 - 20 -8.678 0.912 48 -9.519 <.0001
## 14 - 21 -9.617 0.912 48 -10.550 <.0001
## 14 - 22 -10.495 0.912 48 -11.513 <.0001
## 14 - 23 -11.312 0.912 48 -12.410 <.0001
## 14 - 24 -12.078 0.912 48 -13.249 <.0001
## 15 - 16 -1.742 0.912 48 -1.912 0.9540
## 15 - 17 -3.280 0.912 48 -3.598 0.0954
## 15 - 18 -4.620 0.912 48 -5.068 0.0013
## 15 - 19 -5.785 0.912 48 -6.346 <.0001
## 15 - 20 -6.805 0.912 48 -7.465 <.0001
## 15 - 21 -7.745 0.912 48 -8.496 <.0001
## 15 - 22 -8.623 0.912 48 -9.459 <.0001
## 15 - 23 -9.440 0.912 48 -10.356 <.0001
## 15 - 24 -10.205 0.912 48 -11.195 <.0001
## 16 - 17 -1.538 0.912 48 -1.687 0.9876
## 16 - 18 -2.877 0.912 48 -3.157 0.2506
## 16 - 19 -4.043 0.912 48 -4.435 0.0097
## 16 - 20 -5.062 0.912 48 -5.554 0.0003
## 16 - 21 -6.003 0.912 48 -6.585 <.0001
## 16 - 22 -6.880 0.912 48 -7.547 <.0001
## 16 - 23 -7.697 0.912 48 -8.444 <.0001
## 16 - 24 -8.463 0.912 48 -9.283 <.0001
## 17 - 18 -1.340 0.912 48 -1.470 0.9978
## 17 - 19 -2.505 0.912 48 -2.748 0.4971
## 17 - 20 -3.525 0.912 48 -3.867 0.0484
## 17 - 21 -4.465 0.912 48 -4.898 0.0023
## 17 - 22 -5.343 0.912 48 -5.861 0.0001
## 17 - 23 -6.160 0.912 48 -6.757 <.0001
## 17 - 24 -6.925 0.912 48 -7.597 <.0001
## 18 - 19 -1.165 0.912 48 -1.278 0.9997
## 18 - 20 -2.185 0.912 48 -2.397 0.7388
## 18 - 21 -3.125 0.912 48 -3.428 0.1418
## 18 - 22 -4.003 0.912 48 -4.391 0.0111
## 18 - 23 -4.820 0.912 48 -5.288 0.0007
## 18 - 24 -5.585 0.912 48 -6.127 <.0001
## 19 - 20 -1.020 0.912 48 -1.119 1.0000
## 19 - 21 -1.960 0.912 48 -2.150 0.8741
## 19 - 22 -2.837 0.912 48 -3.113 0.2726
## 19 - 23 -3.655 0.912 48 -4.010 0.0330
## 19 - 24 -4.420 0.912 48 -4.849 0.0027
## 20 - 21 -0.940 0.912 48 -1.031 1.0000
## 20 - 22 -1.817 0.912 48 -1.994 0.9323
## 20 - 23 -2.635 0.912 48 -2.891 0.4017
## 20 - 24 -3.400 0.912 48 -3.730 0.0689
## 21 - 22 -0.877 0.912 48 -0.963 1.0000
## 21 - 23 -1.695 0.912 48 -1.859 0.9649
## 21 - 24 -2.460 0.912 48 -2.699 0.5315
## 22 - 23 -0.818 0.912 48 -0.897 1.0000
## 22 - 24 -1.583 0.912 48 -1.736 0.9829
## 23 - 24 -0.765 0.912 48 -0.839 1.0000
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 24 estimates
emmeans(scatmodel1, list(pairwise~ScatF|TimeF), adjust="bonferroni")
## $`emmeans of ScatF | TimeF`
## TimeF = 1:
## ScatF emmean SE df lower.CL upper.CL
## A 4.93 1.15 120 2.664 7.21
## P 3.09 1.15 120 0.819 5.36
##
## TimeF = 2:
## ScatF emmean SE df lower.CL upper.CL
## A 6.61 1.15 120 4.334 8.88
## P 3.44 1.15 120 1.174 5.72
##
## TimeF = 3:
## ScatF emmean SE df lower.CL upper.CL
## A 8.80 1.15 120 6.534 11.08
## P 4.39 1.15 120 2.124 6.67
##
## TimeF = 4:
## ScatF emmean SE df lower.CL upper.CL
## A 11.34 1.15 120 9.069 13.61
## P 5.70 1.15 120 3.429 7.97
##
## TimeF = 5:
## ScatF emmean SE df lower.CL upper.CL
## A 12.98 1.15 120 10.714 15.26
## P 7.17 1.15 120 4.899 9.44
##
## TimeF = 6:
## ScatF emmean SE df lower.CL upper.CL
## A 13.91 1.15 120 11.639 16.18
## P 8.82 1.15 120 6.549 11.09
##
## TimeF = 7:
## ScatF emmean SE df lower.CL upper.CL
## A 14.42 1.15 120 12.149 16.69
## P 10.57 1.15 120 8.299 12.84
##
## TimeF = 8:
## ScatF emmean SE df lower.CL upper.CL
## A 14.68 1.15 120 12.409 16.95
## P 12.38 1.15 120 10.109 14.65
##
## TimeF = 9:
## ScatF emmean SE df lower.CL upper.CL
## A 15.33 1.15 120 13.059 17.60
## P 14.62 1.15 120 12.344 16.89
##
## TimeF = 10:
## ScatF emmean SE df lower.CL upper.CL
## A 16.43 1.15 120 14.164 18.71
## P 17.21 1.15 120 14.944 19.49
##
## TimeF = 11:
## ScatF emmean SE df lower.CL upper.CL
## A 17.86 1.15 120 15.594 20.14
## P 20.08 1.15 120 17.809 22.35
##
## TimeF = 12:
## ScatF emmean SE df lower.CL upper.CL
## A 19.54 1.15 120 17.264 21.81
## P 23.15 1.15 120 20.879 25.42
##
## TimeF = 13:
## ScatF emmean SE df lower.CL upper.CL
## A 21.04 1.15 120 18.769 23.31
## P 26.11 1.15 120 23.834 28.38
##
## TimeF = 14:
## ScatF emmean SE df lower.CL upper.CL
## A 22.30 1.15 120 20.029 24.57
## P 28.91 1.15 120 26.634 31.18
##
## TimeF = 15:
## ScatF emmean SE df lower.CL upper.CL
## A 23.37 1.15 120 21.099 25.64
## P 31.58 1.15 120 29.309 33.85
##
## TimeF = 16:
## ScatF emmean SE df lower.CL upper.CL
## A 24.29 1.15 120 22.014 26.56
## P 34.15 1.15 120 31.879 36.42
##
## TimeF = 17:
## ScatF emmean SE df lower.CL upper.CL
## A 25.07 1.15 120 22.804 27.35
## P 36.44 1.15 120 34.164 38.71
##
## TimeF = 18:
## ScatF emmean SE df lower.CL upper.CL
## A 25.77 1.15 120 23.494 28.04
## P 38.42 1.15 120 36.154 40.70
##
## TimeF = 19:
## ScatF emmean SE df lower.CL upper.CL
## A 26.36 1.15 120 24.094 28.64
## P 40.16 1.15 120 37.884 42.43
##
## TimeF = 20:
## ScatF emmean SE df lower.CL upper.CL
## A 26.89 1.15 120 24.619 29.16
## P 41.67 1.15 120 39.399 43.94
##
## TimeF = 21:
## ScatF emmean SE df lower.CL upper.CL
## A 27.41 1.15 120 25.134 29.68
## P 43.03 1.15 120 40.764 45.31
##
## TimeF = 22:
## ScatF emmean SE df lower.CL upper.CL
## A 27.93 1.15 120 25.664 30.21
## P 44.26 1.15 120 41.989 46.53
##
## TimeF = 23:
## ScatF emmean SE df lower.CL upper.CL
## A 28.48 1.15 120 26.204 30.75
## P 45.35 1.15 120 43.084 47.63
##
## TimeF = 24:
## ScatF emmean SE df lower.CL upper.CL
## A 29.02 1.15 120 26.754 31.30
## P 46.34 1.15 120 44.064 48.61
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF | TimeF`
## TimeF = 1:
## 2 estimate SE df t.ratio p.value
## A - P 1.845 1.62 120 1.138 0.2575
##
## TimeF = 2:
## 2 estimate SE df t.ratio p.value
## A - P 3.160 1.62 120 1.948 0.0537
##
## TimeF = 3:
## 2 estimate SE df t.ratio p.value
## A - P 4.410 1.62 120 2.719 0.0075
##
## TimeF = 4:
## 2 estimate SE df t.ratio p.value
## A - P 5.640 1.62 120 3.478 0.0007
##
## TimeF = 5:
## 2 estimate SE df t.ratio p.value
## A - P 5.815 1.62 120 3.586 0.0005
##
## TimeF = 6:
## 2 estimate SE df t.ratio p.value
## A - P 5.090 1.62 120 3.138 0.0021
##
## TimeF = 7:
## 2 estimate SE df t.ratio p.value
## A - P 3.850 1.62 120 2.374 0.0192
##
## TimeF = 8:
## 2 estimate SE df t.ratio p.value
## A - P 2.300 1.62 120 1.418 0.1587
##
## TimeF = 9:
## 2 estimate SE df t.ratio p.value
## A - P 0.715 1.62 120 0.441 0.6601
##
## TimeF = 10:
## 2 estimate SE df t.ratio p.value
## A - P -0.780 1.62 120 -0.481 0.6314
##
## TimeF = 11:
## 2 estimate SE df t.ratio p.value
## A - P -2.215 1.62 120 -1.366 0.1746
##
## TimeF = 12:
## 2 estimate SE df t.ratio p.value
## A - P -3.615 1.62 120 -2.229 0.0277
##
## TimeF = 13:
## 2 estimate SE df t.ratio p.value
## A - P -5.065 1.62 120 -3.123 0.0022
##
## TimeF = 14:
## 2 estimate SE df t.ratio p.value
## A - P -6.605 1.62 120 -4.073 0.0001
##
## TimeF = 15:
## 2 estimate SE df t.ratio p.value
## A - P -8.210 1.62 120 -5.062 <.0001
##
## TimeF = 16:
## 2 estimate SE df t.ratio p.value
## A - P -9.865 1.62 120 -6.083 <.0001
##
## TimeF = 17:
## 2 estimate SE df t.ratio p.value
## A - P -11.360 1.62 120 -7.005 <.0001
##
## TimeF = 18:
## 2 estimate SE df t.ratio p.value
## A - P -12.660 1.62 120 -7.806 <.0001
##
## TimeF = 19:
## 2 estimate SE df t.ratio p.value
## A - P -13.790 1.62 120 -8.503 <.0001
##
## TimeF = 20:
## 2 estimate SE df t.ratio p.value
## A - P -14.780 1.62 120 -9.113 <.0001
##
## TimeF = 21:
## 2 estimate SE df t.ratio p.value
## A - P -15.630 1.62 120 -9.637 <.0001
##
## TimeF = 22:
## 2 estimate SE df t.ratio p.value
## A - P -16.325 1.62 120 -10.066 <.0001
##
## TimeF = 23:
## 2 estimate SE df t.ratio p.value
## A - P -16.880 1.62 120 -10.408 <.0001
##
## TimeF = 24:
## 2 estimate SE df t.ratio p.value
## A - P -17.310 1.62 120 -10.673 <.0001
##
## Degrees-of-freedom method: kenward-roger
scatmodel2<-lmer(NNitr~ScatF*TimeF+(1|Subject),data=scatted_soilsdata)
anova(scatmodel2)
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## ScatF 194.0 194.02 1 0 133.246 1
## TimeF 11775.3 511.97 23 48 351.596 <2e-16 ***
## ScatF:TimeF 1561.5 67.89 23 48 46.624 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
eta_sq(scatmodel2,partial=TRUE)
## term partial.etasq
## 1 ScatF 1.000
## 2 TimeF 0.994
## 3 ScatF:TimeF 0.957
emmeans(scatmodel2, list(pairwise~ScatF), adjust="bonferroni")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of ScatF`
## ScatF emmean SE df lower.CL upper.CL
## A 18.1 0.313 498 17.5 18.7
## P 23.2 0.313 498 22.6 23.8
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF`
## 1 estimate SE df t.ratio p.value
## A - P -5.1 0.442 498 -11.543 <.0001
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
emmeans(scatmodel2, list(pairwise~ScatF), adjust="tukey")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of ScatF`
## ScatF emmean SE df lower.CL upper.CL
## A 18.1 0.313 498 17.5 18.7
## P 23.2 0.313 498 22.6 23.8
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF`
## 1 estimate SE df t.ratio p.value
## A - P -5.1 0.442 498 -11.543 <.0001
##
## Results are averaged over the levels of: TimeF
## Degrees-of-freedom method: kenward-roger
emmeans(scatmodel2, list(pairwise~TimeF), adjust="bonferroni")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of TimeF`
## TimeF emmean SE df lower.CL upper.CL
## 1 2.30 0.631 57.3 1.04 3.56
## 2 3.54 0.631 57.3 2.28 4.81
## 3 5.28 0.631 57.3 4.02 6.55
## 4 7.25 0.631 57.3 5.99 8.52
## 5 8.84 0.631 57.3 7.58 10.11
## 6 10.15 0.631 57.3 8.89 11.41
## 7 11.30 0.631 57.3 10.03 12.56
## 8 12.35 0.631 57.3 11.09 13.62
## 9 13.81 0.631 57.3 12.55 15.08
## 10 15.68 0.631 57.3 14.41 16.94
## 11 17.84 0.631 57.3 16.58 19.10
## 12 20.23 0.631 57.3 18.96 21.49
## 13 22.45 0.631 57.3 21.19 23.72
## 14 24.47 0.631 57.3 23.21 25.74
## 15 26.32 0.631 57.3 25.06 27.59
## 16 28.04 0.631 57.3 26.77 29.30
## 17 29.55 0.631 57.3 28.28 30.81
## 18 30.86 0.631 57.3 29.59 32.12
## 19 32.00 0.631 57.3 30.73 33.26
## 20 32.99 0.631 57.3 31.73 34.26
## 21 33.92 0.631 57.3 32.65 35.18
## 22 34.78 0.631 57.3 33.51 36.04
## 23 35.58 0.631 57.3 34.32 36.84
## 24 36.33 0.631 57.3 35.07 37.60
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of TimeF`
## 1 estimate SE df t.ratio p.value
## 1 - 2 -1.245 0.853 48 -1.459 1.0000
## 1 - 3 -2.982 0.853 48 -3.495 0.2841
## 1 - 4 -4.955 0.853 48 -5.807 0.0001
## 1 - 5 -6.543 0.853 48 -7.668 <.0001
## 1 - 6 -7.850 0.853 48 -9.200 <.0001
## 1 - 7 -8.998 0.853 48 -10.545 <.0001
## 1 - 8 -10.053 0.853 48 -11.781 <.0001
## 1 - 9 -11.515 0.853 48 -13.495 <.0001
## 1 - 10 -13.377 0.853 48 -15.678 <.0001
## 1 - 11 -15.540 0.853 48 -18.212 <.0001
## 1 - 12 -17.925 0.853 48 -21.007 <.0001
## 1 - 13 -20.155 0.853 48 -23.621 <.0001
## 1 - 14 -22.172 0.853 48 -25.985 <.0001
## 1 - 15 -24.025 0.853 48 -28.156 <.0001
## 1 - 16 -25.735 0.853 48 -30.160 <.0001
## 1 - 17 -27.245 0.853 48 -31.930 <.0001
## 1 - 18 -28.555 0.853 48 -33.465 <.0001
## 1 - 19 -29.695 0.853 48 -34.801 <.0001
## 1 - 20 -30.695 0.853 48 -35.973 <.0001
## 1 - 21 -31.617 0.853 48 -37.055 <.0001
## 1 - 22 -32.477 0.853 48 -38.062 <.0001
## 1 - 23 -33.280 0.853 48 -39.003 <.0001
## 1 - 24 -34.032 0.853 48 -39.885 <.0001
## 2 - 3 -1.738 0.853 48 -2.036 1.0000
## 2 - 4 -3.710 0.853 48 -4.348 0.0196
## 2 - 5 -5.298 0.853 48 -6.208 <.0001
## 2 - 6 -6.605 0.853 48 -7.741 <.0001
## 2 - 7 -7.753 0.853 48 -9.086 <.0001
## 2 - 8 -8.807 0.853 48 -10.322 <.0001
## 2 - 9 -10.270 0.853 48 -12.036 <.0001
## 2 - 10 -12.133 0.853 48 -14.219 <.0001
## 2 - 11 -14.295 0.853 48 -16.753 <.0001
## 2 - 12 -16.680 0.853 48 -19.548 <.0001
## 2 - 13 -18.910 0.853 48 -22.162 <.0001
## 2 - 14 -20.927 0.853 48 -24.526 <.0001
## 2 - 15 -22.780 0.853 48 -26.697 <.0001
## 2 - 16 -24.490 0.853 48 -28.701 <.0001
## 2 - 17 -26.000 0.853 48 -30.471 <.0001
## 2 - 18 -27.310 0.853 48 -32.006 <.0001
## 2 - 19 -28.450 0.853 48 -33.342 <.0001
## 2 - 20 -29.450 0.853 48 -34.514 <.0001
## 2 - 21 -30.372 0.853 48 -35.595 <.0001
## 2 - 22 -31.233 0.853 48 -36.603 <.0001
## 2 - 23 -32.035 0.853 48 -37.544 <.0001
## 2 - 24 -32.788 0.853 48 -38.426 <.0001
## 3 - 4 -1.972 0.853 48 -2.312 1.0000
## 3 - 5 -3.560 0.853 48 -4.172 0.0347
## 3 - 6 -4.867 0.853 48 -5.705 0.0002
## 3 - 7 -6.015 0.853 48 -7.049 <.0001
## 3 - 8 -7.070 0.853 48 -8.286 <.0001
## 3 - 9 -8.533 0.853 48 -10.000 <.0001
## 3 - 10 -10.395 0.853 48 -12.183 <.0001
## 3 - 11 -12.557 0.853 48 -14.717 <.0001
## 3 - 12 -14.943 0.853 48 -17.512 <.0001
## 3 - 13 -17.172 0.853 48 -20.126 <.0001
## 3 - 14 -19.190 0.853 48 -22.490 <.0001
## 3 - 15 -21.043 0.853 48 -24.661 <.0001
## 3 - 16 -22.753 0.853 48 -26.665 <.0001
## 3 - 17 -24.262 0.853 48 -28.435 <.0001
## 3 - 18 -25.573 0.853 48 -29.970 <.0001
## 3 - 19 -26.712 0.853 48 -31.306 <.0001
## 3 - 20 -27.712 0.853 48 -32.478 <.0001
## 3 - 21 -28.635 0.853 48 -33.559 <.0001
## 3 - 22 -29.495 0.853 48 -34.567 <.0001
## 3 - 23 -30.297 0.853 48 -35.508 <.0001
## 3 - 24 -31.050 0.853 48 -36.389 <.0001
## 4 - 5 -1.587 0.853 48 -1.860 1.0000
## 4 - 6 -2.895 0.853 48 -3.393 0.3848
## 4 - 7 -4.043 0.853 48 -4.738 0.0054
## 4 - 8 -5.098 0.853 48 -5.974 0.0001
## 4 - 9 -6.560 0.853 48 -7.688 <.0001
## 4 - 10 -8.422 0.853 48 -9.871 <.0001
## 4 - 11 -10.585 0.853 48 -12.405 <.0001
## 4 - 12 -12.970 0.853 48 -15.200 <.0001
## 4 - 13 -15.200 0.853 48 -17.814 <.0001
## 4 - 14 -17.218 0.853 48 -20.178 <.0001
## 4 - 15 -19.070 0.853 48 -22.349 <.0001
## 4 - 16 -20.780 0.853 48 -24.353 <.0001
## 4 - 17 -22.290 0.853 48 -26.123 <.0001
## 4 - 18 -23.600 0.853 48 -27.658 <.0001
## 4 - 19 -24.740 0.853 48 -28.994 <.0001
## 4 - 20 -25.740 0.853 48 -30.166 <.0001
## 4 - 21 -26.663 0.853 48 -31.247 <.0001
## 4 - 22 -27.523 0.853 48 -32.255 <.0001
## 4 - 23 -28.325 0.853 48 -33.196 <.0001
## 4 - 24 -29.078 0.853 48 -34.078 <.0001
## 5 - 6 -1.308 0.853 48 -1.532 1.0000
## 5 - 7 -2.455 0.853 48 -2.877 1.0000
## 5 - 8 -3.510 0.853 48 -4.114 0.0419
## 5 - 9 -4.973 0.853 48 -5.828 0.0001
## 5 - 10 -6.835 0.853 48 -8.010 <.0001
## 5 - 11 -8.998 0.853 48 -10.545 <.0001
## 5 - 12 -11.383 0.853 48 -13.340 <.0001
## 5 - 13 -13.613 0.853 48 -15.953 <.0001
## 5 - 14 -15.630 0.853 48 -18.318 <.0001
## 5 - 15 -17.483 0.853 48 -20.489 <.0001
## 5 - 16 -19.192 0.853 48 -22.493 <.0001
## 5 - 17 -20.703 0.853 48 -24.263 <.0001
## 5 - 18 -22.012 0.853 48 -25.798 <.0001
## 5 - 19 -23.152 0.853 48 -27.134 <.0001
## 5 - 20 -24.152 0.853 48 -28.306 <.0001
## 5 - 21 -25.075 0.853 48 -29.387 <.0001
## 5 - 22 -25.935 0.853 48 -30.395 <.0001
## 5 - 23 -26.738 0.853 48 -31.335 <.0001
## 5 - 24 -27.490 0.853 48 -32.217 <.0001
## 6 - 7 -1.147 0.853 48 -1.345 1.0000
## 6 - 8 -2.203 0.853 48 -2.581 1.0000
## 6 - 9 -3.665 0.853 48 -4.295 0.0233
## 6 - 10 -5.527 0.853 48 -6.478 <.0001
## 6 - 11 -7.690 0.853 48 -9.012 <.0001
## 6 - 12 -10.075 0.853 48 -11.808 <.0001
## 6 - 13 -12.305 0.853 48 -14.421 <.0001
## 6 - 14 -14.322 0.853 48 -16.785 <.0001
## 6 - 15 -16.175 0.853 48 -18.956 <.0001
## 6 - 16 -17.885 0.853 48 -20.961 <.0001
## 6 - 17 -19.395 0.853 48 -22.730 <.0001
## 6 - 18 -20.705 0.853 48 -24.266 <.0001
## 6 - 19 -21.845 0.853 48 -25.602 <.0001
## 6 - 20 -22.845 0.853 48 -26.773 <.0001
## 6 - 21 -23.767 0.853 48 -27.855 <.0001
## 6 - 22 -24.628 0.853 48 -28.863 <.0001
## 6 - 23 -25.430 0.853 48 -29.803 <.0001
## 6 - 24 -26.183 0.853 48 -30.685 <.0001
## 7 - 8 -1.055 0.853 48 -1.236 1.0000
## 7 - 9 -2.518 0.853 48 -2.950 1.0000
## 7 - 10 -4.380 0.853 48 -5.133 0.0014
## 7 - 11 -6.543 0.853 48 -7.668 <.0001
## 7 - 12 -8.928 0.853 48 -10.463 <.0001
## 7 - 13 -11.158 0.853 48 -13.076 <.0001
## 7 - 14 -13.175 0.853 48 -15.441 <.0001
## 7 - 15 -15.027 0.853 48 -17.612 <.0001
## 7 - 16 -16.738 0.853 48 -19.616 <.0001
## 7 - 17 -18.247 0.853 48 -21.385 <.0001
## 7 - 18 -19.558 0.853 48 -22.921 <.0001
## 7 - 19 -20.698 0.853 48 -24.257 <.0001
## 7 - 20 -21.698 0.853 48 -25.429 <.0001
## 7 - 21 -22.620 0.853 48 -26.510 <.0001
## 7 - 22 -23.480 0.853 48 -27.518 <.0001
## 7 - 23 -24.282 0.853 48 -28.458 <.0001
## 7 - 24 -25.035 0.853 48 -29.340 <.0001
## 8 - 9 -1.462 0.853 48 -1.714 1.0000
## 8 - 10 -3.325 0.853 48 -3.897 0.0833
## 8 - 11 -5.487 0.853 48 -6.431 <.0001
## 8 - 12 -7.872 0.853 48 -9.226 <.0001
## 8 - 13 -10.102 0.853 48 -11.840 <.0001
## 8 - 14 -12.120 0.853 48 -14.204 <.0001
## 8 - 15 -13.973 0.853 48 -16.375 <.0001
## 8 - 16 -15.682 0.853 48 -18.379 <.0001
## 8 - 17 -17.192 0.853 48 -20.149 <.0001
## 8 - 18 -18.503 0.853 48 -21.684 <.0001
## 8 - 19 -19.642 0.853 48 -23.020 <.0001
## 8 - 20 -20.642 0.853 48 -24.192 <.0001
## 8 - 21 -21.565 0.853 48 -25.273 <.0001
## 8 - 22 -22.425 0.853 48 -26.281 <.0001
## 8 - 23 -23.227 0.853 48 -27.222 <.0001
## 8 - 24 -23.980 0.853 48 -28.104 <.0001
## 9 - 10 -1.863 0.853 48 -2.183 1.0000
## 9 - 11 -4.025 0.853 48 -4.717 0.0058
## 9 - 12 -6.410 0.853 48 -7.512 <.0001
## 9 - 13 -8.640 0.853 48 -10.126 <.0001
## 9 - 14 -10.658 0.853 48 -12.490 <.0001
## 9 - 15 -12.510 0.853 48 -14.661 <.0001
## 9 - 16 -14.220 0.853 48 -16.665 <.0001
## 9 - 17 -15.730 0.853 48 -18.435 <.0001
## 9 - 18 -17.040 0.853 48 -19.970 <.0001
## 9 - 19 -18.180 0.853 48 -21.306 <.0001
## 9 - 20 -19.180 0.853 48 -22.478 <.0001
## 9 - 21 -20.102 0.853 48 -23.559 <.0001
## 9 - 22 -20.962 0.853 48 -24.567 <.0001
## 9 - 23 -21.765 0.853 48 -25.508 <.0001
## 9 - 24 -22.517 0.853 48 -26.390 <.0001
## 10 - 11 -2.163 0.853 48 -2.534 1.0000
## 10 - 12 -4.548 0.853 48 -5.330 0.0007
## 10 - 13 -6.777 0.853 48 -7.943 <.0001
## 10 - 14 -8.795 0.853 48 -10.307 <.0001
## 10 - 15 -10.648 0.853 48 -12.478 <.0001
## 10 - 16 -12.357 0.853 48 -14.483 <.0001
## 10 - 17 -13.867 0.853 48 -16.252 <.0001
## 10 - 18 -15.178 0.853 48 -17.787 <.0001
## 10 - 19 -16.317 0.853 48 -19.124 <.0001
## 10 - 20 -17.317 0.853 48 -20.295 <.0001
## 10 - 21 -18.240 0.853 48 -21.377 <.0001
## 10 - 22 -19.100 0.853 48 -22.384 <.0001
## 10 - 23 -19.902 0.853 48 -23.325 <.0001
## 10 - 24 -20.655 0.853 48 -24.207 <.0001
## 11 - 12 -2.385 0.853 48 -2.795 1.0000
## 11 - 13 -4.615 0.853 48 -5.409 0.0005
## 11 - 14 -6.633 0.853 48 -7.773 <.0001
## 11 - 15 -8.485 0.853 48 -9.944 <.0001
## 11 - 16 -10.195 0.853 48 -11.948 <.0001
## 11 - 17 -11.705 0.853 48 -13.718 <.0001
## 11 - 18 -13.015 0.853 48 -15.253 <.0001
## 11 - 19 -14.155 0.853 48 -16.589 <.0001
## 11 - 20 -15.155 0.853 48 -17.761 <.0001
## 11 - 21 -16.078 0.853 48 -18.842 <.0001
## 11 - 22 -16.938 0.853 48 -19.850 <.0001
## 11 - 23 -17.740 0.853 48 -20.791 <.0001
## 11 - 24 -18.492 0.853 48 -21.673 <.0001
## 12 - 13 -2.230 0.853 48 -2.613 1.0000
## 12 - 14 -4.247 0.853 48 -4.978 0.0024
## 12 - 15 -6.100 0.853 48 -7.149 <.0001
## 12 - 16 -7.810 0.853 48 -9.153 <.0001
## 12 - 17 -9.320 0.853 48 -10.923 <.0001
## 12 - 18 -10.630 0.853 48 -12.458 <.0001
## 12 - 19 -11.770 0.853 48 -13.794 <.0001
## 12 - 20 -12.770 0.853 48 -14.966 <.0001
## 12 - 21 -13.693 0.853 48 -16.047 <.0001
## 12 - 22 -14.553 0.853 48 -17.055 <.0001
## 12 - 23 -15.355 0.853 48 -17.995 <.0001
## 12 - 24 -16.108 0.853 48 -18.877 <.0001
## 13 - 14 -2.018 0.853 48 -2.364 1.0000
## 13 - 15 -3.870 0.853 48 -4.536 0.0106
## 13 - 16 -5.580 0.853 48 -6.540 <.0001
## 13 - 17 -7.090 0.853 48 -8.309 <.0001
## 13 - 18 -8.400 0.853 48 -9.844 <.0001
## 13 - 19 -9.540 0.853 48 -11.181 <.0001
## 13 - 20 -10.540 0.853 48 -12.352 <.0001
## 13 - 21 -11.463 0.853 48 -13.434 <.0001
## 13 - 22 -12.322 0.853 48 -14.442 <.0001
## 13 - 23 -13.125 0.853 48 -15.382 <.0001
## 13 - 24 -13.877 0.853 48 -16.264 <.0001
## 14 - 15 -1.853 0.853 48 -2.171 1.0000
## 14 - 16 -3.562 0.853 48 -4.175 0.0344
## 14 - 17 -5.072 0.853 48 -5.945 0.0001
## 14 - 18 -6.383 0.853 48 -7.480 <.0001
## 14 - 19 -7.522 0.853 48 -8.816 <.0001
## 14 - 20 -8.523 0.853 48 -9.988 <.0001
## 14 - 21 -9.445 0.853 48 -11.069 <.0001
## 14 - 22 -10.305 0.853 48 -12.077 <.0001
## 14 - 23 -11.107 0.853 48 -13.018 <.0001
## 14 - 24 -11.860 0.853 48 -13.899 <.0001
## 15 - 16 -1.710 0.853 48 -2.004 1.0000
## 15 - 17 -3.220 0.853 48 -3.774 0.1221
## 15 - 18 -4.530 0.853 48 -5.309 0.0008
## 15 - 19 -5.670 0.853 48 -6.645 <.0001
## 15 - 20 -6.670 0.853 48 -7.817 <.0001
## 15 - 21 -7.593 0.853 48 -8.898 <.0001
## 15 - 22 -8.453 0.853 48 -9.906 <.0001
## 15 - 23 -9.255 0.853 48 -10.847 <.0001
## 15 - 24 -10.008 0.853 48 -11.728 <.0001
## 16 - 17 -1.510 0.853 48 -1.770 1.0000
## 16 - 18 -2.820 0.853 48 -3.305 0.4974
## 16 - 19 -3.960 0.853 48 -4.641 0.0075
## 16 - 20 -4.960 0.853 48 -5.813 0.0001
## 16 - 21 -5.883 0.853 48 -6.894 <.0001
## 16 - 22 -6.742 0.853 48 -7.902 <.0001
## 16 - 23 -7.545 0.853 48 -8.842 <.0001
## 16 - 24 -8.297 0.853 48 -9.724 <.0001
## 17 - 18 -1.310 0.853 48 -1.535 1.0000
## 17 - 19 -2.450 0.853 48 -2.871 1.0000
## 17 - 20 -3.450 0.853 48 -4.043 0.0525
## 17 - 21 -4.372 0.853 48 -5.124 0.0015
## 17 - 22 -5.232 0.853 48 -6.132 <.0001
## 17 - 23 -6.035 0.853 48 -7.073 <.0001
## 17 - 24 -6.787 0.853 48 -7.955 <.0001
## 18 - 19 -1.140 0.853 48 -1.336 1.0000
## 18 - 20 -2.140 0.853 48 -2.508 1.0000
## 18 - 21 -3.062 0.853 48 -3.589 0.2145
## 18 - 22 -3.922 0.853 48 -4.597 0.0086
## 18 - 23 -4.725 0.853 48 -5.538 0.0003
## 18 - 24 -5.478 0.853 48 -6.419 <.0001
## 19 - 20 -1.000 0.853 48 -1.172 1.0000
## 19 - 21 -1.923 0.853 48 -2.253 1.0000
## 19 - 22 -2.783 0.853 48 -3.261 0.5648
## 19 - 23 -3.585 0.853 48 -4.201 0.0316
## 19 - 24 -4.338 0.853 48 -5.083 0.0017
## 20 - 21 -0.922 0.853 48 -1.081 1.0000
## 20 - 22 -1.782 0.853 48 -2.089 1.0000
## 20 - 23 -2.585 0.853 48 -3.030 1.0000
## 20 - 24 -3.337 0.853 48 -3.911 0.0795
## 21 - 22 -0.860 0.853 48 -1.008 1.0000
## 21 - 23 -1.663 0.853 48 -1.948 1.0000
## 21 - 24 -2.415 0.853 48 -2.830 1.0000
## 22 - 23 -0.802 0.853 48 -0.941 1.0000
## 22 - 24 -1.555 0.853 48 -1.822 1.0000
## 23 - 24 -0.752 0.853 48 -0.882 1.0000
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## P value adjustment: bonferroni method for 276 tests
emmeans(scatmodel2, list(pairwise~TimeF), adjust="tukey")
## NOTE: Results may be misleading due to involvement in interactions
## $`emmeans of TimeF`
## TimeF emmean SE df lower.CL upper.CL
## 1 2.30 0.631 57.3 1.04 3.56
## 2 3.54 0.631 57.3 2.28 4.81
## 3 5.28 0.631 57.3 4.02 6.55
## 4 7.25 0.631 57.3 5.99 8.52
## 5 8.84 0.631 57.3 7.58 10.11
## 6 10.15 0.631 57.3 8.89 11.41
## 7 11.30 0.631 57.3 10.03 12.56
## 8 12.35 0.631 57.3 11.09 13.62
## 9 13.81 0.631 57.3 12.55 15.08
## 10 15.68 0.631 57.3 14.41 16.94
## 11 17.84 0.631 57.3 16.58 19.10
## 12 20.23 0.631 57.3 18.96 21.49
## 13 22.45 0.631 57.3 21.19 23.72
## 14 24.47 0.631 57.3 23.21 25.74
## 15 26.32 0.631 57.3 25.06 27.59
## 16 28.04 0.631 57.3 26.77 29.30
## 17 29.55 0.631 57.3 28.28 30.81
## 18 30.86 0.631 57.3 29.59 32.12
## 19 32.00 0.631 57.3 30.73 33.26
## 20 32.99 0.631 57.3 31.73 34.26
## 21 33.92 0.631 57.3 32.65 35.18
## 22 34.78 0.631 57.3 33.51 36.04
## 23 35.58 0.631 57.3 34.32 36.84
## 24 36.33 0.631 57.3 35.07 37.60
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of TimeF`
## 1 estimate SE df t.ratio p.value
## 1 - 2 -1.245 0.853 48 -1.459 0.9980
## 1 - 3 -2.982 0.853 48 -3.495 0.1216
## 1 - 4 -4.955 0.853 48 -5.807 0.0001
## 1 - 5 -6.543 0.853 48 -7.668 <.0001
## 1 - 6 -7.850 0.853 48 -9.200 <.0001
## 1 - 7 -8.998 0.853 48 -10.545 <.0001
## 1 - 8 -10.053 0.853 48 -11.781 <.0001
## 1 - 9 -11.515 0.853 48 -13.495 <.0001
## 1 - 10 -13.377 0.853 48 -15.678 <.0001
## 1 - 11 -15.540 0.853 48 -18.212 <.0001
## 1 - 12 -17.925 0.853 48 -21.007 <.0001
## 1 - 13 -20.155 0.853 48 -23.621 <.0001
## 1 - 14 -22.172 0.853 48 -25.985 <.0001
## 1 - 15 -24.025 0.853 48 -28.156 <.0001
## 1 - 16 -25.735 0.853 48 -30.160 <.0001
## 1 - 17 -27.245 0.853 48 -31.930 <.0001
## 1 - 18 -28.555 0.853 48 -33.465 <.0001
## 1 - 19 -29.695 0.853 48 -34.801 <.0001
## 1 - 20 -30.695 0.853 48 -35.973 <.0001
## 1 - 21 -31.617 0.853 48 -37.055 <.0001
## 1 - 22 -32.477 0.853 48 -38.062 <.0001
## 1 - 23 -33.280 0.853 48 -39.003 <.0001
## 1 - 24 -34.032 0.853 48 -39.885 <.0001
## 2 - 3 -1.738 0.853 48 -2.036 0.9187
## 2 - 4 -3.710 0.853 48 -4.348 0.0126
## 2 - 5 -5.298 0.853 48 -6.208 <.0001
## 2 - 6 -6.605 0.853 48 -7.741 <.0001
## 2 - 7 -7.753 0.853 48 -9.086 <.0001
## 2 - 8 -8.807 0.853 48 -10.322 <.0001
## 2 - 9 -10.270 0.853 48 -12.036 <.0001
## 2 - 10 -12.133 0.853 48 -14.219 <.0001
## 2 - 11 -14.295 0.853 48 -16.753 <.0001
## 2 - 12 -16.680 0.853 48 -19.548 <.0001
## 2 - 13 -18.910 0.853 48 -22.162 <.0001
## 2 - 14 -20.927 0.853 48 -24.526 <.0001
## 2 - 15 -22.780 0.853 48 -26.697 <.0001
## 2 - 16 -24.490 0.853 48 -28.701 <.0001
## 2 - 17 -26.000 0.853 48 -30.471 <.0001
## 2 - 18 -27.310 0.853 48 -32.006 <.0001
## 2 - 19 -28.450 0.853 48 -33.342 <.0001
## 2 - 20 -29.450 0.853 48 -34.514 <.0001
## 2 - 21 -30.372 0.853 48 -35.595 <.0001
## 2 - 22 -31.233 0.853 48 -36.603 <.0001
## 2 - 23 -32.035 0.853 48 -37.544 <.0001
## 2 - 24 -32.788 0.853 48 -38.426 <.0001
## 3 - 4 -1.972 0.853 48 -2.312 0.7908
## 3 - 5 -3.560 0.853 48 -4.172 0.0209
## 3 - 6 -4.867 0.853 48 -5.705 0.0002
## 3 - 7 -6.015 0.853 48 -7.049 <.0001
## 3 - 8 -7.070 0.853 48 -8.286 <.0001
## 3 - 9 -8.533 0.853 48 -10.000 <.0001
## 3 - 10 -10.395 0.853 48 -12.183 <.0001
## 3 - 11 -12.557 0.853 48 -14.717 <.0001
## 3 - 12 -14.943 0.853 48 -17.512 <.0001
## 3 - 13 -17.172 0.853 48 -20.126 <.0001
## 3 - 14 -19.190 0.853 48 -22.490 <.0001
## 3 - 15 -21.043 0.853 48 -24.661 <.0001
## 3 - 16 -22.753 0.853 48 -26.665 <.0001
## 3 - 17 -24.262 0.853 48 -28.435 <.0001
## 3 - 18 -25.573 0.853 48 -29.970 <.0001
## 3 - 19 -26.712 0.853 48 -31.306 <.0001
## 3 - 20 -27.712 0.853 48 -32.478 <.0001
## 3 - 21 -28.635 0.853 48 -33.559 <.0001
## 3 - 22 -29.495 0.853 48 -34.567 <.0001
## 3 - 23 -30.297 0.853 48 -35.508 <.0001
## 3 - 24 -31.050 0.853 48 -36.389 <.0001
## 4 - 5 -1.587 0.853 48 -1.860 0.9646
## 4 - 6 -2.895 0.853 48 -3.393 0.1534
## 4 - 7 -4.043 0.853 48 -4.738 0.0039
## 4 - 8 -5.098 0.853 48 -5.974 0.0001
## 4 - 9 -6.560 0.853 48 -7.688 <.0001
## 4 - 10 -8.422 0.853 48 -9.871 <.0001
## 4 - 11 -10.585 0.853 48 -12.405 <.0001
## 4 - 12 -12.970 0.853 48 -15.200 <.0001
## 4 - 13 -15.200 0.853 48 -17.814 <.0001
## 4 - 14 -17.218 0.853 48 -20.178 <.0001
## 4 - 15 -19.070 0.853 48 -22.349 <.0001
## 4 - 16 -20.780 0.853 48 -24.353 <.0001
## 4 - 17 -22.290 0.853 48 -26.123 <.0001
## 4 - 18 -23.600 0.853 48 -27.658 <.0001
## 4 - 19 -24.740 0.853 48 -28.994 <.0001
## 4 - 20 -25.740 0.853 48 -30.166 <.0001
## 4 - 21 -26.663 0.853 48 -31.247 <.0001
## 4 - 22 -27.523 0.853 48 -32.255 <.0001
## 4 - 23 -28.325 0.853 48 -33.196 <.0001
## 4 - 24 -29.078 0.853 48 -34.078 <.0001
## 5 - 6 -1.308 0.853 48 -1.532 0.9962
## 5 - 7 -2.455 0.853 48 -2.877 0.4103
## 5 - 8 -3.510 0.853 48 -4.114 0.0247
## 5 - 9 -4.973 0.853 48 -5.828 0.0001
## 5 - 10 -6.835 0.853 48 -8.010 <.0001
## 5 - 11 -8.998 0.853 48 -10.545 <.0001
## 5 - 12 -11.383 0.853 48 -13.340 <.0001
## 5 - 13 -13.613 0.853 48 -15.953 <.0001
## 5 - 14 -15.630 0.853 48 -18.318 <.0001
## 5 - 15 -17.483 0.853 48 -20.489 <.0001
## 5 - 16 -19.192 0.853 48 -22.493 <.0001
## 5 - 17 -20.703 0.853 48 -24.263 <.0001
## 5 - 18 -22.012 0.853 48 -25.798 <.0001
## 5 - 19 -23.152 0.853 48 -27.134 <.0001
## 5 - 20 -24.152 0.853 48 -28.306 <.0001
## 5 - 21 -25.075 0.853 48 -29.387 <.0001
## 5 - 22 -25.935 0.853 48 -30.395 <.0001
## 5 - 23 -26.738 0.853 48 -31.335 <.0001
## 5 - 24 -27.490 0.853 48 -32.217 <.0001
## 6 - 7 -1.147 0.853 48 -1.345 0.9994
## 6 - 8 -2.203 0.853 48 -2.581 0.6142
## 6 - 9 -3.665 0.853 48 -4.295 0.0147
## 6 - 10 -5.527 0.853 48 -6.478 <.0001
## 6 - 11 -7.690 0.853 48 -9.012 <.0001
## 6 - 12 -10.075 0.853 48 -11.808 <.0001
## 6 - 13 -12.305 0.853 48 -14.421 <.0001
## 6 - 14 -14.322 0.853 48 -16.785 <.0001
## 6 - 15 -16.175 0.853 48 -18.956 <.0001
## 6 - 16 -17.885 0.853 48 -20.961 <.0001
## 6 - 17 -19.395 0.853 48 -22.730 <.0001
## 6 - 18 -20.705 0.853 48 -24.266 <.0001
## 6 - 19 -21.845 0.853 48 -25.602 <.0001
## 6 - 20 -22.845 0.853 48 -26.773 <.0001
## 6 - 21 -23.767 0.853 48 -27.855 <.0001
## 6 - 22 -24.628 0.853 48 -28.863 <.0001
## 6 - 23 -25.430 0.853 48 -29.803 <.0001
## 6 - 24 -26.183 0.853 48 -30.685 <.0001
## 7 - 8 -1.055 0.853 48 -1.236 0.9998
## 7 - 9 -2.518 0.853 48 -2.950 0.3642
## 7 - 10 -4.380 0.853 48 -5.133 0.0011
## 7 - 11 -6.543 0.853 48 -7.668 <.0001
## 7 - 12 -8.928 0.853 48 -10.463 <.0001
## 7 - 13 -11.158 0.853 48 -13.076 <.0001
## 7 - 14 -13.175 0.853 48 -15.441 <.0001
## 7 - 15 -15.027 0.853 48 -17.612 <.0001
## 7 - 16 -16.738 0.853 48 -19.616 <.0001
## 7 - 17 -18.247 0.853 48 -21.385 <.0001
## 7 - 18 -19.558 0.853 48 -22.921 <.0001
## 7 - 19 -20.698 0.853 48 -24.257 <.0001
## 7 - 20 -21.698 0.853 48 -25.429 <.0001
## 7 - 21 -22.620 0.853 48 -26.510 <.0001
## 7 - 22 -23.480 0.853 48 -27.518 <.0001
## 7 - 23 -24.282 0.853 48 -28.458 <.0001
## 7 - 24 -25.035 0.853 48 -29.340 <.0001
## 8 - 9 -1.462 0.853 48 -1.714 0.9852
## 8 - 10 -3.325 0.853 48 -3.897 0.0447
## 8 - 11 -5.487 0.853 48 -6.431 <.0001
## 8 - 12 -7.872 0.853 48 -9.226 <.0001
## 8 - 13 -10.102 0.853 48 -11.840 <.0001
## 8 - 14 -12.120 0.853 48 -14.204 <.0001
## 8 - 15 -13.973 0.853 48 -16.375 <.0001
## 8 - 16 -15.682 0.853 48 -18.379 <.0001
## 8 - 17 -17.192 0.853 48 -20.149 <.0001
## 8 - 18 -18.503 0.853 48 -21.684 <.0001
## 8 - 19 -19.642 0.853 48 -23.020 <.0001
## 8 - 20 -20.642 0.853 48 -24.192 <.0001
## 8 - 21 -21.565 0.853 48 -25.273 <.0001
## 8 - 22 -22.425 0.853 48 -26.281 <.0001
## 8 - 23 -23.227 0.853 48 -27.222 <.0001
## 8 - 24 -23.980 0.853 48 -28.104 <.0001
## 9 - 10 -1.863 0.853 48 -2.183 0.8591
## 9 - 11 -4.025 0.853 48 -4.717 0.0041
## 9 - 12 -6.410 0.853 48 -7.512 <.0001
## 9 - 13 -8.640 0.853 48 -10.126 <.0001
## 9 - 14 -10.658 0.853 48 -12.490 <.0001
## 9 - 15 -12.510 0.853 48 -14.661 <.0001
## 9 - 16 -14.220 0.853 48 -16.665 <.0001
## 9 - 17 -15.730 0.853 48 -18.435 <.0001
## 9 - 18 -17.040 0.853 48 -19.970 <.0001
## 9 - 19 -18.180 0.853 48 -21.306 <.0001
## 9 - 20 -19.180 0.853 48 -22.478 <.0001
## 9 - 21 -20.102 0.853 48 -23.559 <.0001
## 9 - 22 -20.962 0.853 48 -24.567 <.0001
## 9 - 23 -21.765 0.853 48 -25.508 <.0001
## 9 - 24 -22.517 0.853 48 -26.390 <.0001
## 10 - 11 -2.163 0.853 48 -2.534 0.6469
## 10 - 12 -4.548 0.853 48 -5.330 0.0006
## 10 - 13 -6.777 0.853 48 -7.943 <.0001
## 10 - 14 -8.795 0.853 48 -10.307 <.0001
## 10 - 15 -10.648 0.853 48 -12.478 <.0001
## 10 - 16 -12.357 0.853 48 -14.483 <.0001
## 10 - 17 -13.867 0.853 48 -16.252 <.0001
## 10 - 18 -15.178 0.853 48 -17.787 <.0001
## 10 - 19 -16.317 0.853 48 -19.124 <.0001
## 10 - 20 -17.317 0.853 48 -20.295 <.0001
## 10 - 21 -18.240 0.853 48 -21.377 <.0001
## 10 - 22 -19.100 0.853 48 -22.384 <.0001
## 10 - 23 -19.902 0.853 48 -23.325 <.0001
## 10 - 24 -20.655 0.853 48 -24.207 <.0001
## 11 - 12 -2.385 0.853 48 -2.795 0.4648
## 11 - 13 -4.615 0.853 48 -5.409 0.0004
## 11 - 14 -6.633 0.853 48 -7.773 <.0001
## 11 - 15 -8.485 0.853 48 -9.944 <.0001
## 11 - 16 -10.195 0.853 48 -11.948 <.0001
## 11 - 17 -11.705 0.853 48 -13.718 <.0001
## 11 - 18 -13.015 0.853 48 -15.253 <.0001
## 11 - 19 -14.155 0.853 48 -16.589 <.0001
## 11 - 20 -15.155 0.853 48 -17.761 <.0001
## 11 - 21 -16.078 0.853 48 -18.842 <.0001
## 11 - 22 -16.938 0.853 48 -19.850 <.0001
## 11 - 23 -17.740 0.853 48 -20.791 <.0001
## 11 - 24 -18.492 0.853 48 -21.673 <.0001
## 12 - 13 -2.230 0.853 48 -2.613 0.5915
## 12 - 14 -4.247 0.853 48 -4.978 0.0018
## 12 - 15 -6.100 0.853 48 -7.149 <.0001
## 12 - 16 -7.810 0.853 48 -9.153 <.0001
## 12 - 17 -9.320 0.853 48 -10.923 <.0001
## 12 - 18 -10.630 0.853 48 -12.458 <.0001
## 12 - 19 -11.770 0.853 48 -13.794 <.0001
## 12 - 20 -12.770 0.853 48 -14.966 <.0001
## 12 - 21 -13.693 0.853 48 -16.047 <.0001
## 12 - 22 -14.553 0.853 48 -17.055 <.0001
## 12 - 23 -15.355 0.853 48 -17.995 <.0001
## 12 - 24 -16.108 0.853 48 -18.877 <.0001
## 13 - 14 -2.018 0.853 48 -2.364 0.7592
## 13 - 15 -3.870 0.853 48 -4.536 0.0072
## 13 - 16 -5.580 0.853 48 -6.540 <.0001
## 13 - 17 -7.090 0.853 48 -8.309 <.0001
## 13 - 18 -8.400 0.853 48 -9.844 <.0001
## 13 - 19 -9.540 0.853 48 -11.181 <.0001
## 13 - 20 -10.540 0.853 48 -12.352 <.0001
## 13 - 21 -11.463 0.853 48 -13.434 <.0001
## 13 - 22 -12.322 0.853 48 -14.442 <.0001
## 13 - 23 -13.125 0.853 48 -15.382 <.0001
## 13 - 24 -13.877 0.853 48 -16.264 <.0001
## 14 - 15 -1.853 0.853 48 -2.171 0.8646
## 14 - 16 -3.562 0.853 48 -4.175 0.0208
## 14 - 17 -5.072 0.853 48 -5.945 0.0001
## 14 - 18 -6.383 0.853 48 -7.480 <.0001
## 14 - 19 -7.522 0.853 48 -8.816 <.0001
## 14 - 20 -8.523 0.853 48 -9.988 <.0001
## 14 - 21 -9.445 0.853 48 -11.069 <.0001
## 14 - 22 -10.305 0.853 48 -12.077 <.0001
## 14 - 23 -11.107 0.853 48 -13.018 <.0001
## 14 - 24 -11.860 0.853 48 -13.899 <.0001
## 15 - 16 -1.710 0.853 48 -2.004 0.9292
## 15 - 17 -3.220 0.853 48 -3.774 0.0616
## 15 - 18 -4.530 0.853 48 -5.309 0.0006
## 15 - 19 -5.670 0.853 48 -6.645 <.0001
## 15 - 20 -6.670 0.853 48 -7.817 <.0001
## 15 - 21 -7.593 0.853 48 -8.898 <.0001
## 15 - 22 -8.453 0.853 48 -9.906 <.0001
## 15 - 23 -9.255 0.853 48 -10.847 <.0001
## 15 - 24 -10.008 0.853 48 -11.728 <.0001
## 16 - 17 -1.510 0.853 48 -1.770 0.9789
## 16 - 18 -2.820 0.853 48 -3.305 0.1855
## 16 - 19 -3.960 0.853 48 -4.641 0.0052
## 16 - 20 -4.960 0.853 48 -5.813 0.0001
## 16 - 21 -5.883 0.853 48 -6.894 <.0001
## 16 - 22 -6.742 0.853 48 -7.902 <.0001
## 16 - 23 -7.545 0.853 48 -8.842 <.0001
## 16 - 24 -8.297 0.853 48 -9.724 <.0001
## 17 - 18 -1.310 0.853 48 -1.535 0.9961
## 17 - 19 -2.450 0.853 48 -2.871 0.4141
## 17 - 20 -3.450 0.853 48 -4.043 0.0301
## 17 - 21 -4.372 0.853 48 -5.124 0.0011
## 17 - 22 -5.232 0.853 48 -6.132 <.0001
## 17 - 23 -6.035 0.853 48 -7.073 <.0001
## 17 - 24 -6.787 0.853 48 -7.955 <.0001
## 18 - 19 -1.140 0.853 48 -1.336 0.9994
## 18 - 20 -2.140 0.853 48 -2.508 0.6651
## 18 - 21 -3.062 0.853 48 -3.589 0.0975
## 18 - 22 -3.922 0.853 48 -4.597 0.0060
## 18 - 23 -4.725 0.853 48 -5.538 0.0003
## 18 - 24 -5.478 0.853 48 -6.419 <.0001
## 19 - 20 -1.000 0.853 48 -1.172 0.9999
## 19 - 21 -1.923 0.853 48 -2.253 0.8235
## 19 - 22 -2.783 0.853 48 -3.261 0.2033
## 19 - 23 -3.585 0.853 48 -4.201 0.0193
## 19 - 24 -4.338 0.853 48 -5.083 0.0013
## 20 - 21 -0.922 0.853 48 -1.081 1.0000
## 20 - 22 -1.782 0.853 48 -2.089 0.8996
## 20 - 23 -2.585 0.853 48 -3.030 0.3176
## 20 - 24 -3.337 0.853 48 -3.911 0.0430
## 21 - 22 -0.860 0.853 48 -1.008 1.0000
## 21 - 23 -1.663 0.853 48 -1.948 0.9450
## 21 - 24 -2.415 0.853 48 -2.830 0.4411
## 22 - 23 -0.802 0.853 48 -0.941 1.0000
## 22 - 24 -1.555 0.853 48 -1.822 0.9713
## 23 - 24 -0.752 0.853 48 -0.882 1.0000
##
## Results are averaged over the levels of: ScatF
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 24 estimates
emmeans(scatmodel2, list(pairwise~ScatF|TimeF), adjust="bonferroni")
## $`emmeans of ScatF | TimeF`
## TimeF = 1:
## ScatF emmean SE df lower.CL upper.CL
## A 3.22 0.892 57.3 1.434 5.01
## P 1.38 0.892 57.3 -0.406 3.17
##
## TimeF = 2:
## ScatF emmean SE df lower.CL upper.CL
## A 5.12 0.892 57.3 3.334 6.91
## P 1.97 0.892 57.3 0.184 3.76
##
## TimeF = 3:
## ScatF emmean SE df lower.CL upper.CL
## A 7.49 0.892 57.3 5.704 9.28
## P 3.08 0.892 57.3 1.289 4.86
##
## TimeF = 4:
## ScatF emmean SE df lower.CL upper.CL
## A 10.07 0.892 57.3 8.289 11.86
## P 4.43 0.892 57.3 2.649 6.22
##
## TimeF = 5:
## ScatF emmean SE df lower.CL upper.CL
## A 11.74 0.892 57.3 9.959 13.53
## P 5.94 0.892 57.3 4.154 7.73
##
## TimeF = 6:
## ScatF emmean SE df lower.CL upper.CL
## A 12.69 0.892 57.3 10.899 14.47
## P 7.62 0.892 57.3 5.829 9.40
##
## TimeF = 7:
## ScatF emmean SE df lower.CL upper.CL
## A 13.21 0.892 57.3 11.424 15.00
## P 9.38 0.892 57.3 7.599 11.17
##
## TimeF = 8:
## ScatF emmean SE df lower.CL upper.CL
## A 13.48 0.892 57.3 11.694 15.27
## P 11.22 0.892 57.3 9.439 13.01
##
## TimeF = 9:
## ScatF emmean SE df lower.CL upper.CL
## A 14.15 0.892 57.3 12.364 15.94
## P 13.48 0.892 57.3 11.694 15.27
##
## TimeF = 10:
## ScatF emmean SE df lower.CL upper.CL
## A 15.26 0.892 57.3 13.474 17.05
## P 16.09 0.892 57.3 14.309 17.88
##
## TimeF = 11:
## ScatF emmean SE df lower.CL upper.CL
## A 16.70 0.892 57.3 14.919 18.49
## P 18.98 0.892 57.3 17.189 20.76
##
## TimeF = 12:
## ScatF emmean SE df lower.CL upper.CL
## A 18.40 0.892 57.3 16.614 20.19
## P 22.05 0.892 57.3 20.264 23.84
##
## TimeF = 13:
## ScatF emmean SE df lower.CL upper.CL
## A 19.91 0.892 57.3 18.124 21.70
## P 25.00 0.892 57.3 23.214 26.79
##
## TimeF = 14:
## ScatF emmean SE df lower.CL upper.CL
## A 21.16 0.892 57.3 19.379 22.95
## P 27.78 0.892 57.3 25.994 29.57
##
## TimeF = 15:
## ScatF emmean SE df lower.CL upper.CL
## A 22.21 0.892 57.3 20.429 24.00
## P 30.43 0.892 57.3 28.649 32.22
##
## TimeF = 16:
## ScatF emmean SE df lower.CL upper.CL
## A 23.10 0.892 57.3 21.314 24.89
## P 32.97 0.892 57.3 31.184 34.76
##
## TimeF = 17:
## ScatF emmean SE df lower.CL upper.CL
## A 23.86 0.892 57.3 22.079 25.65
## P 35.23 0.892 57.3 33.439 37.01
##
## TimeF = 18:
## ScatF emmean SE df lower.CL upper.CL
## A 24.53 0.892 57.3 22.744 26.32
## P 37.18 0.892 57.3 35.394 38.97
##
## TimeF = 19:
## ScatF emmean SE df lower.CL upper.CL
## A 25.11 0.892 57.3 23.324 26.90
## P 38.88 0.892 57.3 37.094 40.67
##
## TimeF = 20:
## ScatF emmean SE df lower.CL upper.CL
## A 25.61 0.892 57.3 23.829 27.40
## P 40.38 0.892 57.3 38.589 42.16
##
## TimeF = 21:
## ScatF emmean SE df lower.CL upper.CL
## A 26.11 0.892 57.3 24.329 27.90
## P 41.72 0.892 57.3 39.934 43.51
##
## TimeF = 22:
## ScatF emmean SE df lower.CL upper.CL
## A 26.63 0.892 57.3 24.844 28.42
## P 42.92 0.892 57.3 41.139 44.71
##
## TimeF = 23:
## ScatF emmean SE df lower.CL upper.CL
## A 27.16 0.892 57.3 25.369 28.94
## P 44.01 0.892 57.3 42.219 45.79
##
## TimeF = 24:
## ScatF emmean SE df lower.CL upper.CL
## A 27.69 0.892 57.3 25.904 29.48
## P 44.98 0.892 57.3 43.189 46.76
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $`pairwise differences of ScatF | TimeF`
## TimeF = 1:
## 2 estimate SE df t.ratio p.value
## A - P 1.840 1.26 57.3 1.459 0.1501
##
## TimeF = 2:
## 2 estimate SE df t.ratio p.value
## A - P 3.150 1.26 57.3 2.497 0.0154
##
## TimeF = 3:
## 2 estimate SE df t.ratio p.value
## A - P 4.415 1.26 57.3 3.500 0.0009
##
## TimeF = 4:
## 2 estimate SE df t.ratio p.value
## A - P 5.640 1.26 57.3 4.472 <.0001
##
## TimeF = 5:
## 2 estimate SE df t.ratio p.value
## A - P 5.805 1.26 57.3 4.602 <.0001
##
## TimeF = 6:
## 2 estimate SE df t.ratio p.value
## A - P 5.070 1.26 57.3 4.020 0.0002
##
## TimeF = 7:
## 2 estimate SE df t.ratio p.value
## A - P 3.825 1.26 57.3 3.033 0.0036
##
## TimeF = 8:
## 2 estimate SE df t.ratio p.value
## A - P 2.255 1.26 57.3 1.788 0.0791
##
## TimeF = 9:
## 2 estimate SE df t.ratio p.value
## A - P 0.670 1.26 57.3 0.531 0.5973
##
## TimeF = 10:
## 2 estimate SE df t.ratio p.value
## A - P -0.835 1.26 57.3 -0.662 0.5106
##
## TimeF = 11:
## 2 estimate SE df t.ratio p.value
## A - P -2.270 1.26 57.3 -1.800 0.0772
##
## TimeF = 12:
## 2 estimate SE df t.ratio p.value
## A - P -3.650 1.26 57.3 -2.894 0.0054
##
## TimeF = 13:
## 2 estimate SE df t.ratio p.value
## A - P -5.090 1.26 57.3 -4.036 0.0002
##
## TimeF = 14:
## 2 estimate SE df t.ratio p.value
## A - P -6.615 1.26 57.3 -5.245 <.0001
##
## TimeF = 15:
## 2 estimate SE df t.ratio p.value
## A - P -8.220 1.26 57.3 -6.517 <.0001
##
## TimeF = 16:
## 2 estimate SE df t.ratio p.value
## A - P -9.870 1.26 57.3 -7.825 <.0001
##
## TimeF = 17:
## 2 estimate SE df t.ratio p.value
## A - P -11.360 1.26 57.3 -9.007 <.0001
##
## TimeF = 18:
## 2 estimate SE df t.ratio p.value
## A - P -12.650 1.26 57.3 -10.029 <.0001
##
## TimeF = 19:
## 2 estimate SE df t.ratio p.value
## A - P -13.770 1.26 57.3 -10.917 <.0001
##
## TimeF = 20:
## 2 estimate SE df t.ratio p.value
## A - P -14.760 1.26 57.3 -11.702 <.0001
##
## TimeF = 21:
## 2 estimate SE df t.ratio p.value
## A - P -15.605 1.26 57.3 -12.372 <.0001
##
## TimeF = 22:
## 2 estimate SE df t.ratio p.value
## A - P -16.295 1.26 57.3 -12.919 <.0001
##
## TimeF = 23:
## 2 estimate SE df t.ratio p.value
## A - P -16.850 1.26 57.3 -13.359 <.0001
##
## TimeF = 24:
## 2 estimate SE df t.ratio p.value
## A - P -17.285 1.26 57.3 -13.704 <.0001
##
## Degrees-of-freedom method: kenward-roger