Skip to content

Instantly share code, notes, and snippets.

@LeaveNhA

LeaveNhA/UIST624 - GLM, HW.ipynb

Created June 9, 2022 18:15
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save LeaveNhA/9870bcd4b8b2a6e0957d0c42365d5275 to your computer and use it in GitHub Desktop.
Save LeaveNhA/9870bcd4b8b2a6e0957d0c42365d5275 to your computer and use it in GitHub Desktop.
Download ZIP
Homework, postgrad UIST624, OMU - Data Science
Display the source blob
Display the rendered blob
Raw
UIST624 - GLM, HW.ipynb
{
"cells": [
{
"cell_type": "markdown",
"id": "b7008b40",
"metadata": {},
"source": [
"## UIST624 HOMEWORK (NAIVE) IMPLEMENTATION\n",
"\n",
"### Q1-)\n",
"\n",
"$$\n",
"NORMAL GLM(x)=\n",
"\\begin{cases}\n",
"y \\text{~} N(\\mu, \\phi/\\omega), & \\quad \\text{when} \\quad \\mu = \\eta, \\quad V(\\mu) = 1, \\quad \\phi = \\sigma^{2} \\\\\n",
"\\mu = \\eta = \\beta_{0} + \\beta_{1} x_{1} + \\cdots + \\beta_{p} x_{p}, & \\quad w_{i} = 1, \\quad i = [1..n]\n",
"\\end{cases}\n",
"$$\n",
"\n",
"#### Score function:\n",
"\n",
"$$\n",
"\\upsilon(\\beta_{j}) = \\frac{1}{\\phi} \\sum_{i=1}^{n} W_{i} \\frac{\\partial \\eta_{i}}{\\partial \\mu_{i}} (y_{i} - \\mu_{i}) x_{ji}\n",
"$$\n",
"\n",
"#### Weights:\n",
"$$\n",
"W_{i} = \\frac{w_{i}}{V(\\mu_{i})(\\frac{\\partial \\eta_{i}}{\\partial \\mu_{i}})^2}\n",
"$$\n",
"\n",
"#### Fisher info:\n",
"$$\n",
"I_{jk}(\\beta) = \\frac{1}{\\phi} \\sum_{i=1}^{n} W_{i} x_{ji} x_{ki} \\quad j, k = \\bar{0,p} \\quad i = \\bar{1,n}\n",
"$$\n",
"\n",
"#### Solution:\n",
"\n",
"$$\n",
"\\frac{\\partial \\eta}{\\partial \\mu} = 1 \\quad W_{i} = \\frac{1}{1 * 1^2} = 1 \\quad j = \\bar{0..p}\n",
"$$\n",
"\n",
"$$\n",
"\\upsilon(\\beta_{j}) = \\frac{1}{\\sigma^2} \\sum_{i=1}^{n} 1 (y_{i} - \\mu_{i}) x_{ji}\n",
"$$\n",
"\n",
"#### Fisher info, populated:\n",
"$$\n",
"I_{jk}(\\beta) = \\frac{1}{\\sigma^2} \\sum_{i=1}^{n} 1 x_{ji} \\space x_{ki} \\quad j, k = \\bar{0,p} \\quad i = \\bar{1,n}\n",
"$$\n",
"\n",
"### Q2-)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "68930a73",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"dataframe\">\n",
"<caption>A data.frame: 31 × 8</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>Miners</th><th scope=col>Eucs</th><th scope=col>Area</th><th scope=col>Grazed</th><th scope=col>Shrubs</th><th scope=col>Bulokes</th><th scope=col>Timber</th><th scope=col>Minerab</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td>0</td><td> 2</td><td>22</td><td>0</td><td>1</td><td>120</td><td>16</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td>10</td><td>11</td><td>0</td><td>1</td><td> 67</td><td>25</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>16</td><td>51</td><td>0</td><td>1</td><td> 85</td><td>13</td><td> 3</td></tr>\n",
"\t<tr><td>1</td><td>20</td><td>22</td><td>0</td><td>1</td><td> 45</td><td>12</td><td> 2</td></tr>\n",
"\t<tr><td>1</td><td>19</td><td> 4</td><td>0</td><td>1</td><td>160</td><td>14</td><td> 8</td></tr>\n",
"\t<tr><td>1</td><td>18</td><td>61</td><td>0</td><td>1</td><td> 75</td><td> 6</td><td> 1</td></tr>\n",
"\t<tr><td>1</td><td>12</td><td>16</td><td>0</td><td>1</td><td>100</td><td>12</td><td> 8</td></tr>\n",
"\t<tr><td>1</td><td>16</td><td>14</td><td>0</td><td>1</td><td>321</td><td>15</td><td> 5</td></tr>\n",
"\t<tr><td>0</td><td> 3</td><td> 5</td><td>0</td><td>1</td><td>275</td><td> 8</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>12</td><td> 6</td><td>1</td><td>0</td><td>227</td><td>10</td><td> 4</td></tr>\n",
"\t<tr><td>1</td><td>32</td><td> 3</td><td>0</td><td>0</td><td> 23</td><td>61</td><td>19</td></tr>\n",
"\t<tr><td>0</td><td> 2</td><td>13</td><td>1</td><td>0</td><td>277</td><td>22</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>16</td><td>24</td><td>1</td><td>0</td><td>243</td><td>25</td><td> 2</td></tr>\n",
"\t<tr><td>0</td><td> 7</td><td>32</td><td>0</td><td>0</td><td>117</td><td> 9</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>10</td><td>22</td><td>1</td><td>0</td><td> 97</td><td>35</td><td> 3</td></tr>\n",
"\t<tr><td>1</td><td>15</td><td>27</td><td>1</td><td>0</td><td> 99</td><td>25</td><td> 1</td></tr>\n",
"\t<tr><td>1</td><td>30</td><td>18</td><td>0</td><td>0</td><td> 45</td><td>19</td><td> 7</td></tr>\n",
"\t<tr><td>1</td><td> 4</td><td> 9</td><td>1</td><td>0</td><td>123</td><td>22</td><td> 1</td></tr>\n",
"\t<tr><td>0</td><td> 4</td><td>21</td><td>1</td><td>1</td><td> 96</td><td>22</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>19</td><td>14</td><td>0</td><td>0</td><td>261</td><td>11</td><td> 7</td></tr>\n",
"\t<tr><td>0</td><td>11</td><td> 4</td><td>1</td><td>0</td><td>160</td><td>45</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 0</td><td> 5</td><td>1</td><td>0</td><td>192</td><td>22</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 0</td><td>19</td><td>1</td><td>0</td><td>323</td><td>60</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 0</td><td>12</td><td>1</td><td>0</td><td>190</td><td>48</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 3</td><td>15</td><td>0</td><td>1</td><td>222</td><td>26</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 8</td><td>38</td><td>1</td><td>0</td><td>148</td><td>38</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td> 8</td><td>24</td><td>0</td><td>1</td><td>198</td><td>28</td><td> 0</td></tr>\n",
"\t<tr><td>0</td><td>15</td><td>16</td><td>0</td><td>1</td><td> 93</td><td>16</td><td> 0</td></tr>\n",
"\t<tr><td>1</td><td>21</td><td> 6</td><td>0</td><td>0</td><td> 88</td><td>45</td><td> 3</td></tr>\n",
"\t<tr><td>1</td><td>24</td><td>16</td><td>1</td><td>0</td><td>138</td><td>21</td><td> 4</td></tr>\n",
"\t<tr><td>1</td><td>15</td><td>15</td><td>1</td><td>0</td><td> 70</td><td>31</td><td> 6</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 31 × 8\n",
"\\begin{tabular}{llllllll}\n",
" Miners & Eucs & Area & Grazed & Shrubs & Bulokes & Timber & Minerab\\\\\n",
" <int> & <int> & <int> & <int> & <int> & <int> & <int> & <int>\\\\\n",
"\\hline\n",
"\t 0 & 2 & 22 & 0 & 1 & 120 & 16 & 0\\\\\n",
"\t 0 & 10 & 11 & 0 & 1 & 67 & 25 & 0\\\\\n",
"\t 1 & 16 & 51 & 0 & 1 & 85 & 13 & 3\\\\\n",
"\t 1 & 20 & 22 & 0 & 1 & 45 & 12 & 2\\\\\n",
"\t 1 & 19 & 4 & 0 & 1 & 160 & 14 & 8\\\\\n",
"\t 1 & 18 & 61 & 0 & 1 & 75 & 6 & 1\\\\\n",
"\t 1 & 12 & 16 & 0 & 1 & 100 & 12 & 8\\\\\n",
"\t 1 & 16 & 14 & 0 & 1 & 321 & 15 & 5\\\\\n",
"\t 0 & 3 & 5 & 0 & 1 & 275 & 8 & 0\\\\\n",
"\t 1 & 12 & 6 & 1 & 0 & 227 & 10 & 4\\\\\n",
"\t 1 & 32 & 3 & 0 & 0 & 23 & 61 & 19\\\\\n",
"\t 0 & 2 & 13 & 1 & 0 & 277 & 22 & 0\\\\\n",
"\t 1 & 16 & 24 & 1 & 0 & 243 & 25 & 2\\\\\n",
"\t 0 & 7 & 32 & 0 & 0 & 117 & 9 & 0\\\\\n",
"\t 1 & 10 & 22 & 1 & 0 & 97 & 35 & 3\\\\\n",
"\t 1 & 15 & 27 & 1 & 0 & 99 & 25 & 1\\\\\n",
"\t 1 & 30 & 18 & 0 & 0 & 45 & 19 & 7\\\\\n",
"\t 1 & 4 & 9 & 1 & 0 & 123 & 22 & 1\\\\\n",
"\t 0 & 4 & 21 & 1 & 1 & 96 & 22 & 0\\\\\n",
"\t 1 & 19 & 14 & 0 & 0 & 261 & 11 & 7\\\\\n",
"\t 0 & 11 & 4 & 1 & 0 & 160 & 45 & 0\\\\\n",
"\t 0 & 0 & 5 & 1 & 0 & 192 & 22 & 0\\\\\n",
"\t 0 & 0 & 19 & 1 & 0 & 323 & 60 & 0\\\\\n",
"\t 0 & 0 & 12 & 1 & 0 & 190 & 48 & 0\\\\\n",
"\t 0 & 3 & 15 & 0 & 1 & 222 & 26 & 0\\\\\n",
"\t 0 & 8 & 38 & 1 & 0 & 148 & 38 & 0\\\\\n",
"\t 0 & 8 & 24 & 0 & 1 & 198 & 28 & 0\\\\\n",
"\t 0 & 15 & 16 & 0 & 1 & 93 & 16 & 0\\\\\n",
"\t 1 & 21 & 6 & 0 & 0 & 88 & 45 & 3\\\\\n",
"\t 1 & 24 & 16 & 1 & 0 & 138 & 21 & 4\\\\\n",
"\t 1 & 15 & 15 & 1 & 0 & 70 & 31 & 6\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 31 × 8\n",
"\n",
"| Miners &lt;int&gt; | Eucs &lt;int&gt; | Area &lt;int&gt; | Grazed &lt;int&gt; | Shrubs &lt;int&gt; | Bulokes &lt;int&gt; | Timber &lt;int&gt; | Minerab &lt;int&gt; |\n",
"|---|---|---|---|---|---|---|---|\n",
"| 0 | 2 | 22 | 0 | 1 | 120 | 16 | 0 |\n",
"| 0 | 10 | 11 | 0 | 1 | 67 | 25 | 0 |\n",
"| 1 | 16 | 51 | 0 | 1 | 85 | 13 | 3 |\n",
"| 1 | 20 | 22 | 0 | 1 | 45 | 12 | 2 |\n",
"| 1 | 19 | 4 | 0 | 1 | 160 | 14 | 8 |\n",
"| 1 | 18 | 61 | 0 | 1 | 75 | 6 | 1 |\n",
"| 1 | 12 | 16 | 0 | 1 | 100 | 12 | 8 |\n",
"| 1 | 16 | 14 | 0 | 1 | 321 | 15 | 5 |\n",
"| 0 | 3 | 5 | 0 | 1 | 275 | 8 | 0 |\n",
"| 1 | 12 | 6 | 1 | 0 | 227 | 10 | 4 |\n",
"| 1 | 32 | 3 | 0 | 0 | 23 | 61 | 19 |\n",
"| 0 | 2 | 13 | 1 | 0 | 277 | 22 | 0 |\n",
"| 1 | 16 | 24 | 1 | 0 | 243 | 25 | 2 |\n",
"| 0 | 7 | 32 | 0 | 0 | 117 | 9 | 0 |\n",
"| 1 | 10 | 22 | 1 | 0 | 97 | 35 | 3 |\n",
"| 1 | 15 | 27 | 1 | 0 | 99 | 25 | 1 |\n",
"| 1 | 30 | 18 | 0 | 0 | 45 | 19 | 7 |\n",
"| 1 | 4 | 9 | 1 | 0 | 123 | 22 | 1 |\n",
"| 0 | 4 | 21 | 1 | 1 | 96 | 22 | 0 |\n",
"| 1 | 19 | 14 | 0 | 0 | 261 | 11 | 7 |\n",
"| 0 | 11 | 4 | 1 | 0 | 160 | 45 | 0 |\n",
"| 0 | 0 | 5 | 1 | 0 | 192 | 22 | 0 |\n",
"| 0 | 0 | 19 | 1 | 0 | 323 | 60 | 0 |\n",
"| 0 | 0 | 12 | 1 | 0 | 190 | 48 | 0 |\n",
"| 0 | 3 | 15 | 0 | 1 | 222 | 26 | 0 |\n",
"| 0 | 8 | 38 | 1 | 0 | 148 | 38 | 0 |\n",
"| 0 | 8 | 24 | 0 | 1 | 198 | 28 | 0 |\n",
"| 0 | 15 | 16 | 0 | 1 | 93 | 16 | 0 |\n",
"| 1 | 21 | 6 | 0 | 0 | 88 | 45 | 3 |\n",
"| 1 | 24 | 16 | 1 | 0 | 138 | 21 | 4 |\n",
"| 1 | 15 | 15 | 1 | 0 | 70 | 31 | 6 |\n",
"\n"
],
"text/plain": [
" Miners Eucs Area Grazed Shrubs Bulokes Timber Minerab\n",
"1 0 2 22 0 1 120 16 0 \n",
"2 0 10 11 0 1 67 25 0 \n",
"3 1 16 51 0 1 85 13 3 \n",
"4 1 20 22 0 1 45 12 2 \n",
"5 1 19 4 0 1 160 14 8 \n",
"6 1 18 61 0 1 75 6 1 \n",
"7 1 12 16 0 1 100 12 8 \n",
"8 1 16 14 0 1 321 15 5 \n",
"9 0 3 5 0 1 275 8 0 \n",
"10 1 12 6 1 0 227 10 4 \n",
"11 1 32 3 0 0 23 61 19 \n",
"12 0 2 13 1 0 277 22 0 \n",
"13 1 16 24 1 0 243 25 2 \n",
"14 0 7 32 0 0 117 9 0 \n",
"15 1 10 22 1 0 97 35 3 \n",
"16 1 15 27 1 0 99 25 1 \n",
"17 1 30 18 0 0 45 19 7 \n",
"18 1 4 9 1 0 123 22 1 \n",
"19 0 4 21 1 1 96 22 0 \n",
"20 1 19 14 0 0 261 11 7 \n",
"21 0 11 4 1 0 160 45 0 \n",
"22 0 0 5 1 0 192 22 0 \n",
"23 0 0 19 1 0 323 60 0 \n",
"24 0 0 12 1 0 190 48 0 \n",
"25 0 3 15 0 1 222 26 0 \n",
"26 0 8 38 1 0 148 38 0 \n",
"27 0 8 24 0 1 198 28 0 \n",
"28 0 15 16 0 1 93 16 0 \n",
"29 1 21 6 0 0 88 45 3 \n",
"30 1 24 16 1 0 138 21 4 \n",
"31 1 15 15 1 0 70 31 6 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#install.packages(\"GLMsData\")\n",
"# First, get the data:\n",
"library(GLMsData)\n",
"data(nminer)\n",
"\n",
"nminer"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "13bef513",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
" Miners Eucs Area Grazed \n",
" Min. :0.0000 Min. : 0 Min. : 3.00 Min. :0.0000 \n",
" 1st Qu.:0.0000 1st Qu.: 4 1st Qu.:10.00 1st Qu.:0.0000 \n",
" Median :1.0000 Median :12 Median :16.00 Median :0.0000 \n",
" Mean :0.5484 Mean :12 Mean :18.23 Mean :0.4516 \n",
" 3rd Qu.:1.0000 3rd Qu.:17 3rd Qu.:22.00 3rd Qu.:1.0000 \n",
" Max. :1.0000 Max. :32 Max. :61.00 Max. :1.0000 \n",
" Shrubs Bulokes Timber Minerab \n",
" Min. :0.0000 Min. : 23.0 Min. : 6.00 Min. : 0.00 \n",
" 1st Qu.:0.0000 1st Qu.: 90.5 1st Qu.:13.50 1st Qu.: 0.00 \n",
" Median :0.0000 Median :123.0 Median :22.00 Median : 1.00 \n",
" Mean :0.4194 Mean :150.9 Mean :24.58 Mean : 2.71 \n",
" 3rd Qu.:1.0000 3rd Qu.:210.0 3rd Qu.:29.50 3rd Qu.: 4.00 \n",
" Max. :1.0000 Max. :323.0 Max. :61.00 Max. :19.00 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"summary(nminer)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "8270d0cd",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Deviance = 82.14668 Iterations - 1\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message:\n",
"“glm.fit: algorithm did not converge”\n"
]
}
],
"source": [
"nm.m1 <- glm(Minerab ~ Eucs,\n",
" data = nminer,\n",
" family = poisson(link=\"log\"),\n",
" control = list(maxit=1, trace=TRUE))"
]
},
{
"cell_type": "code",
"execution_count": 79,
"id": "24fdea97",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<style>\n",
".dl-inline {width: auto; margin:0; padding: 0}\n",
".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
"</style><dl class=dl-inline><dt>(Intercept)</dt><dd>0.122335668948972</dd><dt>Eucs</dt><dd>0.0810710013086576</dd></dl>\n"
],
"text/latex": [
"\\begin{description*}\n",
"\\item[(Intercept)] 0.122335668948972\n",
"\\item[Eucs] 0.0810710013086576\n",
"\\end{description*}\n"
],
"text/markdown": [
"(Intercept)\n",
": 0.122335668948972Eucs\n",
": 0.0810710013086576\n",
"\n"
],
"text/plain": [
"(Intercept) Eucs \n",
" 0.1223357 0.0810710 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"coef_iter1 <- coef(nm.m1)\n",
"\n",
"coef_iter1"
]
},
{
"cell_type": "code",
"execution_count": 122,
"id": "fa76fa7a",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Deviance = 82.14668 Iterations - 1\n",
"Deviance = 64.49515 Iterations - 2\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message:\n",
"“glm.fit: algorithm did not converge”\n"
]
}
],
"source": [
"nm.m2 <- glm(Minerab ~ Eucs,\n",
" data = nminer,\n",
" family = poisson(link=\"log\"),\n",
" control = list(maxit=2, trace=TRUE))"
]
},
{
"cell_type": "code",
"execution_count": 123,
"id": "efad2554",
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/html": [
"<style>\n",
".dl-inline {width: auto; margin:0; padding: 0}\n",
".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
"</style><dl class=dl-inline><dt>(Intercept)</dt><dd>-0.589798413763625</dd><dt>Eucs</dt><dd>0.103745243386491</dd></dl>\n"
],
"text/latex": [
"\\begin{description*}\n",
"\\item[(Intercept)] -0.589798413763625\n",
"\\item[Eucs] 0.103745243386491\n",
"\\end{description*}\n"
],
"text/markdown": [
"(Intercept)\n",
": -0.589798413763625Eucs\n",
": 0.103745243386491\n",
"\n"
],
"text/plain": [
"(Intercept) Eucs \n",
" -0.5897984 0.1037452 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"coef_iter2 <- coef(nm.m2)\n",
"\n",
"coef_iter2"
]
},
{
"cell_type": "code",
"execution_count": 128,
"id": "b75f6cb6",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Deviance = 82.14668 Iterations - 1\n",
"Deviance = 64.49515 Iterations - 2\n",
"Deviance = 63.32603 Iterations - 3\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message:\n",
"“glm.fit: algorithm did not converge”\n"
]
}
],
"source": [
"nm.m3 <- glm(Minerab ~ Eucs,\n",
" data = nminer,\n",
" family = poisson(link=\"log\"),\n",
" control = list(maxit=3, trace=TRUE))"
]
},
{
"cell_type": "code",
"execution_count": 129,
"id": "48c5b2dc",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"\n",
"Call: glm(formula = Minerab ~ Eucs, family = poisson(link = \"log\"), \n",
" data = nminer, control = list(maxit = 3, trace = TRUE))\n",
"\n",
"Coefficients:\n",
"(Intercept) Eucs \n",
" -0.8520 0.1131 \n",
"\n",
"Degrees of Freedom: 30 Total (i.e. Null); 29 Residual\n",
"Null Deviance:\t 150.5 \n",
"Residual Deviance: 63.33 \tAIC: 121.5"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"coef_iter3 <- coef(nm.m3)\n",
"\n",
"nm.m3"
]
},
{
"cell_type": "code",
"execution_count": 139,
"id": "09731f4e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Deviance = 82.14668 Iterations - 1\n",
"Deviance = 64.49515 Iterations - 2\n",
"Deviance = 63.32603 Iterations - 3\n",
"Deviance = 63.31798 Iterations - 4\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message:\n",
"“glm.fit: algorithm did not converge”\n"
]
}
],
"source": [
"nm.m4 <- glm(Minerab ~ Eucs,\n",
" data = nminer,\n",
" family = poisson(link=\"log\"),\n",
" control = list(maxit=4, trace=TRUE))"
]
},
{
"cell_type": "code",
"execution_count": 141,
"id": "b662924f",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<style>\n",
".dl-inline {width: auto; margin:0; padding: 0}\n",
".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
"</style><dl class=dl-inline><dt>(Intercept)</dt><dd>-0.876030743734818</dd><dt>Eucs</dt><dd>0.113975000660348</dd></dl>\n"
],
"text/latex": [
"\\begin{description*}\n",
"\\item[(Intercept)] -0.876030743734818\n",
"\\item[Eucs] 0.113975000660348\n",
"\\end{description*}\n"
],
"text/markdown": [
"(Intercept)\n",
": -0.876030743734818Eucs\n",
": 0.113975000660348\n",
"\n"
],
"text/plain": [
"(Intercept) Eucs \n",
" -0.8760307 0.1139750 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"coef_iter4 <- coef(nm.m4)\n",
"\n",
"coef_iter4"
]
},
{
"cell_type": "code",
"execution_count": 103,
"id": "c3745835",
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"dataframe\">\n",
"<caption>A data.frame: 31 × 8</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>row</th><th scope=col>yi</th><th scope=col>uyumdegeri</th><th scope=col>calismadegeri</th><th scope=col>agirliklar</th><th scope=col>x</th><th scope=col>mu1i</th><th scope=col>unit_dev</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 1</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 2</td><td> 1.33</td><td>2.66</td></tr>\n",
"\t<tr><td> 2</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td>10</td><td> 2.54</td><td>5.08</td></tr>\n",
"\t<tr><td> 3</td><td> 3</td><td> 3.1</td><td> 1.099</td><td> 3.1</td><td>16</td><td> 4.13</td><td>0.54</td></tr>\n",
"\t<tr><td> 4</td><td> 2</td><td> 2.1</td><td> 0.694</td><td> 2.1</td><td>20</td><td> 5.72</td><td>3.43</td></tr>\n",
"\t<tr><td> 5</td><td> 8</td><td> 8.1</td><td> 2.080</td><td> 8.1</td><td>19</td><td> 5.27</td><td>1.41</td></tr>\n",
"\t<tr><td> 6</td><td> 1</td><td> 1.1</td><td> 0.004</td><td> 1.1</td><td>18</td><td> 4.86</td><td>4.75</td></tr>\n",
"\t<tr><td> 7</td><td> 8</td><td> 8.1</td><td> 2.080</td><td> 8.1</td><td>12</td><td> 2.99</td><td>5.93</td></tr>\n",
"\t<tr><td> 8</td><td> 5</td><td> 5.1</td><td> 1.610</td><td> 5.1</td><td>16</td><td> 4.13</td><td>0.37</td></tr>\n",
"\t<tr><td> 9</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 3</td><td> 1.44</td><td>2.88</td></tr>\n",
"\t<tr><td>10</td><td> 4</td><td> 4.1</td><td> 1.387</td><td> 4.1</td><td>12</td><td> 2.99</td><td>0.51</td></tr>\n",
"\t<tr><td>11</td><td>19</td><td>19.1</td><td> 2.944</td><td>19.1</td><td>32</td><td>15.13</td><td>1.11</td></tr>\n",
"\t<tr><td>12</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 2</td><td> 1.33</td><td>2.66</td></tr>\n",
"\t<tr><td>13</td><td> 2</td><td> 2.1</td><td> 0.694</td><td> 2.1</td><td>16</td><td> 4.13</td><td>1.56</td></tr>\n",
"\t<tr><td>14</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 7</td><td> 1.99</td><td>3.99</td></tr>\n",
"\t<tr><td>15</td><td> 3</td><td> 3.1</td><td> 1.099</td><td> 3.1</td><td>10</td><td> 2.54</td><td>0.27</td></tr>\n",
"\t<tr><td>16</td><td> 1</td><td> 1.1</td><td> 0.004</td><td> 1.1</td><td>15</td><td> 3.81</td><td>3.14</td></tr>\n",
"\t<tr><td>17</td><td> 7</td><td> 7.1</td><td> 1.946</td><td> 7.1</td><td>30</td><td>12.86</td><td>3.41</td></tr>\n",
"\t<tr><td>18</td><td> 1</td><td> 1.1</td><td> 0.004</td><td> 1.1</td><td> 4</td><td> 1.56</td><td>0.42</td></tr>\n",
"\t<tr><td>19</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 4</td><td> 1.56</td><td>3.13</td></tr>\n",
"\t<tr><td>20</td><td> 7</td><td> 7.1</td><td> 1.946</td><td> 7.1</td><td>19</td><td> 5.27</td><td>0.71</td></tr>\n",
"\t<tr><td>21</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td>11</td><td> 2.76</td><td>5.51</td></tr>\n",
"\t<tr><td>22</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 0</td><td> 1.13</td><td>2.26</td></tr>\n",
"\t<tr><td>23</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 0</td><td> 1.13</td><td>2.26</td></tr>\n",
"\t<tr><td>24</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 0</td><td> 1.13</td><td>2.26</td></tr>\n",
"\t<tr><td>25</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 3</td><td> 1.44</td><td>2.88</td></tr>\n",
"\t<tr><td>26</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 8</td><td> 2.16</td><td>4.32</td></tr>\n",
"\t<tr><td>27</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td> 8</td><td> 2.16</td><td>4.32</td></tr>\n",
"\t<tr><td>28</td><td> 0</td><td> 0.1</td><td>-3.303</td><td> 0.1</td><td>15</td><td> 3.81</td><td>7.63</td></tr>\n",
"\t<tr><td>29</td><td> 3</td><td> 3.1</td><td> 1.099</td><td> 3.1</td><td>21</td><td> 6.20</td><td>2.24</td></tr>\n",
"\t<tr><td>30</td><td> 4</td><td> 4.1</td><td> 1.387</td><td> 4.1</td><td>24</td><td> 7.91</td><td>2.56</td></tr>\n",
"\t<tr><td>31</td><td> 6</td><td> 6.1</td><td> 1.792</td><td> 6.1</td><td>15</td><td> 3.81</td><td>1.26</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 31 × 8\n",
"\\begin{tabular}{llllllll}\n",
" row & yi & uyumdegeri & calismadegeri & agirliklar & x & mu1i & unit\\_dev\\\\\n",
" <int> & <int> & <dbl> & <dbl> & <dbl> & <int> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 1 & 0 & 0.1 & -3.303 & 0.1 & 2 & 1.33 & 2.66\\\\\n",
"\t 2 & 0 & 0.1 & -3.303 & 0.1 & 10 & 2.54 & 5.08\\\\\n",
"\t 3 & 3 & 3.1 & 1.099 & 3.1 & 16 & 4.13 & 0.54\\\\\n",
"\t 4 & 2 & 2.1 & 0.694 & 2.1 & 20 & 5.72 & 3.43\\\\\n",
"\t 5 & 8 & 8.1 & 2.080 & 8.1 & 19 & 5.27 & 1.41\\\\\n",
"\t 6 & 1 & 1.1 & 0.004 & 1.1 & 18 & 4.86 & 4.75\\\\\n",
"\t 7 & 8 & 8.1 & 2.080 & 8.1 & 12 & 2.99 & 5.93\\\\\n",
"\t 8 & 5 & 5.1 & 1.610 & 5.1 & 16 & 4.13 & 0.37\\\\\n",
"\t 9 & 0 & 0.1 & -3.303 & 0.1 & 3 & 1.44 & 2.88\\\\\n",
"\t 10 & 4 & 4.1 & 1.387 & 4.1 & 12 & 2.99 & 0.51\\\\\n",
"\t 11 & 19 & 19.1 & 2.944 & 19.1 & 32 & 15.13 & 1.11\\\\\n",
"\t 12 & 0 & 0.1 & -3.303 & 0.1 & 2 & 1.33 & 2.66\\\\\n",
"\t 13 & 2 & 2.1 & 0.694 & 2.1 & 16 & 4.13 & 1.56\\\\\n",
"\t 14 & 0 & 0.1 & -3.303 & 0.1 & 7 & 1.99 & 3.99\\\\\n",
"\t 15 & 3 & 3.1 & 1.099 & 3.1 & 10 & 2.54 & 0.27\\\\\n",
"\t 16 & 1 & 1.1 & 0.004 & 1.1 & 15 & 3.81 & 3.14\\\\\n",
"\t 17 & 7 & 7.1 & 1.946 & 7.1 & 30 & 12.86 & 3.41\\\\\n",
"\t 18 & 1 & 1.1 & 0.004 & 1.1 & 4 & 1.56 & 0.42\\\\\n",
"\t 19 & 0 & 0.1 & -3.303 & 0.1 & 4 & 1.56 & 3.13\\\\\n",
"\t 20 & 7 & 7.1 & 1.946 & 7.1 & 19 & 5.27 & 0.71\\\\\n",
"\t 21 & 0 & 0.1 & -3.303 & 0.1 & 11 & 2.76 & 5.51\\\\\n",
"\t 22 & 0 & 0.1 & -3.303 & 0.1 & 0 & 1.13 & 2.26\\\\\n",
"\t 23 & 0 & 0.1 & -3.303 & 0.1 & 0 & 1.13 & 2.26\\\\\n",
"\t 24 & 0 & 0.1 & -3.303 & 0.1 & 0 & 1.13 & 2.26\\\\\n",
"\t 25 & 0 & 0.1 & -3.303 & 0.1 & 3 & 1.44 & 2.88\\\\\n",
"\t 26 & 0 & 0.1 & -3.303 & 0.1 & 8 & 2.16 & 4.32\\\\\n",
"\t 27 & 0 & 0.1 & -3.303 & 0.1 & 8 & 2.16 & 4.32\\\\\n",
"\t 28 & 0 & 0.1 & -3.303 & 0.1 & 15 & 3.81 & 7.63\\\\\n",
"\t 29 & 3 & 3.1 & 1.099 & 3.1 & 21 & 6.20 & 2.24\\\\\n",
"\t 30 & 4 & 4.1 & 1.387 & 4.1 & 24 & 7.91 & 2.56\\\\\n",
"\t 31 & 6 & 6.1 & 1.792 & 6.1 & 15 & 3.81 & 1.26\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 31 × 8\n",
"\n",
"| row &lt;int&gt; | yi &lt;int&gt; | uyumdegeri &lt;dbl&gt; | calismadegeri &lt;dbl&gt; | agirliklar &lt;dbl&gt; | x &lt;int&gt; | mu1i &lt;dbl&gt; | unit_dev &lt;dbl&gt; |\n",
"|---|---|---|---|---|---|---|---|\n",
"| 1 | 0 | 0.1 | -3.303 | 0.1 | 2 | 1.33 | 2.66 |\n",
"| 2 | 0 | 0.1 | -3.303 | 0.1 | 10 | 2.54 | 5.08 |\n",
"| 3 | 3 | 3.1 | 1.099 | 3.1 | 16 | 4.13 | 0.54 |\n",
"| 4 | 2 | 2.1 | 0.694 | 2.1 | 20 | 5.72 | 3.43 |\n",
"| 5 | 8 | 8.1 | 2.080 | 8.1 | 19 | 5.27 | 1.41 |\n",
"| 6 | 1 | 1.1 | 0.004 | 1.1 | 18 | 4.86 | 4.75 |\n",
"| 7 | 8 | 8.1 | 2.080 | 8.1 | 12 | 2.99 | 5.93 |\n",
"| 8 | 5 | 5.1 | 1.610 | 5.1 | 16 | 4.13 | 0.37 |\n",
"| 9 | 0 | 0.1 | -3.303 | 0.1 | 3 | 1.44 | 2.88 |\n",
"| 10 | 4 | 4.1 | 1.387 | 4.1 | 12 | 2.99 | 0.51 |\n",
"| 11 | 19 | 19.1 | 2.944 | 19.1 | 32 | 15.13 | 1.11 |\n",
"| 12 | 0 | 0.1 | -3.303 | 0.1 | 2 | 1.33 | 2.66 |\n",
"| 13 | 2 | 2.1 | 0.694 | 2.1 | 16 | 4.13 | 1.56 |\n",
"| 14 | 0 | 0.1 | -3.303 | 0.1 | 7 | 1.99 | 3.99 |\n",
"| 15 | 3 | 3.1 | 1.099 | 3.1 | 10 | 2.54 | 0.27 |\n",
"| 16 | 1 | 1.1 | 0.004 | 1.1 | 15 | 3.81 | 3.14 |\n",
"| 17 | 7 | 7.1 | 1.946 | 7.1 | 30 | 12.86 | 3.41 |\n",
"| 18 | 1 | 1.1 | 0.004 | 1.1 | 4 | 1.56 | 0.42 |\n",
"| 19 | 0 | 0.1 | -3.303 | 0.1 | 4 | 1.56 | 3.13 |\n",
"| 20 | 7 | 7.1 | 1.946 | 7.1 | 19 | 5.27 | 0.71 |\n",
"| 21 | 0 | 0.1 | -3.303 | 0.1 | 11 | 2.76 | 5.51 |\n",
"| 22 | 0 | 0.1 | -3.303 | 0.1 | 0 | 1.13 | 2.26 |\n",
"| 23 | 0 | 0.1 | -3.303 | 0.1 | 0 | 1.13 | 2.26 |\n",
"| 24 | 0 | 0.1 | -3.303 | 0.1 | 0 | 1.13 | 2.26 |\n",
"| 25 | 0 | 0.1 | -3.303 | 0.1 | 3 | 1.44 | 2.88 |\n",
"| 26 | 0 | 0.1 | -3.303 | 0.1 | 8 | 2.16 | 4.32 |\n",
"| 27 | 0 | 0.1 | -3.303 | 0.1 | 8 | 2.16 | 4.32 |\n",
"| 28 | 0 | 0.1 | -3.303 | 0.1 | 15 | 3.81 | 7.63 |\n",
"| 29 | 3 | 3.1 | 1.099 | 3.1 | 21 | 6.20 | 2.24 |\n",
"| 30 | 4 | 4.1 | 1.387 | 4.1 | 24 | 7.91 | 2.56 |\n",
"| 31 | 6 | 6.1 | 1.792 | 6.1 | 15 | 3.81 | 1.26 |\n",
"\n"
],
"text/plain": [
" row yi uyumdegeri calismadegeri agirliklar x mu1i unit_dev\n",
"1 1 0 0.1 -3.303 0.1 2 1.33 2.66 \n",
"2 2 0 0.1 -3.303 0.1 10 2.54 5.08 \n",
"3 3 3 3.1 1.099 3.1 16 4.13 0.54 \n",
"4 4 2 2.1 0.694 2.1 20 5.72 3.43 \n",
"5 5 8 8.1 2.080 8.1 19 5.27 1.41 \n",
"6 6 1 1.1 0.004 1.1 18 4.86 4.75 \n",
"7 7 8 8.1 2.080 8.1 12 2.99 5.93 \n",
"8 8 5 5.1 1.610 5.1 16 4.13 0.37 \n",
"9 9 0 0.1 -3.303 0.1 3 1.44 2.88 \n",
"10 10 4 4.1 1.387 4.1 12 2.99 0.51 \n",
"11 11 19 19.1 2.944 19.1 32 15.13 1.11 \n",
"12 12 0 0.1 -3.303 0.1 2 1.33 2.66 \n",
"13 13 2 2.1 0.694 2.1 16 4.13 1.56 \n",
"14 14 0 0.1 -3.303 0.1 7 1.99 3.99 \n",
"15 15 3 3.1 1.099 3.1 10 2.54 0.27 \n",
"16 16 1 1.1 0.004 1.1 15 3.81 3.14 \n",
"17 17 7 7.1 1.946 7.1 30 12.86 3.41 \n",
"18 18 1 1.1 0.004 1.1 4 1.56 0.42 \n",
"19 19 0 0.1 -3.303 0.1 4 1.56 3.13 \n",
"20 20 7 7.1 1.946 7.1 19 5.27 0.71 \n",
"21 21 0 0.1 -3.303 0.1 11 2.76 5.51 \n",
"22 22 0 0.1 -3.303 0.1 0 1.13 2.26 \n",
"23 23 0 0.1 -3.303 0.1 0 1.13 2.26 \n",
"24 24 0 0.1 -3.303 0.1 0 1.13 2.26 \n",
"25 25 0 0.1 -3.303 0.1 3 1.44 2.88 \n",
"26 26 0 0.1 -3.303 0.1 8 2.16 4.32 \n",
"27 27 0 0.1 -3.303 0.1 8 2.16 4.32 \n",
"28 28 0 0.1 -3.303 0.1 15 3.81 7.63 \n",
"29 29 3 3.1 1.099 3.1 21 6.20 2.24 \n",
"30 30 4 4.1 1.387 4.1 24 7.91 2.56 \n",
"31 31 6 6.1 1.792 6.1 15 3.81 1.26 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"yi = (nminer$Minerab)\n",
"etai <- (nminer$Minerab + 0.1)\n",
"calculate_eta <- function(etai_){\n",
" log(etai_)\n",
"}\n",
"all_etai <- calculate_eta(etai)\n",
"\n",
"calculate_z1i <- function(a, b, c){\n",
" a + (b - c) / c\n",
"}\n",
"\n",
"all_z1i <- calculate_z1i(all_etai, yi, etai)\n",
"\n",
"calculate_w1i <- function(a){\n",
" a\n",
"}\n",
"\n",
"all_wieghts <- calculate_w1i(etai)\n",
"\n",
"x <- (nminer$Eucs)\n",
"\n",
"calculate_mu1i <- function(coef_, x_){\n",
" exp(coef_[1] + (coef_[2] * x_))\n",
"}\n",
"\n",
"all_mu1i <- calculate_mu1i(coef_iter1, x)\n",
"\n",
"calculate_unit_dev <- function(y_, etai_, mu1_){\n",
" l_ <- etai_/mu1_\n",
" 2 * (y_ * log(l_) - (y_ - mu1_))\n",
"}\n",
"\n",
"all_unit_dev <- calculate_unit_dev(yi, etai, all_mu1i)\n",
"\n",
"tablo_iter1 = data.frame(row = seq.int(yi), yi = yi, uyumdegeri = round(etai, 3), calismadegeri = round(all_z1i, 3), agirliklar = all_wieghts, x = x, mu1i = round(all_mu1i, 2), unit_dev = round(all_unit_dev, 2))\n",
"tablo_iter1"
]
},
{
"cell_type": "code",
"execution_count": 127,
"id": "1ea5ca99",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"dataframe\">\n",
"<caption>A data.frame: 31 × 9</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>row</th><th scope=col>yi</th><th scope=col>uyumdegeri</th><th scope=col>mu1</th><th scope=col>calismadegeri</th><th scope=col>agirliklar</th><th scope=col>x</th><th scope=col>mu2i</th><th scope=col>unit_dev</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 1</td><td> 0</td><td> 0.1</td><td> 1.329068</td><td>-0.72</td><td> 1.329068</td><td> 2</td><td> 0.68</td><td> 1.36</td></tr>\n",
"\t<tr><td> 2</td><td> 0</td><td> 0.1</td><td> 2.542240</td><td>-0.07</td><td> 2.542240</td><td>10</td><td> 1.56</td><td> 3.13</td></tr>\n",
"\t<tr><td> 3</td><td> 3</td><td> 3.1</td><td> 4.134935</td><td> 1.14</td><td> 4.134935</td><td>16</td><td> 2.92</td><td> 0.20</td></tr>\n",
"\t<tr><td> 4</td><td> 2</td><td> 2.1</td><td> 5.718781</td><td> 1.09</td><td> 5.718781</td><td>20</td><td> 4.42</td><td> 1.86</td></tr>\n",
"\t<tr><td> 5</td><td> 8</td><td> 8.1</td><td> 5.273449</td><td> 2.18</td><td> 5.273449</td><td>19</td><td> 3.98</td><td> 3.33</td></tr>\n",
"\t<tr><td> 6</td><td> 1</td><td> 1.1</td><td> 4.862797</td><td> 0.79</td><td> 4.862797</td><td>18</td><td> 3.59</td><td> 2.81</td></tr>\n",
"\t<tr><td> 7</td><td> 8</td><td> 8.1</td><td> 2.989744</td><td> 2.77</td><td> 2.989744</td><td>12</td><td> 1.93</td><td>10.84</td></tr>\n",
"\t<tr><td> 8</td><td> 5</td><td> 5.1</td><td> 4.134935</td><td> 1.63</td><td> 4.134935</td><td>16</td><td> 2.92</td><td> 1.42</td></tr>\n",
"\t<tr><td> 9</td><td> 0</td><td> 0.1</td><td> 1.441305</td><td>-0.63</td><td> 1.441305</td><td> 3</td><td> 0.76</td><td> 1.51</td></tr>\n",
"\t<tr><td>10</td><td> 4</td><td> 4.1</td><td> 2.989744</td><td> 1.43</td><td> 2.989744</td><td>12</td><td> 1.93</td><td> 1.90</td></tr>\n",
"\t<tr><td>11</td><td>19</td><td>19.1</td><td>15.128913</td><td> 2.97</td><td>15.128913</td><td>32</td><td>15.33</td><td> 1.01</td></tr>\n",
"\t<tr><td>12</td><td> 0</td><td> 0.1</td><td> 1.329068</td><td>-0.72</td><td> 1.329068</td><td> 2</td><td> 0.68</td><td> 1.36</td></tr>\n",
"\t<tr><td>13</td><td> 2</td><td> 2.1</td><td> 4.134935</td><td> 0.90</td><td> 4.134935</td><td>16</td><td> 2.92</td><td> 0.52</td></tr>\n",
"\t<tr><td>14</td><td> 0</td><td> 0.1</td><td> 1.993382</td><td>-0.31</td><td> 1.993382</td><td> 7</td><td> 1.15</td><td> 2.29</td></tr>\n",
"\t<tr><td>15</td><td> 3</td><td> 3.1</td><td> 2.542240</td><td> 1.11</td><td> 2.542240</td><td>10</td><td> 1.56</td><td> 1.23</td></tr>\n",
"\t<tr><td>16</td><td> 1</td><td> 1.1</td><td> 3.812941</td><td> 0.60</td><td> 3.812941</td><td>15</td><td> 2.63</td><td> 1.51</td></tr>\n",
"\t<tr><td>17</td><td> 7</td><td> 7.1</td><td>12.864424</td><td> 2.10</td><td>12.864424</td><td>30</td><td>12.46</td><td> 3.05</td></tr>\n",
"\t<tr><td>18</td><td> 1</td><td> 1.1</td><td> 1.563020</td><td> 0.09</td><td> 1.563020</td><td> 4</td><td> 0.84</td><td> 0.22</td></tr>\n",
"\t<tr><td>19</td><td> 0</td><td> 0.1</td><td> 1.563020</td><td>-0.55</td><td> 1.563020</td><td> 4</td><td> 0.84</td><td> 1.68</td></tr>\n",
"\t<tr><td>20</td><td> 7</td><td> 7.1</td><td> 5.273449</td><td> 1.99</td><td> 5.273449</td><td>19</td><td> 3.98</td><td> 2.06</td></tr>\n",
"\t<tr><td>21</td><td> 0</td><td> 0.1</td><td> 2.756927</td><td> 0.01</td><td> 2.756927</td><td>11</td><td> 1.74</td><td> 3.47</td></tr>\n",
"\t<tr><td>22</td><td> 0</td><td> 0.1</td><td> 1.130133</td><td>-0.88</td><td> 1.130133</td><td> 0</td><td> 0.55</td><td> 1.11</td></tr>\n",
"\t<tr><td>23</td><td> 0</td><td> 0.1</td><td> 1.130133</td><td>-0.88</td><td> 1.130133</td><td> 0</td><td> 0.55</td><td> 1.11</td></tr>\n",
"\t<tr><td>24</td><td> 0</td><td> 0.1</td><td> 1.130133</td><td>-0.88</td><td> 1.130133</td><td> 0</td><td> 0.55</td><td> 1.11</td></tr>\n",
"\t<tr><td>25</td><td> 0</td><td> 0.1</td><td> 1.441305</td><td>-0.63</td><td> 1.441305</td><td> 3</td><td> 0.76</td><td> 1.51</td></tr>\n",
"\t<tr><td>26</td><td> 0</td><td> 0.1</td><td> 2.161719</td><td>-0.23</td><td> 2.161719</td><td> 8</td><td> 1.27</td><td> 2.54</td></tr>\n",
"\t<tr><td>27</td><td> 0</td><td> 0.1</td><td> 2.161719</td><td>-0.23</td><td> 2.161719</td><td> 8</td><td> 1.27</td><td> 2.54</td></tr>\n",
"\t<tr><td>28</td><td> 0</td><td> 0.1</td><td> 3.812941</td><td> 0.34</td><td> 3.812941</td><td>15</td><td> 2.63</td><td> 5.26</td></tr>\n",
"\t<tr><td>29</td><td> 3</td><td> 3.1</td><td> 6.201720</td><td> 1.31</td><td> 6.201720</td><td>21</td><td> 4.90</td><td> 1.05</td></tr>\n",
"\t<tr><td>30</td><td> 4</td><td> 4.1</td><td> 7.909303</td><td> 1.57</td><td> 7.909303</td><td>24</td><td> 6.69</td><td> 1.46</td></tr>\n",
"\t<tr><td>31</td><td> 6</td><td> 6.1</td><td> 3.812941</td><td> 1.91</td><td> 3.812941</td><td>15</td><td> 2.63</td><td> 3.36</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 31 × 9\n",
"\\begin{tabular}{lllllllll}\n",
" row & yi & uyumdegeri & mu1 & calismadegeri & agirliklar & x & mu2i & unit\\_dev\\\\\n",
" <int> & <int> & <dbl> & <dbl> & <dbl> & <dbl> & <int> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 1 & 0 & 0.1 & 1.329068 & -0.72 & 1.329068 & 2 & 0.68 & 1.36\\\\\n",
"\t 2 & 0 & 0.1 & 2.542240 & -0.07 & 2.542240 & 10 & 1.56 & 3.13\\\\\n",
"\t 3 & 3 & 3.1 & 4.134935 & 1.14 & 4.134935 & 16 & 2.92 & 0.20\\\\\n",
"\t 4 & 2 & 2.1 & 5.718781 & 1.09 & 5.718781 & 20 & 4.42 & 1.86\\\\\n",
"\t 5 & 8 & 8.1 & 5.273449 & 2.18 & 5.273449 & 19 & 3.98 & 3.33\\\\\n",
"\t 6 & 1 & 1.1 & 4.862797 & 0.79 & 4.862797 & 18 & 3.59 & 2.81\\\\\n",
"\t 7 & 8 & 8.1 & 2.989744 & 2.77 & 2.989744 & 12 & 1.93 & 10.84\\\\\n",
"\t 8 & 5 & 5.1 & 4.134935 & 1.63 & 4.134935 & 16 & 2.92 & 1.42\\\\\n",
"\t 9 & 0 & 0.1 & 1.441305 & -0.63 & 1.441305 & 3 & 0.76 & 1.51\\\\\n",
"\t 10 & 4 & 4.1 & 2.989744 & 1.43 & 2.989744 & 12 & 1.93 & 1.90\\\\\n",
"\t 11 & 19 & 19.1 & 15.128913 & 2.97 & 15.128913 & 32 & 15.33 & 1.01\\\\\n",
"\t 12 & 0 & 0.1 & 1.329068 & -0.72 & 1.329068 & 2 & 0.68 & 1.36\\\\\n",
"\t 13 & 2 & 2.1 & 4.134935 & 0.90 & 4.134935 & 16 & 2.92 & 0.52\\\\\n",
"\t 14 & 0 & 0.1 & 1.993382 & -0.31 & 1.993382 & 7 & 1.15 & 2.29\\\\\n",
"\t 15 & 3 & 3.1 & 2.542240 & 1.11 & 2.542240 & 10 & 1.56 & 1.23\\\\\n",
"\t 16 & 1 & 1.1 & 3.812941 & 0.60 & 3.812941 & 15 & 2.63 & 1.51\\\\\n",
"\t 17 & 7 & 7.1 & 12.864424 & 2.10 & 12.864424 & 30 & 12.46 & 3.05\\\\\n",
"\t 18 & 1 & 1.1 & 1.563020 & 0.09 & 1.563020 & 4 & 0.84 & 0.22\\\\\n",
"\t 19 & 0 & 0.1 & 1.563020 & -0.55 & 1.563020 & 4 & 0.84 & 1.68\\\\\n",
"\t 20 & 7 & 7.1 & 5.273449 & 1.99 & 5.273449 & 19 & 3.98 & 2.06\\\\\n",
"\t 21 & 0 & 0.1 & 2.756927 & 0.01 & 2.756927 & 11 & 1.74 & 3.47\\\\\n",
"\t 22 & 0 & 0.1 & 1.130133 & -0.88 & 1.130133 & 0 & 0.55 & 1.11\\\\\n",
"\t 23 & 0 & 0.1 & 1.130133 & -0.88 & 1.130133 & 0 & 0.55 & 1.11\\\\\n",
"\t 24 & 0 & 0.1 & 1.130133 & -0.88 & 1.130133 & 0 & 0.55 & 1.11\\\\\n",
"\t 25 & 0 & 0.1 & 1.441305 & -0.63 & 1.441305 & 3 & 0.76 & 1.51\\\\\n",
"\t 26 & 0 & 0.1 & 2.161719 & -0.23 & 2.161719 & 8 & 1.27 & 2.54\\\\\n",
"\t 27 & 0 & 0.1 & 2.161719 & -0.23 & 2.161719 & 8 & 1.27 & 2.54\\\\\n",
"\t 28 & 0 & 0.1 & 3.812941 & 0.34 & 3.812941 & 15 & 2.63 & 5.26\\\\\n",
"\t 29 & 3 & 3.1 & 6.201720 & 1.31 & 6.201720 & 21 & 4.90 & 1.05\\\\\n",
"\t 30 & 4 & 4.1 & 7.909303 & 1.57 & 7.909303 & 24 & 6.69 & 1.46\\\\\n",
"\t 31 & 6 & 6.1 & 3.812941 & 1.91 & 3.812941 & 15 & 2.63 & 3.36\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 31 × 9\n",
"\n",
"| row &lt;int&gt; | yi &lt;int&gt; | uyumdegeri &lt;dbl&gt; | mu1 &lt;dbl&gt; | calismadegeri &lt;dbl&gt; | agirliklar &lt;dbl&gt; | x &lt;int&gt; | mu2i &lt;dbl&gt; | unit_dev &lt;dbl&gt; |\n",
"|---|---|---|---|---|---|---|---|---|\n",
"| 1 | 0 | 0.1 | 1.329068 | -0.72 | 1.329068 | 2 | 0.68 | 1.36 |\n",
"| 2 | 0 | 0.1 | 2.542240 | -0.07 | 2.542240 | 10 | 1.56 | 3.13 |\n",
"| 3 | 3 | 3.1 | 4.134935 | 1.14 | 4.134935 | 16 | 2.92 | 0.20 |\n",
"| 4 | 2 | 2.1 | 5.718781 | 1.09 | 5.718781 | 20 | 4.42 | 1.86 |\n",
"| 5 | 8 | 8.1 | 5.273449 | 2.18 | 5.273449 | 19 | 3.98 | 3.33 |\n",
"| 6 | 1 | 1.1 | 4.862797 | 0.79 | 4.862797 | 18 | 3.59 | 2.81 |\n",
"| 7 | 8 | 8.1 | 2.989744 | 2.77 | 2.989744 | 12 | 1.93 | 10.84 |\n",
"| 8 | 5 | 5.1 | 4.134935 | 1.63 | 4.134935 | 16 | 2.92 | 1.42 |\n",
"| 9 | 0 | 0.1 | 1.441305 | -0.63 | 1.441305 | 3 | 0.76 | 1.51 |\n",
"| 10 | 4 | 4.1 | 2.989744 | 1.43 | 2.989744 | 12 | 1.93 | 1.90 |\n",
"| 11 | 19 | 19.1 | 15.128913 | 2.97 | 15.128913 | 32 | 15.33 | 1.01 |\n",
"| 12 | 0 | 0.1 | 1.329068 | -0.72 | 1.329068 | 2 | 0.68 | 1.36 |\n",
"| 13 | 2 | 2.1 | 4.134935 | 0.90 | 4.134935 | 16 | 2.92 | 0.52 |\n",
"| 14 | 0 | 0.1 | 1.993382 | -0.31 | 1.993382 | 7 | 1.15 | 2.29 |\n",
"| 15 | 3 | 3.1 | 2.542240 | 1.11 | 2.542240 | 10 | 1.56 | 1.23 |\n",
"| 16 | 1 | 1.1 | 3.812941 | 0.60 | 3.812941 | 15 | 2.63 | 1.51 |\n",
"| 17 | 7 | 7.1 | 12.864424 | 2.10 | 12.864424 | 30 | 12.46 | 3.05 |\n",
"| 18 | 1 | 1.1 | 1.563020 | 0.09 | 1.563020 | 4 | 0.84 | 0.22 |\n",
"| 19 | 0 | 0.1 | 1.563020 | -0.55 | 1.563020 | 4 | 0.84 | 1.68 |\n",
"| 20 | 7 | 7.1 | 5.273449 | 1.99 | 5.273449 | 19 | 3.98 | 2.06 |\n",
"| 21 | 0 | 0.1 | 2.756927 | 0.01 | 2.756927 | 11 | 1.74 | 3.47 |\n",
"| 22 | 0 | 0.1 | 1.130133 | -0.88 | 1.130133 | 0 | 0.55 | 1.11 |\n",
"| 23 | 0 | 0.1 | 1.130133 | -0.88 | 1.130133 | 0 | 0.55 | 1.11 |\n",
"| 24 | 0 | 0.1 | 1.130133 | -0.88 | 1.130133 | 0 | 0.55 | 1.11 |\n",
"| 25 | 0 | 0.1 | 1.441305 | -0.63 | 1.441305 | 3 | 0.76 | 1.51 |\n",
"| 26 | 0 | 0.1 | 2.161719 | -0.23 | 2.161719 | 8 | 1.27 | 2.54 |\n",
"| 27 | 0 | 0.1 | 2.161719 | -0.23 | 2.161719 | 8 | 1.27 | 2.54 |\n",
"| 28 | 0 | 0.1 | 3.812941 | 0.34 | 3.812941 | 15 | 2.63 | 5.26 |\n",
"| 29 | 3 | 3.1 | 6.201720 | 1.31 | 6.201720 | 21 | 4.90 | 1.05 |\n",
"| 30 | 4 | 4.1 | 7.909303 | 1.57 | 7.909303 | 24 | 6.69 | 1.46 |\n",
"| 31 | 6 | 6.1 | 3.812941 | 1.91 | 3.812941 | 15 | 2.63 | 3.36 |\n",
"\n"
],
"text/plain": [
" row yi uyumdegeri mu1 calismadegeri agirliklar x mu2i unit_dev\n",
"1 1 0 0.1 1.329068 -0.72 1.329068 2 0.68 1.36 \n",
"2 2 0 0.1 2.542240 -0.07 2.542240 10 1.56 3.13 \n",
"3 3 3 3.1 4.134935 1.14 4.134935 16 2.92 0.20 \n",
"4 4 2 2.1 5.718781 1.09 5.718781 20 4.42 1.86 \n",
"5 5 8 8.1 5.273449 2.18 5.273449 19 3.98 3.33 \n",
"6 6 1 1.1 4.862797 0.79 4.862797 18 3.59 2.81 \n",
"7 7 8 8.1 2.989744 2.77 2.989744 12 1.93 10.84 \n",
"8 8 5 5.1 4.134935 1.63 4.134935 16 2.92 1.42 \n",
"9 9 0 0.1 1.441305 -0.63 1.441305 3 0.76 1.51 \n",
"10 10 4 4.1 2.989744 1.43 2.989744 12 1.93 1.90 \n",
"11 11 19 19.1 15.128913 2.97 15.128913 32 15.33 1.01 \n",
"12 12 0 0.1 1.329068 -0.72 1.329068 2 0.68 1.36 \n",
"13 13 2 2.1 4.134935 0.90 4.134935 16 2.92 0.52 \n",
"14 14 0 0.1 1.993382 -0.31 1.993382 7 1.15 2.29 \n",
"15 15 3 3.1 2.542240 1.11 2.542240 10 1.56 1.23 \n",
"16 16 1 1.1 3.812941 0.60 3.812941 15 2.63 1.51 \n",
"17 17 7 7.1 12.864424 2.10 12.864424 30 12.46 3.05 \n",
"18 18 1 1.1 1.563020 0.09 1.563020 4 0.84 0.22 \n",
"19 19 0 0.1 1.563020 -0.55 1.563020 4 0.84 1.68 \n",
"20 20 7 7.1 5.273449 1.99 5.273449 19 3.98 2.06 \n",
"21 21 0 0.1 2.756927 0.01 2.756927 11 1.74 3.47 \n",
"22 22 0 0.1 1.130133 -0.88 1.130133 0 0.55 1.11 \n",
"23 23 0 0.1 1.130133 -0.88 1.130133 0 0.55 1.11 \n",
"24 24 0 0.1 1.130133 -0.88 1.130133 0 0.55 1.11 \n",
"25 25 0 0.1 1.441305 -0.63 1.441305 3 0.76 1.51 \n",
"26 26 0 0.1 2.161719 -0.23 2.161719 8 1.27 2.54 \n",
"27 27 0 0.1 2.161719 -0.23 2.161719 8 1.27 2.54 \n",
"28 28 0 0.1 3.812941 0.34 3.812941 15 2.63 5.26 \n",
"29 29 3 3.1 6.201720 1.31 6.201720 21 4.90 1.05 \n",
"30 30 4 4.1 7.909303 1.57 7.909303 24 6.69 1.46 \n",
"31 31 6 6.1 3.812941 1.91 3.812941 15 2.63 3.36 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"yi = (nminer$Minerab)\n",
"etai <- (nminer$Minerab + 0.1)\n",
"\n",
"calculate_eta <- function(etai_){\n",
" log(etai_)\n",
"}\n",
"all_etai <- calculate_eta(etai)\n",
"\n",
"calculate_z1i <- function(a, b, c){\n",
" a + (b - c) / c\n",
"}\n",
"\n",
"calculate_w1i <- function(a){\n",
" a\n",
"}\n",
"\n",
"all_wieghts <- calculate_w1i(etai)\n",
"\n",
"x <- (nminer$Eucs)\n",
"\n",
"calculate_mu1i <- function(coef_, x_){\n",
" exp(coef_[1] + (coef_[2] * x_))\n",
"}\n",
"\n",
"all_mu2i <- calculate_mu1i(coef_iter2, x)\n",
"\n",
"calculate_unit_dev <- function(y_, etai_, mu1_){\n",
" l_ <- etai_/mu1_\n",
" 2 * (y_ * log(l_) - (y_ - mu1_))\n",
"}\n",
"\n",
"all_z2i <- calculate_z1i(log(all_mu1i), yi, all_mu1i)\n",
"\n",
"all_unit_dev2 <- calculate_unit_dev(yi, etai, all_mu2i)\n",
"\n",
"tablo_iter2 = data.frame(row = seq.int(yi), yi = yi, uyumdegeri = round(etai, 3), mu1 = all_mu1i, calismadegeri = round(all_z2i, 2), agirliklar = all_mu1i, x = x, mu2i = round(all_mu2i, 2), unit_dev = round(all_unit_dev2, 2))\n",
"\n",
"tablo_iter2"
]
},
{
"cell_type": "code",
"execution_count": 138,
"id": "f8aa3fd2",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"dataframe\">\n",
"<caption>A data.frame: 31 × 9</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>row</th><th scope=col>yi</th><th scope=col>uyumdegeri</th><th scope=col>mu2</th><th scope=col>calismadegeri</th><th scope=col>agirliklar</th><th scope=col>x</th><th scope=col>mu3i</th><th scope=col>unit_dev</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 1</td><td> 0</td><td> 0.1</td><td> 0.6822849</td><td>-1.38</td><td> 0.6822849</td><td> 2</td><td> 0.53</td><td> 1.07</td></tr>\n",
"\t<tr><td> 2</td><td> 0</td><td> 0.1</td><td> 1.5646373</td><td>-0.55</td><td> 1.5646373</td><td>10</td><td> 1.32</td><td> 2.64</td></tr>\n",
"\t<tr><td> 3</td><td> 3</td><td> 3.1</td><td> 2.9157453</td><td> 0.80</td><td> 2.9157453</td><td>16</td><td> 2.61</td><td> 0.25</td></tr>\n",
"\t<tr><td> 4</td><td> 2</td><td> 2.1</td><td> 4.4154354</td><td> 0.83</td><td> 4.4154354</td><td>20</td><td> 4.10</td><td> 1.52</td></tr>\n",
"\t<tr><td> 5</td><td> 8</td><td> 8.1</td><td> 3.9803160</td><td> 1.90</td><td> 3.9803160</td><td>19</td><td> 3.66</td><td> 4.03</td></tr>\n",
"\t<tr><td> 6</td><td> 1</td><td> 1.1</td><td> 3.5880754</td><td> 0.48</td><td> 3.5880754</td><td>18</td><td> 3.27</td><td> 2.36</td></tr>\n",
"\t<tr><td> 7</td><td> 8</td><td> 8.1</td><td> 1.9254207</td><td> 2.33</td><td> 1.9254207</td><td>12</td><td> 1.66</td><td>12.70</td></tr>\n",
"\t<tr><td> 8</td><td> 5</td><td> 5.1</td><td> 2.9157453</td><td> 1.28</td><td> 2.9157453</td><td>16</td><td> 2.61</td><td> 1.93</td></tr>\n",
"\t<tr><td> 9</td><td> 0</td><td> 0.1</td><td> 0.7568708</td><td>-1.28</td><td> 0.7568708</td><td> 3</td><td> 0.60</td><td> 1.20</td></tr>\n",
"\t<tr><td>10</td><td> 4</td><td> 4.1</td><td> 1.9254207</td><td> 0.99</td><td> 1.9254207</td><td>12</td><td> 1.66</td><td> 2.56</td></tr>\n",
"\t<tr><td>11</td><td>19</td><td>19.1</td><td>15.3336441</td><td> 2.99</td><td>15.3336441</td><td>32</td><td>15.93</td><td> 0.76</td></tr>\n",
"\t<tr><td>12</td><td> 0</td><td> 0.1</td><td> 0.6822849</td><td>-1.38</td><td> 0.6822849</td><td> 2</td><td> 0.53</td><td> 1.07</td></tr>\n",
"\t<tr><td>13</td><td> 2</td><td> 2.1</td><td> 2.9157453</td><td> 0.55</td><td> 2.9157453</td><td>16</td><td> 2.61</td><td> 0.35</td></tr>\n",
"\t<tr><td>14</td><td> 0</td><td> 0.1</td><td> 1.1461612</td><td>-0.86</td><td> 1.1461612</td><td> 7</td><td> 0.94</td><td> 1.88</td></tr>\n",
"\t<tr><td>15</td><td> 3</td><td> 3.1</td><td> 1.5646373</td><td> 0.63</td><td> 1.5646373</td><td>10</td><td> 1.32</td><td> 1.76</td></tr>\n",
"\t<tr><td>16</td><td> 1</td><td> 1.1</td><td> 2.6284130</td><td> 0.23</td><td> 2.6284130</td><td>15</td><td> 2.33</td><td> 1.16</td></tr>\n",
"\t<tr><td>17</td><td> 7</td><td> 7.1</td><td>12.4604408</td><td> 2.07</td><td>12.4604408</td><td>30</td><td>12.70</td><td> 3.26</td></tr>\n",
"\t<tr><td>18</td><td> 1</td><td> 1.1</td><td> 0.8396103</td><td>-0.54</td><td> 0.8396103</td><td> 4</td><td> 0.67</td><td> 0.33</td></tr>\n",
"\t<tr><td>19</td><td> 0</td><td> 0.1</td><td> 0.8396103</td><td>-1.17</td><td> 0.8396103</td><td> 4</td><td> 0.67</td><td> 1.34</td></tr>\n",
"\t<tr><td>20</td><td> 7</td><td> 7.1</td><td> 3.9803160</td><td> 1.71</td><td> 3.9803160</td><td>19</td><td> 3.66</td><td> 2.60</td></tr>\n",
"\t<tr><td>21</td><td> 0</td><td> 0.1</td><td> 1.7356800</td><td>-0.45</td><td> 1.7356800</td><td>11</td><td> 1.48</td><td> 2.96</td></tr>\n",
"\t<tr><td>22</td><td> 0</td><td> 0.1</td><td> 0.5544390</td><td>-1.59</td><td> 0.5544390</td><td> 0</td><td> 0.43</td><td> 0.85</td></tr>\n",
"\t<tr><td>23</td><td> 0</td><td> 0.1</td><td> 0.5544390</td><td>-1.59</td><td> 0.5544390</td><td> 0</td><td> 0.43</td><td> 0.85</td></tr>\n",
"\t<tr><td>24</td><td> 0</td><td> 0.1</td><td> 0.5544390</td><td>-1.59</td><td> 0.5544390</td><td> 0</td><td> 0.43</td><td> 0.85</td></tr>\n",
"\t<tr><td>25</td><td> 0</td><td> 0.1</td><td> 0.7568708</td><td>-1.28</td><td> 0.7568708</td><td> 3</td><td> 0.60</td><td> 1.20</td></tr>\n",
"\t<tr><td>26</td><td> 0</td><td> 0.1</td><td> 1.2714571</td><td>-0.76</td><td> 1.2714571</td><td> 8</td><td> 1.05</td><td> 2.11</td></tr>\n",
"\t<tr><td>27</td><td> 0</td><td> 0.1</td><td> 1.2714571</td><td>-0.76</td><td> 1.2714571</td><td> 8</td><td> 1.05</td><td> 2.11</td></tr>\n",
"\t<tr><td>28</td><td> 0</td><td> 0.1</td><td> 2.6284130</td><td>-0.03</td><td> 2.6284130</td><td>15</td><td> 2.33</td><td> 4.66</td></tr>\n",
"\t<tr><td>29</td><td> 3</td><td> 3.1</td><td> 4.8981212</td><td> 1.07</td><td> 4.8981212</td><td>21</td><td> 4.59</td><td> 0.82</td></tr>\n",
"\t<tr><td>30</td><td> 4</td><td> 4.1</td><td> 6.6864790</td><td> 1.41</td><td> 6.6864790</td><td>24</td><td> 6.44</td><td> 1.27</td></tr>\n",
"\t<tr><td>31</td><td> 6</td><td> 6.1</td><td> 2.6284130</td><td> 1.54</td><td> 2.6284130</td><td>15</td><td> 2.33</td><td> 4.22</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 31 × 9\n",
"\\begin{tabular}{lllllllll}\n",
" row & yi & uyumdegeri & mu2 & calismadegeri & agirliklar & x & mu3i & unit\\_dev\\\\\n",
" <int> & <int> & <dbl> & <dbl> & <dbl> & <dbl> & <int> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 1 & 0 & 0.1 & 0.6822849 & -1.38 & 0.6822849 & 2 & 0.53 & 1.07\\\\\n",
"\t 2 & 0 & 0.1 & 1.5646373 & -0.55 & 1.5646373 & 10 & 1.32 & 2.64\\\\\n",
"\t 3 & 3 & 3.1 & 2.9157453 & 0.80 & 2.9157453 & 16 & 2.61 & 0.25\\\\\n",
"\t 4 & 2 & 2.1 & 4.4154354 & 0.83 & 4.4154354 & 20 & 4.10 & 1.52\\\\\n",
"\t 5 & 8 & 8.1 & 3.9803160 & 1.90 & 3.9803160 & 19 & 3.66 & 4.03\\\\\n",
"\t 6 & 1 & 1.1 & 3.5880754 & 0.48 & 3.5880754 & 18 & 3.27 & 2.36\\\\\n",
"\t 7 & 8 & 8.1 & 1.9254207 & 2.33 & 1.9254207 & 12 & 1.66 & 12.70\\\\\n",
"\t 8 & 5 & 5.1 & 2.9157453 & 1.28 & 2.9157453 & 16 & 2.61 & 1.93\\\\\n",
"\t 9 & 0 & 0.1 & 0.7568708 & -1.28 & 0.7568708 & 3 & 0.60 & 1.20\\\\\n",
"\t 10 & 4 & 4.1 & 1.9254207 & 0.99 & 1.9254207 & 12 & 1.66 & 2.56\\\\\n",
"\t 11 & 19 & 19.1 & 15.3336441 & 2.99 & 15.3336441 & 32 & 15.93 & 0.76\\\\\n",
"\t 12 & 0 & 0.1 & 0.6822849 & -1.38 & 0.6822849 & 2 & 0.53 & 1.07\\\\\n",
"\t 13 & 2 & 2.1 & 2.9157453 & 0.55 & 2.9157453 & 16 & 2.61 & 0.35\\\\\n",
"\t 14 & 0 & 0.1 & 1.1461612 & -0.86 & 1.1461612 & 7 & 0.94 & 1.88\\\\\n",
"\t 15 & 3 & 3.1 & 1.5646373 & 0.63 & 1.5646373 & 10 & 1.32 & 1.76\\\\\n",
"\t 16 & 1 & 1.1 & 2.6284130 & 0.23 & 2.6284130 & 15 & 2.33 & 1.16\\\\\n",
"\t 17 & 7 & 7.1 & 12.4604408 & 2.07 & 12.4604408 & 30 & 12.70 & 3.26\\\\\n",
"\t 18 & 1 & 1.1 & 0.8396103 & -0.54 & 0.8396103 & 4 & 0.67 & 0.33\\\\\n",
"\t 19 & 0 & 0.1 & 0.8396103 & -1.17 & 0.8396103 & 4 & 0.67 & 1.34\\\\\n",
"\t 20 & 7 & 7.1 & 3.9803160 & 1.71 & 3.9803160 & 19 & 3.66 & 2.60\\\\\n",
"\t 21 & 0 & 0.1 & 1.7356800 & -0.45 & 1.7356800 & 11 & 1.48 & 2.96\\\\\n",
"\t 22 & 0 & 0.1 & 0.5544390 & -1.59 & 0.5544390 & 0 & 0.43 & 0.85\\\\\n",
"\t 23 & 0 & 0.1 & 0.5544390 & -1.59 & 0.5544390 & 0 & 0.43 & 0.85\\\\\n",
"\t 24 & 0 & 0.1 & 0.5544390 & -1.59 & 0.5544390 & 0 & 0.43 & 0.85\\\\\n",
"\t 25 & 0 & 0.1 & 0.7568708 & -1.28 & 0.7568708 & 3 & 0.60 & 1.20\\\\\n",
"\t 26 & 0 & 0.1 & 1.2714571 & -0.76 & 1.2714571 & 8 & 1.05 & 2.11\\\\\n",
"\t 27 & 0 & 0.1 & 1.2714571 & -0.76 & 1.2714571 & 8 & 1.05 & 2.11\\\\\n",
"\t 28 & 0 & 0.1 & 2.6284130 & -0.03 & 2.6284130 & 15 & 2.33 & 4.66\\\\\n",
"\t 29 & 3 & 3.1 & 4.8981212 & 1.07 & 4.8981212 & 21 & 4.59 & 0.82\\\\\n",
"\t 30 & 4 & 4.1 & 6.6864790 & 1.41 & 6.6864790 & 24 & 6.44 & 1.27\\\\\n",
"\t 31 & 6 & 6.1 & 2.6284130 & 1.54 & 2.6284130 & 15 & 2.33 & 4.22\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 31 × 9\n",
"\n",
"| row &lt;int&gt; | yi &lt;int&gt; | uyumdegeri &lt;dbl&gt; | mu2 &lt;dbl&gt; | calismadegeri &lt;dbl&gt; | agirliklar &lt;dbl&gt; | x &lt;int&gt; | mu3i &lt;dbl&gt; | unit_dev &lt;dbl&gt; |\n",
"|---|---|---|---|---|---|---|---|---|\n",
"| 1 | 0 | 0.1 | 0.6822849 | -1.38 | 0.6822849 | 2 | 0.53 | 1.07 |\n",
"| 2 | 0 | 0.1 | 1.5646373 | -0.55 | 1.5646373 | 10 | 1.32 | 2.64 |\n",
"| 3 | 3 | 3.1 | 2.9157453 | 0.80 | 2.9157453 | 16 | 2.61 | 0.25 |\n",
"| 4 | 2 | 2.1 | 4.4154354 | 0.83 | 4.4154354 | 20 | 4.10 | 1.52 |\n",
"| 5 | 8 | 8.1 | 3.9803160 | 1.90 | 3.9803160 | 19 | 3.66 | 4.03 |\n",
"| 6 | 1 | 1.1 | 3.5880754 | 0.48 | 3.5880754 | 18 | 3.27 | 2.36 |\n",
"| 7 | 8 | 8.1 | 1.9254207 | 2.33 | 1.9254207 | 12 | 1.66 | 12.70 |\n",
"| 8 | 5 | 5.1 | 2.9157453 | 1.28 | 2.9157453 | 16 | 2.61 | 1.93 |\n",
"| 9 | 0 | 0.1 | 0.7568708 | -1.28 | 0.7568708 | 3 | 0.60 | 1.20 |\n",
"| 10 | 4 | 4.1 | 1.9254207 | 0.99 | 1.9254207 | 12 | 1.66 | 2.56 |\n",
"| 11 | 19 | 19.1 | 15.3336441 | 2.99 | 15.3336441 | 32 | 15.93 | 0.76 |\n",
"| 12 | 0 | 0.1 | 0.6822849 | -1.38 | 0.6822849 | 2 | 0.53 | 1.07 |\n",
"| 13 | 2 | 2.1 | 2.9157453 | 0.55 | 2.9157453 | 16 | 2.61 | 0.35 |\n",
"| 14 | 0 | 0.1 | 1.1461612 | -0.86 | 1.1461612 | 7 | 0.94 | 1.88 |\n",
"| 15 | 3 | 3.1 | 1.5646373 | 0.63 | 1.5646373 | 10 | 1.32 | 1.76 |\n",
"| 16 | 1 | 1.1 | 2.6284130 | 0.23 | 2.6284130 | 15 | 2.33 | 1.16 |\n",
"| 17 | 7 | 7.1 | 12.4604408 | 2.07 | 12.4604408 | 30 | 12.70 | 3.26 |\n",
"| 18 | 1 | 1.1 | 0.8396103 | -0.54 | 0.8396103 | 4 | 0.67 | 0.33 |\n",
"| 19 | 0 | 0.1 | 0.8396103 | -1.17 | 0.8396103 | 4 | 0.67 | 1.34 |\n",
"| 20 | 7 | 7.1 | 3.9803160 | 1.71 | 3.9803160 | 19 | 3.66 | 2.60 |\n",
"| 21 | 0 | 0.1 | 1.7356800 | -0.45 | 1.7356800 | 11 | 1.48 | 2.96 |\n",
"| 22 | 0 | 0.1 | 0.5544390 | -1.59 | 0.5544390 | 0 | 0.43 | 0.85 |\n",
"| 23 | 0 | 0.1 | 0.5544390 | -1.59 | 0.5544390 | 0 | 0.43 | 0.85 |\n",
"| 24 | 0 | 0.1 | 0.5544390 | -1.59 | 0.5544390 | 0 | 0.43 | 0.85 |\n",
"| 25 | 0 | 0.1 | 0.7568708 | -1.28 | 0.7568708 | 3 | 0.60 | 1.20 |\n",
"| 26 | 0 | 0.1 | 1.2714571 | -0.76 | 1.2714571 | 8 | 1.05 | 2.11 |\n",
"| 27 | 0 | 0.1 | 1.2714571 | -0.76 | 1.2714571 | 8 | 1.05 | 2.11 |\n",
"| 28 | 0 | 0.1 | 2.6284130 | -0.03 | 2.6284130 | 15 | 2.33 | 4.66 |\n",
"| 29 | 3 | 3.1 | 4.8981212 | 1.07 | 4.8981212 | 21 | 4.59 | 0.82 |\n",
"| 30 | 4 | 4.1 | 6.6864790 | 1.41 | 6.6864790 | 24 | 6.44 | 1.27 |\n",
"| 31 | 6 | 6.1 | 2.6284130 | 1.54 | 2.6284130 | 15 | 2.33 | 4.22 |\n",
"\n"
],
"text/plain": [
" row yi uyumdegeri mu2 calismadegeri agirliklar x mu3i unit_dev\n",
"1 1 0 0.1 0.6822849 -1.38 0.6822849 2 0.53 1.07 \n",
"2 2 0 0.1 1.5646373 -0.55 1.5646373 10 1.32 2.64 \n",
"3 3 3 3.1 2.9157453 0.80 2.9157453 16 2.61 0.25 \n",
"4 4 2 2.1 4.4154354 0.83 4.4154354 20 4.10 1.52 \n",
"5 5 8 8.1 3.9803160 1.90 3.9803160 19 3.66 4.03 \n",
"6 6 1 1.1 3.5880754 0.48 3.5880754 18 3.27 2.36 \n",
"7 7 8 8.1 1.9254207 2.33 1.9254207 12 1.66 12.70 \n",
"8 8 5 5.1 2.9157453 1.28 2.9157453 16 2.61 1.93 \n",
"9 9 0 0.1 0.7568708 -1.28 0.7568708 3 0.60 1.20 \n",
"10 10 4 4.1 1.9254207 0.99 1.9254207 12 1.66 2.56 \n",
"11 11 19 19.1 15.3336441 2.99 15.3336441 32 15.93 0.76 \n",
"12 12 0 0.1 0.6822849 -1.38 0.6822849 2 0.53 1.07 \n",
"13 13 2 2.1 2.9157453 0.55 2.9157453 16 2.61 0.35 \n",
"14 14 0 0.1 1.1461612 -0.86 1.1461612 7 0.94 1.88 \n",
"15 15 3 3.1 1.5646373 0.63 1.5646373 10 1.32 1.76 \n",
"16 16 1 1.1 2.6284130 0.23 2.6284130 15 2.33 1.16 \n",
"17 17 7 7.1 12.4604408 2.07 12.4604408 30 12.70 3.26 \n",
"18 18 1 1.1 0.8396103 -0.54 0.8396103 4 0.67 0.33 \n",
"19 19 0 0.1 0.8396103 -1.17 0.8396103 4 0.67 1.34 \n",
"20 20 7 7.1 3.9803160 1.71 3.9803160 19 3.66 2.60 \n",
"21 21 0 0.1 1.7356800 -0.45 1.7356800 11 1.48 2.96 \n",
"22 22 0 0.1 0.5544390 -1.59 0.5544390 0 0.43 0.85 \n",
"23 23 0 0.1 0.5544390 -1.59 0.5544390 0 0.43 0.85 \n",
"24 24 0 0.1 0.5544390 -1.59 0.5544390 0 0.43 0.85 \n",
"25 25 0 0.1 0.7568708 -1.28 0.7568708 3 0.60 1.20 \n",
"26 26 0 0.1 1.2714571 -0.76 1.2714571 8 1.05 2.11 \n",
"27 27 0 0.1 1.2714571 -0.76 1.2714571 8 1.05 2.11 \n",
"28 28 0 0.1 2.6284130 -0.03 2.6284130 15 2.33 4.66 \n",
"29 29 3 3.1 4.8981212 1.07 4.8981212 21 4.59 0.82 \n",
"30 30 4 4.1 6.6864790 1.41 6.6864790 24 6.44 1.27 \n",
"31 31 6 6.1 2.6284130 1.54 2.6284130 15 2.33 4.22 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"all_z3i <- calculate_z1i(log(all_mu2i), yi, all_mu1i)\n",
"\n",
"all_mu3i <- calculate_mu1i(coef_iter3, x)\n",
"\n",
"all_unit_dev3 <- calculate_unit_dev(yi, etai, all_mu3i)\n",
"\n",
"tablo_iter3 = data.frame(row = seq.int(yi), yi = yi, uyumdegeri = round(etai, 3), mu2 = all_mu2i, calismadegeri = round(all_z3i, 2), agirliklar = all_mu2i, x = x, mu3i = round(all_mu3i, 2), unit_dev = round(all_unit_dev3, 2))\n",
"\n",
"tablo_iter3"
]
},
{
"cell_type": "code",
"execution_count": 147,
"id": "0c3d08ac",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table class=\"dataframe\">\n",
"<caption>A data.frame: 31 × 9</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>row</th><th scope=col>yi</th><th scope=col>uyumdegeri</th><th scope=col>mu3</th><th scope=col>calismadegeri</th><th scope=col>agirliklar</th><th scope=col>x</th><th scope=col>mu4i</th><th scope=col>unit_dev</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 1</td><td> 0</td><td> 0.1</td><td> 0.5348672</td><td>-1.63</td><td> 0.5348672</td><td> 2</td><td> 0.52</td><td> 1.05</td></tr>\n",
"\t<tr><td> 2</td><td> 0</td><td> 0.1</td><td> 1.3221327</td><td>-0.72</td><td> 1.3221327</td><td>10</td><td> 1.30</td><td> 2.60</td></tr>\n",
"\t<tr><td> 3</td><td> 3</td><td> 3.1</td><td> 2.6064348</td><td> 1.11</td><td> 2.6064348</td><td>16</td><td> 2.58</td><td> 0.26</td></tr>\n",
"\t<tr><td> 4</td><td> 2</td><td> 2.1</td><td> 4.0979005</td><td> 0.90</td><td> 4.0979005</td><td>20</td><td> 4.07</td><td> 1.49</td></tr>\n",
"\t<tr><td> 5</td><td> 8</td><td> 8.1</td><td> 3.6595929</td><td> 2.48</td><td> 3.6595929</td><td>19</td><td> 3.63</td><td> 4.10</td></tr>\n",
"\t<tr><td> 6</td><td> 1</td><td> 1.1</td><td> 3.2681662</td><td> 0.49</td><td> 3.2681662</td><td>18</td><td> 3.24</td><td> 2.32</td></tr>\n",
"\t<tr><td> 7</td><td> 8</td><td> 8.1</td><td> 1.6578007</td><td> 4.33</td><td> 1.6578007</td><td>12</td><td> 1.64</td><td>12.87</td></tr>\n",
"\t<tr><td> 8</td><td> 5</td><td> 5.1</td><td> 2.6064348</td><td> 1.88</td><td> 2.6064348</td><td>16</td><td> 2.58</td><td> 1.98</td></tr>\n",
"\t<tr><td> 9</td><td> 0</td><td> 0.1</td><td> 0.5989280</td><td>-1.51</td><td> 0.5989280</td><td> 3</td><td> 0.59</td><td> 1.17</td></tr>\n",
"\t<tr><td>10</td><td> 4</td><td> 4.1</td><td> 1.6578007</td><td> 1.92</td><td> 1.6578007</td><td>12</td><td> 1.64</td><td> 2.62</td></tr>\n",
"\t<tr><td>11</td><td>19</td><td>19.1</td><td>15.9259380</td><td> 2.96</td><td>15.9259380</td><td>32</td><td>15.98</td><td> 0.74</td></tr>\n",
"\t<tr><td>12</td><td> 0</td><td> 0.1</td><td> 0.5348672</td><td>-1.63</td><td> 0.5348672</td><td> 2</td><td> 0.52</td><td> 1.05</td></tr>\n",
"\t<tr><td>13</td><td> 2</td><td> 2.1</td><td> 2.6064348</td><td> 0.73</td><td> 2.6064348</td><td>16</td><td> 2.58</td><td> 0.34</td></tr>\n",
"\t<tr><td>14</td><td> 0</td><td> 0.1</td><td> 0.9416492</td><td>-1.06</td><td> 0.9416492</td><td> 7</td><td> 0.92</td><td> 1.85</td></tr>\n",
"\t<tr><td>15</td><td> 3</td><td> 3.1</td><td> 1.3221327</td><td> 1.55</td><td> 1.3221327</td><td>10</td><td> 1.30</td><td> 1.81</td></tr>\n",
"\t<tr><td>16</td><td> 1</td><td> 1.1</td><td> 2.3276530</td><td> 0.27</td><td> 2.3276530</td><td>15</td><td> 2.30</td><td> 1.13</td></tr>\n",
"\t<tr><td>17</td><td> 7</td><td> 7.1</td><td>12.7012876</td><td> 2.09</td><td>12.7012876</td><td>30</td><td>12.72</td><td> 3.28</td></tr>\n",
"\t<tr><td>18</td><td> 1</td><td> 1.1</td><td> 0.6706613</td><td> 0.09</td><td> 0.6706613</td><td> 4</td><td> 0.66</td><td> 0.34</td></tr>\n",
"\t<tr><td>19</td><td> 0</td><td> 0.1</td><td> 0.6706613</td><td>-1.40</td><td> 0.6706613</td><td> 4</td><td> 0.66</td><td> 1.31</td></tr>\n",
"\t<tr><td>20</td><td> 7</td><td> 7.1</td><td> 3.6595929</td><td> 2.21</td><td> 3.6595929</td><td>19</td><td> 3.63</td><td> 2.65</td></tr>\n",
"\t<tr><td>21</td><td> 0</td><td> 0.1</td><td> 1.4804839</td><td>-0.61</td><td> 1.4804839</td><td>11</td><td> 1.46</td><td> 2.92</td></tr>\n",
"\t<tr><td>22</td><td> 0</td><td> 0.1</td><td> 0.4265684</td><td>-1.85</td><td> 0.4265684</td><td> 0</td><td> 0.42</td><td> 0.83</td></tr>\n",
"\t<tr><td>23</td><td> 0</td><td> 0.1</td><td> 0.4265684</td><td>-1.85</td><td> 0.4265684</td><td> 0</td><td> 0.42</td><td> 0.83</td></tr>\n",
"\t<tr><td>24</td><td> 0</td><td> 0.1</td><td> 0.4265684</td><td>-1.85</td><td> 0.4265684</td><td> 0</td><td> 0.42</td><td> 0.83</td></tr>\n",
"\t<tr><td>25</td><td> 0</td><td> 0.1</td><td> 0.5989280</td><td>-1.51</td><td> 0.5989280</td><td> 3</td><td> 0.59</td><td> 1.17</td></tr>\n",
"\t<tr><td>26</td><td> 0</td><td> 0.1</td><td> 1.0544301</td><td>-0.95</td><td> 1.0544301</td><td> 8</td><td> 1.04</td><td> 2.07</td></tr>\n",
"\t<tr><td>27</td><td> 0</td><td> 0.1</td><td> 1.0544301</td><td>-0.95</td><td> 1.0544301</td><td> 8</td><td> 1.04</td><td> 2.07</td></tr>\n",
"\t<tr><td>28</td><td> 0</td><td> 0.1</td><td> 2.3276530</td><td>-0.16</td><td> 2.3276530</td><td>15</td><td> 2.30</td><td> 4.60</td></tr>\n",
"\t<tr><td>29</td><td> 3</td><td> 3.1</td><td> 4.5887040</td><td> 1.18</td><td> 4.5887040</td><td>21</td><td> 4.56</td><td> 0.80</td></tr>\n",
"\t<tr><td>30</td><td> 4</td><td> 4.1</td><td> 6.4428192</td><td> 1.48</td><td> 6.4428192</td><td>24</td><td> 6.42</td><td> 1.25</td></tr>\n",
"\t<tr><td>31</td><td> 6</td><td> 6.1</td><td> 2.3276530</td><td> 2.42</td><td> 2.3276530</td><td>15</td><td> 2.30</td><td> 4.30</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 31 × 9\n",
"\\begin{tabular}{lllllllll}\n",
" row & yi & uyumdegeri & mu3 & calismadegeri & agirliklar & x & mu4i & unit\\_dev\\\\\n",
" <int> & <int> & <dbl> & <dbl> & <dbl> & <dbl> & <int> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 1 & 0 & 0.1 & 0.5348672 & -1.63 & 0.5348672 & 2 & 0.52 & 1.05\\\\\n",
"\t 2 & 0 & 0.1 & 1.3221327 & -0.72 & 1.3221327 & 10 & 1.30 & 2.60\\\\\n",
"\t 3 & 3 & 3.1 & 2.6064348 & 1.11 & 2.6064348 & 16 & 2.58 & 0.26\\\\\n",
"\t 4 & 2 & 2.1 & 4.0979005 & 0.90 & 4.0979005 & 20 & 4.07 & 1.49\\\\\n",
"\t 5 & 8 & 8.1 & 3.6595929 & 2.48 & 3.6595929 & 19 & 3.63 & 4.10\\\\\n",
"\t 6 & 1 & 1.1 & 3.2681662 & 0.49 & 3.2681662 & 18 & 3.24 & 2.32\\\\\n",
"\t 7 & 8 & 8.1 & 1.6578007 & 4.33 & 1.6578007 & 12 & 1.64 & 12.87\\\\\n",
"\t 8 & 5 & 5.1 & 2.6064348 & 1.88 & 2.6064348 & 16 & 2.58 & 1.98\\\\\n",
"\t 9 & 0 & 0.1 & 0.5989280 & -1.51 & 0.5989280 & 3 & 0.59 & 1.17\\\\\n",
"\t 10 & 4 & 4.1 & 1.6578007 & 1.92 & 1.6578007 & 12 & 1.64 & 2.62\\\\\n",
"\t 11 & 19 & 19.1 & 15.9259380 & 2.96 & 15.9259380 & 32 & 15.98 & 0.74\\\\\n",
"\t 12 & 0 & 0.1 & 0.5348672 & -1.63 & 0.5348672 & 2 & 0.52 & 1.05\\\\\n",
"\t 13 & 2 & 2.1 & 2.6064348 & 0.73 & 2.6064348 & 16 & 2.58 & 0.34\\\\\n",
"\t 14 & 0 & 0.1 & 0.9416492 & -1.06 & 0.9416492 & 7 & 0.92 & 1.85\\\\\n",
"\t 15 & 3 & 3.1 & 1.3221327 & 1.55 & 1.3221327 & 10 & 1.30 & 1.81\\\\\n",
"\t 16 & 1 & 1.1 & 2.3276530 & 0.27 & 2.3276530 & 15 & 2.30 & 1.13\\\\\n",
"\t 17 & 7 & 7.1 & 12.7012876 & 2.09 & 12.7012876 & 30 & 12.72 & 3.28\\\\\n",
"\t 18 & 1 & 1.1 & 0.6706613 & 0.09 & 0.6706613 & 4 & 0.66 & 0.34\\\\\n",
"\t 19 & 0 & 0.1 & 0.6706613 & -1.40 & 0.6706613 & 4 & 0.66 & 1.31\\\\\n",
"\t 20 & 7 & 7.1 & 3.6595929 & 2.21 & 3.6595929 & 19 & 3.63 & 2.65\\\\\n",
"\t 21 & 0 & 0.1 & 1.4804839 & -0.61 & 1.4804839 & 11 & 1.46 & 2.92\\\\\n",
"\t 22 & 0 & 0.1 & 0.4265684 & -1.85 & 0.4265684 & 0 & 0.42 & 0.83\\\\\n",
"\t 23 & 0 & 0.1 & 0.4265684 & -1.85 & 0.4265684 & 0 & 0.42 & 0.83\\\\\n",
"\t 24 & 0 & 0.1 & 0.4265684 & -1.85 & 0.4265684 & 0 & 0.42 & 0.83\\\\\n",
"\t 25 & 0 & 0.1 & 0.5989280 & -1.51 & 0.5989280 & 3 & 0.59 & 1.17\\\\\n",
"\t 26 & 0 & 0.1 & 1.0544301 & -0.95 & 1.0544301 & 8 & 1.04 & 2.07\\\\\n",
"\t 27 & 0 & 0.1 & 1.0544301 & -0.95 & 1.0544301 & 8 & 1.04 & 2.07\\\\\n",
"\t 28 & 0 & 0.1 & 2.3276530 & -0.16 & 2.3276530 & 15 & 2.30 & 4.60\\\\\n",
"\t 29 & 3 & 3.1 & 4.5887040 & 1.18 & 4.5887040 & 21 & 4.56 & 0.80\\\\\n",
"\t 30 & 4 & 4.1 & 6.4428192 & 1.48 & 6.4428192 & 24 & 6.42 & 1.25\\\\\n",
"\t 31 & 6 & 6.1 & 2.3276530 & 2.42 & 2.3276530 & 15 & 2.30 & 4.30\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 31 × 9\n",
"\n",
"| row &lt;int&gt; | yi &lt;int&gt; | uyumdegeri &lt;dbl&gt; | mu3 &lt;dbl&gt; | calismadegeri &lt;dbl&gt; | agirliklar &lt;dbl&gt; | x &lt;int&gt; | mu4i &lt;dbl&gt; | unit_dev &lt;dbl&gt; |\n",
"|---|---|---|---|---|---|---|---|---|\n",
"| 1 | 0 | 0.1 | 0.5348672 | -1.63 | 0.5348672 | 2 | 0.52 | 1.05 |\n",
"| 2 | 0 | 0.1 | 1.3221327 | -0.72 | 1.3221327 | 10 | 1.30 | 2.60 |\n",
"| 3 | 3 | 3.1 | 2.6064348 | 1.11 | 2.6064348 | 16 | 2.58 | 0.26 |\n",
"| 4 | 2 | 2.1 | 4.0979005 | 0.90 | 4.0979005 | 20 | 4.07 | 1.49 |\n",
"| 5 | 8 | 8.1 | 3.6595929 | 2.48 | 3.6595929 | 19 | 3.63 | 4.10 |\n",
"| 6 | 1 | 1.1 | 3.2681662 | 0.49 | 3.2681662 | 18 | 3.24 | 2.32 |\n",
"| 7 | 8 | 8.1 | 1.6578007 | 4.33 | 1.6578007 | 12 | 1.64 | 12.87 |\n",
"| 8 | 5 | 5.1 | 2.6064348 | 1.88 | 2.6064348 | 16 | 2.58 | 1.98 |\n",
"| 9 | 0 | 0.1 | 0.5989280 | -1.51 | 0.5989280 | 3 | 0.59 | 1.17 |\n",
"| 10 | 4 | 4.1 | 1.6578007 | 1.92 | 1.6578007 | 12 | 1.64 | 2.62 |\n",
"| 11 | 19 | 19.1 | 15.9259380 | 2.96 | 15.9259380 | 32 | 15.98 | 0.74 |\n",
"| 12 | 0 | 0.1 | 0.5348672 | -1.63 | 0.5348672 | 2 | 0.52 | 1.05 |\n",
"| 13 | 2 | 2.1 | 2.6064348 | 0.73 | 2.6064348 | 16 | 2.58 | 0.34 |\n",
"| 14 | 0 | 0.1 | 0.9416492 | -1.06 | 0.9416492 | 7 | 0.92 | 1.85 |\n",
"| 15 | 3 | 3.1 | 1.3221327 | 1.55 | 1.3221327 | 10 | 1.30 | 1.81 |\n",
"| 16 | 1 | 1.1 | 2.3276530 | 0.27 | 2.3276530 | 15 | 2.30 | 1.13 |\n",
"| 17 | 7 | 7.1 | 12.7012876 | 2.09 | 12.7012876 | 30 | 12.72 | 3.28 |\n",
"| 18 | 1 | 1.1 | 0.6706613 | 0.09 | 0.6706613 | 4 | 0.66 | 0.34 |\n",
"| 19 | 0 | 0.1 | 0.6706613 | -1.40 | 0.6706613 | 4 | 0.66 | 1.31 |\n",
"| 20 | 7 | 7.1 | 3.6595929 | 2.21 | 3.6595929 | 19 | 3.63 | 2.65 |\n",
"| 21 | 0 | 0.1 | 1.4804839 | -0.61 | 1.4804839 | 11 | 1.46 | 2.92 |\n",
"| 22 | 0 | 0.1 | 0.4265684 | -1.85 | 0.4265684 | 0 | 0.42 | 0.83 |\n",
"| 23 | 0 | 0.1 | 0.4265684 | -1.85 | 0.4265684 | 0 | 0.42 | 0.83 |\n",
"| 24 | 0 | 0.1 | 0.4265684 | -1.85 | 0.4265684 | 0 | 0.42 | 0.83 |\n",
"| 25 | 0 | 0.1 | 0.5989280 | -1.51 | 0.5989280 | 3 | 0.59 | 1.17 |\n",
"| 26 | 0 | 0.1 | 1.0544301 | -0.95 | 1.0544301 | 8 | 1.04 | 2.07 |\n",
"| 27 | 0 | 0.1 | 1.0544301 | -0.95 | 1.0544301 | 8 | 1.04 | 2.07 |\n",
"| 28 | 0 | 0.1 | 2.3276530 | -0.16 | 2.3276530 | 15 | 2.30 | 4.60 |\n",
"| 29 | 3 | 3.1 | 4.5887040 | 1.18 | 4.5887040 | 21 | 4.56 | 0.80 |\n",
"| 30 | 4 | 4.1 | 6.4428192 | 1.48 | 6.4428192 | 24 | 6.42 | 1.25 |\n",
"| 31 | 6 | 6.1 | 2.3276530 | 2.42 | 2.3276530 | 15 | 2.30 | 4.30 |\n",
"\n"
],
"text/plain": [
" row yi uyumdegeri mu3 calismadegeri agirliklar x mu4i unit_dev\n",
"1 1 0 0.1 0.5348672 -1.63 0.5348672 2 0.52 1.05 \n",
"2 2 0 0.1 1.3221327 -0.72 1.3221327 10 1.30 2.60 \n",
"3 3 3 3.1 2.6064348 1.11 2.6064348 16 2.58 0.26 \n",
"4 4 2 2.1 4.0979005 0.90 4.0979005 20 4.07 1.49 \n",
"5 5 8 8.1 3.6595929 2.48 3.6595929 19 3.63 4.10 \n",
"6 6 1 1.1 3.2681662 0.49 3.2681662 18 3.24 2.32 \n",
"7 7 8 8.1 1.6578007 4.33 1.6578007 12 1.64 12.87 \n",
"8 8 5 5.1 2.6064348 1.88 2.6064348 16 2.58 1.98 \n",
"9 9 0 0.1 0.5989280 -1.51 0.5989280 3 0.59 1.17 \n",
"10 10 4 4.1 1.6578007 1.92 1.6578007 12 1.64 2.62 \n",
"11 11 19 19.1 15.9259380 2.96 15.9259380 32 15.98 0.74 \n",
"12 12 0 0.1 0.5348672 -1.63 0.5348672 2 0.52 1.05 \n",
"13 13 2 2.1 2.6064348 0.73 2.6064348 16 2.58 0.34 \n",
"14 14 0 0.1 0.9416492 -1.06 0.9416492 7 0.92 1.85 \n",
"15 15 3 3.1 1.3221327 1.55 1.3221327 10 1.30 1.81 \n",
"16 16 1 1.1 2.3276530 0.27 2.3276530 15 2.30 1.13 \n",
"17 17 7 7.1 12.7012876 2.09 12.7012876 30 12.72 3.28 \n",
"18 18 1 1.1 0.6706613 0.09 0.6706613 4 0.66 0.34 \n",
"19 19 0 0.1 0.6706613 -1.40 0.6706613 4 0.66 1.31 \n",
"20 20 7 7.1 3.6595929 2.21 3.6595929 19 3.63 2.65 \n",
"21 21 0 0.1 1.4804839 -0.61 1.4804839 11 1.46 2.92 \n",
"22 22 0 0.1 0.4265684 -1.85 0.4265684 0 0.42 0.83 \n",
"23 23 0 0.1 0.4265684 -1.85 0.4265684 0 0.42 0.83 \n",
"24 24 0 0.1 0.4265684 -1.85 0.4265684 0 0.42 0.83 \n",
"25 25 0 0.1 0.5989280 -1.51 0.5989280 3 0.59 1.17 \n",
"26 26 0 0.1 1.0544301 -0.95 1.0544301 8 1.04 2.07 \n",
"27 27 0 0.1 1.0544301 -0.95 1.0544301 8 1.04 2.07 \n",
"28 28 0 0.1 2.3276530 -0.16 2.3276530 15 2.30 4.60 \n",
"29 29 3 3.1 4.5887040 1.18 4.5887040 21 4.56 0.80 \n",
"30 30 4 4.1 6.4428192 1.48 6.4428192 24 6.42 1.25 \n",
"31 31 6 6.1 2.3276530 2.42 2.3276530 15 2.30 4.30 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"all_z4i <- calculate_z1i(log(all_mu3i), yi, all_mu3i)\n",
"\n",
"all_mu4i <- calculate_mu1i(coef_iter4, x)\n",
"\n",
"all_unit_dev4 <- calculate_unit_dev(yi, etai, all_mu4i)\n",
"\n",
"tablo_iter4 = data.frame(row = seq.int(yi), yi = yi, uyumdegeri = round(etai, 3), mu3 = all_mu3i, calismadegeri = round(all_z4i, 2), agirliklar = all_mu3i, x = x, mu4i = round(all_mu4i, 2), unit_dev = round(all_unit_dev4, 2))\n",
"\n",
"tablo_iter4"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "R",
"language": "R",
"name": "ir"
},
"language_info": {
"codemirror_mode": "r",
"file_extension": ".r",
"mimetype": "text/x-r-source",
"name": "R",
"pygments_lexer": "r",
"version": "4.2.0"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment