Heterogeneity Index
Mohammed Fakher Taher
2024-10-20
Dykstra-Parsons coefficient
Hello everyone,I would like to share with you How I am going to determine the heterogeneity index by the Dykstra-Parsons coefficient.
firstly I am goning to set a working directory:
setwd("C:/Users/dell/Desktop/data")
getwd()## [1] "C:/Users/dell/Desktop/data"Now, calling data
data = read.csv("C:/Users/dell/Desktop/data/data.csv")head(data)## depth caliper ind.deep ind.med gamma phi.N R.deep R.med SP
## 1 5667.0 8.685 618.005 569.781 98.823 0.410 1.618 1.755 -56.587
## 2 5667.5 8.686 497.547 419.494 90.640 0.307 2.010 2.384 -61.916
## 3 5668.0 8.686 384.935 300.155 78.087 0.203 2.598 3.332 -55.861
## 4 5668.5 8.686 278.324 205.224 66.232 0.119 3.593 4.873 -41.860
## 5 5669.0 8.686 183.743 131.155 59.807 0.069 5.442 7.625 -34.934
## 6 5669.5 8.686 109.512 75.633 57.109 0.048 9.131 13.222 -39.769
## density.corr density phi.core k.core Facies phi.core.frac
## 1 -0.033 2.205 33.9000 2442.590 F1 0.339000
## 2 -0.067 2.040 33.4131 3006.989 F1 0.334131
## 3 -0.064 1.888 33.1000 3370.000 F1 0.331000
## 4 -0.053 1.794 34.9000 2270.000 F1 0.349000
## 5 -0.054 1.758 35.0644 2530.758 F1 0.350644
## 6 -0.058 1.759 35.3152 2928.314 F1 0.353152data$k.core## [1] 2442.5901 3006.9888 3370.0000 2270.0000 2530.7581 2928.3137
## [7] 3380.0000 3053.2048 2730.0000 2802.1687 2865.1624 2928.1562
## [13] 2960.0000 2368.7905 2550.4099 2759.4919 3040.0000 2558.1506
## [19] 1931.9360 1291.5699 2155.4353 3000.0000 3066.8655 3110.0000
## [25] 2686.7124 2367.5498 2150.3301 2163.2437 2173.0312 2182.8186
## [31] 2192.6060 2202.3936 2214.3201 2613.3105 2990.0000 2832.4385
## [37] 2707.6108 2563.3101 3090.0657 3630.0000 2531.0762 1975.0400
## [43] 2147.4971 2273.6499 1712.4135 1271.3800 418.6400 1362.4678
## [49] 2760.0000 2602.4993 2335.4399 2495.1360 2648.6399 2463.9275
## [55] 2299.0537 2155.4399 2322.1902 2450.0000 2371.0400 2049.1858
## [61] 1550.0200 1722.9451 2056.2000 2023.9753 1990.0415 1956.1077
## [67] 1907.3199 2196.3706 2830.0000 2697.3630 2519.8000 2274.8657
## [73] 2023.0345 1866.9200 2200.3501 2067.3738 1928.4800 1789.1849
## [79] 1648.4301 1866.0350 2302.9199 2225.8699 2149.6299 2276.3159
## [85] 2421.8501 2427.4851 2436.6628 2445.8406 2455.0183 2464.1960
## [91] 2473.3740 2485.9099 2483.5156 2478.3799 2523.2039 2593.2590
## [97] 2660.0000 2638.9934 2620.0000 2423.1160 2284.4299 2011.7695
## [103] 1852.2700 1467.6908 1306.1700 50.1500 822.9766 1969.0728
## [109] 3250.0000 2448.0698 1437.5200 1641.8264 1804.5400 1040.1243
## [115] 271.4900 185.0957 0.8100 92.4887 273.4974 454.5060
## [121] 635.5146 794.9600 193.4600 126.0047 76.8100 85.0097
## [127] 91.1881 97.3665 103.5448 109.7232 115.9015 122.0799
## [133] 128.2582 134.4366 140.6149 146.7933 152.9717 158.1300
## [139] 116.5507 86.9600 63.8231 39.6700 282.9963 601.2100
## [145] 599.3371 597.4316 595.5261 593.6207 591.7152 589.8097
## [151] 588.2900 522.3776 467.5622 412.7467 357.9313 303.1158
## [157] 248.3004 190.5800 181.7682 168.2400 149.2137 114.3747
## [163] 74.1400 146.0790 231.2015 316.3239 406.2400 308.4382
## [169] 204.7999 140.0700 173.4400 73.2072 0.4200 2.2869
## [175] 3.5809 4.8748 5.8300 35.8600 44.5030 55.5300
## [181] 1599.2813 2570.0000 2294.7261 2089.7009 1884.6757 1679.6505
## [187] 1474.6252 1269.6000 1064.5748 859.5496 654.5244 449.4991
## [193] 244.4739 48.6200 43.9438 39.4679 34.9919 32.5300
## [199] 54.2100 906.6724 2434.4473 4390.0000 4804.8022 5480.0000
## [205] 4900.0000 4786.5928 4561.4106 4336.2285 4111.0464 3885.8645
## [211] 3660.6824 3435.5005 3210.3184 2985.1362 2759.9543 2534.7722
## [217] 2309.5901 2084.4082 1859.2261 1634.0441 1408.8621 1183.6799
## [223] 958.4979 733.3159 508.1339 282.9518 162.9400 1109.5569
## [229] 1481.5601 534.7041 30.6700 70.1081 101.8555 133.6028
## [235] 165.3502 197.0975 228.8449 260.5923 292.3396 324.0869
## [241] 355.8343 387.5817 419.3290 451.0764 482.8237 514.5711
## [247] 546.3184 578.0658 609.8132 641.5605 673.3079 705.0552
## [253] 736.8026 768.5499 800.2973 832.0446 863.7920 895.5394
## [259] 927.2867 959.0341 990.7814 1022.5287 1054.2761 1086.0234
## [265] 1117.7709 1149.5182 1181.2655 1213.0129 1244.7603 1276.5076
## [271] 1308.2550 1340.0023 1371.7496 1403.4971 1435.2444 1466.9917
## [277] 1498.7391 1530.4865 1562.2338 1584.9600 1066.3792 662.5491
## [283] 386.1200 1821.4403 2468.4099 2346.7214 2285.0200 4135.5439
## [289] 5640.0000 3580.1563 1061.7000 2019.7881 3169.1899 4318.5918
## [295] 5467.9937 6180.0000 6844.0718 7200.0000 8390.0000 7956.4258
## [301] 7517.4834 7180.0000 562.5400 1348.7238 2818.8896 3680.0000
## [307] 2216.6179 827.6600 1538.3466 2890.0000 2383.2690 2029.4900
## [313] 3820.7578 5630.0000 5202.5576 4560.0000 6232.8345 7430.0000
## [319] 6186.6069 4830.0000 5338.7207 6070.0000 6447.5181 6879.2051
## [325] 7260.0000 5067.7798 2270.0000 1917.3143 1302.1428 965.5400
## [331] 2610.9341 3743.4775 4876.0210 6008.5645 7141.1079 7930.0000
## [337] 8742.7529 9120.0000 9898.4785 10540.0000 9533.3623 8030.0000
## [343] 6911.2437 5430.0000 5048.1919 4678.7793 4490.0000 3840.0000
## [349] 2968.5474 1508.5300 2100.5769 2981.4578 3750.0000 7119.8560
## [355] 8820.0000 2032.8934 1201.6453 226.0300 131.3576 7.3800
## [361] 397.6351 999.0900 841.6011 572.5554 303.5096 143.3300
## [367] 119.7091 102.8927 86.0763 69.2599 52.4436 35.6272
## [373] 22.0900 759.7528 1377.0441 1994.3353 2640.0000 3427.3506
## [379] 4390.0000 3882.9058 3100.0000 3340.3401 3677.3889 4014.4380
## [385] 4243.5298 3936.2375 3750.0000 3700.0000 3181.6086 2457.0735
## [391] 2010.0000 2805.6685 3230.0000 3302.7673 3340.0000 2963.2549
## [397] 2750.0000 2947.7891 3090.0000 3110.4392 3128.8020 3147.1648
## [403] 3170.0000 3033.0540 2844.1299 2137.3313 1366.1600 1479.3800
## [409] 1291.2828 1042.6234 683.6700 815.7291 1053.0554 1290.3817
## [415] 1585.5100 1510.6265 1411.6324 1312.6384 1213.6443 1114.6503
## [421] 1015.6563 916.6622 817.6682 718.6742 619.6802 520.6862
## [427] 421.6921 322.6981 223.7040 124.7100 3.2300 33.7340
## [433] 73.2033 112.6725 152.1418 191.6111 231.0803 270.5496
## [439] 310.0188 349.4881 388.9573 428.4266 467.8958 507.3651
## [445] 546.8344 586.3036 625.7729 665.2421 704.7114 744.1807
## [451] 783.6499 823.1191 862.5884 902.0577 941.5269 980.9962
## [457] 1020.4655 1059.9347 1099.4039 1178.3425 1217.8118 1257.2810
## [463] 1280.9800 1380.1622 1479.8101 1574.8315 1790.6000 2417.5913
## [469] 3660.0000 4203.3662 4969.6675 5810.0000 4703.0840 3430.0000
## [475] 4832.6045 5950.0000 6389.9697 6950.0000 7862.0825 8190.0000
## [481] 7020.0000 6917.3823 6680.0000 6388.7993 5750.0000 5412.3223
## [487] 4954.9937 4630.0000 4451.9521 4320.0000 4047.0544 3840.0000
## [493] 3872.8757 3890.0000 3997.4805 4075.8799 4154.2788 4210.0000
## [499] 2520.0000 4328.6606 6027.3105 7725.9600 8760.0000 7740.0000
## [505] 8306.0459 8820.0000 7918.8984 6860.0000 6171.7534 5200.0000
## [511] 8820.1025 10860.0000 6697.7793 5190.0000 5363.0371 5479.2744
## [517] 5595.5122 5963.4404 6280.0000 6138.8628 6046.8418 5954.8203
## [523] 5880.0000 5266.6743 4749.9390 4420.0000 4119.3125 3940.0000
## [529] 4011.0640 4070.0000 4267.1182 3970.0000 3982.3992 3998.8315
## [535] 4020.0000 3520.5999 2495.1101 2126.4495 1441.6909 1006.8500
## [541] 884.5511 831.5900 2790.4949 4760.0000 3950.7856 2047.1400
## [547] 2114.7256 2234.4399 3956.0886 4810.0000 3983.7849 3290.0000
## [553] 2250.9197 1454.2400 2423.3625 3650.0000 3328.5642 2492.1201
## [559] 2698.5298 3180.0000 3293.4705 3473.3303 3610.0000 1181.8900
## [565] 1652.2738 2590.4043 3320.0000 2433.3545 1868.6801 1899.3845
## [571] 1913.9399 3405.0000 3385.5239 3361.7825 3338.0410 3310.0000
## [577] 2997.9072 2700.0000 2913.8325 2970.0000 3150.3203 3250.0000
## [583] 2283.9041 1820.0200 1922.0161 1978.5601 2853.9758 3420.0000
## [589] 2846.9758 2570.3201 2393.9656 2295.7800 2918.3340 3220.0000
## [595] 2573.2031 2303.7000 2736.5945 3040.0000 3988.5752 4730.0000
## [601] 3795.9231 3340.0000 2749.5947 2531.0000 2348.8201 2895.5957
## [607] 3720.0000 3561.8772 3350.0000 3459.5876 3630.0000 3526.3140
## [613] 3460.0000 2262.0811 1371.3700 1591.2185 1736.2400 1828.3900
## [619] 1952.0089 2091.1111 2183.5400 2331.9900 1619.7656 809.4100
## [625] 603.8804 458.8995 368.3500 383.8092 392.7800 369.4849
## [631] 350.6500 1862.8700 1181.2034 258.1300 465.0500 371.8761
## [637] 287.3767 279.4883 271.5998 263.7113 255.8228 247.9343
## [643] 240.0458 228.5700 242.9852 282.0900 332.4404 458.2000
## [649] 437.7937 399.1000 431.1000 485.5200 459.3311 417.1300
## [655] 480.6394 582.1041 643.0700 499.3700 447.7827 358.4117
## [661] 300.3600 240.2611 210.2700 216.6000 257.7916 339.4923
## [667] 398.5300 669.3860 827.7700 896.8858 939.4300 780.4501
## [673] 677.5000 710.9836 739.4100 729.7252 721.3600 917.3568
## [679] 1064.6899 963.2323 883.1400 999.1555 1074.0900 1209.9685
## [685] 1309.5500 1053.1267 868.0800 865.2035 863.3700 1137.1615
## [691] 1274.9500 1259.5554 1248.4808 1237.4061 1227.8700 1102.0452
## [697] 991.5684 881.0916 770.6147 660.1379 512.2100 537.4426
## [703] 575.6158 613.7889 651.9620 690.1351 728.3082 767.2200
## [709] 729.5110 691.0581 652.6051 614.1521 564.8300 609.7765
## [715] 672.4337 735.0908 797.7479 860.4050 945.3500 770.6907
## [721] 664.4689 558.2472 452.0254 300.2600 321.5066 358.7004
## [727] 395.8943 433.0881 470.2819 507.4758 526.1400 893.3821
## [733] 1138.5061 1383.6300 1628.7540 1873.8779 2014.7900 1743.1687
## [739] 1552.5758 1361.9829 1171.3900 980.7972 857.3100 1089.8239
## [745] 1261.7902 1433.7565 1605.7228 1777.6891 1708.2007 1556.9391
## [751] 1405.6776 1254.4161 1103.1545 921.8100 910.1433 895.5803
## [757] 881.0173 866.4542 851.8912 830.4200 817.8306 793.8796
## [763] 769.9286 745.9776 722.0266 699.5100 831.6595 956.3419
## [769] 1205.7067 1330.3890 1439.0200 1876.8239 2264.6938 2652.5637
## [775] 3040.4338 3428.3037 3650.0000 2797.1609 2200.1130 1603.0649
## [781] 1006.0170 246.6400 298.5033 369.7330 440.9626 512.1923
## [787] 583.4219 679.7400 712.3636 762.7256 813.0876 863.4496
## [793] 913.8117 957.2500 1528.5027 2030.7120 2532.9214 3035.1306
## [799] 3820.0000 5020.0947 5855.1367 6690.1787 7525.2207 9458.1729
## [805] 10649.9580 11841.7432 13033.5283 14225.3135 15600.0000 13544.9785
## [811] 11117.4023 8689.8252 6262.2485 1300.7900 1249.7059 1196.2823
## [817] 1142.8588 1089.4352 1036.0116Ordering data of permeability core(Decreasing)
data = data[order(data$k.core, decreasing=TRUE),]
data$k.core## [1] 15600.0000 14225.3135 13544.9785 13033.5283 11841.7432 11117.4023
## [7] 10860.0000 10649.9580 10540.0000 9898.4785 9533.3623 9458.1729
## [13] 9120.0000 8820.1025 8820.0000 8820.0000 8760.0000 8742.7529
## [19] 8689.8252 8390.0000 8306.0459 8190.0000 8030.0000 7956.4258
## [25] 7930.0000 7918.8984 7862.0825 7740.0000 7725.9600 7525.2207
## [31] 7517.4834 7430.0000 7260.0000 7200.0000 7180.0000 7141.1079
## [37] 7119.8560 7020.0000 6950.0000 6917.3823 6911.2437 6879.2051
## [43] 6860.0000 6844.0718 6697.7793 6690.1787 6680.0000 6447.5181
## [49] 6389.9697 6388.7993 6280.0000 6262.2485 6232.8345 6186.6069
## [55] 6180.0000 6171.7534 6138.8628 6070.0000 6046.8418 6027.3105
## [61] 6008.5645 5963.4404 5954.8203 5950.0000 5880.0000 5855.1367
## [67] 5810.0000 5750.0000 5640.0000 5630.0000 5595.5122 5480.0000
## [73] 5479.2744 5467.9937 5430.0000 5412.3223 5363.0371 5338.7207
## [79] 5266.6743 5202.5576 5200.0000 5190.0000 5067.7798 5048.1919
## [85] 5020.0947 4969.6675 4954.9937 4900.0000 4876.0210 4832.6045
## [91] 4830.0000 4810.0000 4804.8022 4786.5928 4760.0000 4749.9390
## [97] 4730.0000 4703.0840 4678.7793 4630.0000 4561.4106 4560.0000
## [103] 4490.0000 4451.9521 4420.0000 4390.0000 4390.0000 4336.2285
## [109] 4328.6606 4320.0000 4318.5918 4267.1182 4243.5298 4210.0000
## [115] 4203.3662 4154.2788 4135.5439 4119.3125 4111.0464 4075.8799
## [121] 4070.0000 4047.0544 4020.0000 4014.4380 4011.0640 3998.8315
## [127] 3997.4805 3988.5752 3983.7849 3982.3992 3970.0000 3956.0886
## [133] 3950.7856 3940.0000 3936.2375 3890.0000 3885.8645 3882.9058
## [139] 3872.8757 3840.0000 3840.0000 3820.7578 3820.0000 3795.9231
## [145] 3750.0000 3750.0000 3743.4775 3720.0000 3700.0000 3680.0000
## [151] 3677.3889 3660.6824 3660.0000 3650.0000 3650.0000 3630.0000
## [157] 3630.0000 3610.0000 3580.1563 3561.8772 3526.3140 3520.5999
## [163] 3473.3303 3460.0000 3459.5876 3435.5005 3430.0000 3428.3037
## [169] 3427.3506 3420.0000 3405.0000 3385.5239 3380.0000 3370.0000
## [175] 3361.7825 3350.0000 3340.3401 3340.0000 3340.0000 3338.0410
## [181] 3328.5642 3320.0000 3310.0000 3302.7673 3293.4705 3290.0000
## [187] 3250.0000 3250.0000 3230.0000 3220.0000 3210.3184 3181.6086
## [193] 3180.0000 3170.0000 3169.1899 3150.3203 3147.1648 3128.8020
## [199] 3110.4392 3110.0000 3100.0000 3090.0657 3090.0000 3066.8655
## [205] 3053.2048 3040.4338 3040.0000 3040.0000 3035.1306 3033.0540
## [211] 3006.9888 3000.0000 2997.9072 2990.0000 2985.1362 2981.4578
## [217] 2970.0000 2968.5474 2963.2549 2960.0000 2947.7891 2928.3137
## [223] 2928.1562 2918.3340 2913.8325 2895.5957 2890.0000 2865.1624
## [229] 2853.9758 2846.9758 2844.1299 2832.4385 2830.0000 2818.8896
## [235] 2805.6685 2802.1687 2797.1609 2790.4949 2760.0000 2759.9543
## [241] 2759.4919 2750.0000 2749.5947 2736.5945 2730.0000 2707.6108
## [247] 2700.0000 2698.5298 2697.3630 2686.7124 2660.0000 2652.5637
## [253] 2648.6399 2640.0000 2638.9934 2620.0000 2613.3105 2610.9341
## [259] 2602.4993 2593.2590 2590.4043 2573.2031 2570.3201 2570.0000
## [265] 2563.3101 2558.1506 2550.4099 2534.7722 2532.9214 2531.0762
## [271] 2531.0000 2530.7581 2523.2039 2520.0000 2519.8000 2495.1360
## [277] 2495.1101 2492.1201 2485.9099 2483.5156 2478.3799 2473.3740
## [283] 2468.4099 2464.1960 2463.9275 2457.0735 2455.0183 2450.0000
## [289] 2448.0698 2445.8406 2442.5901 2436.6628 2434.4473 2433.3545
## [295] 2427.4851 2423.3625 2423.1160 2421.8501 2417.5913 2393.9656
## [301] 2383.2690 2371.0400 2368.7905 2367.5498 2348.8201 2346.7214
## [307] 2335.4399 2331.9900 2322.1902 2309.5901 2303.7000 2302.9199
## [313] 2299.0537 2295.7800 2294.7261 2285.0200 2284.4299 2283.9041
## [319] 2276.3159 2274.8657 2273.6499 2270.0000 2270.0000 2264.6938
## [325] 2262.0811 2250.9197 2234.4399 2225.8699 2216.6179 2214.3201
## [331] 2202.3936 2200.3501 2200.1130 2196.3706 2192.6060 2183.5400
## [337] 2182.8186 2173.0312 2163.2437 2155.4399 2155.4353 2150.3301
## [343] 2149.6299 2147.4971 2137.3313 2126.4495 2114.7256 2100.5769
## [349] 2091.1111 2089.7009 2084.4082 2067.3738 2056.2000 2049.1858
## [355] 2047.1400 2032.8934 2030.7120 2029.4900 2023.9753 2023.0345
## [361] 2019.7881 2014.7900 2011.7695 2010.0000 1994.3353 1990.0415
## [367] 1978.5601 1975.0400 1969.0728 1956.1077 1952.0089 1931.9360
## [373] 1928.4800 1922.0161 1917.3143 1913.9399 1907.3199 1899.3845
## [379] 1884.6757 1876.8239 1873.8779 1868.6801 1866.9200 1866.0350
## [385] 1862.8700 1859.2261 1852.2700 1828.3900 1821.4403 1820.0200
## [391] 1804.5400 1790.6000 1789.1849 1777.6891 1743.1687 1736.2400
## [397] 1722.9451 1712.4135 1708.2007 1679.6505 1652.2738 1648.4301
## [403] 1641.8264 1634.0441 1628.7540 1619.7656 1605.7228 1603.0649
## [409] 1599.2813 1591.2185 1585.5100 1584.9600 1574.8315 1562.2338
## [415] 1556.9391 1552.5758 1550.0200 1538.3466 1530.4865 1528.5027
## [421] 1510.6265 1508.5300 1498.7391 1481.5601 1479.8101 1479.3800
## [427] 1474.6252 1467.6908 1466.9917 1454.2400 1441.6909 1439.0200
## [433] 1437.5200 1435.2444 1433.7565 1411.6324 1408.8621 1405.6776
## [439] 1403.4971 1383.6300 1380.1622 1377.0441 1371.7496 1371.3700
## [445] 1366.1600 1362.4678 1361.9829 1348.7238 1340.0023 1330.3890
## [451] 1312.6384 1309.5500 1308.2550 1306.1700 1302.1428 1300.7900
## [457] 1291.5699 1291.2828 1290.3817 1280.9800 1276.5076 1274.9500
## [463] 1271.3800 1269.6000 1261.7902 1259.5554 1257.2810 1254.4161
## [469] 1249.7059 1248.4808 1244.7603 1237.4061 1227.8700 1217.8118
## [475] 1213.6443 1213.0129 1209.9685 1205.7067 1201.6453 1196.2823
## [481] 1183.6799 1181.8900 1181.2655 1181.2034 1178.3425 1171.3900
## [487] 1149.5182 1142.8588 1138.5061 1137.1615 1117.7709 1114.6503
## [493] 1109.5569 1103.1545 1102.0452 1099.4039 1089.8239 1089.4352
## [499] 1086.0234 1074.0900 1066.3792 1064.6899 1064.5748 1061.7000
## [505] 1059.9347 1054.2761 1053.1267 1053.0554 1042.6234 1040.1243
## [511] 1036.0116 1022.5287 1020.4655 1015.6563 1006.8500 1006.0170
## [517] 999.1555 999.0900 991.5684 990.7814 980.9962 980.7972
## [523] 965.5400 963.2323 959.0341 958.4979 957.2500 956.3419
## [529] 945.3500 941.5269 939.4300 927.2867 921.8100 917.3568
## [535] 916.6622 913.8117 910.1433 906.6724 902.0577 896.8858
## [541] 895.5803 895.5394 893.3821 884.5511 883.1400 881.0916
## [547] 881.0173 868.0800 866.4542 865.2035 863.7920 863.4496
## [553] 863.3700 862.5884 860.4050 859.5496 857.3100 851.8912
## [559] 841.6011 832.0446 831.6595 831.5900 830.4200 827.7700
## [565] 827.6600 823.1191 822.9766 817.8306 817.6682 815.7291
## [571] 813.0876 809.4100 800.2973 797.7479 794.9600 793.8796
## [577] 783.6499 780.4501 770.6907 770.6147 769.9286 768.5499
## [583] 767.2200 762.7256 759.7528 745.9776 744.1807 739.4100
## [589] 736.8026 735.0908 733.3159 729.7252 729.5110 728.3082
## [595] 722.0266 721.3600 718.6742 712.3636 710.9836 705.0552
## [601] 704.7114 699.5100 691.0581 690.1351 683.6700 679.7400
## [607] 677.5000 673.3079 672.4337 669.3860 665.2421 664.4689
## [613] 662.5491 660.1379 654.5244 652.6051 651.9620 643.0700
## [619] 641.5605 635.5146 625.7729 619.6802 614.1521 613.7889
## [625] 609.8132 609.7765 603.8804 601.2100 599.3371 597.4316
## [631] 595.5261 593.6207 591.7152 589.8097 588.2900 586.3036
## [637] 583.4219 582.1041 578.0658 575.6158 572.5554 564.8300
## [643] 562.5400 558.2472 546.8344 546.3184 537.4426 534.7041
## [649] 526.1400 522.3776 520.6862 514.5711 512.2100 512.1923
## [655] 508.1339 507.4758 507.3651 499.3700 485.5200 482.8237
## [661] 480.6394 470.2819 467.8958 467.5622 465.0500 459.3311
## [667] 458.8995 458.2000 454.5060 452.0254 451.0764 449.4991
## [673] 447.7827 440.9626 437.7937 433.0881 431.1000 428.4266
## [679] 421.6921 419.3290 418.6400 417.1300 412.7467 406.2400
## [685] 399.1000 398.5300 397.6351 395.8943 392.7800 388.9573
## [691] 387.5817 386.1200 383.8092 371.8761 369.7330 369.4849
## [697] 368.3500 358.7004 358.4117 357.9313 355.8343 350.6500
## [703] 349.4881 339.4923 332.4404 324.0869 322.6981 321.5066
## [709] 316.3239 310.0188 308.4382 303.5096 303.1158 300.3600
## [715] 300.2600 298.5033 292.3396 287.3767 282.9963 282.9518
## [721] 282.0900 279.4883 273.4974 271.5998 271.4900 270.5496
## [727] 263.7113 260.5923 258.1300 257.7916 255.8228 248.3004
## [733] 247.9343 246.6400 244.4739 242.9852 240.2611 240.0458
## [739] 231.2015 231.0803 228.8449 228.5700 226.0300 223.7040
## [745] 216.6000 210.2700 204.7999 197.0975 193.4600 191.6111
## [751] 190.5800 185.0957 181.7682 173.4400 168.2400 165.3502
## [757] 162.9400 158.1300 152.9717 152.1418 149.2137 146.7933
## [763] 146.0790 143.3300 140.6149 140.0700 134.4366 133.6028
## [769] 131.3576 128.2582 126.0047 124.7100 122.0799 119.7091
## [775] 116.5507 115.9015 114.3747 112.6725 109.7232 103.5448
## [781] 102.8927 101.8555 97.3665 92.4887 91.1881 86.9600
## [787] 86.0763 85.0097 76.8100 74.1400 73.2072 73.2033
## [793] 70.1081 69.2599 63.8231 55.5300 54.2100 52.4436
## [799] 50.1500 48.6200 44.5030 43.9438 39.6700 39.4679
## [805] 35.8600 35.6272 34.9919 33.7340 32.5300 30.6700
## [811] 22.0900 7.3800 5.8300 4.8748 3.5809 3.2300
## [817] 2.2869 0.8100 0.4200Finding percentage of permeability:
permeability =data$k.core
sample=c(1:length(permeability))
permeability_percent=(sample*100)/length(permeability) Plotting samples versus permeability
xlab = "Portion of Total Samples Having Larger or Equal Permeability "
ylab = "Permeability(md)"
plot(permeability_percent,permeability,log = 'y',xlab = xlab,ylab=ylab,pch=10,cex=0.5,col ="#091c90")
log_permeability = log(permeability)
model = lm(log_permeability~permeability_percent)
plot(permeability_percent,log_permeability,xlab = xlab,ylab=ylab,pch=10,cex=0.5,col ="#091c90")
abline(model,col = 'green',lwd = 2)
summary(model)##
## Call:
## lm(formula = log_permeability ~ permeability_percent)
##
## Residuals:
## Min 1Q Median 3Q Max
## -5.8697 -0.2047 0.1235 0.3150 0.4280
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 9.2584172 0.0377994 244.94 <2e-16 ***
## permeability_percent -0.0425617 0.0006541 -65.07 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5404 on 817 degrees of freedom
## Multiple R-squared: 0.8382, Adjusted R-squared: 0.838
## F-statistic: 4234 on 1 and 817 DF, p-value: < 2.2e-16Now, I am going to calculate the heterogeneity index:
new_data = data.frame(permeability_percent = c(50,84.1))
predicted_values = predict(model,new_data)
heterogeneity_index = (predicted_values[1]-predicted_values[2])/predicted_values[1]
heterogeneity_index## 1
## 0.2035464