Antoine Jardin, dans un post de 2015, compare les performances du package R lme4 et MixedModels écrit en Julia. La version Julia peut être plus beaucoup rapide.

Douglas Bates, ancien développeur de lme4, et créateur du
package Julia MixedModels montre dans ce post comment appeler ce package depuis R.

Pour apprendre à utiliser MixedModels il y a la documentation et une collection de notebooks très détaillés.

Pour exécuter la suite vous devez installer Julia.

Bibliothèque lme4 pour les modèles mixtes en R

## Loading required package: Matrix
## 
## Attaching package: 'Matrix'
## The following object is masked from 'package:tidyr':
## 
##     expand

Données sur les rats fournies par Pierre Barbillon

##    user  system elapsed 
##  13.712   0.014  13.729
## Linear mixed model fit by REML ['lmerMod']
## Formula: gdp ~ tpsF * gest * ps + (1 | bbF)
##    Data: don
## 
## REML criterion at convergence: 23617.7
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -4.0114 -0.4950  0.0112  0.4703  3.9049 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev.
##  bbF      (Intercept) 210.90   14.522  
##  Residual              71.64    8.464  
## Number of obs: 3619, groups:  bbF, 47
## 
## Fixed effects:
##                     Estimate Std. Error t value
## (Intercept)          3.86042    5.94281   0.650
## tpsF2                5.80625    4.23200   1.372
## tpsF3               11.71458    4.23200   2.768
## tpsF4               17.21458    4.23200   4.068
## tpsF5               23.32708    4.23200   5.512
## tpsF6               28.63958    4.23200   6.767
## tpsF7               33.78958    4.23200   7.984
## tpsF8               38.46667    4.23200   9.089
## tpsF9               42.89292    4.23200  10.135
## tpsF10              47.00208    4.23200  11.106
## tpsF11              50.53750    4.23200  11.942
## tpsF12              50.86875    4.23200  12.020
## tpsF13              55.15208    4.23200  13.032
## tpsF14              61.28958    4.23200  14.482
## tpsF15              65.87812    4.23200  15.567
## tpsF16              70.90083    4.23200  16.754
## tpsF17              72.87708    4.23200  17.220
## tpsF18              75.85208    4.23200  17.923
## tpsF19              81.50208    4.23200  19.259
## tpsF20              83.91458    4.23200  19.829
## tpsF21              89.80208    4.23200  21.220
## tpsF22              93.54375    4.23200  22.104
## tpsF23              97.79625    4.23200  23.109
## tpsF24              99.36458    4.23200  23.479
## tpsF25             101.21458    4.23200  23.916
## tpsF26             106.08958    4.23200  25.068
## tpsF27             106.70208    4.23200  25.213
## tpsF28             111.03958    4.23200  26.238
## tpsF29             113.41042    4.23200  26.798
## tpsF30             115.90167    4.23200  27.387
## tpsF31             119.85208    4.23200  28.320
## tpsF32             119.78958    4.23200  28.306
## tpsF33             122.73958    4.23200  29.003
## tpsF34             123.45208    4.23200  29.171
## tpsF35             128.93958    4.23200  30.468
## tpsF36             132.84792    4.23200  31.391
## tpsF37             135.46667    4.23200  32.010
## tpsF38             132.62708    4.23200  31.339
## tpsF39             133.83958    4.23200  31.626
## tpsF40             136.05208    4.23200  32.148
## tpsF41             136.12708    4.23200  32.166
## tpsF42             137.50208    4.23200  32.491
## tpsF43             142.59583    4.23200  33.695
## tpsF44             144.06458    4.23200  34.042
## tpsF45             145.92708    4.23200  34.482
## tpsF46             144.68958    4.23200  34.189
## tpsF47             148.25208    4.23200  35.031
## tpsF48             145.50208    4.23200  34.381
## tpsF49             149.75208    4.23200  35.386
## tpsF50             146.96458    4.23200  34.727
## tpsF51             149.46458    4.23200  35.318
## tpsF52             153.06458    4.23200  36.168
## tpsF53             152.26458    4.23200  35.979
## tpsF54             158.06458    4.23200  37.350
## tpsF55             156.62708    4.23200  37.010
## tpsF56             160.22708    4.23200  37.861
## tpsF57             164.57500    4.23200  38.888
## tpsF58             164.14375    4.23200  38.786
## tpsF59             164.18958    4.23200  38.797
## tpsF60             162.32708    4.23200  38.357
## tpsF61             163.32708    4.23200  38.593
## tpsF62             162.86458    4.23200  38.484
## tpsF63             166.87708    4.23200  39.432
## tpsF64             167.09792    4.23200  39.484
## tpsF65             168.13542    4.23200  39.730
## tpsF66             170.02708    4.23200  40.177
## tpsF67             168.07708    4.23200  39.716
## tpsF68             171.32708    4.23200  40.484
## tpsF69             169.03958    4.23200  39.943
## tpsF70             172.93958    4.23200  40.865
## tpsF71             171.63958    4.23200  40.558
## tpsF72             172.76458    4.23200  40.823
## tpsF73             175.13958    4.23200  41.385
## tpsF74             174.82708    4.23200  41.311
## tpsF75             180.36458    4.23200  42.619
## tpsF76             177.12708    4.23200  41.854
## tpsF77             179.71458    4.23200  42.466
## gestNP               0.47917    8.40440   0.057
## psC2                 1.82500    8.40440   0.217
## psNP                 2.78006    8.69938   0.320
## tpsF2:gestNP        -0.08750    5.98495  -0.015
## tpsF3:gestNP        -0.10417    5.98495  -0.017
## tpsF4:gestNP        -0.40417    5.98495  -0.068
## tpsF5:gestNP        -0.27917    5.98495  -0.047
## tpsF6:gestNP        -0.72917    5.98495  -0.122
## tpsF7:gestNP        -0.47917    5.98495  -0.080
## tpsF8:gestNP        -0.10000    5.98495  -0.017
## tpsF9:gestNP         0.61375    5.98495   0.103
## tpsF10:gestNP        0.77083    5.98495   0.129
## tpsF11:gestNP       -0.43542    5.98495  -0.073
## tpsF12:gestNP        1.94167    5.98495   0.324
## tpsF13:gestNP        0.79583    5.98495   0.133
## tpsF14:gestNP        1.30833    5.98495   0.219
## tpsF15:gestNP        0.56354    5.98495   0.094
## tpsF16:gestNP        0.15583    5.98495   0.026
## tpsF17:gestNP        0.74583    5.98495   0.125
## tpsF18:gestNP       -0.02917    5.98495  -0.005
## tpsF19:gestNP        0.25833    5.98495   0.043
## tpsF20:gestNP        0.44583    5.98495   0.074
## tpsF21:gestNP        0.80833    5.98495   0.135
## tpsF22:gestNP        1.16667    5.98495   0.195
## tpsF23:gestNP        1.50917    5.98495   0.252
## tpsF24:gestNP        0.59583    5.98495   0.100
## tpsF25:gestNP        0.05833    5.98495   0.010
## tpsF26:gestNP        1.65833    5.98495   0.277
## tpsF27:gestNP        0.97083    5.98495   0.162
## tpsF28:gestNP        1.50833    5.98495   0.252
## tpsF29:gestNP        2.06250    5.98495   0.345
## tpsF30:gestNP        2.83625    5.98495   0.474
## tpsF31:gestNP        1.45833    5.98495   0.244
## tpsF32:gestNP        1.65833    5.98495   0.277
## tpsF33:gestNP        2.98333    5.98495   0.498
## tpsF34:gestNP        1.65833    5.98495   0.277
## tpsF35:gestNP        0.95833    5.98495   0.160
## tpsF36:gestNP       -0.77292    5.98495  -0.129
## tpsF37:gestNP       -0.91667    5.98495  -0.153
## tpsF38:gestNP        4.78333    5.98495   0.799
## tpsF39:gestNP        2.42083    5.98495   0.404
## tpsF40:gestNP        3.45833    5.98495   0.578
## tpsF41:gestNP        4.65833    5.98495   0.778
## tpsF42:gestNP        6.23333    5.98495   1.042
## tpsF43:gestNP        1.92708    5.98495   0.322
## tpsF44:gestNP        2.79583    5.98495   0.467
## tpsF45:gestNP        3.53333    5.98495   0.590
## tpsF46:gestNP        5.63333    5.98495   0.941
## tpsF47:gestNP        5.99583    5.98495   1.002
## tpsF48:gestNP        6.09583    5.98495   1.019
## tpsF49:gestNP        4.25833    5.98495   0.712
## tpsF50:gestNP       10.95417    5.98495   1.830
## tpsF51:gestNP        9.25417    5.98495   1.546
## tpsF52:gestNP        7.48333    5.98495   1.250
## tpsF53:gestNP        5.99583    5.98495   1.002
## tpsF54:gestNP        4.08333    5.98495   0.682
## tpsF55:gestNP        3.69583    5.98495   0.618
## tpsF56:gestNP        6.25833    5.98495   1.046
## tpsF57:gestNP        0.23333    5.98495   0.039
## tpsF58:gestNP        2.32083    5.98495   0.388
## tpsF59:gestNP        4.92083    5.98495   0.822
## tpsF60:gestNP        5.47083    5.98495   0.914
## tpsF61:gestNP        9.09583    5.98495   1.520
## tpsF62:gestNP        6.63333    5.98495   1.108
## tpsF63:gestNP        6.37083    5.98495   1.064
## tpsF64:gestNP        6.75000    5.98495   1.128
## tpsF65:gestNP        6.51250    5.98495   1.088
## tpsF66:gestNP        6.53333    5.98495   1.092
## tpsF67:gestNP        7.17083    5.98495   1.198
## tpsF68:gestNP        7.97083    5.98495   1.332
## tpsF69:gestNP       10.95833    5.98495   1.831
## tpsF70:gestNP       10.58333    5.98495   1.768
## tpsF71:gestNP       11.47542    5.98495   1.917
## tpsF72:gestNP       11.43417    5.98495   1.910
## tpsF73:gestNP       11.05833    5.98495   1.848
## tpsF74:gestNP        9.97083    5.98495   1.666
## tpsF75:gestNP        7.74583    5.98495   1.294
## tpsF76:gestNP        9.32083    5.98495   1.557
## tpsF77:gestNP        7.52083    5.98495   1.257
## tpsF2:psC2          -0.20000    5.98495  -0.033
## tpsF3:psC2          -0.58750    5.98495  -0.098
## tpsF4:psC2          -0.22500    5.98495  -0.038
## tpsF5:psC2          -0.98750    5.98495  -0.165
## tpsF6:psC2          -1.05000    5.98495  -0.175
## tpsF7:psC2          -0.15000    5.98495  -0.025
## tpsF8:psC2           0.67500    5.98495   0.113
## tpsF9:psC2           1.89000    5.98495   0.316
## tpsF10:psC2          2.80000    5.98495   0.468
## tpsF11:psC2         -0.28542    5.98495  -0.048
## tpsF12:psC2          0.42917    5.98495   0.072
## tpsF13:psC2          0.57500    5.98495   0.096
## tpsF14:psC2          2.78750    5.98495   0.466
## tpsF15:psC2          4.01354    5.98495   0.671
## tpsF16:psC2          5.42750    5.98495   0.907
## tpsF17:psC2          2.56250    5.98495   0.428
## tpsF18:psC2          2.06250    5.98495   0.345
## tpsF19:psC2          3.00000    5.98495   0.501
## tpsF20:psC2          1.95000    5.98495   0.326
## tpsF21:psC2          4.51250    5.98495   0.754
## tpsF22:psC2          4.48542    5.98495   0.749
## tpsF23:psC2          5.13292    5.98495   0.858
## tpsF24:psC2          4.73750    5.98495   0.792
## tpsF25:psC2          5.01250    5.98495   0.838
## tpsF26:psC2          5.87500    5.98495   0.982
## tpsF27:psC2          6.05000    5.98495   1.011
## tpsF28:psC2          9.36250    5.98495   1.564
## tpsF29:psC2          9.98542    5.98495   1.668
## tpsF30:psC2         11.26292    5.98495   1.882
## tpsF31:psC2         10.57500    5.98495   1.767
## tpsF32:psC2         10.38750    5.98495   1.736
## tpsF33:psC2         12.73750    5.98495   2.128
## tpsF34:psC2         11.66250    5.98495   1.949
## tpsF35:psC2         12.56250    5.98495   2.099
## tpsF36:psC2         10.90833    5.98495   1.823
## tpsF37:psC2         10.00833    5.98495   1.672
## tpsF38:psC2         15.96250    5.98495   2.667
## tpsF39:psC2         13.68750    5.98495   2.287
## tpsF40:psC2         15.28750    5.98495   2.554
## tpsF41:psC2         15.07500    5.98495   2.519
## tpsF42:psC2         16.65000    5.98495   2.782
## tpsF43:psC2         14.16042    5.98495   2.366
## tpsF44:psC2         15.10417    5.98495   2.524
## tpsF45:psC2         15.98750    5.98495   2.671
## tpsF46:psC2         16.55000    5.98495   2.765
## tpsF47:psC2         17.42500    5.98495   2.911
## tpsF48:psC2         15.72500    5.98495   2.627
## tpsF49:psC2         16.03750    5.98495   2.680
## tpsF50:psC2         20.84583    5.98495   3.483
## tpsF51:psC2         19.88333    5.98495   3.322
## tpsF52:psC2         18.50000    5.98495   3.091
## tpsF53:psC2         18.80000    5.98495   3.141
## tpsF54:psC2         16.57500    5.98495   2.769
## tpsF55:psC2         16.83750    5.98495   2.813
## tpsF56:psC2         18.81250    5.98495   3.143
## tpsF57:psC2         10.62083    5.98495   1.775
## tpsF58:psC2         13.89583    5.98495   2.322
## tpsF59:psC2         17.46250    5.98495   2.918
## tpsF60:psC2         19.86250    5.98495   3.319
## tpsF61:psC2         24.01250    5.98495   4.012
## tpsF62:psC2         21.90000    5.98495   3.659
## tpsF63:psC2         22.68750    5.98495   3.791
## tpsF64:psC2         21.61250    5.98495   3.611
## tpsF65:psC2         22.42500    5.98495   3.747
## tpsF66:psC2         23.78750    5.98495   3.975
## tpsF67:psC2         23.33750    5.98495   3.899
## tpsF68:psC2         26.08750    5.98495   4.359
## tpsF69:psC2         27.21250    5.98495   4.547
## tpsF70:psC2         25.78750    5.98495   4.309
## tpsF71:psC2         29.87125    5.98495   4.991
## tpsF72:psC2         29.28500    5.98495   4.893
## tpsF73:psC2         28.01250    5.98495   4.680
## tpsF74:psC2         26.43750    5.98495   4.417
## tpsF75:psC2         23.82500    5.98495   3.981
## tpsF76:psC2         25.23750    5.98495   4.217
## tpsF77:psC2         22.68750    5.98495   3.791
## tpsF2:psNP          -0.29911    6.19501  -0.048
## tpsF3:psNP          -0.54077    6.19501  -0.087
## tpsF4:psNP          -1.01220    6.19501  -0.163
## tpsF5:psNP          -1.13899    6.19501  -0.184
## tpsF6:psNP          -1.66577    6.19501  -0.269
## tpsF7:psNP          -1.34435    6.19501  -0.217
## tpsF8:psNP          -0.99524    6.19501  -0.161
## tpsF9:psNP          -0.28577    6.19501  -0.046
## tpsF10:psNP         -0.38542    6.19501  -0.062
## tpsF11:psNP         -0.45417    6.19501  -0.073
## tpsF12:psNP         -0.98065    6.19501  -0.158
## tpsF13:psNP         -0.67827    6.19501  -0.109
## tpsF14:psNP         -0.25863    6.19501  -0.042
## tpsF15:psNP         -0.08765    6.19501  -0.014
## tpsF16:psNP          0.17107    6.19501   0.028
## tpsF17:psNP         -0.04613    6.19501  -0.007
## tpsF18:psNP          0.30744    6.19501   0.050
## tpsF19:psNP         -0.05685    6.19501  -0.009
## tpsF20:psNP          1.03065    6.19501   0.166
## tpsF21:psNP          0.84315    6.19501   0.136
## tpsF22:psNP          1.66815    6.19501   0.269
## tpsF23:psNP          2.21565    6.19501   0.358
## tpsF24:psNP          2.36637    6.19501   0.382
## tpsF25:psNP          3.87351    6.19501   0.625
## tpsF26:psNP          2.88423    6.19501   0.466
## tpsF27:psNP          2.22887    6.19501   0.360
## tpsF28:psNP          4.54851    6.19501   0.734
## tpsF29:psNP          4.69435    6.19501   0.758
## tpsF30:psNP          5.62881    6.19501   0.909
## tpsF31:psNP          2.96458    6.19501   0.479
## tpsF32:psNP          4.31280    6.19501   0.696
## tpsF33:psNP          7.14851    6.19501   1.154
## tpsF34:psNP          6.82173    6.19501   1.101
## tpsF35:psNP          4.44851    6.19501   0.718
## tpsF36:psNP          3.11875    6.19501   0.503
## tpsF37:psNP          2.63571    6.19501   0.425
## tpsF38:psNP          8.60387    6.19501   1.389
## tpsF39:psNP          6.97708    6.19501   1.126
## tpsF40:psNP          8.63601    6.19501   1.394
## tpsF41:psNP          8.67530    6.19501   1.400
## tpsF42:psNP          9.70030    6.19501   1.566
## tpsF43:psNP          7.39940    6.19501   1.194
## tpsF44:psNP          8.45208    6.19501   1.364
## tpsF45:psNP          9.47530    6.19501   1.530
## tpsF46:psNP          9.94137    6.19501   1.605
## tpsF47:psNP          8.00744    6.19501   1.293
## tpsF48:psNP          9.02887    6.19501   1.457
## tpsF49:psNP          9.16458    6.19501   1.479
## tpsF50:psNP         11.76875    6.19501   1.900
## tpsF51:psNP         11.09018    6.19501   1.790
## tpsF52:psNP         10.09494    6.19501   1.630
## tpsF53:psNP         10.36637    6.19501   1.673
## tpsF54:psNP          8.73780    6.19501   1.410
## tpsF55:psNP          9.41815    6.19501   1.520
## tpsF56:psNP          8.21815    6.19501   1.327
## tpsF57:psNP          5.84643    6.19501   0.944
## tpsF58:psNP          7.90625    6.19501   1.276
## tpsF59:psNP         10.24137    6.19501   1.653
## tpsF60:psNP         11.47530    6.19501   1.852
## tpsF61:psNP         14.36101    6.19501   2.318
## tpsF62:psNP         13.26637    6.19501   2.141
## tpsF63:psNP         13.63958    6.19501   2.202
## tpsF64:psNP         14.77351    6.19501   2.385
## tpsF65:psNP         14.92887    6.19501   2.410
## tpsF66:psNP         14.98958    6.19501   2.420
## tpsF67:psNP         13.88244    6.19501   2.241
## tpsF68:psNP         15.00387    6.19501   2.422
## tpsF69:psNP         16.20565    6.19501   2.616
## tpsF70:psNP         17.70565    6.19501   2.858
## tpsF71:psNP         21.89613    6.19501   3.534
## tpsF72:psNP         21.51399    6.19501   3.473
## tpsF73:psNP         20.31994    6.19501   3.280
## tpsF74:psNP         20.06101    6.19501   3.238
## tpsF75:psNP         16.39494    6.19501   2.646
## tpsF76:psNP         18.37530    6.19501   2.966
## tpsF77:psNP         13.84494    6.19501   2.235
## gestNP:psC2         -1.75208   11.88561  -0.147
## gestNP:psNP         -2.02173   12.09599  -0.167
## tpsF2:gestNP:psC2   -0.51875    8.46400  -0.061
## tpsF3:gestNP:psC2   -0.94792    8.46400  -0.112
## tpsF4:gestNP:psC2   -1.73542    8.46400  -0.205
## tpsF5:gestNP:psC2   -1.98542    8.46400  -0.235
## tpsF6:gestNP:psC2   -1.69792    8.46400  -0.201
## tpsF7:gestNP:psC2   -2.17292    8.46400  -0.257
## tpsF8:gestNP:psC2   -2.63542    8.46400  -0.311
## tpsF9:gestNP:psC2   -3.48292    8.46400  -0.411
## tpsF10:gestNP:psC2  -5.68542    8.46400  -0.672
## tpsF11:gestNP:psC2  -2.03958    8.46400  -0.241
## tpsF12:gestNP:psC2  -4.31319    8.46400  -0.510
## tpsF13:gestNP:psC2  -4.91042    8.46400  -0.580
## tpsF14:gestNP:psC2  -6.09792    8.46400  -0.720
## tpsF15:gestNP:psC2  -7.29583    8.46400  -0.862
## tpsF16:gestNP:psC2  -8.43056    8.46400  -0.996
## tpsF17:gestNP:psC2  -6.17292    8.46400  -0.729
## tpsF18:gestNP:psC2  -6.46042    8.46400  -0.763
## tpsF19:gestNP:psC2  -6.67292    8.46400  -0.788
## tpsF20:gestNP:psC2  -7.18542    8.46400  -0.849
## tpsF21:gestNP:psC2  -8.77292    8.46400  -1.036
## tpsF22:gestNP:psC2  -9.41458    8.46400  -1.112
## tpsF23:gestNP:psC2 -10.32708    8.46400  -1.220
## tpsF24:gestNP:psC2  -7.02292    8.46400  -0.830
## tpsF25:gestNP:psC2 -10.54792    8.46400  -1.246
## tpsF26:gestNP:psC2 -12.69792    8.46400  -1.500
## tpsF27:gestNP:psC2 -14.53542    8.46400  -1.717
## tpsF28:gestNP:psC2 -16.28542    8.46400  -1.924
## tpsF29:gestNP:psC2 -18.62917    8.46400  -2.201
## tpsF30:gestNP:psC2 -21.32042    8.46400  -2.519
## tpsF31:gestNP:psC2 -16.46042    8.46400  -1.945
## tpsF32:gestNP:psC2 -17.06042    8.46400  -2.016
## tpsF33:gestNP:psC2 -17.84792    8.46400  -2.109
## tpsF34:gestNP:psC2 -16.93542    8.46400  -2.001
## tpsF35:gestNP:psC2 -16.61042    8.46400  -1.962
## tpsF36:gestNP:psC2 -12.60417    8.46400  -1.489
## tpsF37:gestNP:psC2 -11.53542    8.46400  -1.363
## tpsF38:gestNP:psC2 -19.16042    8.46400  -2.264
## tpsF39:gestNP:psC2 -17.92292    8.46400  -2.118
## tpsF40:gestNP:psC2 -18.81042    8.46400  -2.222
## tpsF41:gestNP:psC2 -20.71042    8.46400  -2.447
## tpsF42:gestNP:psC2 -21.53542    8.46400  -2.544
## tpsF43:gestNP:psC2 -19.43750    8.46400  -2.296
## tpsF44:gestNP:psC2 -19.96875    8.46400  -2.359
## tpsF45:gestNP:psC2 -20.31042    8.46400  -2.400
## tpsF46:gestNP:psC2 -21.44792    8.46400  -2.534
## tpsF47:gestNP:psC2 -22.49792    8.46400  -2.658
## tpsF48:gestNP:psC2 -22.18542    8.46400  -2.621
## tpsF49:gestNP:psC2 -20.64792    8.46400  -2.439
## tpsF50:gestNP:psC2 -25.81875    8.46400  -3.050
## tpsF51:gestNP:psC2 -24.96875    8.46400  -2.950
## tpsF52:gestNP:psC2 -23.98542    8.46400  -2.834
## tpsF53:gestNP:psC2 -23.53542    8.46400  -2.781
## tpsF54:gestNP:psC2 -22.28542    8.46400  -2.633
## tpsF55:gestNP:psC2 -22.47292    8.46400  -2.655
## tpsF56:gestNP:psC2 -26.07292    8.46400  -3.080
## tpsF57:gestNP:psC2 -15.82708    8.46400  -1.870
## tpsF58:gestNP:psC2 -19.61458    8.46400  -2.317
## tpsF59:gestNP:psC2 -23.52292    8.46400  -2.779
## tpsF60:gestNP:psC2 -26.94792    8.46400  -3.184
## tpsF61:gestNP:psC2 -31.09792    8.46400  -3.674
## tpsF62:gestNP:psC2 -27.26042    8.46400  -3.221
## tpsF63:gestNP:psC2 -28.72292    8.46400  -3.394
## tpsF64:gestNP:psC2 -30.60208    8.46400  -3.616
## tpsF65:gestNP:psC2 -30.22708    8.46400  -3.571
## tpsF66:gestNP:psC2 -30.54792    8.46400  -3.609
## tpsF67:gestNP:psC2 -30.62292    8.46400  -3.618
## tpsF68:gestNP:psC2 -32.09792    8.46400  -3.792
## tpsF69:gestNP:psC2 -35.34792    8.46400  -4.176
## tpsF70:gestNP:psC2 -31.77292    8.46400  -3.754
## tpsF71:gestNP:psC2 -39.35333    8.46400  -4.649
## tpsF72:gestNP:psC2 -38.00958    8.46400  -4.491
## tpsF73:gestNP:psC2 -36.91042    8.46400  -4.361
## tpsF74:gestNP:psC2 -33.92292    8.46400  -4.008
## tpsF75:gestNP:psC2 -30.38542    8.46400  -3.590
## tpsF76:gestNP:psC2 -33.84792    8.46400  -3.999
## tpsF77:gestNP:psC2 -31.01042    8.46400  -3.664
## tpsF2:gestNP:psNP    0.09911    8.61382   0.012
## tpsF3:gestNP:psNP    0.20744    8.61382   0.024
## tpsF4:gestNP:psNP    0.27887    8.61382   0.032
## tpsF5:gestNP:psNP    0.40565    8.61382   0.047
## tpsF6:gestNP:psNP    0.65744    8.61382   0.076
## tpsF7:gestNP:psNP    1.23601    8.61382   0.143
## tpsF8:gestNP:psNP    0.91607    8.61382   0.106
## tpsF9:gestNP:psNP    0.46661    8.61382   0.054
## tpsF10:gestNP:psNP  -0.99792    8.61382  -0.116
## tpsF11:gestNP:psNP   1.50208    8.61382   0.174
## tpsF12:gestNP:psNP   0.06121    8.61382   0.007
## tpsF13:gestNP:psNP   0.21994    8.61382   0.026
## tpsF14:gestNP:psNP  -0.71220    8.61382  -0.083
## tpsF15:gestNP:psNP  -0.56964    8.61382  -0.066
## tpsF16:gestNP:psNP  -0.74788    8.61382  -0.087
## tpsF17:gestNP:psNP   0.10030    8.61382   0.012
## tpsF18:gestNP:psNP  -0.04077    8.61382  -0.005
## tpsF19:gestNP:psNP  -0.56399    8.61382  -0.065
## tpsF20:gestNP:psNP  -1.03899    8.61382  -0.121
## tpsF21:gestNP:psNP  -1.22649    8.61382  -0.142
## tpsF22:gestNP:psNP  -2.34940    8.61382  -0.273
## tpsF23:gestNP:psNP  -3.12690    8.61382  -0.363
## tpsF24:gestNP:psNP  -2.03720    8.61382  -0.237
## tpsF25:gestNP:psNP  -2.90685    8.61382  -0.337
## tpsF26:gestNP:psNP  -4.45506    8.61382  -0.517
## tpsF27:gestNP:psNP  -2.36220    8.61382  -0.274
## tpsF28:gestNP:psNP  -5.26935    8.61382  -0.612
## tpsF29:gestNP:psNP  -5.71310    8.61382  -0.663
## tpsF30:gestNP:psNP  -7.03631    8.61382  -0.817
## tpsF31:gestNP:psNP  -2.94792    8.61382  -0.342
## tpsF32:gestNP:psNP  -3.82113    8.61382  -0.444
## tpsF33:gestNP:psNP  -7.59435    8.61382  -0.882
## tpsF34:gestNP:psNP  -6.30506    8.61382  -0.732
## tpsF35:gestNP:psNP  -3.23185    8.61382  -0.375
## tpsF36:gestNP:psNP  -1.95833    8.61382  -0.227
## tpsF37:gestNP:psNP  -1.03780    8.61382  -0.120
## tpsF38:gestNP:psNP  -7.42470    8.61382  -0.862
## tpsF39:gestNP:psNP  -3.12292    8.61382  -0.363
## tpsF40:gestNP:psNP  -3.98185    8.61382  -0.462
## tpsF41:gestNP:psNP  -5.08363    8.61382  -0.590
## tpsF42:gestNP:psNP  -6.93363    8.61382  -0.805
## tpsF43:gestNP:psNP  -0.09315    8.61382  -0.011
## tpsF44:gestNP:psNP  -1.80208    8.61382  -0.209
## tpsF45:gestNP:psNP  -3.68363    8.61382  -0.428
## tpsF46:gestNP:psNP  -6.39970    8.61382  -0.743
## tpsF47:gestNP:psNP  -5.95327    8.61382  -0.691
## tpsF48:gestNP:psNP  -4.49970    8.61382  -0.522
## tpsF49:gestNP:psNP  -3.23542    8.61382  -0.376
## tpsF50:gestNP:psNP  -9.45208    8.61382  -1.097
## tpsF51:gestNP:psNP  -7.73601    8.61382  -0.898
## tpsF52:gestNP:psNP  -6.19077    8.61382  -0.719
## tpsF53:gestNP:psNP  -3.96220    8.61382  -0.460
## tpsF54:gestNP:psNP  -2.75863    8.61382  -0.320
## tpsF55:gestNP:psNP  -2.46399    8.61382  -0.286
## tpsF56:gestNP:psNP  -2.55149    8.61382  -0.296
## tpsF57:gestNP:psNP   2.90565    8.61382   0.337
## tpsF58:gestNP:psNP   0.76458    8.61382   0.089
## tpsF59:gestNP:psNP  -1.66220    8.61382  -0.193
## tpsF60:gestNP:psNP  -2.94613    8.61382  -0.342
## tpsF61:gestNP:psNP  -5.94435    8.61382  -0.690
## tpsF62:gestNP:psNP  -3.16220    8.61382  -0.367
## tpsF63:gestNP:psNP  -4.39792    8.61382  -0.511
## tpsF64:gestNP:psNP  -6.97351    8.61382  -0.810
## tpsF65:gestNP:psNP  -6.30387    8.61382  -0.732
## tpsF66:gestNP:psNP  -5.38542    8.61382  -0.625
## tpsF67:gestNP:psNP  -4.26577    8.61382  -0.495
## tpsF68:gestNP:psNP  -5.33720    8.61382  -0.620
## tpsF69:gestNP:psNP  -7.71399    8.61382  -0.896
## tpsF70:gestNP:psNP  -8.53899    8.61382  -0.991
## tpsF71:gestNP:psNP  -7.64655    8.61382  -0.888
## tpsF72:gestNP:psNP  -8.57565    8.61382  -0.996
## tpsF73:gestNP:psNP  -8.56577    8.61382  -0.994
## tpsF74:gestNP:psNP -10.28185    8.61382  -1.194
## tpsF75:gestNP:psNP  -8.92827    8.61382  -1.037
## tpsF76:gestNP:psNP  -8.03363    8.61382  -0.933
## tpsF77:gestNP:psNP  -4.41577    8.61382  -0.513
## 
## Correlation matrix not shown by default, as p = 462 > 12.
## Use print(x, correlation=TRUE)  or
##     vcov(x)        if you need it

Démarrer la session julia

Cette étape nécessite d’avoir installé Julia sur son système et d’avoir spécifié JULIA_HOME dans le .Renviron. Il faut également installer les packages Rcall et MixedModels dans Julia (comme indiqué sur la cette page). Parfois il faut également exécuter la commande Pkg.build("Arpack") dans Julia.

## Julia version 1.2.0 at location /home/metienne/install/julia-1.2.0/bin will be used.
## Loading setup script for JuliaCall...
## Finish loading setup script for JuliaCall.

Transférer les données et la formule du modèle vers la session Julia

Chargement du package MixedModels

Exécuter la fonction julia

## [1] 12.79208
## Julia Object of type LinearMixedModel{Float64}.
## Linear mixed model fit by maximum likelihood
##  gdp ~ 1 + tpsF + gest + ps + tpsF & gest + tpsF & ps + gest & ps + tpsF & gest & ps + (1 | bbF)
##      logLik        -2 logLik          AIC             BIC       
##  -1.27451008×10⁴  2.54902016×10⁴  2.64182016×10⁴  2.92921958×10⁴
## 
## Variance components:
##               Column    Variance  Std.Dev.  
##  bbF      (Intercept)  183.97351 13.5636837
##  Residual               62.49385  7.9053053
##  Number of obs: 3619; levels of grouping factors: 47
## 
##   Fixed-effects parameters:
## ──────────────────────────────────────────────────────────────────────────
##                                  Estimate  Std.Error      z value  P(>|z|)
## ──────────────────────────────────────────────────────────────────────────
## (Intercept)                     3.86042      5.55053   0.695504     0.4867
## tpsF: 2                         5.80625      3.95265   1.46895      0.1418
## tpsF: 3                        11.7146       3.95265   2.96373      0.0030
## tpsF: 4                        17.2146       3.95265   4.3552       <1e-4 
## tpsF: 5                        23.3271       3.95265   5.90163      <1e-8 
## tpsF: 6                        28.6396       3.95265   7.24566      <1e-12
## tpsF: 7                        33.7896       3.95265   8.54858      <1e-16
## tpsF: 8                        38.4667       3.95265   9.73186      <1e-21
## tpsF: 9                        42.8929       3.95265  10.8517       <1e-26
## tpsF: 10                       47.0021       3.95265  11.8913       <1e-31
## tpsF: 11                       50.5375       3.95265  12.7857       <1e-36
## tpsF: 12                       50.8687       3.95265  12.8695       <1e-37
## tpsF: 13                       55.1521       3.95265  13.9532       <1e-43
## tpsF: 14                       61.2896       3.95265  15.5059       <1e-53
## tpsF: 15                       65.8781       3.95265  16.6668       <1e-61
## tpsF: 16                       70.9008       3.95265  17.9375       <1e-71
## tpsF: 17                       72.8771       3.95265  18.4375       <1e-75
## tpsF: 18                       75.8521       3.95265  19.1902       <1e-81
## tpsF: 19                       81.5021       3.95265  20.6196       <1e-93
## tpsF: 20                       83.9146       3.95265  21.2299       <1e-99
## tpsF: 21                       89.8021       3.95265  22.7194       <1e-99
## tpsF: 22                       93.5438       3.95265  23.6661       <1e-99
## tpsF: 23                       97.7963       3.95265  24.7419       <1e-99
## tpsF: 24                       99.3646       3.95265  25.1387       <1e-99
## tpsF: 25                      101.215        3.95265  25.6067       <1e-99
## tpsF: 26                      106.09         3.95265  26.8401       <1e-99
## tpsF: 27                      106.702        3.95265  26.9951       <1e-99
## tpsF: 28                      111.04         3.95265  28.0924       <1e-99
## tpsF: 29                      113.41         3.95265  28.6922       <1e-99
## tpsF: 30                      115.902        3.95265  29.3225       <1e-99
## tpsF: 31                      119.852        3.95265  30.3219       <1e-99
## tpsF: 32                      119.79         3.95265  30.3061       <1e-99
## tpsF: 33                      122.74         3.95265  31.0525       <1e-99
## tpsF: 34                      123.452        3.95265  31.2327       <1e-99
## tpsF: 35                      128.94         3.95265  32.621        <1e-99
## tpsF: 36                      132.848        3.95265  33.6098       <1e-99
## tpsF: 37                      135.467        3.95265  34.2723       <1e-99
## tpsF: 38                      132.627        3.95265  33.5539       <1e-99
## tpsF: 39                      133.84         3.95265  33.8607       <1e-99
## tpsF: 40                      136.052        3.95265  34.4205       <1e-99
## tpsF: 41                      136.127        3.95265  34.4394       <1e-99
## tpsF: 42                      137.502        3.95265  34.7873       <1e-99
## tpsF: 43                      142.596        3.95265  36.076        <1e-99
## tpsF: 44                      144.065        3.95265  36.4476       <1e-99
## tpsF: 45                      145.927        3.95265  36.9188       <1e-99
## tpsF: 46                      144.69         3.95265  36.6057       <1e-99
## tpsF: 47                      148.252        3.95265  37.507        <1e-99
## tpsF: 48                      145.502        3.95265  36.8112       <1e-99
## tpsF: 49                      149.752        3.95265  37.8865       <1e-99
## tpsF: 50                      146.965        3.95265  37.1813       <1e-99
## tpsF: 51                      149.465        3.95265  37.8137       <1e-99
## tpsF: 52                      153.065        3.95265  38.7245       <1e-99
## tpsF: 53                      152.265        3.95265  38.5221       <1e-99
## tpsF: 54                      158.065        3.95265  39.9895       <1e-99
## tpsF: 55                      156.627        3.95265  39.6258       <1e-99
## tpsF: 56                      160.227        3.95265  40.5366       <1e-99
## tpsF: 57                      164.575        3.95265  41.6366       <1e-99
## tpsF: 58                      164.144        3.95265  41.5275       <1e-99
## tpsF: 59                      164.19         3.95265  41.5391       <1e-99
## tpsF: 60                      162.327        3.95265  41.0679       <1e-99
## tpsF: 61                      163.327        3.95265  41.3209       <1e-99
## tpsF: 62                      162.865        3.95265  41.2039       <1e-99
## tpsF: 63                      166.877        3.95265  42.219        <1e-99
## tpsF: 64                      167.098        3.95265  42.2749       <1e-99
## tpsF: 65                      168.135        3.95265  42.5374       <1e-99
## tpsF: 66                      170.027        3.95265  43.0159       <1e-99
## tpsF: 67                      168.077        3.95265  42.5226       <1e-99
## tpsF: 68                      171.327        3.95265  43.3448       <1e-99
## tpsF: 69                      169.04         3.95265  42.7661       <1e-99
## tpsF: 70                      172.94         3.95265  43.7528       <1e-99
## tpsF: 71                      171.64         3.95265  43.4239       <1e-99
## tpsF: 72                      172.765        3.95265  43.7085       <1e-99
## tpsF: 73                      175.14         3.95265  44.3094       <1e-99
## tpsF: 74                      174.827        3.95265  44.2303       <1e-99
## tpsF: 75                      180.365        3.95265  45.6313       <1e-99
## tpsF: 76                      177.127        3.95265  44.8122       <1e-99
## tpsF: 77                      179.715        3.95265  45.4668       <1e-99
## gest: NP                        0.479167     7.84964   0.0610431    0.9513
## ps: C2                          1.825        7.84964   0.232495     0.8162
## ps: NP                          2.78006      8.12515   0.342155     0.7322
## tpsF: 2 & gest: NP             -0.0875       5.58989  -0.0156532    0.9875
## tpsF: 3 & gest: NP             -0.104167     5.58989  -0.0186348    0.9851
## tpsF: 4 & gest: NP             -0.404167     5.58989  -0.0723031    0.9424
## tpsF: 5 & gest: NP             -0.279167     5.58989  -0.0499413    0.9602
## tpsF: 6 & gest: NP             -0.729167     5.58989  -0.130444     0.8962
## tpsF: 7 & gest: NP             -0.479167     5.58989  -0.0857202    0.9317
## tpsF: 8 & gest: NP             -0.1          5.58989  -0.0178894    0.9857
## tpsF: 9 & gest: NP              0.61375      5.58989   0.109796     0.9126
## tpsF: 10 & gest: NP             0.770833     5.58989   0.137898     0.8903
## tpsF: 11 & gest: NP            -0.435417     5.58989  -0.0778935    0.9379
## tpsF: 12 & gest: NP             1.94167      5.58989   0.347353     0.7283
## tpsF: 13 & gest: NP             0.795833     5.58989   0.14237      0.8868
## tpsF: 14 & gest: NP             1.30833      5.58989   0.234053     0.8149
## tpsF: 15 & gest: NP             0.563542     5.58989   0.100814     0.9197
## tpsF: 16 & gest: NP             0.155833     5.58989   0.0278777    0.9778
## tpsF: 17 & gest: NP             0.745833     5.58989   0.133425     0.8939
## tpsF: 18 & gest: NP            -0.0291667    5.58989  -0.00521775   0.9958
## tpsF: 19 & gest: NP             0.258333     5.58989   0.0462143    0.9631
## tpsF: 20 & gest: NP             0.445833     5.58989   0.079757     0.9364
## tpsF: 21 & gest: NP             0.808333     5.58989   0.144606     0.8850
## tpsF: 22 & gest: NP             1.16667      5.58989   0.20871      0.8347
## tpsF: 23 & gest: NP             1.50917      5.58989   0.269981     0.7872
## tpsF: 24 & gest: NP             0.595833     5.58989   0.106591     0.9151
## tpsF: 25 & gest: NP             0.0583333    5.58989   0.0104355    0.9917
## tpsF: 26 & gest: NP             1.65833      5.58989   0.296666     0.7667
## tpsF: 27 & gest: NP             0.970833     5.58989   0.173676     0.8621
## tpsF: 28 & gest: NP             1.50833      5.58989   0.269832     0.7873
## tpsF: 29 & gest: NP             2.0625       5.58989   0.368969     0.7122
## tpsF: 30 & gest: NP             2.83625      5.58989   0.507389     0.6119
## tpsF: 31 & gest: NP             1.45833      5.58989   0.260887     0.7942
## tpsF: 32 & gest: NP             1.65833      5.58989   0.296666     0.7667
## tpsF: 33 & gest: NP             2.98333      5.58989   0.533701     0.5935
## tpsF: 34 & gest: NP             1.65833      5.58989   0.296666     0.7667
## tpsF: 35 & gest: NP             0.958333     5.58989   0.17144      0.8639
## tpsF: 36 & gest: NP            -0.772917     5.58989  -0.13827      0.8900
## tpsF: 37 & gest: NP            -0.916667     5.58989  -0.163986     0.8697
## tpsF: 38 & gest: NP             4.78333      5.58989   0.855711     0.3922
## tpsF: 39 & gest: NP             2.42083      5.58989   0.433073     0.6650
## tpsF: 40 & gest: NP             3.45833      5.58989   0.618676     0.5361
## tpsF: 41 & gest: NP             4.65833      5.58989   0.833349     0.4046
## tpsF: 42 & gest: NP             6.23333      5.58989   1.11511      0.2648
## tpsF: 43 & gest: NP             1.92708      5.58989   0.344744     0.7303
## tpsF: 44 & gest: NP             2.79583      5.58989   0.500158     0.6170
## tpsF: 45 & gest: NP             3.53333      5.58989   0.632093     0.5273
## tpsF: 46 & gest: NP             5.63333      5.58989   1.00777      0.3136
## tpsF: 47 & gest: NP             5.99583      5.58989   1.07262      0.2834
## tpsF: 48 & gest: NP             6.09583      5.58989   1.09051      0.2755
## tpsF: 49 & gest: NP             4.25833      5.58989   0.761791     0.4462
## tpsF: 50 & gest: NP            10.9542       5.58989   1.95964      0.0500
## tpsF: 51 & gest: NP             9.25417      5.58989   1.65552      0.0978
## tpsF: 52 & gest: NP             7.48333      5.58989   1.33873      0.1807
## tpsF: 53 & gest: NP             5.99583      5.58989   1.07262      0.2834
## tpsF: 54 & gest: NP             4.08333      5.58989   0.730485     0.4651
## tpsF: 55 & gest: NP             3.69583      5.58989   0.661163     0.5085
## tpsF: 56 & gest: NP             6.25833      5.58989   1.11958      0.2629
## tpsF: 57 & gest: NP             0.233333     5.58989   0.041742     0.9667
## tpsF: 58 & gest: NP             2.32083      5.58989   0.415184     0.6780
## tpsF: 59 & gest: NP             4.92083      5.58989   0.880309     0.3787
## tpsF: 60 & gest: NP             5.47083      5.58989   0.978701     0.3277
## tpsF: 61 & gest: NP             9.09583      5.58989   1.62719      0.1037
## tpsF: 62 & gest: NP             6.63333      5.58989   1.18667      0.2354
## tpsF: 63 & gest: NP             6.37083      5.58989   1.13971      0.2544
## tpsF: 64 & gest: NP             6.75         5.58989   1.20754      0.2272
## tpsF: 65 & gest: NP             6.5125       5.58989   1.16505      0.2440
## tpsF: 66 & gest: NP             6.53333      5.58989   1.16878      0.2425
## tpsF: 67 & gest: NP             7.17083      5.58989   1.28282      0.1996
## tpsF: 68 & gest: NP             7.97083      5.58989   1.42594      0.1539
## tpsF: 69 & gest: NP            10.9583       5.58989   1.96038      0.0500
## tpsF: 70 & gest: NP            10.5833       5.58989   1.8933       0.0583
## tpsF: 71 & gest: NP            11.4754       5.58989   2.05289      0.0401
## tpsF: 72 & gest: NP            11.4342       5.58989   2.04551      0.0408
## tpsF: 73 & gest: NP            11.0583       5.58989   1.97827      0.0479
## tpsF: 74 & gest: NP             9.97083      5.58989   1.78372      0.0745
## tpsF: 75 & gest: NP             7.74583      5.58989   1.38568      0.1658
## tpsF: 76 & gest: NP             9.32083      5.58989   1.66744      0.0954
## tpsF: 77 & gest: NP             7.52083      5.58989   1.34543      0.1785
## tpsF: 2 & ps: C2               -0.2          5.58989  -0.0357788    0.9715
## tpsF: 3 & ps: C2               -0.5875       5.58989  -0.1051       0.9163
## tpsF: 4 & ps: C2               -0.225        5.58989  -0.0402512    0.9679
## tpsF: 5 & ps: C2               -0.9875       5.58989  -0.176658     0.8598
## tpsF: 6 & ps: C2               -1.05         5.58989  -0.187839     0.8510
## tpsF: 7 & ps: C2               -0.15         5.58989  -0.0268341    0.9786
## tpsF: 8 & ps: C2                0.675        5.58989   0.120754     0.9039
## tpsF: 9 & ps: C2                1.89         5.58989   0.33811      0.7353
## tpsF: 10 & ps: C2               2.8          5.58989   0.500904     0.6164
## tpsF: 11 & ps: C2              -0.285417     5.58989  -0.0510594    0.9593
## tpsF: 12 & ps: C2               0.429167     5.58989   0.0767754    0.9388
## tpsF: 13 & ps: C2               0.575        5.58989   0.102864     0.9181
## tpsF: 14 & ps: C2               2.7875       5.58989   0.498668     0.6180
## tpsF: 15 & ps: C2               4.01354      5.58989   0.717999     0.4728
## tpsF: 16 & ps: C2               5.4275       5.58989   0.970948     0.3316
## tpsF: 17 & ps: C2               2.5625       5.58989   0.458416     0.6467
## tpsF: 18 & ps: C2               2.0625       5.58989   0.368969     0.7122
## tpsF: 19 & ps: C2               3.0          5.58989   0.536683     0.5915
## tpsF: 20 & ps: C2               1.95         5.58989   0.348844     0.7272
## tpsF: 21 & ps: C2               4.5125       5.58989   0.80726      0.4195
## tpsF: 22 & ps: C2               4.48542      5.58989   0.802415     0.4223
## tpsF: 23 & ps: C2               5.13292      5.58989   0.918249     0.3585
## tpsF: 24 & ps: C2               4.7375       5.58989   0.847511     0.3967
## tpsF: 25 & ps: C2               5.0125       5.58989   0.896707     0.3699
## tpsF: 26 & ps: C2               5.875        5.58989   1.051        0.2933
## tpsF: 27 & ps: C2               6.05         5.58989   1.08231      0.2791
## tpsF: 28 & ps: C2               9.3625       5.58989   1.6749       0.0940
## tpsF: 29 & ps: C2               9.98542      5.58989   1.78633      0.0740
## tpsF: 30 & ps: C2              11.2629       5.58989   2.01487      0.0439
## tpsF: 31 & ps: C2              10.575        5.58989   1.89181      0.0585
## tpsF: 32 & ps: C2              10.3875       5.58989   1.85826      0.0631
## tpsF: 33 & ps: C2              12.7375       5.58989   2.27867      0.0227
## tpsF: 34 & ps: C2              11.6625       5.58989   2.08635      0.0369
## tpsF: 35 & ps: C2              12.5625       5.58989   2.24736      0.0246
## tpsF: 36 & ps: C2              10.9083       5.58989   1.95144      0.0510
## tpsF: 37 & ps: C2              10.0083       5.58989   1.79043      0.0734
## tpsF: 38 & ps: C2              15.9625       5.58989   2.8556       0.0043
## tpsF: 39 & ps: C2              13.6875       5.58989   2.44861      0.0143
## tpsF: 40 & ps: C2              15.2875       5.58989   2.73485      0.0062
## tpsF: 41 & ps: C2              15.075        5.58989   2.69683      0.0070
## tpsF: 42 & ps: C2              16.65         5.58989   2.97859      0.0029
## tpsF: 43 & ps: C2              14.1604       5.58989   2.53322      0.0113
## tpsF: 44 & ps: C2              15.1042       5.58989   2.70205      0.0069
## tpsF: 45 & ps: C2              15.9875       5.58989   2.86007      0.0042
## tpsF: 46 & ps: C2              16.55         5.58989   2.9607       0.0031
## tpsF: 47 & ps: C2              17.425        5.58989   3.11723      0.0018
## tpsF: 48 & ps: C2              15.725        5.58989   2.81311      0.0049
## tpsF: 49 & ps: C2              16.0375       5.58989   2.86902      0.0041
## tpsF: 50 & ps: C2              20.8458       5.58989   3.7292       0.0002
## tpsF: 51 & ps: C2              19.8833       5.58989   3.55701      0.0004
## tpsF: 52 & ps: C2              18.5          5.58989   3.30954      0.0009
## tpsF: 53 & ps: C2              18.8          5.58989   3.36321      0.0008
## tpsF: 54 & ps: C2              16.575        5.58989   2.96517      0.0030
## tpsF: 55 & ps: C2              16.8375       5.58989   3.01213      0.0026
## tpsF: 56 & ps: C2              18.8125       5.58989   3.36545      0.0008
## tpsF: 57 & ps: C2              10.6208       5.58989   1.90001      0.0574
## tpsF: 58 & ps: C2              13.8958       5.58989   2.48588      0.0129
## tpsF: 59 & ps: C2              17.4625       5.58989   3.12394      0.0018
## tpsF: 60 & ps: C2              19.8625       5.58989   3.55329      0.0004
## tpsF: 61 & ps: C2              24.0125       5.58989   4.2957       <1e-4 
## tpsF: 62 & ps: C2              21.9          5.58989   3.91778      <1e-4 
## tpsF: 63 & ps: C2              22.6875       5.58989   4.05866      <1e-4 
## tpsF: 64 & ps: C2              21.6125       5.58989   3.86635      0.0001
## tpsF: 65 & ps: C2              22.425        5.58989   4.0117       <1e-4 
## tpsF: 66 & ps: C2              23.7875       5.58989   4.25545      <1e-4 
## tpsF: 67 & ps: C2              23.3375       5.58989   4.17494      <1e-4 
## tpsF: 68 & ps: C2              26.0875       5.58989   4.6669       <1e-5 
## tpsF: 69 & ps: C2              27.2125       5.58989   4.86816      <1e-5 
## tpsF: 70 & ps: C2              25.7875       5.58989   4.61324      <1e-5 
## tpsF: 71 & ps: C2              29.8713       5.58989   5.34379      <1e-7 
## tpsF: 72 & ps: C2              29.285        5.58989   5.23892      <1e-6 
## tpsF: 73 & ps: C2              28.0125       5.58989   5.01127      <1e-6 
## tpsF: 74 & ps: C2              26.4375       5.58989   4.72952      <1e-5 
## tpsF: 75 & ps: C2              23.825        5.58989   4.26216      <1e-4 
## tpsF: 76 & ps: C2              25.2375       5.58989   4.51484      <1e-5 
## tpsF: 77 & ps: C2              22.6875       5.58989   4.05866      <1e-4 
## tpsF: 2 & ps: NP               -0.299107     5.78609  -0.0516942    0.9588
## tpsF: 3 & ps: NP               -0.540774     5.78609  -0.093461     0.9255
## tpsF: 4 & ps: NP               -1.0122       5.78609  -0.174937     0.8611
## tpsF: 5 & ps: NP               -1.13899      5.78609  -0.196849     0.8439
## tpsF: 6 & ps: NP               -1.66577      5.78609  -0.287893     0.7734
## tpsF: 7 & ps: NP               -1.34435      5.78609  -0.232341     0.8163
## tpsF: 8 & ps: NP               -0.995238     5.78609  -0.172005     0.8634
## tpsF: 9 & ps: NP               -0.285774     5.78609  -0.0493898    0.9606
## tpsF: 10 & ps: NP              -0.385417     5.78609  -0.0666109    0.9469
## tpsF: 11 & ps: NP              -0.454167     5.78609  -0.0784928    0.9374
## tpsF: 12 & ps: NP              -0.980655     5.78609  -0.169485     0.8654
## tpsF: 13 & ps: NP              -0.678274     5.78609  -0.117225     0.9067
## tpsF: 14 & ps: NP              -0.258631     5.78609  -0.0446987    0.9643
## tpsF: 15 & ps: NP              -0.0876488    5.78609  -0.0151482    0.9879
## tpsF: 16 & ps: NP               0.171071     5.78609   0.029566     0.9764
## tpsF: 17 & ps: NP              -0.046131     5.78609  -0.00797273   0.9936
## tpsF: 18 & ps: NP               0.30744      5.78609   0.0531344    0.9576
## tpsF: 19 & ps: NP              -0.0568452    5.78609  -0.00982446   0.9922
## tpsF: 20 & ps: NP               1.03065      5.78609   0.178126     0.8586
## tpsF: 21 & ps: NP               0.843155     5.78609   0.145721     0.8841
## tpsF: 22 & ps: NP               1.66815      5.78609   0.288304     0.7731
## tpsF: 23 & ps: NP               2.21565      5.78609   0.382928     0.7018
## tpsF: 24 & ps: NP               2.36637      5.78609   0.408975     0.6826
## tpsF: 25 & ps: NP               3.87351      5.78609   0.669452     0.5032
## tpsF: 26 & ps: NP               2.88423      5.78609   0.498476     0.6181
## tpsF: 27 & ps: NP               2.22887      5.78609   0.385212     0.7001
## tpsF: 28 & ps: NP               4.54851      5.78609   0.786111     0.4318
## tpsF: 29 & ps: NP               4.69435      5.78609   0.811315     0.4172
## tpsF: 30 & ps: NP               5.62881      5.78609   0.972817     0.3306
## tpsF: 31 & ps: NP               2.96458      5.78609   0.512364     0.6084
## tpsF: 32 & ps: NP               4.3128       5.78609   0.745373     0.4560
## tpsF: 33 & ps: NP               7.14851      5.78609   1.23546      0.2167
## tpsF: 34 & ps: NP               6.82173      5.78609   1.17899      0.2384
## tpsF: 35 & ps: NP               4.44851      5.78609   0.768829     0.4420
## tpsF: 36 & ps: NP               3.11875      5.78609   0.539008     0.5899
## tpsF: 37 & ps: NP               2.63571      5.78609   0.455526     0.6487
## tpsF: 38 & ps: NP               8.60387      5.78609   1.48699      0.1370
## tpsF: 39 & ps: NP               6.97708      5.78609   1.20584      0.2279
## tpsF: 40 & ps: NP               8.63601      5.78609   1.49255      0.1356
## tpsF: 41 & ps: NP               8.6753       5.78609   1.49934      0.1338
## tpsF: 42 & ps: NP               9.7003       5.78609   1.67649      0.0936
## tpsF: 43 & ps: NP               7.3994       5.78609   1.27883      0.2010
## tpsF: 44 & ps: NP               8.45208      5.78609   1.46076      0.1441
## tpsF: 45 & ps: NP               9.4753       5.78609   1.6376       0.1015
## tpsF: 46 & ps: NP               9.94137      5.78609   1.71815      0.0858
## tpsF: 47 & ps: NP               8.00744      5.78609   1.38391      0.1664
## tpsF: 48 & ps: NP               9.02887      5.78609   1.56044      0.1187
## tpsF: 49 & ps: NP               9.16458      5.78609   1.5839       0.1132
## tpsF: 50 & ps: NP              11.7687       5.78609   2.03397      0.0420
## tpsF: 51 & ps: NP              11.0902       5.78609   1.9167       0.0553
## tpsF: 52 & ps: NP              10.0949       5.78609   1.74469      0.0810
## tpsF: 53 & ps: NP              10.3664       5.78609   1.7916       0.0732
## tpsF: 54 & ps: NP               8.7378       5.78609   1.51014      0.1310
## tpsF: 55 & ps: NP               9.41815      5.78609   1.62772      0.1036
## tpsF: 56 & ps: NP               8.21815      5.78609   1.42033      0.1555
## tpsF: 57 & ps: NP               5.84643      5.78609   1.01043      0.3123
## tpsF: 58 & ps: NP               7.90625      5.78609   1.36642      0.1718
## tpsF: 59 & ps: NP              10.2414       5.78609   1.77         0.0767
## tpsF: 60 & ps: NP              11.4753       5.78609   1.98326      0.0473
## tpsF: 61 & ps: NP              14.361        5.78609   2.48199      0.0131
## tpsF: 62 & ps: NP              13.2664       5.78609   2.2928       0.0219
## tpsF: 63 & ps: NP              13.6396       5.78609   2.35731      0.0184
## tpsF: 64 & ps: NP              14.7735       5.78609   2.55328      0.0107
## tpsF: 65 & ps: NP              14.9289       5.78609   2.58013      0.0099
## tpsF: 66 & ps: NP              14.9896       5.78609   2.59062      0.0096
## tpsF: 67 & ps: NP              13.8824       5.78609   2.39928      0.0164
## tpsF: 68 & ps: NP              15.0039       5.78609   2.59309      0.0095
## tpsF: 69 & ps: NP              16.2057       5.78609   2.8008       0.0051
## tpsF: 70 & ps: NP              17.7057       5.78609   3.06004      0.0022
## tpsF: 71 & ps: NP              21.8961       5.78609   3.78427      0.0002
## tpsF: 72 & ps: NP              21.514        5.78609   3.71822      0.0002
## tpsF: 73 & ps: NP              20.3199       5.78609   3.51186      0.0004
## tpsF: 74 & ps: NP              20.061        5.78609   3.46711      0.0005
## tpsF: 75 & ps: NP              16.3949       5.78609   2.83351      0.0046
## tpsF: 76 & ps: NP              18.3753       5.78609   3.17577      0.0015
## tpsF: 77 & ps: NP              13.8449       5.78609   2.3928       0.0167
## gest: NP & ps: C2              -1.75208     11.1011   -0.15783      0.8746
## gest: NP & ps: NP              -2.02173     11.2976   -0.178952     0.8580
## tpsF: 2 & gest: NP & ps: C2    -0.51875      7.90531  -0.0656205    0.9477
## tpsF: 3 & gest: NP & ps: C2    -0.947917     7.90531  -0.119909     0.9046
## tpsF: 4 & gest: NP & ps: C2    -1.73542      7.90531  -0.219526     0.8262
## tpsF: 5 & gest: NP & ps: C2    -1.98542      7.90531  -0.25115      0.8017
## tpsF: 6 & gest: NP & ps: C2    -1.69792      7.90531  -0.214782     0.8299
## tpsF: 7 & gest: NP & ps: C2    -2.17292      7.90531  -0.274868     0.7834
## tpsF: 8 & gest: NP & ps: C2    -2.63542      7.90531  -0.333373     0.7389
## tpsF: 9 & gest: NP & ps: C2    -3.48292      7.90531  -0.44058      0.6595
## tpsF: 10 & gest: NP & ps: C2   -5.68542      7.90531  -0.71919      0.4720
## tpsF: 11 & gest: NP & ps: C2   -2.03958      7.90531  -0.258002     0.7964
## tpsF: 12 & gest: NP & ps: C2   -4.31319      7.90531  -0.545608     0.5853
## tpsF: 13 & gest: NP & ps: C2   -4.91042      7.90531  -0.621155     0.5345
## tpsF: 14 & gest: NP & ps: C2   -6.09792      7.90531  -0.77137      0.4405
## tpsF: 15 & gest: NP & ps: C2   -7.29583      7.90531  -0.922903     0.3561
## tpsF: 16 & gest: NP & ps: C2   -8.43056      7.90531  -1.06644      0.2862
## tpsF: 17 & gest: NP & ps: C2   -6.17292      7.90531  -0.780857     0.4349
## tpsF: 18 & gest: NP & ps: C2   -6.46042      7.90531  -0.817225     0.4138
## tpsF: 19 & gest: NP & ps: C2   -6.67292      7.90531  -0.844106     0.3986
## tpsF: 20 & gest: NP & ps: C2   -7.18542      7.90531  -0.908936     0.3634
## tpsF: 21 & gest: NP & ps: C2   -8.77292      7.90531  -1.10975      0.2671
## tpsF: 22 & gest: NP & ps: C2   -9.41458      7.90531  -1.19092      0.2337
## tpsF: 23 & gest: NP & ps: C2  -10.3271       7.90531  -1.30635      0.1914
## tpsF: 24 & gest: NP & ps: C2   -7.02292      7.90531  -0.88838      0.3743
## tpsF: 25 & gest: NP & ps: C2  -10.5479       7.90531  -1.33428      0.1821
## tpsF: 26 & gest: NP & ps: C2  -12.6979       7.90531  -1.60625      0.1082
## tpsF: 27 & gest: NP & ps: C2  -14.5354       7.90531  -1.83869      0.0660
## tpsF: 28 & gest: NP & ps: C2  -16.2854       7.90531  -2.06006      0.0394
## tpsF: 29 & gest: NP & ps: C2  -18.6292       7.90531  -2.35654      0.0184
## tpsF: 30 & gest: NP & ps: C2  -21.3204       7.90531  -2.69698      0.0070
## tpsF: 31 & gest: NP & ps: C2  -16.4604       7.90531  -2.0822       0.0373
## tpsF: 32 & gest: NP & ps: C2  -17.0604       7.90531  -2.1581       0.0309
## tpsF: 33 & gest: NP & ps: C2  -17.8479       7.90531  -2.25771      0.0240
## tpsF: 34 & gest: NP & ps: C2  -16.9354       7.90531  -2.14228      0.0322
## tpsF: 35 & gest: NP & ps: C2  -16.6104       7.90531  -2.10117      0.0356
## tpsF: 36 & gest: NP & ps: C2  -12.6042       7.90531  -1.59439      0.1108
## tpsF: 37 & gest: NP & ps: C2  -11.5354       7.90531  -1.4592       0.1445
## tpsF: 38 & gest: NP & ps: C2  -19.1604       7.90531  -2.42374      0.0154
## tpsF: 39 & gest: NP & ps: C2  -17.9229       7.90531  -2.2672       0.0234
## tpsF: 40 & gest: NP & ps: C2  -18.8104       7.90531  -2.37947      0.0173
## tpsF: 41 & gest: NP & ps: C2  -20.7104       7.90531  -2.61981      0.0088
## tpsF: 42 & gest: NP & ps: C2  -21.5354       7.90531  -2.72417      0.0064
## tpsF: 43 & gest: NP & ps: C2  -19.4375       7.90531  -2.45879      0.0139
## tpsF: 44 & gest: NP & ps: C2  -19.9687       7.90531  -2.52599      0.0115
## tpsF: 45 & gest: NP & ps: C2  -20.3104       7.90531  -2.56921      0.0102
## tpsF: 46 & gest: NP & ps: C2  -21.4479       7.90531  -2.7131       0.0067
## tpsF: 47 & gest: NP & ps: C2  -22.4979       7.90531  -2.84593      0.0044
## tpsF: 48 & gest: NP & ps: C2  -22.1854       7.90531  -2.8064       0.0050
## tpsF: 49 & gest: NP & ps: C2  -20.6479       7.90531  -2.61191      0.0090
## tpsF: 50 & gest: NP & ps: C2  -25.8187       7.90531  -3.266        0.0011
## tpsF: 51 & gest: NP & ps: C2  -24.9687       7.90531  -3.15848      0.0016
## tpsF: 52 & gest: NP & ps: C2  -23.9854       7.90531  -3.03409      0.0024
## tpsF: 53 & gest: NP & ps: C2  -23.5354       7.90531  -2.97717      0.0029
## tpsF: 54 & gest: NP & ps: C2  -22.2854       7.90531  -2.81905      0.0048
## tpsF: 55 & gest: NP & ps: C2  -22.4729       7.90531  -2.84276      0.0045
## tpsF: 56 & gest: NP & ps: C2  -26.0729       7.90531  -3.29815      0.0010
## tpsF: 57 & gest: NP & ps: C2  -15.8271       7.90531  -2.00208      0.0453
## tpsF: 58 & gest: NP & ps: C2  -19.6146       7.90531  -2.48119      0.0131
## tpsF: 59 & gest: NP & ps: C2  -23.5229       7.90531  -2.97559      0.0029
## tpsF: 60 & gest: NP & ps: C2  -26.9479       7.90531  -3.40884      0.0007
## tpsF: 61 & gest: NP & ps: C2  -31.0979       7.90531  -3.9338       <1e-4 
## tpsF: 62 & gest: NP & ps: C2  -27.2604       7.90531  -3.44837      0.0006
## tpsF: 63 & gest: NP & ps: C2  -28.7229       7.90531  -3.63337      0.0003
## tpsF: 64 & gest: NP & ps: C2  -30.6021       7.90531  -3.87108      0.0001
## tpsF: 65 & gest: NP & ps: C2  -30.2271       7.90531  -3.82365      0.0001
## tpsF: 66 & gest: NP & ps: C2  -30.5479       7.90531  -3.86423      0.0001
## tpsF: 67 & gest: NP & ps: C2  -30.6229       7.90531  -3.87372      0.0001
## tpsF: 68 & gest: NP & ps: C2  -32.0979       7.90531  -4.0603       <1e-4 
## tpsF: 69 & gest: NP & ps: C2  -35.3479       7.90531  -4.47142      <1e-5 
## tpsF: 70 & gest: NP & ps: C2  -31.7729       7.90531  -4.01919      <1e-4 
## tpsF: 71 & gest: NP & ps: C2  -39.3533       7.90531  -4.97809      <1e-6 
## tpsF: 72 & gest: NP & ps: C2  -38.0096       7.90531  -4.80811      <1e-5 
## tpsF: 73 & gest: NP & ps: C2  -36.9104       7.90531  -4.66907      <1e-5 
## tpsF: 74 & gest: NP & ps: C2  -33.9229       7.90531  -4.29116      <1e-4 
## tpsF: 75 & gest: NP & ps: C2  -30.3854       7.90531  -3.84367      0.0001
## tpsF: 76 & gest: NP & ps: C2  -33.8479       7.90531  -4.28167      <1e-4 
## tpsF: 77 & gest: NP & ps: C2  -31.0104       7.90531  -3.92273      <1e-4 
## tpsF: 2 & gest: NP & ps: NP     0.0991071    8.04523   0.0123187    0.9902
## tpsF: 3 & gest: NP & ps: NP     0.20744      8.04523   0.0257843    0.9794
## tpsF: 4 & gest: NP & ps: NP     0.278869     8.04523   0.0346626    0.9723
## tpsF: 5 & gest: NP & ps: NP     0.405655     8.04523   0.0504218    0.9598
## tpsF: 6 & gest: NP & ps: NP     0.65744      8.04523   0.081718     0.9349
## tpsF: 7 & gest: NP & ps: NP     1.23601      8.04523   0.153633     0.8779
## tpsF: 8 & gest: NP & ps: NP     0.916071     8.04523   0.113865     0.9093
## tpsF: 9 & gest: NP & ps: NP     0.466607     8.04523   0.057998     0.9538
## tpsF: 10 & gest: NP & ps: NP   -0.997917     8.04523  -0.124038     0.9013
## tpsF: 11 & gest: NP & ps: NP    1.50208      8.04523   0.186705     0.8519
## tpsF: 12 & gest: NP & ps: NP    0.0612103    8.04523   0.00760827   0.9939
## tpsF: 13 & gest: NP & ps: NP    0.21994      8.04523   0.027338     0.9782
## tpsF: 14 & gest: NP & ps: NP   -0.712202     8.04523  -0.0885248    0.9295
## tpsF: 15 & gest: NP & ps: NP   -0.569643     8.04523  -0.070805     0.9436
## tpsF: 16 & gest: NP & ps: NP   -0.747877     8.04523  -0.092959     0.9259
## tpsF: 17 & gest: NP & ps: NP    0.100298     8.04523   0.0124667    0.9901
## tpsF: 18 & gest: NP & ps: NP   -0.0407738    8.04523  -0.00506807   0.9960
## tpsF: 19 & gest: NP & ps: NP   -0.563988     8.04523  -0.0701021    0.9441
## tpsF: 20 & gest: NP & ps: NP   -1.03899      8.04523  -0.129143     0.8972
## tpsF: 21 & gest: NP & ps: NP   -1.22649      8.04523  -0.152449     0.8788
## tpsF: 22 & gest: NP & ps: NP   -2.3494       8.04523  -0.292024     0.7703
## tpsF: 23 & gest: NP & ps: NP   -3.1269       8.04523  -0.388666     0.6975
## tpsF: 24 & gest: NP & ps: NP   -2.0372       8.04523  -0.253219     0.8001
## tpsF: 25 & gest: NP & ps: NP   -2.90685      8.04523  -0.361313     0.7179
## tpsF: 26 & gest: NP & ps: NP   -4.45506      8.04523  -0.553751     0.5797
## tpsF: 27 & gest: NP & ps: NP   -2.3622       8.04523  -0.293615     0.7691
## tpsF: 28 & gest: NP & ps: NP   -5.26935      8.04523  -0.654965     0.5125
## tpsF: 29 & gest: NP & ps: NP   -5.7131       8.04523  -0.710122     0.4776
## tpsF: 30 & gest: NP & ps: NP   -7.03631      8.04523  -0.874594     0.3818
## tpsF: 31 & gest: NP & ps: NP   -2.94792      8.04523  -0.366418     0.7141
## tpsF: 32 & gest: NP & ps: NP   -3.82113      8.04523  -0.474956     0.6348
## tpsF: 33 & gest: NP & ps: NP   -7.59435      8.04523  -0.943956     0.3452
## tpsF: 34 & gest: NP & ps: NP   -6.30506      8.04523  -0.783701     0.4332
## tpsF: 35 & gest: NP & ps: NP   -3.23185      8.04523  -0.401709     0.6879
## tpsF: 36 & gest: NP & ps: NP   -1.95833      8.04523  -0.243415     0.8077
## tpsF: 37 & gest: NP & ps: NP   -1.0378       8.04523  -0.128995     0.8974
## tpsF: 38 & gest: NP & ps: NP   -7.4247       8.04523  -0.92287      0.3561
## tpsF: 39 & gest: NP & ps: NP   -3.12292      8.04523  -0.38817      0.6979
## tpsF: 40 & gest: NP & ps: NP   -3.98185      8.04523  -0.494932     0.6206
## tpsF: 41 & gest: NP & ps: NP   -5.08363      8.04523  -0.631881     0.5275
## tpsF: 42 & gest: NP & ps: NP   -6.93363      8.04523  -0.861831     0.3888
## tpsF: 43 & gest: NP & ps: NP   -0.0931548    8.04523  -0.0115789    0.9908
## tpsF: 44 & gest: NP & ps: NP   -1.80208      8.04523  -0.223994     0.8228
## tpsF: 45 & gest: NP & ps: NP   -3.68363      8.04523  -0.457865     0.6470
## tpsF: 46 & gest: NP & ps: NP   -6.3997       8.04523  -0.795465     0.4263
## tpsF: 47 & gest: NP & ps: NP   -5.95327      8.04523  -0.739975     0.4593
## tpsF: 48 & gest: NP & ps: NP   -4.4997       8.04523  -0.5593       0.5760
## tpsF: 49 & gest: NP & ps: NP   -3.23542      8.04523  -0.402153     0.6876
## tpsF: 50 & gest: NP & ps: NP   -9.45208      8.04523  -1.17487      0.2400
## tpsF: 51 & gest: NP & ps: NP   -7.73601      8.04523  -0.961565     0.3363
## tpsF: 52 & gest: NP & ps: NP   -6.19077      8.04523  -0.769496     0.4416
## tpsF: 53 & gest: NP & ps: NP   -3.9622       8.04523  -0.492491     0.6224
## tpsF: 54 & gest: NP & ps: NP   -2.75863      8.04523  -0.34289      0.7317
## tpsF: 55 & gest: NP & ps: NP   -2.46399      8.04523  -0.306267     0.7594
## tpsF: 56 & gest: NP & ps: NP   -2.55149      8.04523  -0.317143     0.7511
## tpsF: 57 & gest: NP & ps: NP    2.90565      8.04523   0.361165     0.7180
## tpsF: 58 & gest: NP & ps: NP    0.764583     8.04523   0.0950356    0.9243
## tpsF: 59 & gest: NP & ps: NP   -1.6622       8.04523  -0.206607     0.8363
## tpsF: 60 & gest: NP & ps: NP   -2.94613      8.04523  -0.366196     0.7142
## tpsF: 61 & gest: NP & ps: NP   -5.94435      8.04523  -0.738866     0.4600
## tpsF: 62 & gest: NP & ps: NP   -3.1622       8.04523  -0.393053     0.6943
## tpsF: 63 & gest: NP & ps: NP   -4.39792      8.04523  -0.546649     0.5846
## tpsF: 64 & gest: NP & ps: NP   -6.97351      8.04523  -0.866788     0.3861
## tpsF: 65 & gest: NP & ps: NP   -6.30387      8.04523  -0.783553     0.4333
## tpsF: 66 & gest: NP & ps: NP   -5.38542      8.04523  -0.669392     0.5032
## tpsF: 67 & gest: NP & ps: NP   -4.26577      8.04523  -0.530224     0.5960
## tpsF: 68 & gest: NP & ps: NP   -5.3372       8.04523  -0.663399     0.5071
## tpsF: 69 & gest: NP & ps: NP   -7.71399      8.04523  -0.958827     0.3376
## tpsF: 70 & gest: NP & ps: NP   -8.53899      8.04523  -1.06137      0.2885
## tpsF: 71 & gest: NP & ps: NP   -7.64655      8.04523  -0.950445     0.3419
## tpsF: 72 & gest: NP & ps: NP   -8.57565      8.04523  -1.06593      0.2865
## tpsF: 73 & gest: NP & ps: NP   -8.56577      8.04523  -1.0647       0.2870
## tpsF: 74 & gest: NP & ps: NP  -10.2818       8.04523  -1.278        0.2012
## tpsF: 75 & gest: NP & ps: NP   -8.92827      8.04523  -1.10976      0.2671
## tpsF: 76 & gest: NP & ps: NP   -8.03363      8.04523  -0.998558     0.3180
## tpsF: 77 & gest: NP & ps: NP   -4.41577      8.04523  -0.548868     0.5831
## ──────────────────────────────────────────────────────────────────────────