-
Notifications
You must be signed in to change notification settings - Fork 0
/
Graphs.R
4447 lines (4308 loc) · 605 KB
/
Graphs.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
##############################################################/############
###### The implementation of irrigation leads to declines in farmland birds ######
############################### Graphics ################################
###########################################################################
# Load the reshape library
library(reshape)
# Load the ggplot2 library
library(ggplot2)
# Load the cowplot library
library(cowplot)
# Load the ggthemes library
library(ggthemes)
Output.model<- output.Model
################################
########### Figure 2 ###########
###### Effect of irrigation ####
alpha2.val <- cbind(c(Output.model$q2.5$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]),
c(Output.model$q25$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]),
c(Output.model$mean$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]),
c(Output.model$q75$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]),
c(Output.model$q97.5$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]),
c(Output.model$overlap0$effect.a0.sp[c(45,1,15,7,22,29,18,33,12,9,47,30,40,5,20,23,3,27,26,44,42,43,38,24,37,13,21,34,2,32,17,4,8,39,10,19,35,25,28,46,41,16,36,14,11,31,6)]))
sppnames2 <- c("56", "55", "54", "53", "52",
"51", "50", "49", "48", "47",
"46", "45", "44", "43", "42",
"41", "40", "39", "38", "37",
"36", "35", "34", "33", "32",
"31", "30", "29", "28", "27",
"26", "25", "24", "23", "22",
"21", "20", "19", "18", "17",
"16", "15", "14", "13", "12",
"11", "10")
values2 <- data.frame(sppnames2, alpha2.val)
colnames(values2) <- c("species", "lower.alpha", "l25.alpha", "mean.alpha", "u75.alpha", "upper.alpha", "alpha.sig")
values2$species <- factor(values2$species, levels = values2$species)
values2$nodiff <- values2$l25.alpha < 0 & values2$u75.alpha > 0
values2$H95 <- values2$lower.alpha > 0
values2$H50 <- values2$lower.alpha < 0 & values2$l25.alpha > 0
values2$L95 <- values2$upper.alpha < 0
values2$L50 <- values2$upper.alpha > 0 & values2$u75.alpha < 0
# Community mean
alpha2.1.val <- cbind(c(mean(output.Model$q2.5$effect.a0.sp)),
c(mean(output.Model$q25$effect.a0.sp)),
c(mean(output.Model$mean$effect.a0.sp)),
c(mean(output.Model$q75$effect.a0.sp)),
c(mean(output.Model$q97.5$effect.a0.sp)),
c(0))
sppnames2.1 <- c("57")
values2.1 <- data.frame(sppnames2.1, alpha2.1.val)
colnames(values2.1) <- c("species", "lower.alpha", "l25.alpha", "mean.alpha", "u75.alpha", "upper.alpha", "alpha.sig")
values2.1$species <- factor(values2.1$species, levels = values2.1$species)
values2.1$nodiff <- values2.1$l25.alpha < 0 & values2.1$u75.alpha > 0
values2.1$H95 <- values2.1$lower.alpha > 0
values2.1$H50 <- values2.1$lower.alpha < 0 & values2.1$l25.alpha > 0
values2.1$L95 <- values2.1$upper.alpha < 0
values2.1$L50 <- values2.1$upper.alpha > 0 & values2.1$u75.alpha < 0
#Figure
Figure2_C <- ggplot() +
geom_vline(xintercept = 0, alpha = 0.75, linetype = "dotdash") +
geom_hline(yintercept = 47.5, size = 0.35) +
geom_hline(yintercept = 32.5, linetype = "dotted", size = 0.6) +
geom_hline(yintercept = 23.5, linetype = "dotted", size = 0.6) +
geom_hline(yintercept = 22.5, linetype = "dotted", size = 0.6) +
geom_hline(yintercept = 19.5, linetype = "dotted", size = 0.6) +
geom_hline(yintercept = 17.5, linetype = "dotted", size = 0.6) +
geom_hline(yintercept = 15.5, linetype = "dotted", size = 0.6) +
geom_errorbarh(data = subset(values2.1, L50 == TRUE), aes(y = species, xmin = lower.alpha, xmax = upper.alpha, color = "red2"),
size = 0.75, height = 0) +
geom_errorbarh(data = subset(values2.1, L50 == TRUE), aes(y = species, xmin = l25.alpha, xmax = u75.alpha, color = "red2"),
size = 2, height = 0) +
geom_errorbarh(data = subset(values2.1, L50 == TRUE), aes(y = species, xmin = mean.alpha, xmax = mean.alpha, color = "black"),
height = 0.5) +
geom_errorbarh(data = subset(values2, nodiff == TRUE), aes(y = species, xmin = lower.alpha, xmax = upper.alpha, color = "gray80"),
height = 0, size = 0.75) +
geom_errorbarh(data = subset(values2, nodiff == TRUE), aes(y = species, xmin = l25.alpha, xmax = u75.alpha, color = "gray80"),
height = 0, size = 2) +
geom_errorbarh(data = subset(values2, nodiff == TRUE), aes(y = species, xmin = mean.alpha, xmax = mean.alpha, color = "black"),
height = 0.5) +
geom_errorbarh(data = subset(values2, H50 == TRUE), aes(y = species, xmin = lower.alpha, xmax = upper.alpha, color = "blue2"),
size = 0.75, height = 0) +
geom_errorbarh(data = subset(values2, H50 == TRUE), aes(y = species, xmin = l25.alpha, xmax = u75.alpha, color = "blue2"),
size = 2, height = 0) +
geom_errorbarh(data = subset(values2, H50 == TRUE), aes(y = species, xmin = mean.alpha, xmax = mean.alpha, color = "black"),
height = 0.5) +
geom_errorbarh(data = subset(values2, L50 == TRUE), aes(y = species, xmin = lower.alpha, xmax = upper.alpha, color = "red2"),
size = 0.75, height = 0) +
geom_errorbarh(data = subset(values2, L50 == TRUE), aes(y = species, xmin = l25.alpha, xmax = u75.alpha, color = "red2"),
size = 2, height = 0) +
geom_errorbarh(data = subset(values2, L50 == TRUE), aes(y = species, xmin = mean.alpha, xmax = mean.alpha, color = "black"),
height = 0.5) +
geom_errorbarh(data = subset(values2, L95 == TRUE), aes(y = species, xmin = lower.alpha, xmax = upper.alpha, color = "red"),
size = 0.75, height = 0) +
geom_errorbarh(data = subset(values2, L95 == TRUE), aes(y = species, xmin = l25.alpha, xmax = u75.alpha, color = "red"),
size = 2, height = 0) +
geom_errorbarh(data = subset(values2, L95 == TRUE), aes(y = species, xmin = mean.alpha, xmax = mean.alpha, color = "black"),
height = 0.5) +
coord_cartesian(xlim = c(-7.3, 7.3)) +
annotate("text", y = 40, x = 5.65, label = "Farmland", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 28, x = 5.65, label = "Shrubland", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 23, x = 5.65, label = "Rocky", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 21, x = 5.65, label = "Forest", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 18.5, x = 5.65, label = "Wetland", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 16.5, x = 5.65, label = "Urban", size=3.4, color="black", fontface = "bold") +
annotate("text", y = 7.5, x = 5.4, label = "Non-specialist", size=3.4, color="black", fontface = "bold") +
scale_color_manual(name = "", values = c("red" = "red3", "red2" = "chocolate2", "blue2" = "steelblue3", "black" = "black", "gray80" = "gray80")) +
scale_y_discrete(labels = c("Common buzzard", "Black kite", "Western marsh harrier", "Carrion crow", "Magpie",
"Common cuckoo", "Spotless starling", "Common blackbird", "Common nightingale",
"Barn swallow", "House sparrow", "Cirl bunting", "European greenfinch", "European serin",
"Common linnet", "Common swift", "Common house martin", "White wagtail", "Mallard",
"Great tit", "Common chaffinch", "Common wood pigeon", "Red-billed chough", "Melodious warbler",
"European stonechat", "Sardinian warbler", "Subalpine warbler", "Dartford warbler",
"Iberian grey shrike", "Woodchat shrike", "Tawny pipit", "Thekla's lark", "Common kestrel",
"Little owl", "European turtle-dove", "European bee-eater", "Eurasian hoopoe",
"European goldfinch", "Zitting cisticola", "Black-eared wheatear", "Corn bunting",
"Calandra lark", "Crested lark", "Greater short-toed lark", "Common quail",
"Red-legged partridge", "Little bustard", "Community mean"))+
scale_x_continuous(position = 'top')+
theme_few()+
xlab("Change in occurrence") +
ylab("") +
theme(plot.margin = unit(c(0.25, 0.25, 0.15, -0.3), "cm"),
text = element_text(size = 11),
panel.background = element_rect(fill = "transparent", color = NA),
plot.background = element_rect(fill = "transparent", color = NA),
axis.text.y = element_text(size = 10, vjust = 0.4, color = "black"),
axis.title.x = element_text(size=11.5, face="bold"),
axis.ticks.y = element_line(size= 1.2, color = c("gray80", "chocolate2", "gray80", "red3", "red3", "gray80",
"chocolate2", "gray80", "red3", "steelblue3", "chocolate2",
"gray80", "gray80", "chocolate2", "gray80", "gray80", "steelblue3",
"steelblue3", "gray80", "chocolate2", "red3", "gray80", "gray80",
"chocolate2", "gray80", "steelblue3", "chocolate2", "gray80",
"chocolate2", "chocolate2", "chocolate2", "chocolate2",
"chocolate2", "red3", "gray80", "red3", "red3", "steelblue3",
"gray80", "red3", "chocolate2", "gray80", "chocolate2", "red3",
"chocolate2", "chocolate2", "red3", "chocolate2")),
legend.position = "none")
Figure2_C # size = 4.9 x 7.2
#####################################################
###################### Figure 3 ######################
## Species richness by different contest and years ##
# Control sampling locations
alpha3.val.control <- cbind(c(mean(Output.model$q2.5$Nsite[c(1,5),1]),
mean(Output.model$q2.5$Nsite[c(1,5),2]),
mean(Output.model$q2.5$Nsite[c(1,5),3]),
mean(Output.model$q2.5$Nsite[c(1,5),4]),
mean(Output.model$q2.5$Nsite[c(1,5),5]),
mean(Output.model$q2.5$Nsite[c(1,5),6]),
mean(Output.model$q2.5$Nsite[c(1,5),7]),
mean(Output.model$q2.5$Nsite[c(1,5),8]),
mean(Output.model$q2.5$Nsite[c(1,5),9]),
mean(Output.model$q2.5$Nsite[c(1,5),10]),
mean(Output.model$q2.5$Nsite[c(1,5),11]),
mean(Output.model$q2.5$Nsite[c(1,5),12])),
c(mean(Output.model$q25$Nsite[c(1,5),1]),
mean(Output.model$q25$Nsite[c(1,5),2]),
mean(Output.model$q25$Nsite[c(1,5),3]),
mean(Output.model$q25$Nsite[c(1,5),4]),
mean(Output.model$q25$Nsite[c(1,5),5]),
mean(Output.model$q25$Nsite[c(1,5),6]),
mean(Output.model$q25$Nsite[c(1,5),7]),
mean(Output.model$q25$Nsite[c(1,5),8]),
mean(Output.model$q25$Nsite[c(1,5),9]),
mean(Output.model$q25$Nsite[c(1,5),10]),
mean(Output.model$q25$Nsite[c(1,5),11]),
mean(Output.model$q25$Nsite[c(1,5),12])),
c(mean(Output.model$mean$Nsite[c(1,5),1]),
mean(Output.model$mean$Nsite[c(1,5),2]),
mean(Output.model$mean$Nsite[c(1,5),3]),
mean(Output.model$mean$Nsite[c(1,5),4]),
mean(Output.model$mean$Nsite[c(1,5),5]),
mean(Output.model$mean$Nsite[c(1,5),6]),
mean(Output.model$mean$Nsite[c(1,5),7]),
mean(Output.model$mean$Nsite[c(1,5),8]),
mean(Output.model$mean$Nsite[c(1,5),9]),
mean(Output.model$mean$Nsite[c(1,5),10]),
mean(Output.model$mean$Nsite[c(1,5),11]),
mean(Output.model$mean$Nsite[c(1,5),12])),
c(mean(Output.model$q75$Nsite[c(1,5),1]),
mean(Output.model$q75$Nsite[c(1,5),2]),
mean(Output.model$q75$Nsite[c(1,5),3]),
mean(Output.model$q75$Nsite[c(1,5),4]),
mean(Output.model$q75$Nsite[c(1,5),5]),
mean(Output.model$q75$Nsite[c(1,5),6]),
mean(Output.model$q75$Nsite[c(1,5),7]),
mean(Output.model$q75$Nsite[c(1,5),8]),
mean(Output.model$q75$Nsite[c(1,5),9]),
mean(Output.model$q75$Nsite[c(1,5),10]),
mean(Output.model$q75$Nsite[c(1,5),11]),
mean(Output.model$q75$Nsite[c(1,5),12])),
c(mean(Output.model$q97.5$Nsite[c(1,5),1]),
mean(Output.model$q97.5$Nsite[c(1,5),2]),
mean(Output.model$q97.5$Nsite[c(1,5),3]),
mean(Output.model$q97.5$Nsite[c(1,5),4]),
mean(Output.model$q97.5$Nsite[c(1,5),5]),
mean(Output.model$q97.5$Nsite[c(1,5),6]),
mean(Output.model$q97.5$Nsite[c(1,5),7]),
mean(Output.model$q97.5$Nsite[c(1,5),8]),
mean(Output.model$q97.5$Nsite[c(1,5),9]),
mean(Output.model$q97.5$Nsite[c(1,5),10]),
mean(Output.model$q97.5$Nsite[c(1,5),11]),
mean(Output.model$q97.5$Nsite[c(1,5),12])))
# Western Sampling locations
alpha3.val.west <- cbind(c(mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),1]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),2]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),3]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),4]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),5]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),6]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),7]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),8]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),9]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),10]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),11]),
mean(Output.model$q2.5$Nsite[c(2,3,4,6,7,8,9,10),12])),
c(mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),1]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),2]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),3]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),4]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),5]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),6]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),7]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),8]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),9]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),10]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),11]),
mean(Output.model$q25$Nsite[c(2,3,4,6,7,8,9,10),12])),
c(mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),1]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),2]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),3]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),4]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),5]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),6]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),7]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),8]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),9]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),10]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),11]),
mean(Output.model$mean$Nsite[c(2,3,4,6,7,8,9,10),12])),
c(mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),1]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),2]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),3]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),4]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),5]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),6]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),7]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),8]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),9]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),10]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),11]),
mean(Output.model$q75$Nsite[c(2,3,4,6,7,8,9,10),12])),
c(mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),1]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),2]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),3]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),4]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),5]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),6]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),7]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),8]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),9]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),10]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),11]),
mean(Output.model$q97.5$Nsite[c(2,3,4,6,7,8,9,10),12])))
# Eastern Sampling locations
alpha3.val.east <- cbind(c(mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),1]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),2]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),3]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),4]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),5]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),6]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),7]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),8]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),9]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),10]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),11]),
mean(Output.model$q2.5$Nsite[c(11,12,13,14,15,16,17,18,19),12])),
c(mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),1]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),2]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),3]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),4]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),5]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),6]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),7]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),8]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),9]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),10]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),11]),
mean(Output.model$q25$Nsite[c(11,12,13,14,15,16,17,18,19),12])),
c(mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),1]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),2]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),3]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),4]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),5]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),6]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),7]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),8]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),9]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),10]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),11]),
mean(Output.model$mean$Nsite[c(11,12,13,14,15,16,17,18,19),12])),
c(mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),1]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),2]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),3]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),4]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),5]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),6]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),7]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),8]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),9]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),10]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),11]),
mean(Output.model$q75$Nsite[c(11,12,13,14,15,16,17,18,19),12])),
c(mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),1]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),2]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),3]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),4]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),5]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),6]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),7]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),8]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),9]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),10]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),11]),
mean(Output.model$q97.5$Nsite[c(11,12,13,14,15,16,17,18,19),12])))
alpha3.val.control
alpha3.val.west
alpha3.val.east
# sampling years
years <- c(1,2,3,4,5,6,7,8,9,10,11,12)
## Figure3.1
values3.control <- data.frame(years, alpha3.val.control)
colnames(values3.control) <- c("years", "lower.alpha", "l25.alpha", "mean.alpha", "u75.alpha", "upper.alpha")
#Figure
Figure3.1 <- ggplot() +
scale_y_continuous(limits = c(11, 43), breaks = c(15,20,25,30,35,40))+
geom_errorbar(data = subset(values3.control), aes(x = years, ymin = lower.alpha, ymax = upper.alpha, color = "yellow"),
size = 0.75, width = 0) +
geom_errorbar(data = subset(values3.control), aes(x = years, ymin = l25.alpha, ymax = u75.alpha, color = "yellow"),
size = 2, width = 0) +
geom_errorbar(data = subset(values3.control), aes(x = years, ymin = mean.alpha, ymax = mean.alpha, color = "black"),
width = 0.3) +
scale_color_manual(name = "", values = c("yellow" = "#FFF685","black" = "black", "blue" = "#00BFFF")) +
scale_x_continuous(limits=c(0.8, 12.2), breaks = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12))+
theme_few() +
ylab("Species richness") +
xlab("Year") +
theme(plot.margin = unit(c(0.5, -0.1, 0.5, 0.3), "cm"),
text = element_text(size = 11),
panel.background = element_rect(fill = "transparent", color = NA),
plot.background = element_rect(fill = "transparent", color = NA),
axis.text.x = element_text(angle = 45, hjust = 1, vjust = 0.5, color = "black"),
axis.title.y = element_text(size=11, vjust = 2, face="bold"),
axis.title.x = element_text(size=11, vjust = -2),
legend.position = "none")
## Figure3.2
values3.west <- data.frame(years, alpha3.val.west)
colnames(values3.west) <- c("years", "lower.alpha", "l25.alpha", "mean.alpha", "u75.alpha", "upper.alpha")
#Identify irrigated and non-irrigated years
values3.west$sec <- values3.west$years < 10
values3.west$irri <- values3.west$years > 9
#Figure
Figure3.2 <- ggplot() +
scale_y_continuous(limits = c(11, 43), breaks = c(15,20,25,30,35,40))+
geom_errorbar(data = subset(values3.west, irri == TRUE), aes(x = years, ymin = lower.alpha, ymax = upper.alpha, color = "blue"),
size = 0.75, width = 0) +
geom_errorbar(data = subset(values3.west, irri == TRUE), aes(x = years, ymin = l25.alpha, ymax = u75.alpha, color = "blue"),
size = 2, width = 0) +
geom_errorbar(data = subset(values3.west, irri == TRUE), aes(x = years, ymin = mean.alpha, ymax = mean.alpha, color = "black"),
width = 0.3) +
geom_errorbar(data = subset(values3.west, sec == TRUE), aes(x = years, ymin = lower.alpha, ymax = upper.alpha, color = "yellow"),
width = 0, size = 0.75) +
geom_errorbar(data = subset(values3.west, sec == TRUE), aes(x = years, ymin = l25.alpha, ymax = u75.alpha, color = "yellow"),
width = 0, size = 2) +
geom_errorbar(data = subset(values3.west, sec == TRUE), aes(x = years, ymin = mean.alpha, ymax = mean.alpha, color = "black"),
width = 0.3) +
scale_color_manual(name = "", values = c("yellow" = "#FFF685","black" = "black", "blue" = "#00BFFF")) +
scale_x_continuous(limits=c(0.8, 12.2), breaks = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12))+
theme_few() +
ylab("") +
xlab("Year") +
theme(plot.margin = unit(c(0.5, 0.1, 0.5, 0.1), "cm"),
text = element_text(size = 11),
panel.background = element_rect(fill = "transparent", color = NA),
plot.background = element_rect(fill = "transparent", color = NA),
axis.text.x = element_text(angle = 45, hjust = 1, vjust = 0.5, color = "black"),
axis.title.y = element_text(size=11, vjust = 2, face="bold"),
axis.title.x = element_text(size=11, vjust = -2),
legend.position = "none")
Figure3.2
# Figure3.3
values3.east <- data.frame(years, alpha3.val.east)
colnames(values3.east) <- c("years", "lower.alpha", "l25.alpha", "mean.alpha", "u75.alpha", "upper.alpha")
#Identify irrigated and non-irrigated years
values3.east$sec <- values3.east$years < 3
values3.east$irri <- values3.east$years > 2
#Figure
Figure3.3 <- ggplot() +
scale_y_continuous(limits = c(11, 43), breaks = c(15,20,25,30,35,40))+
geom_errorbar(data = subset(values3.east, irri == TRUE), aes(x = years, ymin = lower.alpha, ymax = upper.alpha, color = "blue"),
size = 0.75, width = 0) +
geom_errorbar(data = subset(values3.east, irri == TRUE), aes(x = years, ymin = l25.alpha, ymax = u75.alpha, color = "blue"),
size = 2, width = 0) +
geom_errorbar(data = subset(values3.east, irri == TRUE), aes(x = years, ymin = mean.alpha, ymax = mean.alpha, color = "black"),
width = 0.3) +
geom_errorbar(data = subset(values3.east, sec == TRUE), aes(x = years, ymin = lower.alpha, ymax = upper.alpha, color = "yellow"),
width = 0, size = 0.75) +
geom_errorbar(data = subset(values3.east, sec == TRUE), aes(x = years, ymin = l25.alpha, ymax = u75.alpha, color = "yellow"),
width = 0, size = 2) +
geom_errorbar(data = subset(values3.east, sec == TRUE), aes(x = years, ymin = mean.alpha, ymax = mean.alpha, color = "black"),
width = 0.3) +
scale_color_manual(name = "", values = c("yellow" = "#FFF685","black" = "black", "blue" = "#00BFFF")) +
scale_x_continuous(limits=c(0.8, 12.2), breaks = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12))+
theme_few() +
ylab("") +
xlab("Year") +
theme(plot.margin = unit(c(0.5, 0.3, 0.5, -0.1), "cm"),
text = element_text(size = 11),
panel.background = element_rect(fill = "transparent", color = NA),
plot.background = element_rect(fill = "transparent", color = NA),
axis.text.x = element_text(angle = 45, hjust = 1, vjust = 0.5, color = "black"),
axis.title.y = element_text(size=11, vjust = 2, face="bold"),
axis.title.x = element_text(size=11, vjust = -2),
legend.position = "none")
Figure3.3
#Figure 3 # 6.8 x 2.25
plot_grid(Figure3.1 , Figure3.2, Figure3.3, labels = c("a)", "b)", "c)"), hjust = -0.5,
label_size = 11, ncol = 3, nrow = 1,rel_widths = c(1, 1, 1))
#######################################
############## Figure A.1 ##############
### Effect of occurrence probability ##
alphaS1.val <- cbind(c(mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),45]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),45]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q2.5$psi[c(11),c(1,2),45]),
mean(Output.model$q2.5$psi[c(12),c(1,2),45]),mean(Output.model$q2.5$psi[c(13),c(1,2),45]),mean(Output.model$q2.5$psi[c(14),c(1,2),45]),
mean(Output.model$q2.5$psi[c(15),c(1,2),45]),mean(Output.model$q2.5$psi[c(16),c(1,2),45]),mean(Output.model$q2.5$psi[c(17),c(1,2),45]),
mean(Output.model$q2.5$psi[c(18),c(1,2),45]),mean(Output.model$q2.5$psi[c(19),c(1,2),45]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),1]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),1]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q2.5$psi[c(11),c(1,2),1]),
mean(Output.model$q2.5$psi[c(12),c(1,2),1]),mean(Output.model$q2.5$psi[c(13),c(1,2),1]),mean(Output.model$q2.5$psi[c(14),c(1,2),1]),
mean(Output.model$q2.5$psi[c(15),c(1,2),1]),mean(Output.model$q2.5$psi[c(16),c(1,2),1]),mean(Output.model$q2.5$psi[c(17),c(1,2),1]),
mean(Output.model$q2.5$psi[c(18),c(1,2),1]),mean(Output.model$q2.5$psi[c(19),c(1,2),1]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),15]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),15]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q2.5$psi[c(11),c(1,2),15]),
mean(Output.model$q2.5$psi[c(12),c(1,2),15]),mean(Output.model$q2.5$psi[c(13),c(1,2),15]),mean(Output.model$q2.5$psi[c(14),c(1,2),15]),
mean(Output.model$q2.5$psi[c(15),c(1,2),15]),mean(Output.model$q2.5$psi[c(16),c(1,2),15]),mean(Output.model$q2.5$psi[c(17),c(1,2),15]),
mean(Output.model$q2.5$psi[c(18),c(1,2),15]),mean(Output.model$q2.5$psi[c(19),c(1,2),15]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),7]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),7]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q2.5$psi[c(11),c(1,2),7]),
mean(Output.model$q2.5$psi[c(12),c(1,2),7]),mean(Output.model$q2.5$psi[c(13),c(1,2),7]),mean(Output.model$q2.5$psi[c(14),c(1,2),7]),
mean(Output.model$q2.5$psi[c(15),c(1,2),7]),mean(Output.model$q2.5$psi[c(16),c(1,2),7]),mean(Output.model$q2.5$psi[c(17),c(1,2),7]),
mean(Output.model$q2.5$psi[c(18),c(1,2),7]),mean(Output.model$q2.5$psi[c(19),c(1,2),7]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),22]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),22]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q2.5$psi[c(11),c(1,2),22]),
mean(Output.model$q2.5$psi[c(12),c(1,2),22]),mean(Output.model$q2.5$psi[c(13),c(1,2),22]),mean(Output.model$q2.5$psi[c(14),c(1,2),22]),
mean(Output.model$q2.5$psi[c(15),c(1,2),22]),mean(Output.model$q2.5$psi[c(16),c(1,2),22]),mean(Output.model$q2.5$psi[c(17),c(1,2),22]),
mean(Output.model$q2.5$psi[c(18),c(1,2),22]),mean(Output.model$q2.5$psi[c(19),c(1,2),22]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),29]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),29]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q2.5$psi[c(11),c(1,2),29]),
mean(Output.model$q2.5$psi[c(12),c(1,2),29]),mean(Output.model$q2.5$psi[c(13),c(1,2),29]),mean(Output.model$q2.5$psi[c(14),c(1,2),29]),
mean(Output.model$q2.5$psi[c(15),c(1,2),29]),mean(Output.model$q2.5$psi[c(16),c(1,2),29]),mean(Output.model$q2.5$psi[c(17),c(1,2),29]),
mean(Output.model$q2.5$psi[c(18),c(1,2),29]),mean(Output.model$q2.5$psi[c(19),c(1,2),29]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),18]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),18]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q2.5$psi[c(11),c(1,2),18]),
mean(Output.model$q2.5$psi[c(12),c(1,2),18]),mean(Output.model$q2.5$psi[c(13),c(1,2),18]),mean(Output.model$q2.5$psi[c(14),c(1,2),18]),
mean(Output.model$q2.5$psi[c(15),c(1,2),18]),mean(Output.model$q2.5$psi[c(16),c(1,2),18]),mean(Output.model$q2.5$psi[c(17),c(1,2),18]),
mean(Output.model$q2.5$psi[c(18),c(1,2),18]),mean(Output.model$q2.5$psi[c(19),c(1,2),18]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),33]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),33]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q2.5$psi[c(11),c(1,2),33]),
mean(Output.model$q2.5$psi[c(12),c(1,2),33]),mean(Output.model$q2.5$psi[c(13),c(1,2),33]),mean(Output.model$q2.5$psi[c(14),c(1,2),33]),
mean(Output.model$q2.5$psi[c(15),c(1,2),33]),mean(Output.model$q2.5$psi[c(16),c(1,2),33]),mean(Output.model$q2.5$psi[c(17),c(1,2),33]),
mean(Output.model$q2.5$psi[c(18),c(1,2),33]),mean(Output.model$q2.5$psi[c(19),c(1,2),33]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),12]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),12]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q2.5$psi[c(11),c(1,2),12]),
mean(Output.model$q2.5$psi[c(12),c(1,2),12]),mean(Output.model$q2.5$psi[c(13),c(1,2),12]),mean(Output.model$q2.5$psi[c(14),c(1,2),12]),
mean(Output.model$q2.5$psi[c(15),c(1,2),12]),mean(Output.model$q2.5$psi[c(16),c(1,2),12]),mean(Output.model$q2.5$psi[c(17),c(1,2),12]),
mean(Output.model$q2.5$psi[c(18),c(1,2),12]),mean(Output.model$q2.5$psi[c(19),c(1,2),12]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),9]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),9]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q2.5$psi[c(11),c(1,2),9]),
mean(Output.model$q2.5$psi[c(12),c(1,2),9]),mean(Output.model$q2.5$psi[c(13),c(1,2),9]),mean(Output.model$q2.5$psi[c(14),c(1,2),9]),
mean(Output.model$q2.5$psi[c(15),c(1,2),9]),mean(Output.model$q2.5$psi[c(16),c(1,2),9]),mean(Output.model$q2.5$psi[c(17),c(1,2),9]),
mean(Output.model$q2.5$psi[c(18),c(1,2),9]),mean(Output.model$q2.5$psi[c(19),c(1,2),9]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),47]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),47]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q2.5$psi[c(11),c(1,2),47]),
mean(Output.model$q2.5$psi[c(12),c(1,2),47]),mean(Output.model$q2.5$psi[c(13),c(1,2),47]),mean(Output.model$q2.5$psi[c(14),c(1,2),47]),
mean(Output.model$q2.5$psi[c(15),c(1,2),47]),mean(Output.model$q2.5$psi[c(16),c(1,2),47]),mean(Output.model$q2.5$psi[c(17),c(1,2),47]),
mean(Output.model$q2.5$psi[c(18),c(1,2),47]),mean(Output.model$q2.5$psi[c(19),c(1,2),47]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),30]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),30]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q2.5$psi[c(11),c(1,2),30]),
mean(Output.model$q2.5$psi[c(12),c(1,2),30]),mean(Output.model$q2.5$psi[c(13),c(1,2),30]),mean(Output.model$q2.5$psi[c(14),c(1,2),30]),
mean(Output.model$q2.5$psi[c(15),c(1,2),30]),mean(Output.model$q2.5$psi[c(16),c(1,2),30]),mean(Output.model$q2.5$psi[c(17),c(1,2),30]),
mean(Output.model$q2.5$psi[c(18),c(1,2),30]),mean(Output.model$q2.5$psi[c(19),c(1,2),30]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),40]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),40]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q2.5$psi[c(11),c(1,2),40]),
mean(Output.model$q2.5$psi[c(12),c(1,2),40]),mean(Output.model$q2.5$psi[c(13),c(1,2),40]),mean(Output.model$q2.5$psi[c(14),c(1,2),40]),
mean(Output.model$q2.5$psi[c(15),c(1,2),40]),mean(Output.model$q2.5$psi[c(16),c(1,2),40]),mean(Output.model$q2.5$psi[c(17),c(1,2),40]),
mean(Output.model$q2.5$psi[c(18),c(1,2),40]),mean(Output.model$q2.5$psi[c(19),c(1,2),40]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),5]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),5]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q2.5$psi[c(11),c(1,2),5]),
mean(Output.model$q2.5$psi[c(12),c(1,2),5]),mean(Output.model$q2.5$psi[c(13),c(1,2),5]),mean(Output.model$q2.5$psi[c(14),c(1,2),5]),
mean(Output.model$q2.5$psi[c(15),c(1,2),5]),mean(Output.model$q2.5$psi[c(16),c(1,2),5]),mean(Output.model$q2.5$psi[c(17),c(1,2),5]),
mean(Output.model$q2.5$psi[c(18),c(1,2),5]),mean(Output.model$q2.5$psi[c(19),c(1,2),5]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),20]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),20]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q2.5$psi[c(11),c(1,2),20]),
mean(Output.model$q2.5$psi[c(12),c(1,2),20]),mean(Output.model$q2.5$psi[c(13),c(1,2),20]),mean(Output.model$q2.5$psi[c(14),c(1,2),20]),
mean(Output.model$q2.5$psi[c(15),c(1,2),20]),mean(Output.model$q2.5$psi[c(16),c(1,2),20]),mean(Output.model$q2.5$psi[c(17),c(1,2),20]),
mean(Output.model$q2.5$psi[c(18),c(1,2),20]),mean(Output.model$q2.5$psi[c(19),c(1,2),20]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),23]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),23]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q2.5$psi[c(11),c(1,2),23]),
mean(Output.model$q2.5$psi[c(12),c(1,2),23]),mean(Output.model$q2.5$psi[c(13),c(1,2),23]),mean(Output.model$q2.5$psi[c(14),c(1,2),23]),
mean(Output.model$q2.5$psi[c(15),c(1,2),23]),mean(Output.model$q2.5$psi[c(16),c(1,2),23]),mean(Output.model$q2.5$psi[c(17),c(1,2),23]),
mean(Output.model$q2.5$psi[c(18),c(1,2),23]),mean(Output.model$q2.5$psi[c(19),c(1,2),23]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),3]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),3]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q2.5$psi[c(11),c(1,2),3]),
mean(Output.model$q2.5$psi[c(12),c(1,2),3]),mean(Output.model$q2.5$psi[c(13),c(1,2),3]),mean(Output.model$q2.5$psi[c(14),c(1,2),3]),
mean(Output.model$q2.5$psi[c(15),c(1,2),3]),mean(Output.model$q2.5$psi[c(16),c(1,2),3]),mean(Output.model$q2.5$psi[c(17),c(1,2),3]),
mean(Output.model$q2.5$psi[c(18),c(1,2),3]),mean(Output.model$q2.5$psi[c(19),c(1,2),3]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),27]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),27]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q2.5$psi[c(11),c(1,2),27]),
mean(Output.model$q2.5$psi[c(12),c(1,2),27]),mean(Output.model$q2.5$psi[c(13),c(1,2),27]),mean(Output.model$q2.5$psi[c(14),c(1,2),27]),
mean(Output.model$q2.5$psi[c(15),c(1,2),27]),mean(Output.model$q2.5$psi[c(16),c(1,2),27]),mean(Output.model$q2.5$psi[c(17),c(1,2),27]),
mean(Output.model$q2.5$psi[c(18),c(1,2),27]),mean(Output.model$q2.5$psi[c(19),c(1,2),27]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),26]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),26]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q2.5$psi[c(11),c(1,2),26]),
mean(Output.model$q2.5$psi[c(12),c(1,2),26]),mean(Output.model$q2.5$psi[c(13),c(1,2),26]),mean(Output.model$q2.5$psi[c(14),c(1,2),26]),
mean(Output.model$q2.5$psi[c(15),c(1,2),26]),mean(Output.model$q2.5$psi[c(16),c(1,2),26]),mean(Output.model$q2.5$psi[c(17),c(1,2),26]),
mean(Output.model$q2.5$psi[c(18),c(1,2),26]),mean(Output.model$q2.5$psi[c(19),c(1,2),26]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),44]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),44]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q2.5$psi[c(11),c(1,2),44]),
mean(Output.model$q2.5$psi[c(12),c(1,2),44]),mean(Output.model$q2.5$psi[c(13),c(1,2),44]),mean(Output.model$q2.5$psi[c(14),c(1,2),44]),
mean(Output.model$q2.5$psi[c(15),c(1,2),44]),mean(Output.model$q2.5$psi[c(16),c(1,2),44]),mean(Output.model$q2.5$psi[c(17),c(1,2),44]),
mean(Output.model$q2.5$psi[c(18),c(1,2),44]),mean(Output.model$q2.5$psi[c(19),c(1,2),44]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),42]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),42]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q2.5$psi[c(11),c(1,2),42]),
mean(Output.model$q2.5$psi[c(12),c(1,2),42]),mean(Output.model$q2.5$psi[c(13),c(1,2),42]),mean(Output.model$q2.5$psi[c(14),c(1,2),42]),
mean(Output.model$q2.5$psi[c(15),c(1,2),42]),mean(Output.model$q2.5$psi[c(16),c(1,2),42]),mean(Output.model$q2.5$psi[c(17),c(1,2),42]),
mean(Output.model$q2.5$psi[c(18),c(1,2),42]),mean(Output.model$q2.5$psi[c(19),c(1,2),42]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),43]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),43]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q2.5$psi[c(11),c(1,2),43]),
mean(Output.model$q2.5$psi[c(12),c(1,2),43]),mean(Output.model$q2.5$psi[c(13),c(1,2),43]),mean(Output.model$q2.5$psi[c(14),c(1,2),43]),
mean(Output.model$q2.5$psi[c(15),c(1,2),43]),mean(Output.model$q2.5$psi[c(16),c(1,2),43]),mean(Output.model$q2.5$psi[c(17),c(1,2),43]),
mean(Output.model$q2.5$psi[c(18),c(1,2),43]),mean(Output.model$q2.5$psi[c(19),c(1,2),43]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),38]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),38]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q2.5$psi[c(11),c(1,2),38]),
mean(Output.model$q2.5$psi[c(12),c(1,2),38]),mean(Output.model$q2.5$psi[c(13),c(1,2),38]),mean(Output.model$q2.5$psi[c(14),c(1,2),38]),
mean(Output.model$q2.5$psi[c(15),c(1,2),38]),mean(Output.model$q2.5$psi[c(16),c(1,2),38]),mean(Output.model$q2.5$psi[c(17),c(1,2),38]),
mean(Output.model$q2.5$psi[c(18),c(1,2),38]),mean(Output.model$q2.5$psi[c(19),c(1,2),38]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),24]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),24]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q2.5$psi[c(11),c(1,2),24]),
mean(Output.model$q2.5$psi[c(12),c(1,2),24]),mean(Output.model$q2.5$psi[c(13),c(1,2),24]),mean(Output.model$q2.5$psi[c(14),c(1,2),24]),
mean(Output.model$q2.5$psi[c(15),c(1,2),24]),mean(Output.model$q2.5$psi[c(16),c(1,2),24]),mean(Output.model$q2.5$psi[c(17),c(1,2),24]),
mean(Output.model$q2.5$psi[c(18),c(1,2),24]),mean(Output.model$q2.5$psi[c(19),c(1,2),24]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),37]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),37]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q2.5$psi[c(11),c(1,2),37]),
mean(Output.model$q2.5$psi[c(12),c(1,2),37]),mean(Output.model$q2.5$psi[c(13),c(1,2),37]),mean(Output.model$q2.5$psi[c(14),c(1,2),37]),
mean(Output.model$q2.5$psi[c(15),c(1,2),37]),mean(Output.model$q2.5$psi[c(16),c(1,2),37]),mean(Output.model$q2.5$psi[c(17),c(1,2),37]),
mean(Output.model$q2.5$psi[c(18),c(1,2),37]),mean(Output.model$q2.5$psi[c(19),c(1,2),37]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),13]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),13]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q2.5$psi[c(11),c(1,2),13]),
mean(Output.model$q2.5$psi[c(12),c(1,2),13]),mean(Output.model$q2.5$psi[c(13),c(1,2),13]),mean(Output.model$q2.5$psi[c(14),c(1,2),13]),
mean(Output.model$q2.5$psi[c(15),c(1,2),13]),mean(Output.model$q2.5$psi[c(16),c(1,2),13]),mean(Output.model$q2.5$psi[c(17),c(1,2),13]),
mean(Output.model$q2.5$psi[c(18),c(1,2),13]),mean(Output.model$q2.5$psi[c(19),c(1,2),13]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),21]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),21]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q2.5$psi[c(11),c(1,2),21]),
mean(Output.model$q2.5$psi[c(12),c(1,2),21]),mean(Output.model$q2.5$psi[c(13),c(1,2),21]),mean(Output.model$q2.5$psi[c(14),c(1,2),21]),
mean(Output.model$q2.5$psi[c(15),c(1,2),21]),mean(Output.model$q2.5$psi[c(16),c(1,2),21]),mean(Output.model$q2.5$psi[c(17),c(1,2),21]),
mean(Output.model$q2.5$psi[c(18),c(1,2),21]),mean(Output.model$q2.5$psi[c(19),c(1,2),21]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),34]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),34]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q2.5$psi[c(11),c(1,2),34]),
mean(Output.model$q2.5$psi[c(12),c(1,2),34]),mean(Output.model$q2.5$psi[c(13),c(1,2),34]),mean(Output.model$q2.5$psi[c(14),c(1,2),34]),
mean(Output.model$q2.5$psi[c(15),c(1,2),34]),mean(Output.model$q2.5$psi[c(16),c(1,2),34]),mean(Output.model$q2.5$psi[c(17),c(1,2),34]),
mean(Output.model$q2.5$psi[c(18),c(1,2),34]),mean(Output.model$q2.5$psi[c(19),c(1,2),34]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),2]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),2]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q2.5$psi[c(11),c(1,2),2]),
mean(Output.model$q2.5$psi[c(12),c(1,2),2]),mean(Output.model$q2.5$psi[c(13),c(1,2),2]),mean(Output.model$q2.5$psi[c(14),c(1,2),2]),
mean(Output.model$q2.5$psi[c(15),c(1,2),2]),mean(Output.model$q2.5$psi[c(16),c(1,2),2]),mean(Output.model$q2.5$psi[c(17),c(1,2),2]),
mean(Output.model$q2.5$psi[c(18),c(1,2),2]),mean(Output.model$q2.5$psi[c(19),c(1,2),2]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),32]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),32]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q2.5$psi[c(11),c(1,2),32]),
mean(Output.model$q2.5$psi[c(12),c(1,2),32]),mean(Output.model$q2.5$psi[c(13),c(1,2),32]),mean(Output.model$q2.5$psi[c(14),c(1,2),32]),
mean(Output.model$q2.5$psi[c(15),c(1,2),32]),mean(Output.model$q2.5$psi[c(16),c(1,2),32]),mean(Output.model$q2.5$psi[c(17),c(1,2),32]),
mean(Output.model$q2.5$psi[c(18),c(1,2),32]),mean(Output.model$q2.5$psi[c(19),c(1,2),32]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),17]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),17]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q2.5$psi[c(11),c(1,2),17]),
mean(Output.model$q2.5$psi[c(12),c(1,2),17]),mean(Output.model$q2.5$psi[c(13),c(1,2),17]),mean(Output.model$q2.5$psi[c(14),c(1,2),17]),
mean(Output.model$q2.5$psi[c(15),c(1,2),17]),mean(Output.model$q2.5$psi[c(16),c(1,2),17]),mean(Output.model$q2.5$psi[c(17),c(1,2),17]),
mean(Output.model$q2.5$psi[c(18),c(1,2),17]),mean(Output.model$q2.5$psi[c(19),c(1,2),17]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),4]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),4]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q2.5$psi[c(11),c(1,2),4]),
mean(Output.model$q2.5$psi[c(12),c(1,2),4]),mean(Output.model$q2.5$psi[c(13),c(1,2),4]),mean(Output.model$q2.5$psi[c(14),c(1,2),4]),
mean(Output.model$q2.5$psi[c(15),c(1,2),4]),mean(Output.model$q2.5$psi[c(16),c(1,2),4]),mean(Output.model$q2.5$psi[c(17),c(1,2),4]),
mean(Output.model$q2.5$psi[c(18),c(1,2),4]),mean(Output.model$q2.5$psi[c(19),c(1,2),4]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),8]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),8]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q2.5$psi[c(11),c(1,2),8]),
mean(Output.model$q2.5$psi[c(12),c(1,2),8]),mean(Output.model$q2.5$psi[c(13),c(1,2),8]),mean(Output.model$q2.5$psi[c(14),c(1,2),8]),
mean(Output.model$q2.5$psi[c(15),c(1,2),8]),mean(Output.model$q2.5$psi[c(16),c(1,2),8]),mean(Output.model$q2.5$psi[c(17),c(1,2),8]),
mean(Output.model$q2.5$psi[c(18),c(1,2),8]),mean(Output.model$q2.5$psi[c(19),c(1,2),8]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),39]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),39]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q2.5$psi[c(11),c(1,2),39]),
mean(Output.model$q2.5$psi[c(12),c(1,2),39]),mean(Output.model$q2.5$psi[c(13),c(1,2),39]),mean(Output.model$q2.5$psi[c(14),c(1,2),39]),
mean(Output.model$q2.5$psi[c(15),c(1,2),39]),mean(Output.model$q2.5$psi[c(16),c(1,2),39]),mean(Output.model$q2.5$psi[c(17),c(1,2),39]),
mean(Output.model$q2.5$psi[c(18),c(1,2),39]),mean(Output.model$q2.5$psi[c(19),c(1,2),39]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),10]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),10]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q2.5$psi[c(11),c(1,2),10]),
mean(Output.model$q2.5$psi[c(12),c(1,2),10]),mean(Output.model$q2.5$psi[c(13),c(1,2),10]),mean(Output.model$q2.5$psi[c(14),c(1,2),10]),
mean(Output.model$q2.5$psi[c(15),c(1,2),10]),mean(Output.model$q2.5$psi[c(16),c(1,2),10]),mean(Output.model$q2.5$psi[c(17),c(1,2),10]),
mean(Output.model$q2.5$psi[c(18),c(1,2),10]),mean(Output.model$q2.5$psi[c(19),c(1,2),10]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),19]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),19]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q2.5$psi[c(11),c(1,2),19]),
mean(Output.model$q2.5$psi[c(12),c(1,2),19]),mean(Output.model$q2.5$psi[c(13),c(1,2),19]),mean(Output.model$q2.5$psi[c(14),c(1,2),19]),
mean(Output.model$q2.5$psi[c(15),c(1,2),19]),mean(Output.model$q2.5$psi[c(16),c(1,2),19]),mean(Output.model$q2.5$psi[c(17),c(1,2),19]),
mean(Output.model$q2.5$psi[c(18),c(1,2),19]),mean(Output.model$q2.5$psi[c(19),c(1,2),19]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),35]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),35]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q2.5$psi[c(11),c(1,2),35]),
mean(Output.model$q2.5$psi[c(12),c(1,2),35]),mean(Output.model$q2.5$psi[c(13),c(1,2),35]),mean(Output.model$q2.5$psi[c(14),c(1,2),35]),
mean(Output.model$q2.5$psi[c(15),c(1,2),35]),mean(Output.model$q2.5$psi[c(16),c(1,2),35]),mean(Output.model$q2.5$psi[c(17),c(1,2),35]),
mean(Output.model$q2.5$psi[c(18),c(1,2),35]),mean(Output.model$q2.5$psi[c(19),c(1,2),35]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),25]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),25]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q2.5$psi[c(11),c(1,2),25]),
mean(Output.model$q2.5$psi[c(12),c(1,2),25]),mean(Output.model$q2.5$psi[c(13),c(1,2),25]),mean(Output.model$q2.5$psi[c(14),c(1,2),25]),
mean(Output.model$q2.5$psi[c(15),c(1,2),25]),mean(Output.model$q2.5$psi[c(16),c(1,2),25]),mean(Output.model$q2.5$psi[c(17),c(1,2),25]),
mean(Output.model$q2.5$psi[c(18),c(1,2),25]),mean(Output.model$q2.5$psi[c(19),c(1,2),25]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),28]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),28]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q2.5$psi[c(11),c(1,2),28]),
mean(Output.model$q2.5$psi[c(12),c(1,2),28]),mean(Output.model$q2.5$psi[c(13),c(1,2),28]),mean(Output.model$q2.5$psi[c(14),c(1,2),28]),
mean(Output.model$q2.5$psi[c(15),c(1,2),28]),mean(Output.model$q2.5$psi[c(16),c(1,2),28]),mean(Output.model$q2.5$psi[c(17),c(1,2),28]),
mean(Output.model$q2.5$psi[c(18),c(1,2),28]),mean(Output.model$q2.5$psi[c(19),c(1,2),28]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),46]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),46]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q2.5$psi[c(11),c(1,2),46]),
mean(Output.model$q2.5$psi[c(12),c(1,2),46]),mean(Output.model$q2.5$psi[c(13),c(1,2),46]),mean(Output.model$q2.5$psi[c(14),c(1,2),46]),
mean(Output.model$q2.5$psi[c(15),c(1,2),46]),mean(Output.model$q2.5$psi[c(16),c(1,2),46]),mean(Output.model$q2.5$psi[c(17),c(1,2),46]),
mean(Output.model$q2.5$psi[c(18),c(1,2),46]),mean(Output.model$q2.5$psi[c(19),c(1,2),46]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),41]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),41]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q2.5$psi[c(11),c(1,2),41]),
mean(Output.model$q2.5$psi[c(12),c(1,2),41]),mean(Output.model$q2.5$psi[c(13),c(1,2),41]),mean(Output.model$q2.5$psi[c(14),c(1,2),41]),
mean(Output.model$q2.5$psi[c(15),c(1,2),41]),mean(Output.model$q2.5$psi[c(16),c(1,2),41]),mean(Output.model$q2.5$psi[c(17),c(1,2),41]),
mean(Output.model$q2.5$psi[c(18),c(1,2),41]),mean(Output.model$q2.5$psi[c(19),c(1,2),41]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),16]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),16]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q2.5$psi[c(11),c(1,2),16]),
mean(Output.model$q2.5$psi[c(12),c(1,2),16]),mean(Output.model$q2.5$psi[c(13),c(1,2),16]),mean(Output.model$q2.5$psi[c(14),c(1,2),16]),
mean(Output.model$q2.5$psi[c(15),c(1,2),16]),mean(Output.model$q2.5$psi[c(16),c(1,2),16]),mean(Output.model$q2.5$psi[c(17),c(1,2),16]),
mean(Output.model$q2.5$psi[c(18),c(1,2),16]),mean(Output.model$q2.5$psi[c(19),c(1,2),16]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),36]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),36]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q2.5$psi[c(11),c(1,2),36]),
mean(Output.model$q2.5$psi[c(12),c(1,2),36]),mean(Output.model$q2.5$psi[c(13),c(1,2),36]),mean(Output.model$q2.5$psi[c(14),c(1,2),36]),
mean(Output.model$q2.5$psi[c(15),c(1,2),36]),mean(Output.model$q2.5$psi[c(16),c(1,2),36]),mean(Output.model$q2.5$psi[c(17),c(1,2),36]),
mean(Output.model$q2.5$psi[c(18),c(1,2),36]),mean(Output.model$q2.5$psi[c(19),c(1,2),36]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),14]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),14]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),14]),mean(Output.model$q2.5$psi[c(11),c(1,2),14]),
mean(Output.model$q2.5$psi[c(12),c(1,2),14]),mean(Output.model$q2.5$psi[c(13),c(1,2),14]),mean(Output.model$q2.5$psi[c(14),c(1,2),14]),
mean(Output.model$q2.5$psi[c(15),c(1,2),14]),mean(Output.model$q2.5$psi[c(16),c(1,2),14]),mean(Output.model$q2.5$psi[c(17),c(1,2),14]),
mean(Output.model$q2.5$psi[c(18),c(1,2),14]),mean(Output.model$q2.5$psi[c(19),c(1,2),14]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),11]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),11]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),11]),mean(Output.model$q2.5$psi[c(11),c(1,2),11]),
mean(Output.model$q2.5$psi[c(12),c(1,2),11]),mean(Output.model$q2.5$psi[c(13),c(1,2),11]),mean(Output.model$q2.5$psi[c(14),c(1,2),11]),
mean(Output.model$q2.5$psi[c(15),c(1,2),11]),mean(Output.model$q2.5$psi[c(16),c(1,2),11]),mean(Output.model$q2.5$psi[c(17),c(1,2),11]),
mean(Output.model$q2.5$psi[c(18),c(1,2),11]),mean(Output.model$q2.5$psi[c(19),c(1,2),11]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),31]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),31]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),31]),mean(Output.model$q2.5$psi[c(11),c(1,2),31]),
mean(Output.model$q2.5$psi[c(12),c(1,2),31]),mean(Output.model$q2.5$psi[c(13),c(1,2),31]),mean(Output.model$q2.5$psi[c(14),c(1,2),31]),
mean(Output.model$q2.5$psi[c(15),c(1,2),31]),mean(Output.model$q2.5$psi[c(16),c(1,2),31]),mean(Output.model$q2.5$psi[c(17),c(1,2),31]),
mean(Output.model$q2.5$psi[c(18),c(1,2),31]),mean(Output.model$q2.5$psi[c(19),c(1,2),31]))),
mean(c(mean(Output.model$q2.5$psi[c(2),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(3),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(4),c(1,2,3,4,5,6,7,8,9),6]),
mean(Output.model$q2.5$psi[c(6),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(7),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(8),c(1,2,3,4,5,6,7,8,9),6]),
mean(Output.model$q2.5$psi[c(9),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(10),c(1,2,3,4,5,6,7,8,9),6]),mean(Output.model$q2.5$psi[c(11),c(1,2),6]),
mean(Output.model$q2.5$psi[c(12),c(1,2),6]),mean(Output.model$q2.5$psi[c(13),c(1,2),6]),mean(Output.model$q2.5$psi[c(14),c(1,2),6]),
mean(Output.model$q2.5$psi[c(15),c(1,2),6]),mean(Output.model$q2.5$psi[c(16),c(1,2),6]),mean(Output.model$q2.5$psi[c(17),c(1,2),6]),
mean(Output.model$q2.5$psi[c(18),c(1,2),6]),mean(Output.model$q2.5$psi[c(19),c(1,2),6])))),
c(mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),45]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),45]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),45]),mean(Output.model$q25$psi[c(11),c(1,2),45]),
mean(Output.model$q25$psi[c(12),c(1,2),45]),mean(Output.model$q25$psi[c(13),c(1,2),45]),mean(Output.model$q25$psi[c(14),c(1,2),45]),
mean(Output.model$q25$psi[c(15),c(1,2),45]),mean(Output.model$q25$psi[c(16),c(1,2),45]),mean(Output.model$q25$psi[c(17),c(1,2),45]),
mean(Output.model$q25$psi[c(18),c(1,2),45]),mean(Output.model$q25$psi[c(19),c(1,2),45]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),1]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),1]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),1]),mean(Output.model$q25$psi[c(11),c(1,2),1]),
mean(Output.model$q25$psi[c(12),c(1,2),1]),mean(Output.model$q25$psi[c(13),c(1,2),1]),mean(Output.model$q25$psi[c(14),c(1,2),1]),
mean(Output.model$q25$psi[c(15),c(1,2),1]),mean(Output.model$q25$psi[c(16),c(1,2),1]),mean(Output.model$q25$psi[c(17),c(1,2),1]),
mean(Output.model$q25$psi[c(18),c(1,2),1]),mean(Output.model$q25$psi[c(19),c(1,2),1]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),15]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),15]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),15]),mean(Output.model$q25$psi[c(11),c(1,2),15]),
mean(Output.model$q25$psi[c(12),c(1,2),15]),mean(Output.model$q25$psi[c(13),c(1,2),15]),mean(Output.model$q25$psi[c(14),c(1,2),15]),
mean(Output.model$q25$psi[c(15),c(1,2),15]),mean(Output.model$q25$psi[c(16),c(1,2),15]),mean(Output.model$q25$psi[c(17),c(1,2),15]),
mean(Output.model$q25$psi[c(18),c(1,2),15]),mean(Output.model$q25$psi[c(19),c(1,2),15]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),7]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),7]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),7]),mean(Output.model$q25$psi[c(11),c(1,2),7]),
mean(Output.model$q25$psi[c(12),c(1,2),7]),mean(Output.model$q25$psi[c(13),c(1,2),7]),mean(Output.model$q25$psi[c(14),c(1,2),7]),
mean(Output.model$q25$psi[c(15),c(1,2),7]),mean(Output.model$q25$psi[c(16),c(1,2),7]),mean(Output.model$q25$psi[c(17),c(1,2),7]),
mean(Output.model$q25$psi[c(18),c(1,2),7]),mean(Output.model$q25$psi[c(19),c(1,2),7]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),22]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),22]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),22]),mean(Output.model$q25$psi[c(11),c(1,2),22]),
mean(Output.model$q25$psi[c(12),c(1,2),22]),mean(Output.model$q25$psi[c(13),c(1,2),22]),mean(Output.model$q25$psi[c(14),c(1,2),22]),
mean(Output.model$q25$psi[c(15),c(1,2),22]),mean(Output.model$q25$psi[c(16),c(1,2),22]),mean(Output.model$q25$psi[c(17),c(1,2),22]),
mean(Output.model$q25$psi[c(18),c(1,2),22]),mean(Output.model$q25$psi[c(19),c(1,2),22]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),29]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),29]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),29]),mean(Output.model$q25$psi[c(11),c(1,2),29]),
mean(Output.model$q25$psi[c(12),c(1,2),29]),mean(Output.model$q25$psi[c(13),c(1,2),29]),mean(Output.model$q25$psi[c(14),c(1,2),29]),
mean(Output.model$q25$psi[c(15),c(1,2),29]),mean(Output.model$q25$psi[c(16),c(1,2),29]),mean(Output.model$q25$psi[c(17),c(1,2),29]),
mean(Output.model$q25$psi[c(18),c(1,2),29]),mean(Output.model$q25$psi[c(19),c(1,2),29]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),18]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),18]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),18]),mean(Output.model$q25$psi[c(11),c(1,2),18]),
mean(Output.model$q25$psi[c(12),c(1,2),18]),mean(Output.model$q25$psi[c(13),c(1,2),18]),mean(Output.model$q25$psi[c(14),c(1,2),18]),
mean(Output.model$q25$psi[c(15),c(1,2),18]),mean(Output.model$q25$psi[c(16),c(1,2),18]),mean(Output.model$q25$psi[c(17),c(1,2),18]),
mean(Output.model$q25$psi[c(18),c(1,2),18]),mean(Output.model$q25$psi[c(19),c(1,2),18]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),33]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),33]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),33]),mean(Output.model$q25$psi[c(11),c(1,2),33]),
mean(Output.model$q25$psi[c(12),c(1,2),33]),mean(Output.model$q25$psi[c(13),c(1,2),33]),mean(Output.model$q25$psi[c(14),c(1,2),33]),
mean(Output.model$q25$psi[c(15),c(1,2),33]),mean(Output.model$q25$psi[c(16),c(1,2),33]),mean(Output.model$q25$psi[c(17),c(1,2),33]),
mean(Output.model$q25$psi[c(18),c(1,2),33]),mean(Output.model$q25$psi[c(19),c(1,2),33]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),12]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),12]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),12]),mean(Output.model$q25$psi[c(11),c(1,2),12]),
mean(Output.model$q25$psi[c(12),c(1,2),12]),mean(Output.model$q25$psi[c(13),c(1,2),12]),mean(Output.model$q25$psi[c(14),c(1,2),12]),
mean(Output.model$q25$psi[c(15),c(1,2),12]),mean(Output.model$q25$psi[c(16),c(1,2),12]),mean(Output.model$q25$psi[c(17),c(1,2),12]),
mean(Output.model$q25$psi[c(18),c(1,2),12]),mean(Output.model$q25$psi[c(19),c(1,2),12]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),9]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),9]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),9]),mean(Output.model$q25$psi[c(11),c(1,2),9]),
mean(Output.model$q25$psi[c(12),c(1,2),9]),mean(Output.model$q25$psi[c(13),c(1,2),9]),mean(Output.model$q25$psi[c(14),c(1,2),9]),
mean(Output.model$q25$psi[c(15),c(1,2),9]),mean(Output.model$q25$psi[c(16),c(1,2),9]),mean(Output.model$q25$psi[c(17),c(1,2),9]),
mean(Output.model$q25$psi[c(18),c(1,2),9]),mean(Output.model$q25$psi[c(19),c(1,2),9]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),47]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),47]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),47]),mean(Output.model$q25$psi[c(11),c(1,2),47]),
mean(Output.model$q25$psi[c(12),c(1,2),47]),mean(Output.model$q25$psi[c(13),c(1,2),47]),mean(Output.model$q25$psi[c(14),c(1,2),47]),
mean(Output.model$q25$psi[c(15),c(1,2),47]),mean(Output.model$q25$psi[c(16),c(1,2),47]),mean(Output.model$q25$psi[c(17),c(1,2),47]),
mean(Output.model$q25$psi[c(18),c(1,2),47]),mean(Output.model$q25$psi[c(19),c(1,2),47]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),30]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),30]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),30]),mean(Output.model$q25$psi[c(11),c(1,2),30]),
mean(Output.model$q25$psi[c(12),c(1,2),30]),mean(Output.model$q25$psi[c(13),c(1,2),30]),mean(Output.model$q25$psi[c(14),c(1,2),30]),
mean(Output.model$q25$psi[c(15),c(1,2),30]),mean(Output.model$q25$psi[c(16),c(1,2),30]),mean(Output.model$q25$psi[c(17),c(1,2),30]),
mean(Output.model$q25$psi[c(18),c(1,2),30]),mean(Output.model$q25$psi[c(19),c(1,2),30]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),40]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),40]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),40]),mean(Output.model$q25$psi[c(11),c(1,2),40]),
mean(Output.model$q25$psi[c(12),c(1,2),40]),mean(Output.model$q25$psi[c(13),c(1,2),40]),mean(Output.model$q25$psi[c(14),c(1,2),40]),
mean(Output.model$q25$psi[c(15),c(1,2),40]),mean(Output.model$q25$psi[c(16),c(1,2),40]),mean(Output.model$q25$psi[c(17),c(1,2),40]),
mean(Output.model$q25$psi[c(18),c(1,2),40]),mean(Output.model$q25$psi[c(19),c(1,2),40]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),5]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),5]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),5]),mean(Output.model$q25$psi[c(11),c(1,2),5]),
mean(Output.model$q25$psi[c(12),c(1,2),5]),mean(Output.model$q25$psi[c(13),c(1,2),5]),mean(Output.model$q25$psi[c(14),c(1,2),5]),
mean(Output.model$q25$psi[c(15),c(1,2),5]),mean(Output.model$q25$psi[c(16),c(1,2),5]),mean(Output.model$q25$psi[c(17),c(1,2),5]),
mean(Output.model$q25$psi[c(18),c(1,2),5]),mean(Output.model$q25$psi[c(19),c(1,2),5]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),20]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),20]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),20]),mean(Output.model$q25$psi[c(11),c(1,2),20]),
mean(Output.model$q25$psi[c(12),c(1,2),20]),mean(Output.model$q25$psi[c(13),c(1,2),20]),mean(Output.model$q25$psi[c(14),c(1,2),20]),
mean(Output.model$q25$psi[c(15),c(1,2),20]),mean(Output.model$q25$psi[c(16),c(1,2),20]),mean(Output.model$q25$psi[c(17),c(1,2),20]),
mean(Output.model$q25$psi[c(18),c(1,2),20]),mean(Output.model$q25$psi[c(19),c(1,2),20]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),23]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),23]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),23]),mean(Output.model$q25$psi[c(11),c(1,2),23]),
mean(Output.model$q25$psi[c(12),c(1,2),23]),mean(Output.model$q25$psi[c(13),c(1,2),23]),mean(Output.model$q25$psi[c(14),c(1,2),23]),
mean(Output.model$q25$psi[c(15),c(1,2),23]),mean(Output.model$q25$psi[c(16),c(1,2),23]),mean(Output.model$q25$psi[c(17),c(1,2),23]),
mean(Output.model$q25$psi[c(18),c(1,2),23]),mean(Output.model$q25$psi[c(19),c(1,2),23]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),3]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),3]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),3]),mean(Output.model$q25$psi[c(11),c(1,2),3]),
mean(Output.model$q25$psi[c(12),c(1,2),3]),mean(Output.model$q25$psi[c(13),c(1,2),3]),mean(Output.model$q25$psi[c(14),c(1,2),3]),
mean(Output.model$q25$psi[c(15),c(1,2),3]),mean(Output.model$q25$psi[c(16),c(1,2),3]),mean(Output.model$q25$psi[c(17),c(1,2),3]),
mean(Output.model$q25$psi[c(18),c(1,2),3]),mean(Output.model$q25$psi[c(19),c(1,2),3]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),27]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),27]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),27]),mean(Output.model$q25$psi[c(11),c(1,2),27]),
mean(Output.model$q25$psi[c(12),c(1,2),27]),mean(Output.model$q25$psi[c(13),c(1,2),27]),mean(Output.model$q25$psi[c(14),c(1,2),27]),
mean(Output.model$q25$psi[c(15),c(1,2),27]),mean(Output.model$q25$psi[c(16),c(1,2),27]),mean(Output.model$q25$psi[c(17),c(1,2),27]),
mean(Output.model$q25$psi[c(18),c(1,2),27]),mean(Output.model$q25$psi[c(19),c(1,2),27]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),26]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),26]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),26]),mean(Output.model$q25$psi[c(11),c(1,2),26]),
mean(Output.model$q25$psi[c(12),c(1,2),26]),mean(Output.model$q25$psi[c(13),c(1,2),26]),mean(Output.model$q25$psi[c(14),c(1,2),26]),
mean(Output.model$q25$psi[c(15),c(1,2),26]),mean(Output.model$q25$psi[c(16),c(1,2),26]),mean(Output.model$q25$psi[c(17),c(1,2),26]),
mean(Output.model$q25$psi[c(18),c(1,2),26]),mean(Output.model$q25$psi[c(19),c(1,2),26]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),44]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),44]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),44]),mean(Output.model$q25$psi[c(11),c(1,2),44]),
mean(Output.model$q25$psi[c(12),c(1,2),44]),mean(Output.model$q25$psi[c(13),c(1,2),44]),mean(Output.model$q25$psi[c(14),c(1,2),44]),
mean(Output.model$q25$psi[c(15),c(1,2),44]),mean(Output.model$q25$psi[c(16),c(1,2),44]),mean(Output.model$q25$psi[c(17),c(1,2),44]),
mean(Output.model$q25$psi[c(18),c(1,2),44]),mean(Output.model$q25$psi[c(19),c(1,2),44]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),42]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),42]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),42]),mean(Output.model$q25$psi[c(11),c(1,2),42]),
mean(Output.model$q25$psi[c(12),c(1,2),42]),mean(Output.model$q25$psi[c(13),c(1,2),42]),mean(Output.model$q25$psi[c(14),c(1,2),42]),
mean(Output.model$q25$psi[c(15),c(1,2),42]),mean(Output.model$q25$psi[c(16),c(1,2),42]),mean(Output.model$q25$psi[c(17),c(1,2),42]),
mean(Output.model$q25$psi[c(18),c(1,2),42]),mean(Output.model$q25$psi[c(19),c(1,2),42]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),43]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),43]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),43]),mean(Output.model$q25$psi[c(11),c(1,2),43]),
mean(Output.model$q25$psi[c(12),c(1,2),43]),mean(Output.model$q25$psi[c(13),c(1,2),43]),mean(Output.model$q25$psi[c(14),c(1,2),43]),
mean(Output.model$q25$psi[c(15),c(1,2),43]),mean(Output.model$q25$psi[c(16),c(1,2),43]),mean(Output.model$q25$psi[c(17),c(1,2),43]),
mean(Output.model$q25$psi[c(18),c(1,2),43]),mean(Output.model$q25$psi[c(19),c(1,2),43]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),38]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),38]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),38]),mean(Output.model$q25$psi[c(11),c(1,2),38]),
mean(Output.model$q25$psi[c(12),c(1,2),38]),mean(Output.model$q25$psi[c(13),c(1,2),38]),mean(Output.model$q25$psi[c(14),c(1,2),38]),
mean(Output.model$q25$psi[c(15),c(1,2),38]),mean(Output.model$q25$psi[c(16),c(1,2),38]),mean(Output.model$q25$psi[c(17),c(1,2),38]),
mean(Output.model$q25$psi[c(18),c(1,2),38]),mean(Output.model$q25$psi[c(19),c(1,2),38]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),24]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),24]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),24]),mean(Output.model$q25$psi[c(11),c(1,2),24]),
mean(Output.model$q25$psi[c(12),c(1,2),24]),mean(Output.model$q25$psi[c(13),c(1,2),24]),mean(Output.model$q25$psi[c(14),c(1,2),24]),
mean(Output.model$q25$psi[c(15),c(1,2),24]),mean(Output.model$q25$psi[c(16),c(1,2),24]),mean(Output.model$q25$psi[c(17),c(1,2),24]),
mean(Output.model$q25$psi[c(18),c(1,2),24]),mean(Output.model$q25$psi[c(19),c(1,2),24]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),37]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),37]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),37]),mean(Output.model$q25$psi[c(11),c(1,2),37]),
mean(Output.model$q25$psi[c(12),c(1,2),37]),mean(Output.model$q25$psi[c(13),c(1,2),37]),mean(Output.model$q25$psi[c(14),c(1,2),37]),
mean(Output.model$q25$psi[c(15),c(1,2),37]),mean(Output.model$q25$psi[c(16),c(1,2),37]),mean(Output.model$q25$psi[c(17),c(1,2),37]),
mean(Output.model$q25$psi[c(18),c(1,2),37]),mean(Output.model$q25$psi[c(19),c(1,2),37]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),13]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),13]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),13]),mean(Output.model$q25$psi[c(11),c(1,2),13]),
mean(Output.model$q25$psi[c(12),c(1,2),13]),mean(Output.model$q25$psi[c(13),c(1,2),13]),mean(Output.model$q25$psi[c(14),c(1,2),13]),
mean(Output.model$q25$psi[c(15),c(1,2),13]),mean(Output.model$q25$psi[c(16),c(1,2),13]),mean(Output.model$q25$psi[c(17),c(1,2),13]),
mean(Output.model$q25$psi[c(18),c(1,2),13]),mean(Output.model$q25$psi[c(19),c(1,2),13]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),21]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),21]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),21]),mean(Output.model$q25$psi[c(11),c(1,2),21]),
mean(Output.model$q25$psi[c(12),c(1,2),21]),mean(Output.model$q25$psi[c(13),c(1,2),21]),mean(Output.model$q25$psi[c(14),c(1,2),21]),
mean(Output.model$q25$psi[c(15),c(1,2),21]),mean(Output.model$q25$psi[c(16),c(1,2),21]),mean(Output.model$q25$psi[c(17),c(1,2),21]),
mean(Output.model$q25$psi[c(18),c(1,2),21]),mean(Output.model$q25$psi[c(19),c(1,2),21]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),34]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),34]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),34]),mean(Output.model$q25$psi[c(11),c(1,2),34]),
mean(Output.model$q25$psi[c(12),c(1,2),34]),mean(Output.model$q25$psi[c(13),c(1,2),34]),mean(Output.model$q25$psi[c(14),c(1,2),34]),
mean(Output.model$q25$psi[c(15),c(1,2),34]),mean(Output.model$q25$psi[c(16),c(1,2),34]),mean(Output.model$q25$psi[c(17),c(1,2),34]),
mean(Output.model$q25$psi[c(18),c(1,2),34]),mean(Output.model$q25$psi[c(19),c(1,2),34]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),2]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),2]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),2]),mean(Output.model$q25$psi[c(11),c(1,2),2]),
mean(Output.model$q25$psi[c(12),c(1,2),2]),mean(Output.model$q25$psi[c(13),c(1,2),2]),mean(Output.model$q25$psi[c(14),c(1,2),2]),
mean(Output.model$q25$psi[c(15),c(1,2),2]),mean(Output.model$q25$psi[c(16),c(1,2),2]),mean(Output.model$q25$psi[c(17),c(1,2),2]),
mean(Output.model$q25$psi[c(18),c(1,2),2]),mean(Output.model$q25$psi[c(19),c(1,2),2]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),32]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),32]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),32]),mean(Output.model$q25$psi[c(11),c(1,2),32]),
mean(Output.model$q25$psi[c(12),c(1,2),32]),mean(Output.model$q25$psi[c(13),c(1,2),32]),mean(Output.model$q25$psi[c(14),c(1,2),32]),
mean(Output.model$q25$psi[c(15),c(1,2),32]),mean(Output.model$q25$psi[c(16),c(1,2),32]),mean(Output.model$q25$psi[c(17),c(1,2),32]),
mean(Output.model$q25$psi[c(18),c(1,2),32]),mean(Output.model$q25$psi[c(19),c(1,2),32]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),17]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),17]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),17]),mean(Output.model$q25$psi[c(11),c(1,2),17]),
mean(Output.model$q25$psi[c(12),c(1,2),17]),mean(Output.model$q25$psi[c(13),c(1,2),17]),mean(Output.model$q25$psi[c(14),c(1,2),17]),
mean(Output.model$q25$psi[c(15),c(1,2),17]),mean(Output.model$q25$psi[c(16),c(1,2),17]),mean(Output.model$q25$psi[c(17),c(1,2),17]),
mean(Output.model$q25$psi[c(18),c(1,2),17]),mean(Output.model$q25$psi[c(19),c(1,2),17]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),4]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),4]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),4]),mean(Output.model$q25$psi[c(11),c(1,2),4]),
mean(Output.model$q25$psi[c(12),c(1,2),4]),mean(Output.model$q25$psi[c(13),c(1,2),4]),mean(Output.model$q25$psi[c(14),c(1,2),4]),
mean(Output.model$q25$psi[c(15),c(1,2),4]),mean(Output.model$q25$psi[c(16),c(1,2),4]),mean(Output.model$q25$psi[c(17),c(1,2),4]),
mean(Output.model$q25$psi[c(18),c(1,2),4]),mean(Output.model$q25$psi[c(19),c(1,2),4]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),8]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),8]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),8]),mean(Output.model$q25$psi[c(11),c(1,2),8]),
mean(Output.model$q25$psi[c(12),c(1,2),8]),mean(Output.model$q25$psi[c(13),c(1,2),8]),mean(Output.model$q25$psi[c(14),c(1,2),8]),
mean(Output.model$q25$psi[c(15),c(1,2),8]),mean(Output.model$q25$psi[c(16),c(1,2),8]),mean(Output.model$q25$psi[c(17),c(1,2),8]),
mean(Output.model$q25$psi[c(18),c(1,2),8]),mean(Output.model$q25$psi[c(19),c(1,2),8]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),39]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),39]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),39]),mean(Output.model$q25$psi[c(11),c(1,2),39]),
mean(Output.model$q25$psi[c(12),c(1,2),39]),mean(Output.model$q25$psi[c(13),c(1,2),39]),mean(Output.model$q25$psi[c(14),c(1,2),39]),
mean(Output.model$q25$psi[c(15),c(1,2),39]),mean(Output.model$q25$psi[c(16),c(1,2),39]),mean(Output.model$q25$psi[c(17),c(1,2),39]),
mean(Output.model$q25$psi[c(18),c(1,2),39]),mean(Output.model$q25$psi[c(19),c(1,2),39]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),10]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),10]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),10]),mean(Output.model$q25$psi[c(11),c(1,2),10]),
mean(Output.model$q25$psi[c(12),c(1,2),10]),mean(Output.model$q25$psi[c(13),c(1,2),10]),mean(Output.model$q25$psi[c(14),c(1,2),10]),
mean(Output.model$q25$psi[c(15),c(1,2),10]),mean(Output.model$q25$psi[c(16),c(1,2),10]),mean(Output.model$q25$psi[c(17),c(1,2),10]),
mean(Output.model$q25$psi[c(18),c(1,2),10]),mean(Output.model$q25$psi[c(19),c(1,2),10]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),19]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),19]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),19]),mean(Output.model$q25$psi[c(11),c(1,2),19]),
mean(Output.model$q25$psi[c(12),c(1,2),19]),mean(Output.model$q25$psi[c(13),c(1,2),19]),mean(Output.model$q25$psi[c(14),c(1,2),19]),
mean(Output.model$q25$psi[c(15),c(1,2),19]),mean(Output.model$q25$psi[c(16),c(1,2),19]),mean(Output.model$q25$psi[c(17),c(1,2),19]),
mean(Output.model$q25$psi[c(18),c(1,2),19]),mean(Output.model$q25$psi[c(19),c(1,2),19]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),35]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),35]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),35]),mean(Output.model$q25$psi[c(11),c(1,2),35]),
mean(Output.model$q25$psi[c(12),c(1,2),35]),mean(Output.model$q25$psi[c(13),c(1,2),35]),mean(Output.model$q25$psi[c(14),c(1,2),35]),
mean(Output.model$q25$psi[c(15),c(1,2),35]),mean(Output.model$q25$psi[c(16),c(1,2),35]),mean(Output.model$q25$psi[c(17),c(1,2),35]),
mean(Output.model$q25$psi[c(18),c(1,2),35]),mean(Output.model$q25$psi[c(19),c(1,2),35]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),25]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),25]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),25]),mean(Output.model$q25$psi[c(11),c(1,2),25]),
mean(Output.model$q25$psi[c(12),c(1,2),25]),mean(Output.model$q25$psi[c(13),c(1,2),25]),mean(Output.model$q25$psi[c(14),c(1,2),25]),
mean(Output.model$q25$psi[c(15),c(1,2),25]),mean(Output.model$q25$psi[c(16),c(1,2),25]),mean(Output.model$q25$psi[c(17),c(1,2),25]),
mean(Output.model$q25$psi[c(18),c(1,2),25]),mean(Output.model$q25$psi[c(19),c(1,2),25]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),28]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),28]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),28]),mean(Output.model$q25$psi[c(11),c(1,2),28]),
mean(Output.model$q25$psi[c(12),c(1,2),28]),mean(Output.model$q25$psi[c(13),c(1,2),28]),mean(Output.model$q25$psi[c(14),c(1,2),28]),
mean(Output.model$q25$psi[c(15),c(1,2),28]),mean(Output.model$q25$psi[c(16),c(1,2),28]),mean(Output.model$q25$psi[c(17),c(1,2),28]),
mean(Output.model$q25$psi[c(18),c(1,2),28]),mean(Output.model$q25$psi[c(19),c(1,2),28]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),46]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),46]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),46]),mean(Output.model$q25$psi[c(11),c(1,2),46]),
mean(Output.model$q25$psi[c(12),c(1,2),46]),mean(Output.model$q25$psi[c(13),c(1,2),46]),mean(Output.model$q25$psi[c(14),c(1,2),46]),
mean(Output.model$q25$psi[c(15),c(1,2),46]),mean(Output.model$q25$psi[c(16),c(1,2),46]),mean(Output.model$q25$psi[c(17),c(1,2),46]),
mean(Output.model$q25$psi[c(18),c(1,2),46]),mean(Output.model$q25$psi[c(19),c(1,2),46]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),41]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),41]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),41]),mean(Output.model$q25$psi[c(11),c(1,2),41]),
mean(Output.model$q25$psi[c(12),c(1,2),41]),mean(Output.model$q25$psi[c(13),c(1,2),41]),mean(Output.model$q25$psi[c(14),c(1,2),41]),
mean(Output.model$q25$psi[c(15),c(1,2),41]),mean(Output.model$q25$psi[c(16),c(1,2),41]),mean(Output.model$q25$psi[c(17),c(1,2),41]),
mean(Output.model$q25$psi[c(18),c(1,2),41]),mean(Output.model$q25$psi[c(19),c(1,2),41]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),16]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),16]),
mean(Output.model$q25$psi[c(9),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(10),c(1,2,3,4,5,6,7,8,9),16]),mean(Output.model$q25$psi[c(11),c(1,2),16]),
mean(Output.model$q25$psi[c(12),c(1,2),16]),mean(Output.model$q25$psi[c(13),c(1,2),16]),mean(Output.model$q25$psi[c(14),c(1,2),16]),
mean(Output.model$q25$psi[c(15),c(1,2),16]),mean(Output.model$q25$psi[c(16),c(1,2),16]),mean(Output.model$q25$psi[c(17),c(1,2),16]),
mean(Output.model$q25$psi[c(18),c(1,2),16]),mean(Output.model$q25$psi[c(19),c(1,2),16]))),
mean(c(mean(Output.model$q25$psi[c(2),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q25$psi[c(3),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q25$psi[c(4),c(1,2,3,4,5,6,7,8,9),36]),
mean(Output.model$q25$psi[c(6),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q25$psi[c(7),c(1,2,3,4,5,6,7,8,9),36]),mean(Output.model$q25$psi[c(8),c(1,2,3,4,5,6,7,8,9),36]),