HW01
Zeena Ahmed
2024-10-14
R Markdown
Simple linear regression
setwd("C:/Users/karar alshikh/Desktop/R")
getwd()## [1] "C:/Users/karar alshikh/Desktop/R"data2<-read.csv("C:/Users/karar alshikh/Desktop/R/karpur.csv")
plot(data2$phi.N,data2$phi.core.fraction)
porosity_model <- lm(data2$phi.core.fraction~data2$phi.N,data=data2)
summary(porosity_model)##
## Call:
## lm(formula = data2$phi.core.fraction ~ data2$phi.N, data = data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.135237 -0.030779 0.009432 0.033563 0.104025
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.30962 0.00485 63.846 <2e-16 ***
## data2$phi.N -0.18207 0.02080 -8.753 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.04368 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-16corrected_porosity <- predict(porosity_model)
head(cbind(data2$phi.core.fraction,predict(porosity_model)))## [,1] [,2]
## 1 0.339000 0.2349747
## 2 0.334131 0.2537279
## 3 0.331000 0.2726631
## 4 0.349000 0.2879570
## 5 0.350644 0.2970604
## 6 0.353152 0.3008839permeability_model <- lm(data2$k.core~corrected_porosity,data=data2)
summary(permeability_model)##
## Call:
## lm(formula = data2$k.core ~ corrected_porosity, data = data2)
##
## 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 1529 -5.789 1.01e-08 ***
## corrected_porosity 41228 5670 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# porosity corrected to log scale
print(corrected_porosity)## 1 2 3 4 5 6 7 8
## 0.2349747 0.2537279 0.2726631 0.2879570 0.2970604 0.3008839 0.3010660 0.2996094
## 9 10 11 12 13 14 15 16
## 0.2976066 0.2961501 0.2950577 0.2966963 0.2992453 0.3007018 0.3008839 0.3005198
## 17 18 19 20 21 22 23 24
## 0.2990632 0.2972425 0.2954218 0.2936011 0.2903239 0.2874107 0.2874107 0.2890494
## 25 26 27 28 29 30 31 32
## 0.2895956 0.2886852 0.2883211 0.2879570 0.2852259 0.2812204 0.2801280 0.2835873
## 33 34 35 36 37 38 39 40
## 0.2879570 0.2915983 0.2943294 0.2954218 0.2965142 0.2954218 0.2934190 0.2923266
## 41 42 43 44 45 46 47 48
## 0.2919625 0.2914163 0.2905059 0.2885032 0.2844976 0.2859542 0.2874107 0.2915983
## 49 50 51 52 53 54 55 56
## 0.2941473 0.2952397 0.2956039 0.2963322 0.2965142 0.2965142 0.2963322 0.2956039
## 57 58 59 60 61 62 63 64
## 0.2954218 0.2959680 0.2963322 0.2957859 0.2950577 0.2945115 0.2943294 0.2939653
## 65 66 67 68 69 70 71 72
## 0.2934190 0.2930549 0.2925087 0.2923266 0.2932370 0.2941473 0.2946935 0.2948756
## 73 74 75 76 77 78 79 80
## 0.2948756 0.2950577 0.2952397 0.2954218 0.2956039 0.2957859 0.2957859 0.2956039
## 81 82 83 84 85 86 87 88
## 0.2959680 0.2959680 0.2956039 0.2957859 0.2959680 0.2961501 0.2961501 0.2961501
## 89 90 91 92 93 94 95 96
## 0.2961501 0.2961501 0.2959680 0.2959680 0.2961501 0.2961501 0.2961501 0.2961501
## 97 98 99 100 101 102 103 104
## 0.2959680 0.2957859 0.2957859 0.2961501 0.2965142 0.2968784 0.2966963 0.2956039
## 105 106 107 108 109 110 111 112
## 0.2943294 0.2932370 0.2897777 0.2865004 0.2855901 0.2870466 0.2894135 0.2901418
## 113 114 115 116 117 118 119 120
## 0.2885032 0.2850438 0.2783073 0.2693859 0.2633776 0.2606465 0.2579155 0.2577334
## 121 122 123 124 125 126 127 128
## 0.2606465 0.2633776 0.2648341 0.2661086 0.2662907 0.2664727 0.2666548 0.2662907
## 129 130 131 132 133 134 135 136
## 0.2642879 0.2610107 0.2568231 0.2533637 0.2528175 0.2540920 0.2551844 0.2555486
## 137 138 139 140 141 142 143 144
## 0.2562768 0.2577334 0.2590079 0.2597362 0.2611927 0.2637417 0.2650162 0.2644700
## 145 146 147 148 149 150 151 152
## 0.2637417 0.2608286 0.2570051 0.2539099 0.2522713 0.2535458 0.2573693 0.2617389
## 153 154 155 156 157 158 159 160
## 0.2644700 0.2666548 0.2688397 0.2708424 0.2717528 0.2721169 0.2722990 0.2721169
## 161 162 163 164 165 166 167 168
## 0.2708424 0.2686576 0.2673831 0.2672010 0.2670190 0.2666548 0.2664727 0.2668369
## 169 170 171 172 173 174 175 176
## 0.2670190 0.2661086 0.2642879 0.2622851 0.2599182 0.2571872 0.2560948 0.2584617
## 177 178 179 180 181 182 183 184
## 0.2624672 0.2650162 0.2664727 0.2684755 0.2730273 0.2810383 0.2832231 0.2852259
## 185 186 187 188 189 190 191 192
## 0.2844976 0.2792176 0.2730273 0.2697500 0.2675652 0.2655624 0.2637417 0.2630134
## 193 194 195 196 197 198 199 200
## 0.2641058 0.2637417 0.2644700 0.2682934 0.2688397 0.2688397 0.2710245 0.2775790
## 201 202 203 204 205 206 207 208
## 0.2874107 0.2952397 0.2974246 0.2996094 0.3016122 0.3005198 0.2965142 0.2899597
## 209 210 211 212 213 214 215 216
## 0.2823128 0.2753942 0.2712066 0.2688397 0.2677472 0.2675652 0.2675652 0.2673831
## 217 218 219 220 221 222 223 224
## 0.2673831 0.2672010 0.2666548 0.2655624 0.2641058 0.2626493 0.2624672 0.2621031
## 225 226 227 228 229 230 231 232
## 0.2630134 0.2659265 0.2712066 0.2773969 0.2817666 0.2824949 0.2797638 0.2761224
## 233 234 235 236 237 238 239 240
## 0.2726631 0.2701141 0.2681114 0.2670190 0.2659265 0.2646520 0.2637417 0.2637417
## 241 242 243 244 245 246 247 248
## 0.2641058 0.2642879 0.2642879 0.2637417 0.2630134 0.2626493 0.2622851 0.2619210
## 249 250 251 252 253 254 255 256
## 0.2615569 0.2615569 0.2621031 0.2628314 0.2633776 0.2631955 0.2630134 0.2624672
## 257 258 259 260 261 262 263 264
## 0.2617389 0.2611927 0.2610107 0.2601003 0.2577334 0.2571872 0.2568231 0.2575513
## 265 266 267 268 269 270 271 272
## 0.2573693 0.2570051 0.2571872 0.2577334 0.2584617 0.2593720 0.2601003 0.2610107
## 273 274 275 276 277 278 279 280
## 0.2615569 0.2619210 0.2619210 0.2617389 0.2622851 0.2653803 0.2712066 0.2773969
## 281 282 283 284 285 286 287 288
## 0.2792176 0.2808562 0.2821307 0.2837694 0.2861363 0.2886852 0.2914163 0.2941473
## 289 290 291 292 293 294 295 296
## 0.2966963 0.2988811 0.3005198 0.3017942 0.3028867 0.3039791 0.3048894 0.3057998
## 297 298 299 300 301 302 303 304
## 0.3065281 0.3068922 0.3068922 0.3067101 0.3061639 0.3056177 0.3034329 0.2994273
## 305 306 307 308 309 310 311 312
## 0.2936011 0.2848618 0.2744838 0.2662907 0.2626493 0.2613748 0.2626493 0.2648341
## 313 314 315 316 317 318 319 320
## 0.2664727 0.2690217 0.2715707 0.2728452 0.2733914 0.2735735 0.2735735 0.2741197
## 321 322 323 324 325 326 327 328
## 0.2748479 0.2752121 0.2752121 0.2748479 0.2743017 0.2743017 0.2746659 0.2750300
## 329 330 331 332 333 334 335 336
## 0.2753942 0.2755762 0.2759404 0.2763045 0.2761224 0.2744838 0.2724810 0.2719348
## 337 338 339 340 341 342 343 344
## 0.2719348 0.2704783 0.2693859 0.2701141 0.2712066 0.2719348 0.2721169 0.2715707
## 345 346 347 348 349 350 351 352
## 0.2704783 0.2692038 0.2684755 0.2686576 0.2688397 0.2684755 0.2686576 0.2710245
## 353 354 355 356 357 358 359 360
## 0.2732093 0.2746659 0.2744838 0.2715707 0.2668369 0.2619210 0.2575513 0.2559127
## 361 362 363 364 365 366 367 368
## 0.2542741 0.2531817 0.2542741 0.2571872 0.2597362 0.2608286 0.2611927 0.2615569
## 369 370 371 372 373 374 375 376
## 0.2621031 0.2635596 0.2653803 0.2668369 0.2684755 0.2692038 0.2693859 0.2692038
## 377 378 379 380 381 382 383 384
## 0.2686576 0.2681114 0.2673831 0.2662907 0.2657445 0.2651983 0.2650162 0.2659265
## 385 386 387 388 389 390 391 392
## 0.2661086 0.2653803 0.2650162 0.2642879 0.2635596 0.2631955 0.2633776 0.2637417
## 393 394 395 396 397 398 399 400
## 0.2639238 0.2641058 0.2637417 0.2628314 0.2615569 0.2611927 0.2617389 0.2621031
## 401 402 403 404 405 406 407 408
## 0.2622851 0.2624672 0.2633776 0.2648341 0.2653803 0.2650162 0.2648341 0.2653803
## 409 410 411 412 413 414 415 416
## 0.2666548 0.2672010 0.2673831 0.2679293 0.2688397 0.2695679 0.2706603 0.2713886
## 417 418 419 420 421 422 423 424
## 0.2712066 0.2717528 0.2717528 0.2710245 0.2693859 0.2672010 0.2664727 0.2662907
## 425 426 427 428 429 430 431 432
## 0.2644700 0.2631955 0.2621031 0.2599182 0.2571872 0.2550024 0.2535458 0.2526355
## 433 434 435 436 437 438 439 440
## 0.2519072 0.2515430 0.2515430 0.2515430 0.2504506 0.2488120 0.2480837 0.2480837
## 441 442 443 444 445 446 447 448
## 0.2455347 0.2402547 0.2377058 0.2377058 0.2366133 0.2364313 0.2362492 0.2364313
## 449 450 451 452 453 454 455 456
## 0.2377058 0.2397085 0.2420754 0.2446244 0.2477196 0.2499044 0.2515430 0.2533637
## 457 458 459 460 461 462 463 464
## 0.2559127 0.2582796 0.2588258 0.2588258 0.2584617 0.2586438 0.2582796 0.2573693
## 465 466 467 468 469 470 471 472
## 0.2568231 0.2570051 0.2575513 0.2580975 0.2590079 0.2606465 0.2622851 0.2635596
## 473 474 475 476 477 478 479 480
## 0.2641058 0.2646520 0.2653803 0.2659265 0.2657445 0.2653803 0.2648341 0.2633776
## 481 482 483 484 485 486 487 488
## 0.2617389 0.2610107 0.2608286 0.2611927 0.2615569 0.2615569 0.2615569 0.2615569
## 489 490 491 492 493 494 495 496
## 0.2613748 0.2615569 0.2615569 0.2624672 0.2619210 0.2619210 0.2617389 0.2619210
## 497 498 499 500 501 502 503 504
## 0.2628314 0.2648341 0.2681114 0.2708424 0.2724810 0.2728452 0.2732093 0.2733914
## 505 506 507 508 509 510 511 512
## 0.2739376 0.2746659 0.2753942 0.2763045 0.2770328 0.2770328 0.2759404 0.2737555
## 513 514 515 516 517 518 519 520
## 0.2713886 0.2697500 0.2692038 0.2690217 0.2688397 0.2675652 0.2662907 0.2653803
## 521 522 523 524 525 526 527 528
## 0.2650162 0.2651983 0.2655624 0.2657445 0.2648341 0.2639238 0.2633776 0.2635596
## 529 530 531 532 533 534 535 536
## 0.2641058 0.2641058 0.2639238 0.2641058 0.2641058 0.2639238 0.2639238 0.2633776
## 537 538 539 540 541 542 543 544
## 0.2622851 0.2611927 0.2604644 0.2602824 0.2608286 0.2626493 0.2646520 0.2653803
## 545 546 547 548 549 550 551 552
## 0.2650162 0.2646520 0.2644700 0.2641058 0.2639238 0.2635596 0.2633776 0.2628314
## 553 554 555 556 557 558 559 560
## 0.2624672 0.2628314 0.2626493 0.2615569 0.2602824 0.2588258 0.2580975 0.2579155
## 561 562 563 564 565 566 567 568
## 0.2579155 0.2579155 0.2579155 0.2580975 0.2584617 0.2588258 0.2591900 0.2595541
## 569 570 571 572 573 574 575 576
## 0.2599182 0.2608286 0.2613748 0.2615569 0.2608286 0.2601003 0.2601003 0.2604644
## 577 578 579 580 581 582 583 584
## 0.2602824 0.2599182 0.2593720 0.2586438 0.2582796 0.2582796 0.2584617 0.2580975
## 585 586 587 588 589 590 591 592
## 0.2582796 0.2586438 0.2590079 0.2599182 0.2601003 0.2599182 0.2602824 0.2611927
## 593 594 595 596 597 598 599 600
## 0.2630134 0.2661086 0.2688397 0.2688397 0.2662907 0.2633776 0.2622851 0.2622851
## 601 602 603 604 605 606 607 608
## 0.2617389 0.2604644 0.2593720 0.2588258 0.2590079 0.2595541 0.2597362 0.2601003
## 609 610 611 612 613 614 615 616
## 0.2599182 0.2595541 0.2586438 0.2570051 0.2557306 0.2557306 0.2560948 0.2566410
## 617 618 619 620 621 622 623 624
## 0.2570051 0.2570051 0.2564589 0.2555486 0.2550024 0.2550024 0.2553665 0.2557306
## 625 626 627 628 629 630 631 632
## 0.2564589 0.2573693 0.2586438 0.2602824 0.2619210 0.2628314 0.2635596 0.2635596
## 633 634 635 636 637 638 639 640
## 0.2630134 0.2624672 0.2622851 0.2624672 0.2622851 0.2621031 0.2622851 0.2624672
## 641 642 643 644 645 646 647 648
## 0.2622851 0.2617389 0.2615569 0.2615569 0.2619210 0.2621031 0.2619210 0.2615569
## 649 650 651 652 653 654 655 656
## 0.2610107 0.2602824 0.2595541 0.2593720 0.2595541 0.2602824 0.2608286 0.2610107
## 657 658 659 660 661 662 663 664
## 0.2608286 0.2610107 0.2611927 0.2617389 0.2619210 0.2617389 0.2621031 0.2630134
## 665 666 667 668 669 670 671 672
## 0.2637417 0.2642879 0.2646520 0.2644700 0.2637417 0.2633776 0.2630134 0.2621031
## 673 674 675 676 677 678 679 680
## 0.2617389 0.2617389 0.2617389 0.2617389 0.2615569 0.2611927 0.2608286 0.2604644
## 681 682 683 684 685 686 687 688
## 0.2602824 0.2601003 0.2601003 0.2602824 0.2606465 0.2611927 0.2615569 0.2615569
## 689 690 691 692 693 694 695 696
## 0.2615569 0.2615569 0.2613748 0.2604644 0.2593720 0.2586438 0.2582796 0.2579155
## 697 698 699 700 701 702 703 704
## 0.2575513 0.2573693 0.2577334 0.2580975 0.2580975 0.2580975 0.2584617 0.2593720
## 705 706 707 708 709 710 711 712
## 0.2606465 0.2615569 0.2630134 0.2644700 0.2641058 0.2631955 0.2630134 0.2635596
## 713 714 715 716 717 718 719 720
## 0.2635596 0.2633776 0.2628314 0.2630134 0.2631955 0.2628314 0.2624672 0.2630134
## 721 722 723 724 725 726 727 728
## 0.2641058 0.2651983 0.2657445 0.2657445 0.2650162 0.2639238 0.2624672 0.2619210
## 729 730 731 732 733 734 735 736
## 0.2626493 0.2635596 0.2641058 0.2641058 0.2641058 0.2642879 0.2648341 0.2650162
## 737 738 739 740 741 742 743 744
## 0.2644700 0.2644700 0.2679293 0.2743017 0.2788535 0.2792176 0.2753942 0.2695679
## 745 746 747 748 749 750 751 752
## 0.2677472 0.2672010 0.2679293 0.2692038 0.2686576 0.2684755 0.2686576 0.2690217
## 753 754 755 756 757 758 759 760
## 0.2695679 0.2701141 0.2708424 0.2733914 0.2755762 0.2761224 0.2741197 0.2710245
## 761 762 763 764 765 766 767 768
## 0.2686576 0.2675652 0.2673831 0.2675652 0.2675652 0.2670190 0.2661086 0.2657445
## 769 770 771 772 773 774 775 776
## 0.2672010 0.2677472 0.2679293 0.2682934 0.2686576 0.2686576 0.2688397 0.2684755
## 777 778 779 780 781 782 783 784
## 0.2672010 0.2659265 0.2659265 0.2670190 0.2692038 0.2708424 0.2706603 0.2688397
## 785 786 787 788 789 790 791 792
## 0.2673831 0.2675652 0.2679293 0.2675652 0.2673831 0.2677472 0.2681114 0.2692038
## 793 794 795 796 797 798 799 800
## 0.2692038 0.2677472 0.2659265 0.2644700 0.2637417 0.2642879 0.2686576 0.2695679
## 801 802 803 804 805 806 807 808
## 0.2697500 0.2702962 0.2717528 0.2726631 0.2708424 0.2702962 0.2695679 0.2693859
## 809 810 811 812 813 814 815 816
## 0.2695679 0.2690217 0.2684755 0.2686576 0.2697500 0.2717528 0.2706603 0.2681114
## 817 818 819
## 0.2653803 0.2630134 0.2613748