GRUPO SULFATOS: LISBETH BELIZARIO, JULITZA CUTIPA, RUSSELL QUISPE
library(ggplot2)getwd()[1] "C:/Users/Lisbeth/Documents"file.choose()[1] "C:\\Users\\Lisbeth\\Downloads\\DISEÑO EXPERIMENTAL (TRABAJOS)\\SULFATOS.xlsx"ruta_propulsora <- "C:\\Users\\Lisbeth\\Downloads\\DISEÑO EXPERIMENTAL (TRABAJOS)\\SULFATOS.xlsx"excel_sheets(ruta_propulsora)[1] "Hoja1"Propul<-read_excel(ruta_propulsora)print(head(Propul))View(Propul)attach(Propul)names(Propul)[1] "Año" "Distrito" "sulfatos (mg/L)"summary(Propul) Año Distrito sulfatos (mg/L)
Min. :2012 Length:1613 Length:1613
1st Qu.:2014 Class :character Class :character
Median :2016 Mode :character Mode :character
Mean :2017
3rd Qu.:2021
Max. :2027 str(Propul)tibble [1,613 × 3] (S3: tbl_df/tbl/data.frame)
$ Año : num [1:1613] 2012 2012 2012 2012 2012 ...
$ Distrito : chr [1:1613] "Chupa" "Ilave" "Orurillo" "Limbani" ...
$ sulfatos (mg/L): chr [1:1613] "4" "4" "760" "16" ...plot(`sulfatos (mg/L)`)Warning: NAs introduced by coercion
plot(SULFATOS)
AÑ <- factor(Propul$`Año`)DIS <-as.vector(Propul$`Distrito`)SULFATES <-as.vector(Propul$`sulfatos (mg/L)`)SULFATES1 >-as.numeric(SULFATES)Warning: NAs introduced by coercion [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[10] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[19] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[28] TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[37] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[46] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[55] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[64] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[73] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[82] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[91] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[100] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[109] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[118] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[127] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[136] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[145] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[154] TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE
[163] FALSE FALSE FALSE FALSE FALSE FALSE TRUE TRUE TRUE
[172] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[181] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[190] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[199] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[208] TRUE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE
[217] TRUE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE
[226] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[235] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[244] TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE
[253] FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[262] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[271] FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[280] TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE
[289] FALSE FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[298] FALSE FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[307] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[316] FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[325] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[334] TRUE FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[343] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[352] TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE TRUE
[361] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[370] TRUE TRUE FALSE FALSE TRUE TRUE TRUE TRUE TRUE
[379] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[388] TRUE TRUE TRUE FALSE TRUE TRUE FALSE FALSE FALSE
[397] FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[406] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE
[415] FALSE FALSE FALSE NA NA FALSE FALSE NA FALSE
[424] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[433] FALSE FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE
[442] FALSE FALSE FALSE FALSE FALSE FALSE NA TRUE TRUE
[451] NA NA TRUE NA TRUE TRUE TRUE NA TRUE
[460] TRUE TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[469] FALSE FALSE FALSE FALSE FALSE TRUE NA TRUE TRUE
[478] TRUE TRUE NA NA NA NA TRUE NA FALSE
[487] FALSE FALSE FALSE TRUE NA FALSE FALSE NA TRUE
[496] FALSE TRUE NA TRUE NA TRUE TRUE TRUE TRUE
[505] TRUE NA TRUE NA TRUE TRUE NA TRUE NA
[514] NA NA TRUE TRUE NA TRUE NA TRUE NA
[523] TRUE TRUE TRUE FALSE FALSE FALSE TRUE FALSE FALSE
[532] TRUE NA TRUE TRUE NA TRUE NA NA NA
[541] TRUE NA TRUE FALSE FALSE NA TRUE TRUE TRUE
[550] TRUE TRUE TRUE NA TRUE TRUE TRUE TRUE FALSE
[559] FALSE NA TRUE TRUE FALSE TRUE TRUE TRUE TRUE
[568] TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE FALSE
[577] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[586] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[595] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[604] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[613] FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE FALSE
[622] FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[631] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[640] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[649] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE
[658] FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[667] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[676] TRUE FALSE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[685] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[694] FALSE FALSE FALSE TRUE TRUE FALSE TRUE TRUE TRUE
[703] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[712] FALSE TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE
[721] TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE FALSE
[730] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[739] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[748] TRUE FALSE FALSE FALSE FALSE TRUE TRUE TRUE TRUE
[757] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[766] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[775] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[784] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[793] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[802] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[811] FALSE FALSE FALSE FALSE TRUE TRUE TRUE TRUE TRUE
[820] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[829] TRUE FALSE TRUE FALSE TRUE TRUE TRUE FALSE TRUE
[838] FALSE TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE
[847] TRUE TRUE TRUE FALSE FALSE FALSE FALSE FALSE FALSE
[856] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[865] TRUE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[874] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[883] TRUE TRUE TRUE FALSE TRUE TRUE TRUE TRUE TRUE
[892] FALSE NA NA NA NA NA NA NA NA
[901] NA NA NA NA NA NA NA NA FALSE
[910] FALSE FALSE NA NA NA NA NA NA NA
[919] NA FALSE NA FALSE FALSE NA NA NA NA
[928] NA NA NA NA NA NA FALSE FALSE FALSE
[937] FALSE NA NA NA NA NA FALSE FALSE NA
[946] NA NA NA NA NA NA NA NA NA
[955] NA NA NA NA FALSE FALSE FALSE NA FALSE
[964] FALSE FALSE NA NA NA NA NA NA FALSE
[973] FALSE NA NA NA FALSE NA FALSE FALSE FALSE
[982] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[991] TRUE NA NA TRUE FALSE FALSE FALSE FALSE FALSE
[1000] FALSE
[ reached getOption("max.print") -- omitted 613 entries ]par(mfrow=c(1,1))boxplot(split(SULFATES1,AÑ),xlab="Año", ylab="Sulfatos")
boxplot(split(SULFATES1,DIS),xlab="Distrito", ylab="Sulfatos")
Ra.aov<-aov(SULFATES1 ~ AÑ+DIS)anova(Ra.aov)Analysis of Variance Table
Response: SULFATES1
Df Sum Sq Mean Sq F value Pr(>F)
AÑ 9 36159326 4017703 4.7272 3.298e-06 ***
DIS 160 4489964541 28062278 33.0176 < 2.2e-16 ***
Residuals 1339 1138041877 849919
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1ggplot(Ra.aov, aes(x = AÑ, DIS, fill=AÑ, DIS, color=AÑ, DIS)) + geom_boxplot() + geom_jitter() + theme(legend.position = "none") + geom_point(color = 'red', fill = 'red', size= 5, shape = 18, alpha = 0.5) + geom_jitter(size = 2, color = 'gray', alpha = 0.8) + geom_boxplot() + theme_bw()Warning: Duplicated aesthetics after name standardisation:
Año. <- ggplot(Propul, aes(x = AÑ, y = ALC, fill=AÑ)) + geom_boxplot() + theme(legend.position = "none")Dist. <- ggplot(Propul, aes(x = DIS, y = ALC, fill=DIS)) + geom_boxplot() + theme(legend.position = "none")cv.model(Ra.aov)SULFATO.lm <- lm(SULFATES1 ~ DIS+AÑ)anova(SULFATO.lm , test="F")Analysis of Variance Table
Response: SULFATES1
Df Sum Sq Mean Sq F value Pr(>F)
DIS 160 4501101144 28131882 33.0995 <2e-16 ***
AÑ 9 25022723 2780303 3.2713 6e-04 ***
Residuals 1339 1138041877 849919
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Grafica <- interactionMeans(SULFATO.lm)
plot(Grafica)
par(mfrow=c(1,1))shapiro.test(SULFATO.lm$res)
Shapiro-Wilk normality test
data: SULFATO.lm$res
W = 0.1233, p-value < 2.2e-16qqPlot(Ra.aov)Error in model.frame.default(formula = Y ~ X - 1, drop.unused.levels = TRUE) :
variable lengths differ (found for 'X')
TukeyHSD(Ra.aov,conf.level = 0.95) Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = SULFATES1 ~ AÑ + DIS)
$AÑ
diff lwr upr p adj
2013-2012 -39.798233 -871.58050 791.98403 1.0000000
2014-2012 -41.576351 -884.15473 801.00202 1.0000000
2015-2012 251.537338 -615.92840 1119.00308 0.9958848
2016-2012 -17.794318 -846.89326 811.30463 1.0000000
2017-2012 -85.365705 -943.05118 772.31977 0.9999995
2018-2012 500.808607 -350.91109 1352.52830 0.6939362
2021-2012 2.727528 -841.95871 847.41377 1.0000000
2022-2012 86.569360 -738.85867 911.99739 0.9999992
2027-2012 -94.361538 -3126.23636 2937.51328 1.0000000
2014-2013 -1.778118 -299.46814 295.91190 1.0000000
2015-2013 291.335571 -70.84924 653.52038 0.2434406
2016-2013 22.003914 -235.06824 279.07607 0.9999999
2017-2013 -45.567472 -383.65830 292.52335 0.9999926
2018-2013 540.606839 217.94995 863.26373 0.0000057
2021-2013 42.525761 -261.07905 346.13057 0.9999896
2022-2013 126.367593 -118.60683 371.34202 0.8309869
2027-2013 -54.563306 -2982.18110 2873.05449 1.0000000
2015-2014 293.113689 -93.22091 679.44829 0.3238491
2016-2014 23.782033 -266.32601 313.89007 0.9999999
2017-2014 -43.789354 -407.63271 320.05400 0.9999972
2018-2014 542.384958 192.83660 891.93331 0.0000432
2021-2014 44.303879 -287.73898 376.34674 0.9999932
2022-2014 128.145711 -151.29847 407.58989 0.9098744
2027-2014 -52.785187 -2983.48861 2877.91823 1.0000000
2016-2015 -269.331656 -625.31082 86.64751 0.3278798
2017-2015 -336.903043 -755.16191 81.35582 0.2416647
2018-2015 249.271269 -156.61362 655.15616 0.6371160
2021-2015 -248.809810 -639.72023 142.10061 0.5870298
2022-2015 -164.967978 -512.31153 182.37557 0.8903509
2027-2015 -345.898876 -3283.85412 2592.05637 0.9999977
2017-2016 -67.571387 -399.00574 263.86297 0.9997496
2018-2016 518.602925 202.92778 834.27807 0.0000099
2021-2016 20.521846 -275.65244 316.69613 1.0000000
2022-2016 104.363678 -131.33909 340.06645 0.9264248
2027-2016 -76.567220 -3003.42377 2850.28933 1.0000000
2018-2017 586.174312 201.63549 970.71314 0.0000661
2021-2017 88.093233 -280.60518 456.79165 0.9990892
2022-2017 171.935065 -150.20638 494.07651 0.8004896
2027-2017 -8.995833 -2944.07822 2926.08655 1.0000000
2021-2018 -498.081078 -852.68026 -143.48189 0.0003924
2022-2018 -414.239247 -720.14311 -108.33539 0.0007952
2027-2018 -595.170145 -3528.51476 2338.17447 0.9997594
2022-2021 83.841832 -201.89508 369.57874 0.9955004
2027-2021 -97.089067 -3028.39919 2834.22106 1.0000000
2027-2022 -180.930899 -3106.74970 2744.88790 1.0000000
$DIS
diff lwr
Acora-Achaya 5.061555e+02 -3773.14029
Ajoyani-Achaya 4.407408e+02 -3861.01869
Alto Selva Alegre-Achaya 5.081368e+02 -3909.68323
Amantani-Achaya 7.026052e+02 -4398.65405
Ananea-Achaya 4.736235e+02 -3767.96851
Antauta-Achaya 3.648637e+02 -4052.95640
Arapa-Achaya 1.906446e+02 -4910.61459
Asilla-Achaya 5.111150e+02 -5379.31181
Asillo-Achaya 5.519794e+02 -3838.48549
Atuncolla-Achaya 3.159005e+02 -4136.84356
Atuncollla-Achaya -3.247000e+01 -4841.98331
Ayapata-Achaya 5.161368e+02 -4585.12238
Ayaviri-Achaya 4.175213e+02 -3767.61630
Azangaro-Achaya 4.772825e+02 -3723.63131
Azángaro-Achaya 5.502751e+02 -3761.07584
Barrio Chanu Chanu-Achaya 6.861329e+02 -5204.29384
Cabana-Achaya 4.821191e+02 -4016.76870
Cabanilla-Achaya 5.094695e+02 -3989.41823
Cabanillas-Achaya 4.296189e+02 -3988.20121
Campamento lechemayo-Achaya 4.910379e+02 -4610.22129
Candarave-Achaya 4.932329e+02 -5397.19384
Capachica-Achaya 3.567242e+02 -4300.06699
Capazo-Achaya 4.947996e+02 -4314.71372
Carabaya-Achaya 2.106767e+02 -4017.12152
Caracoto-Achaya 5.602167e+02 -3790.15179
Cavanillas-Achaya 6.161368e+02 -5274.28993
Cayarani-Achaya 4.153074e+02 -3843.46409
Challhuahuacho-Achaya 4.939329e+02 -5396.49384
Chucuito-Achaya 5.185609e+02 -3717.80430
Chupa-Achaya 4.723386e+02 -5418.08816
Coasa-Achaya 3.781312e+02 -4278.66002
Coata-Achaya 5.678343e+02 -3754.55726
Cochauma-Achaya 5.123550e+02 -5378.07181
Cojata-Achaya 5.353612e+02 -3917.38291
Collao-Achaya 2.228013e+02 -4433.98997
Comunidad San Jose de Calala-Achaya 4.908329e+02 -5399.59384
Condensuyo-Achaya -4.700000e-01 -5890.89676
Conduriri-Achaya 6.021368e+02 -4207.37647
Conima-Achaya 5.044689e+02 -4305.04445
Cooperativa Minera "Los Andes"-Achaya 6.071329e+02 -4494.12629
Copani-Achaya 5.060239e+02 -4595.23528
Copani del Rosario-Achaya 5.289329e+02 -5361.49384
Cupi-Achaya 6.552499e+02 -4446.00933
Cusipata-Achaya 8.061329e+02 -5084.29384
Cutucuyo-Achaya 4.901329e+02 -5400.29384
Cuyo Cuyo-Achaya 6.008072e+02 -4208.70613
Cuyocuyo-Achaya 5.026729e+02 -5387.75384
Cuyuraya-Achaya 4.847412e+02 -3933.07889
Desaguadero-Achaya 5.773564e+02 -4232.15693
Espinar-Achaya 3.961926e+02 -3972.26479
Huacullani-Achaya 5.043629e+02 -4596.89629
Huaje-Achaya 5.155550e+02 -4585.70426
Huancane-Achaya 4.961127e+02 -3726.50126
Huancané-Achaya 5.428579e+02 -4558.40129
Huata-Achaya 4.656111e+02 -3869.62564
Huatasani-Achaya 5.641577e+02 -4245.35557
Huatta-Achaya 5.203150e+02 -5370.11181
Huaytapata-Achaya 5.269350e+02 -5363.49181
Huaytará-Achaya 5.056329e+02 -4595.62629
Humachiri-Achaya 4.956011e+02 -4605.65814
Humajalso-Achaya 5.561329e+02 -5334.29384
Ilave-Achaya 4.747529e+02 -3751.21630
Inambiri-Achaya 5.365650e+02 -4120.22628
Inchupalla-Achaya 5.313368e+02 -4031.36811
Ituata-Achaya 5.196125e+02 -4043.09249
J.D. Choquehuanca-Achaya 5.549150e+02 -5335.51181
Jallihuaya-Achaya 5.751743e+02 -5315.25245
Jose Domingo Choquehuanca-Achaya 5.081368e+02 -4593.12238
Juli-Achaya 5.772963e+02 -3633.88973
Juliaca-Achaya 1.137057e+03 -3048.67233
Julica-Achaya 5.417192e+02 -5348.70752
Kellullo-Iayala-Wilacollo-Achaya -3.247000e+01 -5922.89676
Kelluyo-Achaya 5.667172e+02 -3696.46067
Lampa-Achaya 2.738358e+02 -4005.46008
Laraqueri-Achaya 5.121368e+02 -4297.37647
Las lomas-Achaya 1.240300e+02 -5766.39676
Limbani-Achaya 5.191668e+02 -3979.72093
llalli-Achaya 5.537043e+02 -5336.72245
Llalli-Achaya 4.762718e+02 -3941.54824
Llica-Achaya 9.630000e+00 -5091.62922
Macari-Achaya 5.864851e+02 -4070.30611
Macusani-Achaya 5.444552e+02 -3846.00968
Mañaso-Achaya 5.815761e+02 -4519.68314
Mañazo-Achaya 5.008242e+02 -3916.99590
Marcahuasi-Achaya 5.501329e+02 -5340.29384
Mariscal Nieto-Achaya 5.191329e+02 -5371.29384
Melgar-Achaya 1.447176e+02 -4245.74733
Miraflores-Achaya 1.066133e+03 -4824.29384
Moho-Achaya 5.664427e+02 -3768.79406
Muñani-Achaya 5.163629e+02 -4584.89629
Nicasio-Achaya -2.270000e+00 -5892.69676
Nuñoa-Achaya 5.246929e+02 -4576.56629
Occorrachi-Achaya 5.159150e+02 -4293.59836
Ocuviri-Achaya 5.568362e+02 -3626.39468
Ollachea-Achaya 5.486368e+02 -4108.15440
Ollaraya-Achaya 4.881329e+02 -5402.29384
Orurillo-Achaya 9.248695e+02 -3884.64386
pajcha-Achaya 4.547492e+02 -5435.67752
Pajcha-Achaya 4.944257e+02 -3958.31835
Patambuco-Achaya 4.960029e+02 -4160.78832
Paucarcolla-Achaya 5.031182e+02 -4059.58677
Pekosani-Achaya 4.977286e+02 -3955.01549
Phara-Achaya 5.121368e+02 -4297.37647
Pichacani-Achaya 4.659045e+02 -4343.60883
Pilcuyo-Achaya 5.483014e+02 -4014.40350
Pisacoma-Achaya 6.684753e+04 60957.10324
Pisacona-Achaya 4.901329e+02 -5400.29384
Plateria-Achaya 4.528528e+02 -4046.03496
Platería-Achaya 5.541368e+02 -5336.28993
pokopaka-Achaya 4.830092e+02 -5407.41752
Pokopaka-Achaya 5.005587e+02 -3917.26139
Pomata-Achaya 5.922310e+02 -3906.65677
Potoni-Achaya 6.001368e+02 -5290.28993
Pucará-Achaya 5.121368e+02 -5378.28993
Puno-Achaya 4.483070e+02 -3726.52885
Pusi-Achaya 5.745868e+02 -3760.64994
Púsi-Achaya 7.219150e+02 -3840.78999
Putina-Achaya 4.128843e+02 -3977.58055
Quilca Punco-Achaya 5.081368e+02 -4148.65440
Quilcapunco-Achaya 5.185150e+02 -4582.74426
Rio Asillo-Achaya 5.636329e+02 -5326.79384
Rio Lechemayo grande-Achaya 5.012429e+02 -5389.18384
Rosasapata-Achaya 5.119050e+02 -5378.52181
Saman-Achaya 1.127385e+03 -3973.87459
San Anton-Achaya 5.867910e+02 -3912.09678
San Antonio-Achaya 5.680499e+02 -4088.74130
San Antonio de Esquilachi-Achaya 4.023592e+02 -4254.43199
San Gaban-Achaya 5.240049e+02 -3674.61148
San Gabán-Achaya 5.161368e+02 -5374.28993
San Pablo-Achaya 4.277542e+02 -3990.06583
San Pedro-Achaya 7.049150e+02 -4396.34426
San Pedro de Putina Punco-Achaya 4.723386e+02 -5418.08816
Sandia-Achaya 5.119150e+02 -5378.51181
Santa Lucia-Achaya 4.660261e+02 -4190.76516
Santa rosa-Achaya 2.638013e+02 -5626.62550
Santa Rosa-Achaya 5.180489e+02 -3804.34263
Santa Rosa-Mazo Cruz-Achaya 5.281368e+02 -5362.28993
Santiago B-Achaya -3.247000e+01 -5922.89676
Santiago de Pupuja-Achaya 5.232418e+02 -4578.01738
Sicuani-Achaya 4.881329e+02 -4613.12629
Sina-Achaya 4.921329e+02 -5398.29384
Taraco-Achaya 2.289434e+02 -4005.06764
Tarata-Achaya 5.883404e+02 -3746.89629
Thunco-Achaya 6.421329e+02 -5248.29384
Tinicachi-Achaya 5.141368e+02 -4587.12238
Tiquilla-Achaya 5.315329e+02 -5358.89384
Tiquillaca-Achaya 5.183071e+02 -4291.20618
Torata-Achaya -3.247000e+01 -4689.26124
Torococha-Achaya 4.656111e+02 -5424.81569
Umachiri-Achaya 5.338925e+02 -4028.81249
URB. Chanu Chanu-Achaya 5.181329e+02 -5372.29384
Ururillo-Achaya 5.225368e+02 -4578.72238
Usicayo-Achaya 3.817692e+02 -4180.93571
Usicayos-Achaya 5.358195e+02 -3916.92461
Vilque-Achaya 1.417943e+02 -4667.71901
Vilquechico-Achaya 6.383639e+02 -4462.89528
Yanacocha-Achaya 5.055911e+02 -4595.66814
Yanahuaya-Achaya 5.191150e+02 -5371.31181
Yunguyo-Achaya 5.258184e+02 -3926.92569
Zepita-Achaya 6.296814e+02 -3705.55530
Ajoyani-Acora -6.541475e+01 -1521.56722
Alto Selva Alegre-Acora 1.981301e+00 -1767.87168
Amantani-Acora 1.964496e+02 -2908.07787
Ananea-Acora -3.253200e+01 -1299.95004
Antauta-Acora -1.412919e+02 -1911.14484
Arapa-Acora -3.155109e+02 -3420.03841
Asilla-Acora 4.959419e+00 -4274.33641
Asillo-Acora 4.582386e+01 -1654.59588
Atuncolla-Acora -1.902550e+02 -2045.56505
Atuncollla-Acora -5.386255e+02 -3136.05959
Ayapata-Acora 9.981301e+00 -3094.54619
Ayaviri-Acora -8.863422e+01 -1151.94142
Azangaro-Acora -2.887303e+01 -1152.67118
Azángaro-Acora 4.411955e+01 -1440.12848
Barrio Chanu Chanu-Acora 1.799774e+02 -4099.31845
Cabana-Acora -2.403648e+01 -1987.51207
Cabanilla-Acora 3.313981e+00 -1960.16161
Cabanillas-Acora -7.653667e+01 -1846.38965
Campamento lechemayo-Acora -1.511761e+01 -3119.64511
Candarave-Acora -1.292261e+01 -4292.21845
Capachica-Acora -1.494313e+02 -2451.81050
Capazo-Acora -1.135595e+01 -2608.79000
Carabaya-Acora -2.954788e+02 -1515.93906
Caracoto-Acora 5.406115e+01 -1539.97500
Cavanillas-Acora 1.099813e+02 -4169.31453
Cayarani-Acora -9.084811e+01 -1414.62308
Challhuahuacho-Acora -1.222261e+01 -4291.51845
Chucuito-Acora 1.240535e+01 -1237.40881
Chupa-Acora -3.381693e+01 -4313.11277
Coasa-Acora -1.280243e+02 -2430.40353
Coata-Acora 6.167875e+01 -1454.34038
Cochauma-Acora 6.199419e+00 -4273.09641
Cojata-Acora 2.920565e+01 -1826.10439
Collao-Acora -2.833543e+02 -2585.73348
Comunidad San Jose de Calala-Acora -1.532261e+01 -4294.61845
Condensuyo-Acora -5.066255e+02 -4785.92137
Conduriri-Acora 9.598130e+01 -2501.45276
Conima-Acora -1.686671e+00 -2599.12073
Cooperativa Minera "Los Andes"-Acora 1.009774e+02 -3003.55011
Copani-Acora -1.315975e-01 -3104.65909
Copani del Rosario-Acora 2.277739e+01 -4256.51845
Cupi-Acora 1.490943e+02 -2955.43315
Cusipata-Acora 2.999774e+02 -3979.31845
Cutucuyo-Acora -1.602261e+01 -4295.31845
Cuyo Cuyo-Acora 9.465164e+01 -2502.78241
Cuyocuyo-Acora -3.482614e+00 -4282.77845
Cuyuraya-Acora -2.141436e+01 -1791.26733
Desaguadero-Acora 7.120085e+01 -2526.23321
Espinar-Acora -1.099629e+02 -1752.72446
Huacullani-Acora -1.792614e+00 -3106.32011
Huaje-Acora 9.399419e+00 -3095.12808
Huancane-Acora -1.004289e+01 -1212.42122
Huancané-Acora 3.670239e+01 -3067.82511
Huata-Acora -4.054446e+01 -1592.80821
Huatasani-Acora 5.800221e+01 -2539.43185
Huatta-Acora 1.415942e+01 -4265.13641
Huaytapata-Acora 2.077942e+01 -4258.51641
Huaytará-Acora -5.226139e-01 -3105.05011
Humachiri-Acora -1.055446e+01 -3115.08195
Humajalso-Acora 4.997739e+01 -4229.31845
Ilave-Acora -3.140264e+01 -1245.51177
Inambiri-Acora 3.040942e+01 -2271.96979
Inchupalla-Acora 2.518130e+01 -2080.41169
Ituata-Acora 1.345692e+01 -2092.13607
J.D. Choquehuanca-Acora 4.875942e+01 -4230.53641
Jallihuaya-Acora 6.901877e+01 -4210.27706
Jose Domingo Choquehuanca-Acora 1.981301e+00 -3102.54619
Juli-Acora 7.114072e+01 -1090.46711
Juliaca-Acora 6.309017e+02 -434.73292
Julica-Acora 3.556371e+01 -4243.73213
Kellullo-Iayala-Wilacollo-Acora -5.386255e+02 -4817.92137
Kelluyo-Acora 6.056169e+01 -1277.32140
Lampa-Acora -2.323198e+02 -1620.70669
Laraqueri-Acora 5.981301e+00 -2591.45276
Las lomas-Acora -3.821255e+02 -4661.42137
Limbani-Acora 1.301128e+01 -1950.46431
llalli-Acora 4.754877e+01 -4231.74706
Llalli-Acora -2.988371e+01 -1799.73669
Llica-Acora -4.965255e+02 -3601.05303
Macari-Acora 8.032959e+01 -2222.04962
Macusani-Acora 3.829967e+01 -1662.12007
Mañaso-Acora 7.542054e+01 -3029.10695
Mañazo-Acora -5.331370e+00 -1775.18435
Marcahuasi-Acora 4.397739e+01 -4235.31845
Mariscal Nieto-Acora 1.297739e+01 -4266.31845
Melgar-Acora -3.614380e+02 -2061.85772
Miraflores-Acora 5.599774e+02 -3719.31845
Moho-Acora 6.028712e+01 -1491.97663
Muñani-Acora 1.020739e+01 -3094.32011
Nicasio-Acora -5.084255e+02 -4787.72137
upr p adj
Acora-Achaya 4785.45137 1.0000000
Ajoyani-Achaya 4742.50027 1.0000000
Alto Selva Alegre-Achaya 4925.95691 1.0000000
Amantani-Achaya 5803.86438 1.0000000
Ananea-Achaya 4715.21559 1.0000000
Antauta-Achaya 4782.68375 1.0000000
Arapa-Achaya 5291.90384 1.0000000
Asilla-Achaya 6401.54172 1.0000000
Asillo-Achaya 4942.44429 1.0000000
Atuncolla-Achaya 4768.64463 1.0000000
Atuncollla-Achaya 4777.04331 1.0000000
Ayapata-Achaya 5617.39606 1.0000000
Ayaviri-Achaya 4602.65893 1.0000000
Azangaro-Achaya 4678.19632 1.0000000
Azángaro-Achaya 4861.62602 1.0000000
Barrio Chanu Chanu-Achaya 6576.55969 1.0000000
Cabana-Achaya 4981.00681 1.0000000
Cabanilla-Achaya 5008.35727 1.0000000
Cabanillas-Achaya 4847.43894 1.0000000
Campamento lechemayo-Achaya 5592.29714 1.0000000
Candarave-Achaya 6383.65969 1.0000000
Capachica-Achaya 5013.51549 1.0000000
Capazo-Achaya 5304.31290 1.0000000
Carabaya-Achaya 4438.47492 1.0000000
Caracoto-Achaya 4910.58517 1.0000000
Cavanillas-Achaya 6506.56360 1.0000000
Cayarani-Achaya 4674.07894 1.0000000
Challhuahuacho-Achaya 6384.35969 1.0000000
Chucuito-Achaya 4754.92609 1.0000000
Chupa-Achaya 6362.76537 1.0000000
Coasa-Achaya 5034.92246 1.0000000
Coata-Achaya 4890.22583 1.0000000
Cochauma-Achaya 6402.78172 1.0000000
Cojata-Achaya 4988.10528 1.0000000
Collao-Achaya 4879.59251 1.0000000
Comunidad San Jose de Calala-Achaya 6381.25969 1.0000000
Condensuyo-Achaya 5889.95676 1.0000000
Conduriri-Achaya 5411.65015 1.0000000
Conima-Achaya 5313.98218 1.0000000
Cooperativa Minera "Los Andes"-Achaya 5708.39214 1.0000000
Copani-Achaya 5607.28316 1.0000000
Copani del Rosario-Achaya 6419.35969 1.0000000
Cupi-Achaya 5756.50910 1.0000000
Cusipata-Achaya 6696.55969 1.0000000
Cutucuyo-Achaya 6380.55969 1.0000000
Cuyo Cuyo-Achaya 5410.32050 1.0000000
Cuyocuyo-Achaya 6393.09969 1.0000000
Cuyuraya-Achaya 4902.56126 1.0000000
Desaguadero-Achaya 5386.86970 1.0000000
Espinar-Achaya 4764.65002 1.0000000
Huacullani-Achaya 5605.62214 1.0000000
Huaje-Achaya 5616.81417 1.0000000
Huancane-Achaya 4718.72657 1.0000000
Huancané-Achaya 5644.11714 1.0000000
Huata-Achaya 4800.84779 1.0000000
Huatasani-Achaya 5373.67106 1.0000000
Huatta-Achaya 6410.74172 1.0000000
Huaytapata-Achaya 6417.36172 1.0000000
Huaytará-Achaya 5606.89214 1.0000000
Humachiri-Achaya 5596.86030 1.0000000
Humajalso-Achaya 6446.55969 1.0000000
Ilave-Achaya 4700.72209 1.0000000
Inambiri-Achaya 5193.35620 1.0000000
Inchupalla-Achaya 5094.04179 1.0000000
Ituata-Achaya 5082.31742 1.0000000
J.D. Choquehuanca-Achaya 6445.34172 1.0000000
Jallihuaya-Achaya 6465.60108 1.0000000
Jose Domingo Choquehuanca-Achaya 5609.39606 1.0000000
Juli-Achaya 4788.48224 1.0000000
Juliaca-Achaya 5322.78672 1.0000000
Julica-Achaya 6432.14601 1.0000000
Kellullo-Iayala-Wilacollo-Achaya 5857.95676 1.0000000
Kelluyo-Achaya 4829.89513 1.0000000
Lampa-Achaya 4553.13159 1.0000000
Laraqueri-Achaya 5321.65015 1.0000000
Las lomas-Achaya 6014.45676 1.0000000
Limbani-Achaya 5018.05457 1.0000000
llalli-Achaya 6444.13108 1.0000000
Llalli-Achaya 4894.09190 1.0000000
Llica-Achaya 5110.88922 1.0000000
Macari-Achaya 5243.27637 1.0000000
Macusani-Achaya 4934.92010 1.0000000
Mañaso-Achaya 5682.83530 1.0000000
Mañazo-Achaya 4918.64424 1.0000000
Marcahuasi-Achaya 6440.55969 1.0000000
Mariscal Nieto-Achaya 6409.55969 1.0000000
Melgar-Achaya 4535.18244 1.0000000
Miraflores-Achaya 6956.55969 1.0000000
Moho-Achaya 4901.67937 1.0000000
Muñani-Achaya 5617.62214 1.0000000
Nicasio-Achaya 5888.15676 1.0000000
Nuñoa-Achaya 5625.95214 1.0000000
Occorrachi-Achaya 5325.42827 1.0000000
Ocuviri-Achaya 4740.06713 1.0000000
Ollachea-Achaya 5205.42808 1.0000000
Ollaraya-Achaya 6378.55969 1.0000000
Orurillo-Achaya 5734.38277 1.0000000
pajcha-Achaya 6345.17601 1.0000000
Pajcha-Achaya 4947.16984 1.0000000
Patambuco-Achaya 5152.79417 1.0000000
Paucarcolla-Achaya 5065.82314 1.0000000
Pekosani-Achaya 4950.47270 1.0000000
Phara-Achaya 5321.65015 1.0000000
Pichacani-Achaya 5275.41780 1.0000000
Pilcuyo-Achaya 5111.00640 1.0000000
Pisacoma-Achaya 72737.95676 0.0000000
Pisacona-Achaya 6380.55969 1.0000000
Plateria-Achaya 4951.74055 1.0000000
Platería-Achaya 6444.56360 1.0000000
pokopaka-Achaya 6373.43601 1.0000000
Pokopaka-Achaya 4918.37876 1.0000000
Pomata-Achaya 5091.11874 1.0000000
Potoni-Achaya 6490.56360 1.0000000
Pucará-Achaya 6402.56360 1.0000000
Puno-Achaya 4623.14293 1.0000000
Pusi-Achaya 4909.82349 1.0000000
Púsi-Achaya 5284.61991 1.0000000
Putina-Achaya 4803.34923 1.0000000
Quilca Punco-Achaya 5164.92808 1.0000000
Quilcapunco-Achaya 5619.77417 1.0000000
Rio Asillo-Achaya 6454.05969 1.0000000
Rio Lechemayo grande-Achaya 6391.66969 1.0000000
Rosasapata-Achaya 6402.33172 1.0000000
Saman-Achaya 6228.64384 1.0000000
San Anton-Achaya 5085.67873 1.0000000
San Antonio-Achaya 5224.84119 1.0000000
San Antonio de Esquilachi-Achaya 5059.15049 1.0000000
San Gaban-Achaya 4722.62122 1.0000000
San Gabán-Achaya 6406.56360 1.0000000
San Pablo-Achaya 4845.57432 1.0000000
San Pedro-Achaya 5806.17417 1.0000000
San Pedro de Putina Punco-Achaya 6362.76537 1.0000000
Sandia-Achaya 6402.34172 1.0000000
Santa Lucia-Achaya 5122.81732 1.0000000
Santa rosa-Achaya 6154.22803 1.0000000
Santa Rosa-Achaya 4840.44045 1.0000000
Santa Rosa-Mazo Cruz-Achaya 6418.56360 1.0000000
Santiago B-Achaya 5857.95676 1.0000000
Santiago de Pupuja-Achaya 5624.50106 1.0000000
Sicuani-Achaya 5589.39214 1.0000000
Sina-Achaya 6382.55969 1.0000000
Taraco-Achaya 4462.95437 1.0000000
Tarata-Achaya 4923.57714 1.0000000
Thunco-Achaya 6532.55969 1.0000000
Tinicachi-Achaya 5615.39606 1.0000000
Tiquilla-Achaya 6421.95969 1.0000000
Tiquillaca-Achaya 5327.82045 1.0000000
Torata-Achaya 4624.32124 1.0000000
Torococha-Achaya 6356.03784 1.0000000
Umachiri-Achaya 5096.59742 1.0000000
URB. Chanu Chanu-Achaya 6408.55969 1.0000000
Ururillo-Achaya 5623.79606 1.0000000
Usicayo-Achaya 4944.47420 1.0000000
Usicayos-Achaya 4988.56358 1.0000000
Vilque-Achaya 4951.30762 1.0000000
Vilquechico-Achaya 5739.62316 1.0000000
Yanacocha-Achaya 5606.85030 1.0000000
Yanahuaya-Achaya 6409.54172 1.0000000
Yunguyo-Achaya 4978.56250 1.0000000
Zepita-Achaya 4964.91813 1.0000000
Ajoyani-Acora 1390.73772 1.0000000
Alto Selva Alegre-Acora 1771.83428 1.0000000
Amantani-Acora 3300.97712 1.0000000
Ananea-Acora 1234.88604 1.0000000
Antauta-Acora 1628.56111 1.0000000
Arapa-Acora 2789.01658 1.0000000
Asilla-Acora 4284.25525 1.0000000
Asillo-Acora 1746.24360 1.0000000
Atuncolla-Acora 1665.05503 1.0000000
Atuncollla-Acora 2058.80852 1.0000000
Ayapata-Acora 3114.50879 1.0000000
Ayaviri-Acora 974.67298 1.0000000
Azangaro-Acora 1094.92511 1.0000000
Azángaro-Acora 1528.36758 1.0000000
Barrio Chanu Chanu-Acora 4459.27322 1.0000000
Cabana-Acora 1939.43911 1.0000000
Cabanilla-Acora 1966.78957 1.0000000
Cabanillas-Acora 1693.31630 1.0000000
Campamento lechemayo-Acora 3089.40988 1.0000000
Candarave-Acora 4266.37322 1.0000000
Capachica-Acora 2152.94792 1.0000000
Capazo-Acora 2586.07811 1.0000000
Carabaya-Acora 924.98138 1.0000000
Caracoto-Acora 1648.09730 1.0000000
Cavanillas-Acora 4389.27713 1.0000000
Cayarani-Acora 1232.92686 1.0000000
Challhuahuacho-Acora 4267.07322 1.0000000
Chucuito-Acora 1262.21951 1.0000000
Chupa-Acora 4245.47890 1.0000000
Coasa-Acora 2174.35489 1.0000000
Coata-Acora 1577.69787 1.0000000
Cochauma-Acora 4285.49525 1.0000000
Cojata-Acora 1884.51569 1.0000000
Collao-Acora 2019.02494 1.0000000
Comunidad San Jose de Calala-Acora 4263.97322 1.0000000
Condensuyo-Acora 3772.67029 1.0000000
Conduriri-Acora 2693.41536 1.0000000
Conima-Acora 2595.74738 1.0000000
Cooperativa Minera "Los Andes"-Acora 3205.50488 1.0000000
Copani-Acora 3104.39590 1.0000000
Copani del Rosario-Acora 4302.07322 1.0000000
Cupi-Acora 3253.62184 1.0000000
Cusipata-Acora 4579.27322 1.0000000
Cutucuyo-Acora 4263.27322 1.0000000
Cuyo Cuyo-Acora 2692.08570 1.0000000
Cuyocuyo-Acora 4275.81322 1.0000000
Cuyuraya-Acora 1748.43862 1.0000000
Desaguadero-Acora 2668.63490 1.0000000
Espinar-Acora 1532.79862 1.0000000
Huacullani-Acora 3102.73488 1.0000000
Huaje-Acora 3113.92691 1.0000000
Huancane-Acora 1192.33544 1.0000000
Huancané-Acora 3141.22988 1.0000000
Huata-Acora 1511.71929 1.0000000
Huatasani-Acora 2655.43626 1.0000000
Huatta-Acora 4293.45525 1.0000000
Huaytapata-Acora 4300.07525 1.0000000
Huaytará-Acora 3104.00488 1.0000000
Humachiri-Acora 3093.97303 1.0000000
Humajalso-Acora 4329.27322 1.0000000
Ilave-Acora 1182.70648 1.0000000
Inambiri-Acora 2332.78863 1.0000000
Inchupalla-Acora 2130.77429 1.0000000
Ituata-Acora 2119.04992 1.0000000
J.D. Choquehuanca-Acora 4328.05525 1.0000000
Jallihuaya-Acora 4348.31461 1.0000000
Jose Domingo Choquehuanca-Acora 3106.50879 1.0000000
Juli-Acora 1232.74854 1.0000000
Juliaca-Acora 1696.53623 0.9996508
Julica-Acora 4314.85954 1.0000000
Kellullo-Iayala-Wilacollo-Acora 3740.67029 1.0000000
Kelluyo-Acora 1398.44479 1.0000000
Lampa-Acora 1156.06712 1.0000000
Laraqueri-Acora 2603.41536 1.0000000
Las lomas-Acora 3897.17029 1.0000000
Limbani-Acora 1976.48687 1.0000000
llalli-Acora 4326.84461 1.0000000
Llalli-Acora 1739.96927 1.0000000
Llica-Acora 2608.00196 1.0000000
Macari-Acora 2382.70880 1.0000000
Macusani-Acora 1738.71941 1.0000000
Mañaso-Acora 3179.94803 1.0000000
Mañazo-Acora 1764.52161 1.0000000
Marcahuasi-Acora 4323.27322 1.0000000
Mariscal Nieto-Acora 4292.27322 1.0000000
Melgar-Acora 1338.98176 1.0000000
Miraflores-Acora 4839.27322 1.0000000
Moho-Acora 1612.55087 1.0000000
Muñani-Acora 3114.73488 1.0000000
Nicasio-Acora 3770.87029 1.0000000
[ reached getOption("max.print") -- omitted 12630 rows ]plot(TukeyHSD(Ra.aov))
LSD.test(y = Ra.aov, trt = "SULF" ,group = T, console = T)
Study: Ra.aov ~ "SULF"
LSD t Test for SULFATES1
Mean Square Error: 849919.3
SULFATES1, means and individual ( 95 %) CI
Alpha: 0.05 ; DF Error: 1339
Critical Value of t: 1.961737
Groups according to probability of means differences and alpha level( 0.05 )
Treatments with the same letter are not significantly different.LSD.test(y = Ra.aov, trt = "AÑ" ,group = T, console = T)
Study: Ra.aov ~ "AÑ"
LSD t Test for SULFATES1
Mean Square Error: 849919.3
AÑ, means and individual ( 95 %) CI
Alpha: 0.05 ; DF Error: 1339
Critical Value of t: 1.961737
Groups according to probability of means differences and alpha level( 0.05 )
Treatments with the same letter are not significantly different.LSD.test(Ra.aov,"DIS",console=TRUE)
Study: Ra.aov ~ "DIS"
LSD t Test for SULFATES1
Mean Square Error: 849919.3
DIS, means and individual ( 95 %) CI
Alpha: 0.05 ; DF Error: 1339
Critical Value of t: 1.961737
Groups according to probability of means differences and alpha level( 0.05 )
Treatments with the same letter are not significantly different.HSD.test(Ra.aov, "DIS",console=TRUE)
Study: Ra.aov ~ "DIS"
HSD Test for SULFATES1
Mean Square Error: 849919.3
DIS, means
Alpha: 0.05 ; DF Error: 1339
Critical Value of Studentized Range: 6.389368
Groups according to probability of means differences and alpha level( 0.05 )
Treatments with the same letter are not significantly different.SNK.test(Ra.aov, "DIS",console=TRUE)
Study: Ra.aov ~ "DIS"
Student Newman Keuls Test
for SULFATES1
Mean Square Error: 849919.3
DIS, means
Groups according to probability of means differences and alpha level( 0.05 )
Means with the same letter are not significantly different.scheffe.test(Ra.aov, "DIS",console=TRUE)
Study: Ra.aov ~ "DIS"
Scheffe Test for SULFATES1
Mean Square Error : 849919.3
DIS, means
Alpha: 0.05 ; DF Error: 1339
Critical Value of F: 1.204911
Groups according to probability of means differences and alpha level( 0.05 )
Means with the same letter are not significantly different.duncan.test(Ra.aov, "DIS",console=TRUE)
Study: Ra.aov ~ "DIS"
Duncan's new multiple range test
for SULFATES1
Mean Square Error: 849919.3
DIS, means
Groups according to probability of means differences and alpha level( 0.05 )
Means with the same letter are not significantly different.shapiro.test(Ra.aov$res)
Shapiro-Wilk normality test
data: Ra.aov$res
W = 0.1233, p-value < 2.2e-16summary(Ra.aov$residuals) Min. 1st Qu. Median Mean 3rd Qu. Max.
-1429.744 -65.810 -7.045 0.000 28.363 27135.507 LS0tDQp0aXRsZTogIkdSVVBPIFNVTEZBVE9TOiBMSVNCRVRIIEJFTElaQVJJTywgSlVMSVRaQSBDVVRJUEEsIFJVU1NFTEwgUVVJU1BFIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KYGBge3J9DQpsaWJyYXJ5KGdncGxvdDIpDQpgYGANCg0KYGBge3J9DQpnZXR3ZCgpDQpgYGANCg0KYGBge3J9DQpmaWxlLmNob29zZSgpDQpgYGANCg0KYGBge3J9DQpydXRhX3Byb3B1bHNvcmEgPC0gIkM6XFxVc2Vyc1xcTGlzYmV0aFxcRG93bmxvYWRzXFxESVNFw5FPIEVYUEVSSU1FTlRBTCAoVFJBQkFKT1MpXFxTVUxGQVRPUy54bHN4Ig0KYGBgDQoNCg0KYGBge3J9DQpleGNlbF9zaGVldHMocnV0YV9wcm9wdWxzb3JhKQ0KYGBgDQoNCmBgYHtyfQ0KUHJvcHVsPC1yZWFkX2V4Y2VsKHJ1dGFfcHJvcHVsc29yYSkNCmBgYA0KDQpgYGB7cn0NCnByaW50KGhlYWQoUHJvcHVsKSkNCmBgYA0KDQpgYGB7cn0NClZpZXcoUHJvcHVsKQ0KYGBgDQoNCmBgYHtyfQ0KYXR0YWNoKFByb3B1bCkNCmBgYA0KDQpgYGB7cn0NCm5hbWVzKFByb3B1bCkNCmBgYA0KDQpgYGB7cn0NCnN1bW1hcnkoUHJvcHVsKQ0KYGBgDQoNCmBgYHtyfQ0Kc3RyKFByb3B1bCkNCmBgYA0KDQpgYGB7cn0NCnBsb3QoYHN1bGZhdG9zIChtZy9MKWApDQpgYGANCg0KYGBge3J9DQpwbG90KFNVTEZBVE9TKQ0KYGBgDQoNCmBgYHtyfQ0KQcORIDwtIGZhY3RvcihQcm9wdWwkYEHDsW9gKQ0KYGBgDQoNCmBgYHtyfQ0KRElTIDwtYXMudmVjdG9yKFByb3B1bCRgRGlzdHJpdG9gKQ0KYGBgDQoNCmBgYHtyfQ0KU1VMRkFURVMgPC1hcy52ZWN0b3IoUHJvcHVsJGBzdWxmYXRvcyAobWcvTClgKQ0KYGBgDQoNCg0KYGBge3J9DQpTVUxGQVRFUzEgPi1hcy5udW1lcmljKFNVTEZBVEVTKQ0KYGBgDQoNCmBgYHtyfQ0KcGFyKG1mcm93PWMoMSwxKSkNCmBgYA0KDQpgYGB7cn0NCmJveHBsb3Qoc3BsaXQoU1VMRkFURVMxLEHDkSkseGxhYj0iQcOxbyIsIHlsYWI9IlN1bGZhdG9zIikNCmBgYA0KDQpgYGB7cn0NCmJveHBsb3Qoc3BsaXQoU1VMRkFURVMxLERJUykseGxhYj0iRGlzdHJpdG8iLCB5bGFiPSJTdWxmYXRvcyIpDQpgYGANCg0KYGBge3J9DQpSYS5hb3Y8LWFvdihTVUxGQVRFUzEgfiBBw5ErRElTKQ0KYGBgDQoNCmBgYHtyfQ0KYW5vdmEoUmEuYW92KQ0KYGBgDQoNCmBgYHtyfQ0KZ2dwbG90KFJhLmFvdiwgYWVzKHggPSBBw5EsIERJUywgZmlsbD1Bw5EsIERJUywgY29sb3I9QcORLCBESVMpKSArIGdlb21fYm94cGxvdCgpICsgZ2VvbV9qaXR0ZXIoKSArIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgKyBnZW9tX3BvaW50KGNvbG9yID0gJ3JlZCcsIGZpbGwgPSAncmVkJywgc2l6ZT0gNSwgc2hhcGUgPSAxOCwgYWxwaGEgPSAwLjUpICsgZ2VvbV9qaXR0ZXIoc2l6ZSA9IDIsIGNvbG9yID0gJ2dyYXknLCBhbHBoYSA9IDAuOCkgKyBnZW9tX2JveHBsb3QoKSArIHRoZW1lX2J3KCkNCmBgYA0KDQpgYGB7cn0NCkHDsW8uIDwtIGdncGxvdChQcm9wdWwsIGFlcyh4ID0gQcORLCB5ID0gQUxDLCBmaWxsPUHDkSkpICsgZ2VvbV9ib3hwbG90KCkgKyB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibm9uZSIpDQpgYGANCg0KYGBge3J9DQpEaXN0LiA8LSBnZ3Bsb3QoUHJvcHVsLCBhZXMoeCA9IERJUywgeSA9IEFMQywgZmlsbD1ESVMpKSArIGdlb21fYm94cGxvdCgpICsgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKQ0KYGBgDQoNCg0KYGBge3J9DQpjdi5tb2RlbChSYS5hb3YpDQpgYGANCg0KYGBge3J9DQpTVUxGQVRPLmxtIDwtIGxtKFNVTEZBVEVTMSB+IERJUytBw5EpDQpgYGANCg0KYGBge3J9DQphbm92YShTVUxGQVRPLmxtICwgdGVzdD0iRiIpDQpgYGANCg0KYGBge3J9DQpHcmFmaWNhIDwtIGludGVyYWN0aW9uTWVhbnMoU1VMRkFUTy5sbSkNCnBsb3QoR3JhZmljYSkNCnBhcihtZnJvdz1jKDEsMSkpDQpgYGANCg0KYGBge3J9DQpzaGFwaXJvLnRlc3QoU1VMRkFUTy5sbSRyZXMpDQpgYGANCg0KYGBge3J9DQpxcVBsb3QoUmEuYW92KQ0KYGBgDQoNCmBgYHtyfQ0KVHVrZXlIU0QoUmEuYW92LGNvbmYubGV2ZWwgPSAwLjk1KQ0KYGBgDQoNCmBgYHtyfQ0KcGxvdChUdWtleUhTRChSYS5hb3YpKQ0KYGBgDQoNCmBgYHtyfQ0KTFNELnRlc3QoeSA9IFJhLmFvdiwgdHJ0ID0gIlNVTEYiICxncm91cCA9IFQsIGNvbnNvbGUgPSBUKQ0KYGBgDQoNCmBgYHtyfQ0KTFNELnRlc3QoeSA9IFJhLmFvdiwgdHJ0ID0gIkHDkSIgLGdyb3VwID0gVCwgY29uc29sZSA9IFQpDQpgYGANCg0KYGBge3J9DQpMU0QudGVzdChSYS5hb3YsIkRJUyIsY29uc29sZT1UUlVFKQ0KYGBgDQoNCmBgYHtyfQ0KSFNELnRlc3QoUmEuYW92LCAiRElTIixjb25zb2xlPVRSVUUpDQpgYGANCg0KYGBge3J9DQpTTksudGVzdChSYS5hb3YsICJESVMiLGNvbnNvbGU9VFJVRSkNCmBgYA0KDQpgYGB7cn0NCnNjaGVmZmUudGVzdChSYS5hb3YsICJESVMiLGNvbnNvbGU9VFJVRSkNCmBgYA0KDQpgYGB7cn0NCmR1bmNhbi50ZXN0KFJhLmFvdiwgIkRJUyIsY29uc29sZT1UUlVFKQ0KYGBgDQoNCmBgYHtyfQ0Kc2hhcGlyby50ZXN0KFJhLmFvdiRyZXMpIA0KYGBgDQoNCmBgYHtyfQ0Kc3VtbWFyeShSYS5hb3YkcmVzaWR1YWxzKQ0KYGBgDQoNCg==