Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- inputMat <- as.matrix(read.table(textConnection("
- 4.268063 4 5 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 4.835864 4 5 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 6.090391 6 7 1 10
- 6.926898 6 7 1 10
- 2.757738 2 3 1 10
- 6.753170 6 7 1 10
- 6.632133 6 7 1 10
- 7.240152 7 8 1 10
- 3.888131 3 4 1 10
- 4.875474 4 5 1 10
- 4.988239 4 5 1 10
- 3.997805 3 4 1 10
- 5.783397 5 6 1 10
- Inf 20 Inf 4 10
- 6.422350 6 7 1 10
- Inf 20 Inf 4 10
- 6.284683 6 7 1 10
- 5.869865 5 6 1 10
- 4.301658 4 5 1 10
- 6.089626 6 7 1 10
- 4.898094 4 5 1 10
- Inf 20 Inf 4 10
- 6.246641 6 7 1 10
- 6.799316 6 7 1 10
- 5.039809 5 6 1 10
- 8.137007 8 9 1 10
- Inf 20 Inf 4 10
- 7.163037 7 8 1 10
- 5.233179 5 6 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 4.934695 4 5 1 10
- 3.862736 3 4 1 10
- 4.686978 4 5 1 10
- 5.787939 5 6 1 10
- 7.943111 7 8 1 10
- Inf 20 Inf 4 10
- 8.297172 8 9 1 10
- 4.621190 4 5 1 10
- Inf 20 Inf 4 10
- 2.784991 2 3 1 10
- 6.684080 6 7 1 10
- 4.837746 4 5 1 10
- Inf 20 Inf 4 10
- 4.834292 4 5 1 10
- 6.685557 6 7 1 10
- 4.181680 4 5 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 6.825699 6 7 1 10
- 5.549326 5 6 1 10
- 4.365653 4 5 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 6.502134 6 7 1 10
- 4.847652 4 5 1 10
- 8.350025 8 9 1 10
- 6.899775 6 7 1 10
- Inf 20 Inf 4 10
- 7.965559 7 8 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 6.496165 6 7 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 5.844762 5 6 1 10
- Inf 20 Inf 4 10
- 3.651155 3 4 1 10
- 5.086461 5 6 1 10
- 6.046691 6 7 1 10
- 7.058280 7 8 1 10
- Inf 20 Inf 4 10
- 4.375909 4 5 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 3.178064 3 4 1 10
- 6.513465 6 7 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 7.614067 7 8 1 10
- 10.506601 10 11 1 10
- 10.393137 10 11 1 10
- 4.841498 4 5 1 10
- 7.306547 7 8 1 10
- 7.233084 7 8 1 10
- Inf 20 Inf 4 10
- 4.925290 4 5 1 10
- 7.011571 7 8 1 10
- 9.928764 9 10 1 10
- Inf 20 Inf 4 10
- 6.671448 6 7 1 10
- Inf 20 Inf 4 10
- Inf 20 Inf 4 10
- 4.493113 4 5 1 10
- 5.974026 5 6 1 16
- 9.422095 9 10 1 16
- 5.288058 5 6 1 16
- 9.766532 9 10 1 16
- 6.841781 6 7 1 16
- Inf 20 Inf 4 16
- 6.473939 6 7 1 16
- 6.608567 6 7 1 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- 8.712419 8 9 1 16
- 8.492222 8 9 1 16
- 11.626487 11 12 1 16
- 8.239599 8 9 1 16
- 7.076040 7 8 1 16
- 6.826600 6 7 1 16
- 4.417569 4 5 1 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- 6.018623 6 7 1 16
- 8.945612 8 9 1 16
- 7.030062 7 8 1 16
- Inf 20 Inf 4 16
- 7.561805 7 8 1 16
- 17.295557 17 18 1 16
- 7.299832 7 8 1 16
- 8.750735 8 9 1 16
- 5.530297 5 6 1 16
- 4.493255 4 5 1 16
- 7.962182 7 8 1 16
- Inf 20 Inf 4 16
- 4.073075 4 5 1 16
- 12.590208 12 13 1 16
- 5.810693 5 6 1 16
- 9.173796 9 10 1 16
- 12.656002 12 13 1 16
- 10.541377 10 11 1 16
- Inf 20 Inf 4 16
- 5.807266 5 6 1 16
- 5.216672 5 6 1 16
- 8.350314 8 9 1 16
- 2.912009 2 3 1 16
- 11.072793 11 12 1 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- 8.067160 8 9 1 16
- 7.866341 7 8 1 16
- Inf 20 Inf 4 16
- 10.543505 10 11 1 16
- 6.328423 6 7 1 16
- 4.235567 4 5 1 16
- Inf 20 Inf 4 16
- 6.308195 6 7 1 16
- 4.989120 4 5 1 16
- 7.017843 7 8 1 16
- 4.940731 4 5 1 16
- 6.568252 6 7 1 16
- 4.713135 4 5 1 16
- 3.668359 3 4 1 16
- 5.927443 5 6 1 16
- 9.645991 9 10 1 16
- Inf 20 Inf 4 16
- 5.497114 5 6 1 16
- Inf 20 Inf 4 16
- 8.330205 8 9 1 16
- Inf 20 Inf 4 16
- Inf 20 Inf 4 16
- 9.686394 9 10 1 16
- 8.399837 8 9 1 16
- 4.835181 4 5 1 16
- 8.511252 8 9 1 16
- 10.483904 10 11 1 16
- 6.081267 6 7 1 16
- 4.973101 4 5 1 16
- 11.152892 11 12 1 16
- 6.602480 6 7 1 16
- 5.118556 5 6 1 16
- 7.453080 7 8 1 16
- 25.512731 20 Inf 0 16
- Inf 20 Inf 4 16
- 5.214515 5 6 1 16
- Inf 20 Inf 4 16
- 8.809781 8 9 1 16
- 7.113735 7 8 1 16
- Inf 20 Inf 4 16
- 5.626853 5 6 1 16
- Inf 20 Inf 4 16
- 5.494501 5 6 1 16
- 6.816956 6 7 1 16
- Inf 20 Inf 4 16
- 6.924974 6 7 1 16
- 3.589199 3 4 1 16
- 8.248053 8 9 1 16
- 6.077767 6 7 1 16
- Inf 20 Inf 4 16
- 12.749272 12 13 1 16
- Inf 20 Inf 4 22
- 8.145331 8 9 1 22
- 7.018483 7 8 1 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- 4.689269 4 5 1 22
- 16.291933 16 17 1 22
- 6.835840 6 7 1 22
- 5.746658 5 6 1 22
- 5.392693 5 6 1 22
- 5.373915 5 6 1 22
- 7.563032 7 8 1 22
- 6.450659 6 7 1 22
- 7.855748 7 8 1 22
- 5.518934 5 6 1 22
- 7.413918 7 8 1 22
- Inf 20 Inf 4 22
- 4.347408 4 5 1 22
- 6.326859 6 7 1 22
- 4.247252 4 5 1 22
- Inf 20 Inf 4 22
- 4.233826 4 5 1 22
- 5.114650 5 6 1 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- 7.610498 7 8 1 22
- 8.848758 8 9 1 22
- 5.839719 5 6 1 22
- 8.664659 8 9 1 22
- 9.092371 9 10 1 22
- 7.155833 7 8 1 22
- Inf 20 Inf 4 22
- 7.312117 7 8 1 22
- 9.412312 9 10 1 22
- 4.259247 4 5 1 22
- 6.629548 6 7 1 22
- 5.819321 5 6 1 22
- Inf 20 Inf 4 22
- 6.912066 6 7 1 22
- 3.599371 3 4 1 22
- Inf 20 Inf 4 22
- 8.539910 8 9 1 22
- 3.439026 3 4 1 22
- 10.755944 10 11 1 22
- 7.131071 7 8 1 22
- 3.734676 3 4 1 22
- 7.460057 7 8 1 22
- 11.043886 11 12 1 22
- 6.166909 6 7 1 22
- 14.270150 14 15 1 22
- 6.502952 6 7 1 22
- Inf 20 Inf 4 22
- 7.863157 7 8 1 22
- 10.357448 10 11 1 22
- 9.210890 9 10 1 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- 5.389687 5 6 1 22
- 5.298091 5 6 1 22
- 7.573046 7 8 1 22
- 8.900940 8 9 1 22
- 8.006156 8 9 1 22
- Inf 20 Inf 4 22
- 9.686941 9 10 1 22
- 6.884800 6 7 1 22
- 10.327921 10 11 1 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- 7.914521 7 8 1 22
- 3.695012 3 4 1 22
- 5.467132 5 6 1 22
- Inf 20 Inf 4 22
- 7.840289 7 8 1 22
- 6.560040 6 7 1 22
- Inf 20 Inf 4 22
- 4.400657 4 5 1 22
- 6.855404 6 7 1 22
- 10.140989 10 11 1 22
- Inf 20 Inf 4 22
- 7.600298 7 8 1 22
- 6.428420 6 7 1 22
- 8.586318 8 9 1 22
- 9.911352 9 10 1 22
- 10.944423 10 11 1 22
- 6.684860 6 7 1 22
- 9.088183 9 10 1 22
- 9.765874 9 10 1 22
- 5.112833 5 6 1 22
- 6.890751 6 7 1 22
- 6.753094 6 7 1 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- Inf 20 Inf 4 22
- 11.343692 11 12 1 22
- Inf 20 Inf 4 22
- 6.732804 6 7 1 22
- 5.118488 5 6 1 22
- 6.445765 6 7 1 22
- 8.453932 8 9 1 28
- Inf 20 Inf 4 28
- 7.611811 7 8 1 28
- 7.148196 7 8 1 28
- 5.690943 5 6 1 28
- 5.883728 5 6 1 28
- Inf 20 Inf 4 28
- 10.483547 10 11 1 28
- 5.955929 5 6 1 28
- 4.186291 4 5 1 28
- 11.521526 11 12 1 28
- 8.009429 8 9 1 28
- 7.655484 7 8 1 28
- Inf 20 Inf 4 28
- 9.991222 9 10 1 28
- 7.172954 7 8 1 28
- Inf 20 Inf 4 28
- 8.935009 8 9 1 28
- Inf 20 Inf 4 28
- Inf 20 Inf 4 28
- 10.318172 10 11 1 28
- 8.680718 8 9 1 28
- 8.086856 8 9 1 28
- 10.837159 10 11 1 28
- 7.156762 7 8 1 28
- 6.001851 6 7 1 28
- 8.922795 8 9 1 28
- 7.860494 7 8 1 28
- 6.308847 6 7 1 28
- 10.353015 10 11 1 28
- Inf 20 Inf 4 28
- 6.297438 6 7 1 28
- Inf 20 Inf 4 28
- 14.801072 14 15 1 28
- 9.601222 9 10 1 28
- 6.458924 6 7 1 28
- 6.742558 6 7 1 28
- 6.844412 6 7 1 28
- Inf 20 Inf 4 28
- 5.455491 5 6 1 28
- Inf 20 Inf 4 28
- 6.402477 6 7 1 28
- 2.726490 2 3 1 28
- 9.490706 9 10 1 28
- 9.752210 9 10 1 28
- Inf 20 Inf 4 28
- 11.661637 11 12 1 28
- 8.057321 8 9 1 28
- 7.209312 7 8 1 28
- 8.300370 8 9 1 28
- 6.078291 6 7 1 28
- 9.043891 9 10 1 28
- 5.721226 5 6 1 28
- 4.233122 4 5 1 28
- Inf 20 Inf 4 28
- 11.172754 11 12 1 28
- 5.819324 5 6 1 28
- 5.051789 5 6 1 28
- Inf 20 Inf 4 28
- 8.709141 8 9 1 28
- 10.876180 10 11 1 28
- 10.207926 10 11 1 28
- 7.286342 7 8 1 28
- 9.643852 9 10 1 28
- 7.909810 7 8 1 28
- 10.759215 10 11 1 28
- Inf 20 Inf 4 28
- 4.366978 4 5 1 28
- 7.092675 7 8 1 28
- 4.382632 4 5 1 28
- 10.566793 10 11 1 28
- 7.042897 7 8 1 28
- 10.648184 10 11 1 28
- Inf 20 Inf 4 28
- 8.493182 8 9 1 28
- 5.336394 5 6 1 28
- 5.330173 5 6 1 28
- 10.830822 10 11 1 28
- 13.201346 13 14 1 28
- Inf 20 Inf 4 28
- 15.328194 15 16 1 28
- 7.212677 7 8 1 28
- 6.440357 6 7 1 28
- 7.833290 7 8 1 28
- Inf 20 Inf 4 28
- 5.532838 5 6 1 28
- Inf 20 Inf 4 28
- 5.776575 5 6 1 28
- 7.013364 7 8 1 28
- 8.430898 8 9 1 28
- 9.576473 9 10 1 28
- 4.758495 4 5 1 28
- 14.082692 14 15 1 28
- 6.052728 6 7 1 28
- 6.120784 6 7 1 28
- 6.832399 6 7 1 28
- 10.047205 10 11 1 28
- Inf 20 Inf 4 28
- 6.282310 6 7 1 28
- 9.192467 9 10 1 28
- Inf 20 Inf 4 34
- 8.258754 8 9 1 34
- 8.849809 8 9 1 34
- 9.778415 9 10 1 34
- 9.424452 9 10 1 34
- Inf 20 Inf 4 34
- 12.167809 12 13 1 34
- 7.230517 7 8 1 34
- 11.645215 11 12 1 34
- 10.852363 10 11 1 34
- Inf 20 Inf 4 34
- Inf 20 Inf 4 34
- 9.383850 9 10 1 34
- 7.129181 7 8 1 34
- Inf 20 Inf 4 34
- 9.919855 9 10 1 34
- 8.460769 8 9 1 34
- 8.609479 8 9 1 34
- 6.989502 6 7 1 34
- 4.788868 4 5 1 34
- 5.921015 5 6 1 34
- 5.076072 5 6 1 34
- 15.193696 15 16 1 34
- Inf 20 Inf 4 34
- Inf 20 Inf 4 34
- 9.723276 9 10 1 34
- 4.786926 4 5 1 34
- Inf 20 Inf 4 34
- 14.604087 14 15 1 34
- 4.804427 4 5 1 34
- 6.469944 6 7 1 34
- 9.856422 9 10 1 34
- 7.347468 7 8 1 34
- 9.711560 9 10 1 34
- 9.456971 9 10 1 34
- 11.158517 11 12 1 34
- Inf 20 Inf 4 34
- 9.050926 9 10 1 34
- 8.850817 8 9 1 34
- 6.129082 6 7 1 34
- 11.496230 11 12 1 34
- 8.636476 8 9 1 34
- 6.310575 6 7 1 34
- 8.790823 8 9 1 34
- 8.567049 8 9 1 34
- 8.187000 8 9 1 34
- Inf 20 Inf 4 34
- 4.597418 4 5 1 34
- 11.688279 11 12 1 34
- 7.131484 7 8 1 34
- 18.364503 18 19 1 34
- 4.875456 4 5 1 34
- 7.269105 7 8 1 34
- Inf 20 Inf 4 34
- 6.761910 6 7 1 34
- 10.835314 10 11 1 34
- 5.410301 5 6 1 34
- 5.843082 5 6 1 34
- 9.866109 9 10 1 34
- 3.446969 3 4 1 34
- 16.965309 16 17 1 34
- 7.846356 7 8 1 34
- 7.515407 7 8 1 34
- 8.438184 8 9 1 34
- 15.361471 15 16 1 34
- Inf 20 Inf 4 34
- 6.191788 6 7 1 34
- 7.565918 7 8 1 34
- 11.367279 11 12 1 34
- 6.190009 6 7 1 34
- 9.063369 9 10 1 34
- 13.982467 13 14 1 34
- 11.730461 11 12 1 34
- Inf 20 Inf 4 34
- 9.248405 9 10 1 34
- 6.996990 6 7 1 34
- 10.955868 10 11 1 34
- 6.577806 6 7 1 34
- 11.513186 11 12 1 34
- 7.882403 7 8 1 34
- 7.696856 7 8 1 34
- 11.102816 11 12 1 34
- 5.333107 5 6 1 34
- 6.758993 6 7 1 34
- Inf 20 Inf 4 34
- 8.255823 8 9 1 34
- 8.196053 8 9 1 34
- 13.925059 13 14 1 34
- 7.464187 7 8 1 34
- 5.908341 5 6 1 34
- 7.709907 7 8 1 34
- 11.052416 11 12 1 34
- Inf 20 Inf 4 34
- 16.801992 16 17 1 34
- 8.781210 8 9 1 34
- 11.534025 11 12 1 34
- 3.276015 3 4 1 34
- 12.553650 12 13 1 34
- 7.302265 7 8 1 34
- 12.735121 12 13 1 34
- Inf 20 Inf 4 40
- Inf 20 Inf 4 40
- Inf 20 Inf 4 40
- 10.976337 10 11 1 40
- 8.286781 8 9 1 40
- 6.304428 6 7 1 40
- 6.252275 6 7 1 40
- 14.204029 14 15 1 40
- 15.753046 15 16 1 40
- Inf 20 Inf 4 40
- 9.313259 9 10 1 40
- 7.526972 7 8 1 40
- 17.843488 17 18 1 40
- 13.163388 13 14 1 40
- Inf 20 Inf 4 40
- 9.425878 9 10 1 40
- 9.171105 9 10 1 40
- Inf 20 Inf 4 40
- 7.309866 7 8 1 40
- 7.230562 7 8 1 40
- 8.786555 8 9 1 40
- 4.558044 4 5 1 40
- Inf 20 Inf 4 40
- 5.562887 5 6 1 40
- 8.683998 8 9 1 40
- Inf 20 Inf 4 40
- 7.769392 7 8 1 40
- 6.025881 6 7 1 40
- Inf 20 Inf 4 40
- 16.385345 16 17 1 40
- 10.658288 10 11 1 40
- 12.029552 12 13 1 40
- 6.867404 6 7 1 40
- 14.717769 14 15 1 40
- 7.038278 7 8 1 40
- 5.232182 5 6 1 40
- 8.087490 8 9 1 40
- 7.226209 7 8 1 40
- 7.488762 7 8 1 40
- 7.542690 7 8 1 40
- 6.763303 6 7 1 40
- 10.540067 10 11 1 40
- 5.791338 5 6 1 40
- 13.673177 13 14 1 40
- 6.602495 6 7 1 40
- 8.930889 8 9 1 40
- 11.227102 11 12 1 40
- 12.301289 12 13 1 40
- 8.671360 8 9 1 40
- 7.689139 7 8 1 40
- 17.134939 17 18 1 40
- 8.831534 8 9 1 40
- 10.654551 10 11 1 40
- 14.065886 14 15 1 40
- 8.182912 8 9 1 40
- 11.333733 11 12 1 40
- Inf 20 Inf 4 40
- 5.041513 5 6 1 40
- 3.647567 3 4 1 40
- Inf 20 Inf 4 40
- 8.707224 8 9 1 40
- 4.336742 4 5 1 40
- 11.928365 11 12 1 40
- 4.544854 4 5 1 40
- 22.131423 20 Inf 0 40
- 9.164520 9 10 1 40
- 17.255478 17 18 1 40
- 7.684772 7 8 1 40
- 11.322955 11 12 1 40
- 15.313612 15 16 1 40
- 15.352947 15 16 1 40
- 4.502466 4 5 1 40
- 6.812818 6 7 1 40
- 9.947180 9 10 1 40
- Inf 20 Inf 4 40
- 12.027548 12 13 1 40
- 6.352271 6 7 1 40
- 17.328514 17 18 1 40
- 8.712545 8 9 1 40
- 11.901630 11 12 1 40
- 10.504628 10 11 1 40
- Inf 20 Inf 4 40
- 13.647926 13 14 1 40
- 26.793991 20 Inf 0 40
- 15.229589 15 16 1 40
- 16.711487 16 17 1 40
- 8.036608 8 9 1 40
- Inf 20 Inf 4 40
- 7.597567 7 8 1 40
- 7.049394 7 8 1 40
- Inf 20 Inf 4 40
- Inf 20 Inf 4 40
- 8.310225 8 9 1 40
- 6.167071 6 7 1 40
- 13.942419 13 14 1 40
- 19.600884 19 20 1 40
- 7.356598 7 8 1 40
- Inf 20 Inf 4 40
- 7.841882 7 8 1 40
- 8.174693 8 9 1 40")))
- objFunc <- function(x, t) {
- p <- exp(-exp(x[1] + x[2] * (t[, 5] / 100)))
- alpha <- x[3] + x[4] * (t[, 5] / 10)
- beta <- rep(x[5], nrow(t))
- gFunc <- function(t, a, b) 1 / (1 + (t / a)^b)
- funcList <- function(t){
- return(switch(as.character(t),
- "0" = function(p, t, alpha, beta) -sum(log(p + (1-p) * gFunc(t[, 2], alpha, beta))),
- "1" = function(p, t, alpha, beta) -sum(log((1-p) * (gFunc(t[, 2], alpha, beta) - gFunc(t[, 3], alpha, beta)))),
- "2" = function(p, t, alpha, beta) -sum(log(p + (1-p) * (1 - gFunc(t[, 3], alpha, beta)))),
- "3" = function(p, t, alpha, beta) -sum(log((1-p) * (beta / alpha) * (t[, 1] / alpha) ^ (beta - 1) * (gFunc(t[, 1], alpha, beta)) ^ 2)),
- function(p, t, alpha, beta) -sum(log(p + (1-p) * gFunc(t[ , 2], alpha, beta)))))
- }
- sum(sapply(sort(unique(t[,4])), function(i) funcList(i)(p[t[, 4]==i], matrix(t[t[, 4]==i, ], nrow = sum(t[, 4]==i)), alpha[t[, 4]==i], beta[t[, 4]==i])))
- }
- # test
- objFunc(c(0.1, 1, 5, 1, 5), inputMat) # 1399.626
- R <- 1000L
- iter <- 1L
- bootRes <- matrix(NA, R, 5)
- repeat{
- idx <- sample(nrow(inputMat), nrow(inputMat), TRUE)
- optimRes <- optim(c(0.1, 1, 5, 1, 5), function(x) objFunc(x, inputMat[idx, ]),
- control = list(maxit = 5e4L))
- if (optimRes$conve == 0) {
- bootRes[iter, ] <- optimRes$par
- iter <- iter + 1L
- }
- if (iter > R)
- break
- }
Advertisement
Add Comment
Please, Sign In to add comment
