R Markdown

setwd("C:/Users/Intel/Desktop/Yousif")

yousif2<-read.csv("c:/Users/Intel/Desktop/Yousif/Karpur.csv")
plot(yousif2$phi.N,yousif2$phi.core)

porosity_model <- lm(yousif2$phi.core~yousif2$phi.N)

summary(porosity_model)
## 
## Call:
## lm(formula = yousif2$phi.core ~ yousif2$phi.N)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -13.5237  -3.0779   0.9432   3.3563  10.4025 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     30.962      0.485  63.846   <2e-16 ***
## yousif2$phi.N  -18.207      2.080  -8.753   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.368 on 817 degrees of freedom
## Multiple R-squared:  0.08573,    Adjusted R-squared:  0.08462 
## F-statistic: 76.61 on 1 and 817 DF,  p-value: < 2.2e-16
corrected_porosity <- predict(porosity_model)
head(cbind(yousif2$phi.core,predict(porosity_model)))
##      [,1]     [,2]
## 1 33.9000 23.49747
## 2 33.4131 25.37279
## 3 33.1000 27.26631
## 4 34.9000 28.79570
## 5 35.0644 29.70604
## 6 35.3152 30.08839
# porosity corrected to log Scale

print(corrected_porosity)
##        1        2        3        4        5        6        7        8 
## 23.49747 25.37279 27.26631 28.79570 29.70604 30.08839 30.10660 29.96094 
##        9       10       11       12       13       14       15       16 
## 29.76066 29.61501 29.50577 29.66963 29.92453 30.07018 30.08839 30.05198 
##       17       18       19       20       21       22       23       24 
## 29.90632 29.72425 29.54218 29.36011 29.03239 28.74107 28.74107 28.90494 
##       25       26       27       28       29       30       31       32 
## 28.95956 28.86852 28.83211 28.79570 28.52259 28.12204 28.01280 28.35873 
##       33       34       35       36       37       38       39       40 
## 28.79570 29.15983 29.43294 29.54218 29.65142 29.54218 29.34190 29.23266 
##       41       42       43       44       45       46       47       48 
## 29.19625 29.14163 29.05059 28.85032 28.44976 28.59542 28.74107 29.15983 
##       49       50       51       52       53       54       55       56 
## 29.41473 29.52397 29.56039 29.63322 29.65142 29.65142 29.63322 29.56039 
##       57       58       59       60       61       62       63       64 
## 29.54218 29.59680 29.63322 29.57859 29.50577 29.45115 29.43294 29.39653 
##       65       66       67       68       69       70       71       72 
## 29.34190 29.30549 29.25087 29.23266 29.32370 29.41473 29.46935 29.48756 
##       73       74       75       76       77       78       79       80 
## 29.48756 29.50577 29.52397 29.54218 29.56039 29.57859 29.57859 29.56039 
##       81       82       83       84       85       86       87       88 
## 29.59680 29.59680 29.56039 29.57859 29.59680 29.61501 29.61501 29.61501 
##       89       90       91       92       93       94       95       96 
## 29.61501 29.61501 29.59680 29.59680 29.61501 29.61501 29.61501 29.61501 
##       97       98       99      100      101      102      103      104 
## 29.59680 29.57859 29.57859 29.61501 29.65142 29.68784 29.66963 29.56039 
##      105      106      107      108      109      110      111      112 
## 29.43294 29.32370 28.97777 28.65004 28.55901 28.70466 28.94135 29.01418 
##      113      114      115      116      117      118      119      120 
## 28.85032 28.50438 27.83073 26.93859 26.33776 26.06465 25.79155 25.77334 
##      121      122      123      124      125      126      127      128 
## 26.06465 26.33776 26.48341 26.61086 26.62907 26.64727 26.66548 26.62907 
##      129      130      131      132      133      134      135      136 
## 26.42879 26.10107 25.68231 25.33637 25.28175 25.40920 25.51844 25.55486 
##      137      138      139      140      141      142      143      144 
## 25.62768 25.77334 25.90079 25.97362 26.11927 26.37417 26.50162 26.44700 
##      145      146      147      148      149      150      151      152 
## 26.37417 26.08286 25.70051 25.39099 25.22713 25.35458 25.73693 26.17389 
##      153      154      155      156      157      158      159      160 
## 26.44700 26.66548 26.88397 27.08424 27.17528 27.21169 27.22990 27.21169 
##      161      162      163      164      165      166      167      168 
## 27.08424 26.86576 26.73831 26.72010 26.70190 26.66548 26.64727 26.68369 
##      169      170      171      172      173      174      175      176 
## 26.70190 26.61086 26.42879 26.22851 25.99182 25.71872 25.60948 25.84617 
##      177      178      179      180      181      182      183      184 
## 26.24672 26.50162 26.64727 26.84755 27.30273 28.10383 28.32231 28.52259 
##      185      186      187      188      189      190      191      192 
## 28.44976 27.92176 27.30273 26.97500 26.75652 26.55624 26.37417 26.30134 
##      193      194      195      196      197      198      199      200 
## 26.41058 26.37417 26.44700 26.82934 26.88397 26.88397 27.10245 27.75790 
##      201      202      203      204      205      206      207      208 
## 28.74107 29.52397 29.74246 29.96094 30.16122 30.05198 29.65142 28.99597 
##      209      210      211      212      213      214      215      216 
## 28.23128 27.53942 27.12066 26.88397 26.77472 26.75652 26.75652 26.73831 
##      217      218      219      220      221      222      223      224 
## 26.73831 26.72010 26.66548 26.55624 26.41058 26.26493 26.24672 26.21031 
##      225      226      227      228      229      230      231      232 
## 26.30134 26.59265 27.12066 27.73969 28.17666 28.24949 27.97638 27.61224 
##      233      234      235      236      237      238      239      240 
## 27.26631 27.01141 26.81114 26.70190 26.59265 26.46520 26.37417 26.37417 
##      241      242      243      244      245      246      247      248 
## 26.41058 26.42879 26.42879 26.37417 26.30134 26.26493 26.22851 26.19210 
##      249      250      251      252      253      254      255      256 
## 26.15569 26.15569 26.21031 26.28314 26.33776 26.31955 26.30134 26.24672 
##      257      258      259      260      261      262      263      264 
## 26.17389 26.11927 26.10107 26.01003 25.77334 25.71872 25.68231 25.75513 
##      265      266      267      268      269      270      271      272 
## 25.73693 25.70051 25.71872 25.77334 25.84617 25.93720 26.01003 26.10107 
##      273      274      275      276      277      278      279      280 
## 26.15569 26.19210 26.19210 26.17389 26.22851 26.53803 27.12066 27.73969 
##      281      282      283      284      285      286      287      288 
## 27.92176 28.08562 28.21307 28.37694 28.61363 28.86852 29.14163 29.41473 
##      289      290      291      292      293      294      295      296 
## 29.66963 29.88811 30.05198 30.17942 30.28867 30.39791 30.48894 30.57998 
##      297      298      299      300      301      302      303      304 
## 30.65281 30.68922 30.68922 30.67101 30.61639 30.56177 30.34329 29.94273 
##      305      306      307      308      309      310      311      312 
## 29.36011 28.48618 27.44838 26.62907 26.26493 26.13748 26.26493 26.48341 
##      313      314      315      316      317      318      319      320 
## 26.64727 26.90217 27.15707 27.28452 27.33914 27.35735 27.35735 27.41197 
##      321      322      323      324      325      326      327      328 
## 27.48479 27.52121 27.52121 27.48479 27.43017 27.43017 27.46659 27.50300 
##      329      330      331      332      333      334      335      336 
## 27.53942 27.55762 27.59404 27.63045 27.61224 27.44838 27.24810 27.19348 
##      337      338      339      340      341      342      343      344 
## 27.19348 27.04783 26.93859 27.01141 27.12066 27.19348 27.21169 27.15707 
##      345      346      347      348      349      350      351      352 
## 27.04783 26.92038 26.84755 26.86576 26.88397 26.84755 26.86576 27.10245 
##      353      354      355      356      357      358      359      360 
## 27.32093 27.46659 27.44838 27.15707 26.68369 26.19210 25.75513 25.59127 
##      361      362      363      364      365      366      367      368 
## 25.42741 25.31817 25.42741 25.71872 25.97362 26.08286 26.11927 26.15569 
##      369      370      371      372      373      374      375      376 
## 26.21031 26.35596 26.53803 26.68369 26.84755 26.92038 26.93859 26.92038 
##      377      378      379      380      381      382      383      384 
## 26.86576 26.81114 26.73831 26.62907 26.57445 26.51983 26.50162 26.59265 
##      385      386      387      388      389      390      391      392 
## 26.61086 26.53803 26.50162 26.42879 26.35596 26.31955 26.33776 26.37417 
##      393      394      395      396      397      398      399      400 
## 26.39238 26.41058 26.37417 26.28314 26.15569 26.11927 26.17389 26.21031 
##      401      402      403      404      405      406      407      408 
## 26.22851 26.24672 26.33776 26.48341 26.53803 26.50162 26.48341 26.53803 
##      409      410      411      412      413      414      415      416 
## 26.66548 26.72010 26.73831 26.79293 26.88397 26.95679 27.06603 27.13886 
##      417      418      419      420      421      422      423      424 
## 27.12066 27.17528 27.17528 27.10245 26.93859 26.72010 26.64727 26.62907 
##      425      426      427      428      429      430      431      432 
## 26.44700 26.31955 26.21031 25.99182 25.71872 25.50024 25.35458 25.26355 
##      433      434      435      436      437      438      439      440 
## 25.19072 25.15430 25.15430 25.15430 25.04506 24.88120 24.80837 24.80837 
##      441      442      443      444      445      446      447      448 
## 24.55347 24.02547 23.77058 23.77058 23.66133 23.64313 23.62492 23.64313 
##      449      450      451      452      453      454      455      456 
## 23.77058 23.97085 24.20754 24.46244 24.77196 24.99044 25.15430 25.33637 
##      457      458      459      460      461      462      463      464 
## 25.59127 25.82796 25.88258 25.88258 25.84617 25.86438 25.82796 25.73693 
##      465      466      467      468      469      470      471      472 
## 25.68231 25.70051 25.75513 25.80975 25.90079 26.06465 26.22851 26.35596 
##      473      474      475      476      477      478      479      480 
## 26.41058 26.46520 26.53803 26.59265 26.57445 26.53803 26.48341 26.33776 
##      481      482      483      484      485      486      487      488 
## 26.17389 26.10107 26.08286 26.11927 26.15569 26.15569 26.15569 26.15569 
##      489      490      491      492      493      494      495      496 
## 26.13748 26.15569 26.15569 26.24672 26.19210 26.19210 26.17389 26.19210 
##      497      498      499      500      501      502      503      504 
## 26.28314 26.48341 26.81114 27.08424 27.24810 27.28452 27.32093 27.33914 
##      505      506      507      508      509      510      511      512 
## 27.39376 27.46659 27.53942 27.63045 27.70328 27.70328 27.59404 27.37555 
##      513      514      515      516      517      518      519      520 
## 27.13886 26.97500 26.92038 26.90217 26.88397 26.75652 26.62907 26.53803 
##      521      522      523      524      525      526      527      528 
## 26.50162 26.51983 26.55624 26.57445 26.48341 26.39238 26.33776 26.35596 
##      529      530      531      532      533      534      535      536 
## 26.41058 26.41058 26.39238 26.41058 26.41058 26.39238 26.39238 26.33776 
##      537      538      539      540      541      542      543      544 
## 26.22851 26.11927 26.04644 26.02824 26.08286 26.26493 26.46520 26.53803 
##      545      546      547      548      549      550      551      552 
## 26.50162 26.46520 26.44700 26.41058 26.39238 26.35596 26.33776 26.28314 
##      553      554      555      556      557      558      559      560 
## 26.24672 26.28314 26.26493 26.15569 26.02824 25.88258 25.80975 25.79155 
##      561      562      563      564      565      566      567      568 
## 25.79155 25.79155 25.79155 25.80975 25.84617 25.88258 25.91900 25.95541 
##      569      570      571      572      573      574      575      576 
## 25.99182 26.08286 26.13748 26.15569 26.08286 26.01003 26.01003 26.04644 
##      577      578      579      580      581      582      583      584 
## 26.02824 25.99182 25.93720 25.86438 25.82796 25.82796 25.84617 25.80975 
##      585      586      587      588      589      590      591      592 
## 25.82796 25.86438 25.90079 25.99182 26.01003 25.99182 26.02824 26.11927 
##      593      594      595      596      597      598      599      600 
## 26.30134 26.61086 26.88397 26.88397 26.62907 26.33776 26.22851 26.22851 
##      601      602      603      604      605      606      607      608 
## 26.17389 26.04644 25.93720 25.88258 25.90079 25.95541 25.97362 26.01003 
##      609      610      611      612      613      614      615      616 
## 25.99182 25.95541 25.86438 25.70051 25.57306 25.57306 25.60948 25.66410 
##      617      618      619      620      621      622      623      624 
## 25.70051 25.70051 25.64589 25.55486 25.50024 25.50024 25.53665 25.57306 
##      625      626      627      628      629      630      631      632 
## 25.64589 25.73693 25.86438 26.02824 26.19210 26.28314 26.35596 26.35596 
##      633      634      635      636      637      638      639      640 
## 26.30134 26.24672 26.22851 26.24672 26.22851 26.21031 26.22851 26.24672 
##      641      642      643      644      645      646      647      648 
## 26.22851 26.17389 26.15569 26.15569 26.19210 26.21031 26.19210 26.15569 
##      649      650      651      652      653      654      655      656 
## 26.10107 26.02824 25.95541 25.93720 25.95541 26.02824 26.08286 26.10107 
##      657      658      659      660      661      662      663      664 
## 26.08286 26.10107 26.11927 26.17389 26.19210 26.17389 26.21031 26.30134 
##      665      666      667      668      669      670      671      672 
## 26.37417 26.42879 26.46520 26.44700 26.37417 26.33776 26.30134 26.21031 
##      673      674      675      676      677      678      679      680 
## 26.17389 26.17389 26.17389 26.17389 26.15569 26.11927 26.08286 26.04644 
##      681      682      683      684      685      686      687      688 
## 26.02824 26.01003 26.01003 26.02824 26.06465 26.11927 26.15569 26.15569 
##      689      690      691      692      693      694      695      696 
## 26.15569 26.15569 26.13748 26.04644 25.93720 25.86438 25.82796 25.79155 
##      697      698      699      700      701      702      703      704 
## 25.75513 25.73693 25.77334 25.80975 25.80975 25.80975 25.84617 25.93720 
##      705      706      707      708      709      710      711      712 
## 26.06465 26.15569 26.30134 26.44700 26.41058 26.31955 26.30134 26.35596 
##      713      714      715      716      717      718      719      720 
## 26.35596 26.33776 26.28314 26.30134 26.31955 26.28314 26.24672 26.30134 
##      721      722      723      724      725      726      727      728 
## 26.41058 26.51983 26.57445 26.57445 26.50162 26.39238 26.24672 26.19210 
##      729      730      731      732      733      734      735      736 
## 26.26493 26.35596 26.41058 26.41058 26.41058 26.42879 26.48341 26.50162 
##      737      738      739      740      741      742      743      744 
## 26.44700 26.44700 26.79293 27.43017 27.88535 27.92176 27.53942 26.95679 
##      745      746      747      748      749      750      751      752 
## 26.77472 26.72010 26.79293 26.92038 26.86576 26.84755 26.86576 26.90217 
##      753      754      755      756      757      758      759      760 
## 26.95679 27.01141 27.08424 27.33914 27.55762 27.61224 27.41197 27.10245 
##      761      762      763      764      765      766      767      768 
## 26.86576 26.75652 26.73831 26.75652 26.75652 26.70190 26.61086 26.57445 
##      769      770      771      772      773      774      775      776 
## 26.72010 26.77472 26.79293 26.82934 26.86576 26.86576 26.88397 26.84755 
##      777      778      779      780      781      782      783      784 
## 26.72010 26.59265 26.59265 26.70190 26.92038 27.08424 27.06603 26.88397 
##      785      786      787      788      789      790      791      792 
## 26.73831 26.75652 26.79293 26.75652 26.73831 26.77472 26.81114 26.92038 
##      793      794      795      796      797      798      799      800 
## 26.92038 26.77472 26.59265 26.44700 26.37417 26.42879 26.86576 26.95679 
##      801      802      803      804      805      806      807      808 
## 26.97500 27.02962 27.17528 27.26631 27.08424 27.02962 26.95679 26.93859 
##      809      810      811      812      813      814      815      816 
## 26.95679 26.90217 26.84755 26.86576 26.97500 27.17528 27.06603 26.81114 
##      817      818      819 
## 26.53803 26.30134 26.13748
permeability_model <- lm(yousif2$k.core~corrected_porosity)

summary(permeability_model)
## 
## Call:
## lm(formula = yousif2$k.core ~ corrected_porosity)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3187.3 -1406.4  -768.6   971.2 13338.4 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         -8852.2     1529.1  -5.789 1.01e-08 ***
## corrected_porosity    412.3       56.7   7.271 8.38e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2168 on 817 degrees of freedom
## Multiple R-squared:  0.06077,    Adjusted R-squared:  0.05962 
## F-statistic: 52.86 on 1 and 817 DF,  p-value: 8.376e-13
corrected_permeability <- predict(permeability_model)

head(cbind(yousif2$k.core,corrected_permeability))
##            corrected_permeability
## 1 2442.590               835.3875
## 2 3006.989              1608.5507
## 3 3370.000              2389.2202
## 4 2270.000              3019.7610
## 5 2530.758              3395.0829
## 6 2928.314              3552.7181
# pearmeability corrected to log scale

print(corrected_permeability)
##         1         2         3         4         5         6         7         8 
##  835.3875 1608.5507 2389.2202 3019.7610 3395.0829 3552.7181 3560.2246 3500.1731 
##         9        10        11        12        13        14        15        16 
## 3417.6022 3357.5507 3312.5121 3380.0700 3485.1602 3545.2117 3552.7181 3537.7052 
##        17        18        19        20        21        22        23        24 
## 3477.6537 3402.5894 3327.5250 3252.4606 3117.3447 2997.2417 2997.2417 3064.7996 
##        25        26        27        28        29        30        31        32 
## 3087.3190 3049.7868 3034.7739 3019.7610 2907.1644 2742.0228 2696.9842 2839.6065 
##        33        34        35        36        37        38        39        40 
## 3019.7610 3169.8898 3282.4863 3327.5250 3372.5636 3327.5250 3244.9542 3199.9155 
##        41        42        43        44        45        46        47        48 
## 3184.9027 3162.3833 3124.8511 3042.2803 2877.1387 2937.1902 2997.2417 3169.8898 
##        49        50        51        52        53        54        55        56 
## 3274.9799 3320.0185 3335.0314 3365.0572 3372.5636 3372.5636 3365.0572 3335.0314 
##        57        58        59        60        61        62        63        64 
## 3327.5250 3350.0443 3365.0572 3342.5379 3312.5121 3289.9928 3282.4863 3267.4735 
##        65        66        67        68        69        70        71        72 
## 3244.9542 3229.9413 3207.4220 3199.9155 3237.4477 3274.9799 3297.4992 3305.0057 
##        73        74        75        76        77        78        79        80 
## 3305.0057 3312.5121 3320.0185 3327.5250 3335.0314 3342.5379 3342.5379 3335.0314 
##        81        82        83        84        85        86        87        88 
## 3350.0443 3350.0443 3335.0314 3342.5379 3350.0443 3357.5507 3357.5507 3357.5507 
##        89        90        91        92        93        94        95        96 
## 3357.5507 3357.5507 3350.0443 3350.0443 3357.5507 3357.5507 3357.5507 3357.5507 
##        97        98        99       100       101       102       103       104 
## 3350.0443 3342.5379 3342.5379 3357.5507 3372.5636 3387.5765 3380.0700 3335.0314 
##       105       106       107       108       109       110       111       112 
## 3282.4863 3237.4477 3094.8254 2959.7095 2922.1773 2982.2288 3079.8125 3109.8383 
##       113       114       115       116       117       118       119       120 
## 3042.2803 2899.6580 2621.9198 2254.1043 2006.3919 1893.7953 1781.1987 1773.6923 
##       121       122       123       124       125       126       127       128 
## 1893.7953 2006.3919 2066.4434 2118.9884 2126.4949 2134.0013 2141.5078 2126.4949 
##       129       130       131       132       133       134       135       136 
## 2043.9241 1908.8082 1736.1601 1593.5378 1571.0185 1623.5635 1668.6022 1683.6150 
##       137       138       139       140       141       142       143       144 
## 1713.6408 1773.6923 1826.2374 1856.2631 1916.3146 2021.4047 2073.9498 2051.4305 
##       145       146       147       148       149       150       151       152 
## 2021.4047 1901.3017 1743.6665 1616.0571 1548.4991 1601.0442 1758.6794 1938.8339 
##       153       154       155       156       157       158       159       160 
## 2051.4305 2141.5078 2231.5850 2314.1558 2351.6880 2366.7009 2374.2073 2366.7009 
##       161       162       163       164       165       166       167       168 
## 2314.1558 2224.0786 2171.5335 2164.0271 2156.5206 2141.5078 2134.0013 2149.0142 
##       169       170       171       172       173       174       175       176 
## 2156.5206 2118.9884 2043.9241 1961.3532 1863.7695 1751.1730 1706.1343 1803.7180 
##       177       178       179       180       181       182       183       184 
## 1968.8597 2073.9498 2134.0013 2216.5721 2404.2331 2734.5164 2824.5936 2907.1644 
##       185       186       187       188       189       190       191       192 
## 2877.1387 2659.4520 2404.2331 2269.1172 2179.0399 2096.4691 2021.4047 1991.3790 
##       193       194       195       196       197       198       199       200 
## 2036.4176 2021.4047 2051.4305 2209.0657 2231.5850 2231.5850 2321.6623 2591.8940 
##       201       202       203       204       205       206       207       208 
## 2997.2417 3320.0185 3410.0958 3500.1731 3582.7439 3537.7052 3372.5636 3102.3318 
##       209       210       211       212       213       214       215       216 
## 2787.0614 2501.8168 2329.1687 2231.5850 2186.5464 2179.0399 2179.0399 2171.5335 
##       217       218       219       220       221       222       223       224 
## 2171.5335 2164.0271 2141.5078 2096.4691 2036.4176 1976.3661 1968.8597 1953.8468 
##       225       226       227       228       229       230       231       232 
## 1991.3790 2111.4820 2329.1687 2584.3876 2764.5421 2794.5679 2681.9713 2531.8425 
##       233       234       235       236       237       238       239       240 
## 2389.2202 2284.1301 2201.5593 2156.5206 2111.4820 2058.9369 2021.4047 2021.4047 
##       241       242       243       244       245       246       247       248 
## 2036.4176 2043.9241 2043.9241 2021.4047 1991.3790 1976.3661 1961.3532 1946.3404 
##       249       250       251       252       253       254       255       256 
## 1931.3275 1931.3275 1953.8468 1983.8726 2006.3919 1998.8854 1991.3790 1968.8597 
##       257       258       259       260       261       262       263       264 
## 1938.8339 1916.3146 1908.8082 1871.2760 1773.6923 1751.1730 1736.1601 1766.1859 
##       265       266       267       268       269       270       271       272 
## 1758.6794 1743.6665 1751.1730 1773.6923 1803.7180 1841.2502 1871.2760 1908.8082 
##       273       274       275       276       277       278       279       280 
## 1931.3275 1946.3404 1946.3404 1938.8339 1961.3532 2088.9627 2329.1687 2584.3876 
##       281       282       283       284       285       286       287       288 
## 2659.4520 2727.0099 2779.5550 2847.1129 2944.6966 3049.7868 3162.3833 3274.9799 
##       289       290       291       292       293       294       295       296 
## 3380.0700 3470.1473 3537.7052 3590.2503 3635.2889 3680.3276 3717.8598 3755.3919 
##       297       298       299       300       301       302       303       304 
## 3785.4177 3800.4306 3800.4306 3792.9241 3770.4048 3747.8855 3657.8082 3492.6666 
##       305       306       307       308       309       310       311       312 
## 3252.4606 2892.1516 2464.2846 2126.4949 1976.3661 1923.8211 1976.3661 2066.4434 
##       313       314       315       316       317       318       319       320 
## 2134.0013 2239.0915 2344.1816 2396.7267 2419.2460 2426.7524 2426.7524 2449.2717 
##       321       322       323       324       325       326       327       328 
## 2479.2975 2494.3103 2494.3103 2479.2975 2456.7782 2456.7782 2471.7910 2486.8039 
##       329       330       331       332       333       334       335       336 
## 2501.8168 2509.3232 2524.3361 2539.3490 2531.8425 2464.2846 2381.7138 2359.1945 
##       337       338       339       340       341       342       343       344 
## 2359.1945 2299.1430 2254.1043 2284.1301 2329.1687 2359.1945 2366.7009 2344.1816 
##       345       346       347       348       349       350       351       352 
## 2299.1430 2246.5979 2216.5721 2224.0786 2231.5850 2216.5721 2224.0786 2321.6623 
##       353       354       355       356       357       358       359       360 
## 2411.7395 2471.7910 2464.2846 2344.1816 2149.0142 1946.3404 1766.1859 1698.6279 
##       361       362       363       364       365       366       367       368 
## 1631.0700 1586.0313 1631.0700 1751.1730 1856.2631 1901.3017 1916.3146 1931.3275 
##       369       370       371       372       373       374       375       376 
## 1953.8468 2013.8983 2088.9627 2149.0142 2216.5721 2246.5979 2254.1043 2246.5979 
##       377       378       379       380       381       382       383       384 
## 2224.0786 2201.5593 2171.5335 2126.4949 2103.9756 2081.4563 2073.9498 2111.4820 
##       385       386       387       388       389       390       391       392 
## 2118.9884 2088.9627 2073.9498 2043.9241 2013.8983 1998.8854 2006.3919 2021.4047 
##       393       394       395       396       397       398       399       400 
## 2028.9112 2036.4176 2021.4047 1983.8726 1931.3275 1916.3146 1938.8339 1953.8468 
##       401       402       403       404       405       406       407       408 
## 1961.3532 1968.8597 2006.3919 2066.4434 2088.9627 2073.9498 2066.4434 2088.9627 
##       409       410       411       412       413       414       415       416 
## 2141.5078 2164.0271 2171.5335 2194.0528 2231.5850 2261.6108 2306.6494 2336.6751 
##       417       418       419       420       421       422       423       424 
## 2329.1687 2351.6880 2351.6880 2321.6623 2254.1043 2164.0271 2134.0013 2126.4949 
##       425       426       427       428       429       430       431       432 
## 2051.4305 1998.8854 1953.8468 1863.7695 1751.1730 1661.0957 1601.0442 1563.5120 
##       433       434       435       436       437       438       439       440 
## 1533.4863 1518.4734 1518.4734 1518.4734 1473.4348 1405.8768 1375.8511 1375.8511 
##       441       442       443       444       445       446       447       448 
## 1270.7609 1053.0742  947.9841  947.9841  902.9455  895.4390  887.9326  895.4390 
##       449       450       451       452       453       454       455       456 
##  947.9841 1030.5549 1128.1386 1233.2288 1360.8382 1450.9155 1518.4734 1593.5378 
##       457       458       459       460       461       462       463       464 
## 1698.6279 1796.2116 1818.7309 1818.7309 1803.7180 1811.2245 1796.2116 1758.6794 
##       465       466       467       468       469       470       471       472 
## 1736.1601 1743.6665 1766.1859 1788.7052 1826.2374 1893.7953 1961.3532 2013.8983 
##       473       474       475       476       477       478       479       480 
## 2036.4176 2058.9369 2088.9627 2111.4820 2103.9756 2088.9627 2066.4434 2006.3919 
##       481       482       483       484       485       486       487       488 
## 1938.8339 1908.8082 1901.3017 1916.3146 1931.3275 1931.3275 1931.3275 1931.3275 
##       489       490       491       492       493       494       495       496 
## 1923.8211 1931.3275 1931.3275 1968.8597 1946.3404 1946.3404 1938.8339 1946.3404 
##       497       498       499       500       501       502       503       504 
## 1983.8726 2066.4434 2201.5593 2314.1558 2381.7138 2396.7267 2411.7395 2419.2460 
##       505       506       507       508       509       510       511       512 
## 2441.7653 2471.7910 2501.8168 2539.3490 2569.3747 2569.3747 2524.3361 2434.2588 
##       513       514       515       516       517       518       519       520 
## 2336.6751 2269.1172 2246.5979 2239.0915 2231.5850 2179.0399 2126.4949 2088.9627 
##       521       522       523       524       525       526       527       528 
## 2073.9498 2081.4563 2096.4691 2103.9756 2066.4434 2028.9112 2006.3919 2013.8983 
##       529       530       531       532       533       534       535       536 
## 2036.4176 2036.4176 2028.9112 2036.4176 2036.4176 2028.9112 2028.9112 2006.3919 
##       537       538       539       540       541       542       543       544 
## 1961.3532 1916.3146 1886.2889 1878.7824 1901.3017 1976.3661 2058.9369 2088.9627 
##       545       546       547       548       549       550       551       552 
## 2073.9498 2058.9369 2051.4305 2036.4176 2028.9112 2013.8983 2006.3919 1983.8726 
##       553       554       555       556       557       558       559       560 
## 1968.8597 1983.8726 1976.3661 1931.3275 1878.7824 1818.7309 1788.7052 1781.1987 
##       561       562       563       564       565       566       567       568 
## 1781.1987 1781.1987 1781.1987 1788.7052 1803.7180 1818.7309 1833.7438 1848.7567 
##       569       570       571       572       573       574       575       576 
## 1863.7695 1901.3017 1923.8211 1931.3275 1901.3017 1871.2760 1871.2760 1886.2889 
##       577       578       579       580       581       582       583       584 
## 1878.7824 1863.7695 1841.2502 1811.2245 1796.2116 1796.2116 1803.7180 1788.7052 
##       585       586       587       588       589       590       591       592 
## 1796.2116 1811.2245 1826.2374 1863.7695 1871.2760 1863.7695 1878.7824 1916.3146 
##       593       594       595       596       597       598       599       600 
## 1991.3790 2118.9884 2231.5850 2231.5850 2126.4949 2006.3919 1961.3532 1961.3532 
##       601       602       603       604       605       606       607       608 
## 1938.8339 1886.2889 1841.2502 1818.7309 1826.2374 1848.7567 1856.2631 1871.2760 
##       609       610       611       612       613       614       615       616 
## 1863.7695 1848.7567 1811.2245 1743.6665 1691.1215 1691.1215 1706.1343 1728.6537 
##       617       618       619       620       621       622       623       624 
## 1743.6665 1743.6665 1721.1472 1683.6150 1661.0957 1661.0957 1676.1086 1691.1215 
##       625       626       627       628       629       630       631       632 
## 1721.1472 1758.6794 1811.2245 1878.7824 1946.3404 1983.8726 2013.8983 2013.8983 
##       633       634       635       636       637       638       639       640 
## 1991.3790 1968.8597 1961.3532 1968.8597 1961.3532 1953.8468 1961.3532 1968.8597 
##       641       642       643       644       645       646       647       648 
## 1961.3532 1938.8339 1931.3275 1931.3275 1946.3404 1953.8468 1946.3404 1931.3275 
##       649       650       651       652       653       654       655       656 
## 1908.8082 1878.7824 1848.7567 1841.2502 1848.7567 1878.7824 1901.3017 1908.8082 
##       657       658       659       660       661       662       663       664 
## 1901.3017 1908.8082 1916.3146 1938.8339 1946.3404 1938.8339 1953.8468 1991.3790 
##       665       666       667       668       669       670       671       672 
## 2021.4047 2043.9241 2058.9369 2051.4305 2021.4047 2006.3919 1991.3790 1953.8468 
##       673       674       675       676       677       678       679       680 
## 1938.8339 1938.8339 1938.8339 1938.8339 1931.3275 1916.3146 1901.3017 1886.2889 
##       681       682       683       684       685       686       687       688 
## 1878.7824 1871.2760 1871.2760 1878.7824 1893.7953 1916.3146 1931.3275 1931.3275 
##       689       690       691       692       693       694       695       696 
## 1931.3275 1931.3275 1923.8211 1886.2889 1841.2502 1811.2245 1796.2116 1781.1987 
##       697       698       699       700       701       702       703       704 
## 1766.1859 1758.6794 1773.6923 1788.7052 1788.7052 1788.7052 1803.7180 1841.2502 
##       705       706       707       708       709       710       711       712 
## 1893.7953 1931.3275 1991.3790 2051.4305 2036.4176 1998.8854 1991.3790 2013.8983 
##       713       714       715       716       717       718       719       720 
## 2013.8983 2006.3919 1983.8726 1991.3790 1998.8854 1983.8726 1968.8597 1991.3790 
##       721       722       723       724       725       726       727       728 
## 2036.4176 2081.4563 2103.9756 2103.9756 2073.9498 2028.9112 1968.8597 1946.3404 
##       729       730       731       732       733       734       735       736 
## 1976.3661 2013.8983 2036.4176 2036.4176 2036.4176 2043.9241 2066.4434 2073.9498 
##       737       738       739       740       741       742       743       744 
## 2051.4305 2051.4305 2194.0528 2456.7782 2644.4391 2659.4520 2501.8168 2261.6108 
##       745       746       747       748       749       750       751       752 
## 2186.5464 2164.0271 2194.0528 2246.5979 2224.0786 2216.5721 2224.0786 2239.0915 
##       753       754       755       756       757       758       759       760 
## 2261.6108 2284.1301 2314.1558 2419.2460 2509.3232 2531.8425 2449.2717 2321.6623 
##       761       762       763       764       765       766       767       768 
## 2224.0786 2179.0399 2171.5335 2179.0399 2179.0399 2156.5206 2118.9884 2103.9756 
##       769       770       771       772       773       774       775       776 
## 2164.0271 2186.5464 2194.0528 2209.0657 2224.0786 2224.0786 2231.5850 2216.5721 
##       777       778       779       780       781       782       783       784 
## 2164.0271 2111.4820 2111.4820 2156.5206 2246.5979 2314.1558 2306.6494 2231.5850 
##       785       786       787       788       789       790       791       792 
## 2171.5335 2179.0399 2194.0528 2179.0399 2171.5335 2186.5464 2201.5593 2246.5979 
##       793       794       795       796       797       798       799       800 
## 2246.5979 2186.5464 2111.4820 2051.4305 2021.4047 2043.9241 2224.0786 2261.6108 
##       801       802       803       804       805       806       807       808 
## 2269.1172 2291.6365 2351.6880 2389.2202 2314.1558 2291.6365 2261.6108 2254.1043 
##       809       810       811       812       813       814       815       816 
## 2261.6108 2239.0915 2216.5721 2224.0786 2269.1172 2351.6880 2306.6494 2201.5593 
##       817       818       819 
## 2088.9627 1991.3790 1923.8211