# Load Library dan Data
library(readxl)
library(mirt)
## Warning: package 'mirt' was built under R version 4.4.3
## Loading required package: stats4
## Loading required package: lattice
datasampel <- read_excel("C:/Users/ASUS/Documents/UNY/MySta/SEM 4/Statistika Pendidikan/UTS_Statpen/Dataset_Lana.xlsx", sheet = 2)
datasampel
## # A tibble: 150 × 37
## b1 b2 b3 b4 b5 b6 b7 b8 b9 b10 b11 b12 b13
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 0 1 0 0 0 0 0 0 0 1 0 1 0
## 2 1 0 1 0 1 1 0 1 0 1 1 0 1
## 3 1 1 1 1 1 1 1 1 1 1 1 1 1
## 4 0 0 1 1 1 0 1 0 0 0 0 0 0
## 5 0 1 0 1 0 0 0 0 1 0 1 0 0
## 6 0 0 0 0 0 0 0 0 0 1 0 0 0
## 7 0 0 0 0 0 0 0 0 0 1 0 0 0
## 8 1 1 1 1 1 1 1 1 1 1 1 1 1
## 9 1 1 1 1 1 1 1 1 1 1 1 1 1
## 10 1 1 1 1 1 1 1 1 1 1 1 1 1
## # ℹ 140 more rows
## # ℹ 24 more variables: b14 <dbl>, b15 <dbl>, b16 <dbl>, b17 <dbl>, b18 <dbl>,
## # b19 <dbl>, b20 <dbl>, b21 <dbl>, b22 <dbl>, b23 <dbl>, b24 <dbl>,
## # b25 <dbl>, b26 <dbl>, b27 <dbl>, b28 <dbl>, b29 <dbl>, b30 <dbl>,
## # b31 <dbl>, b32 <dbl>, b33 <dbl>, b34 <dbl>, b35 <dbl>, b36 <dbl>, b37 <dbl>
# Estimasi Model Rasch, 1PL, 2PL, 3PL, 4PL
# Model Rasch
mod_rasch <- mirt(datasampel, 1, itemtype = "Rasch", SE = TRUE)
## Iteration: 1, Log-Lik: -2858.002, Max-Change: 0.88713Iteration: 2, Log-Lik: -2827.568, Max-Change: 0.55961Iteration: 3, Log-Lik: -2823.147, Max-Change: 0.23844Iteration: 4, Log-Lik: -2822.559, Max-Change: 0.08298Iteration: 5, Log-Lik: -2822.399, Max-Change: 0.02880Iteration: 6, Log-Lik: -2822.312, Max-Change: 0.00994Iteration: 7, Log-Lik: -2822.134, Max-Change: 0.00883Iteration: 8, Log-Lik: -2822.084, Max-Change: 0.00729Iteration: 9, Log-Lik: -2822.040, Max-Change: 0.00688Iteration: 10, Log-Lik: -2821.838, Max-Change: 0.00458Iteration: 11, Log-Lik: -2821.821, Max-Change: 0.00430Iteration: 12, Log-Lik: -2821.805, Max-Change: 0.00405Iteration: 13, Log-Lik: -2821.736, Max-Change: 0.00267Iteration: 14, Log-Lik: -2821.730, Max-Change: 0.00255Iteration: 15, Log-Lik: -2821.724, Max-Change: 0.00240Iteration: 16, Log-Lik: -2821.700, Max-Change: 0.00210Iteration: 17, Log-Lik: -2821.698, Max-Change: 0.00152Iteration: 18, Log-Lik: -2821.696, Max-Change: 0.00142Iteration: 19, Log-Lik: -2821.687, Max-Change: 0.00149Iteration: 20, Log-Lik: -2821.687, Max-Change: 0.00090Iteration: 21, Log-Lik: -2821.686, Max-Change: 0.00085Iteration: 22, Log-Lik: -2821.683, Max-Change: 0.00095Iteration: 23, Log-Lik: -2821.682, Max-Change: 0.00054Iteration: 24, Log-Lik: -2821.682, Max-Change: 0.00050Iteration: 25, Log-Lik: -2821.681, Max-Change: 0.00060Iteration: 26, Log-Lik: -2821.681, Max-Change: 0.00032Iteration: 27, Log-Lik: -2821.681, Max-Change: 0.00030Iteration: 28, Log-Lik: -2821.681, Max-Change: 0.00037Iteration: 29, Log-Lik: -2821.680, Max-Change: 0.00019Iteration: 30, Log-Lik: -2821.680, Max-Change: 0.00017Iteration: 31, Log-Lik: -2821.680, Max-Change: 0.00021Iteration: 32, Log-Lik: -2821.680, Max-Change: 0.00011Iteration: 33, Log-Lik: -2821.680, Max-Change: 0.00010Iteration: 34, Log-Lik: -2821.680, Max-Change: 0.00013Iteration: 35, Log-Lik: -2821.680, Max-Change: 0.00007
##
## Calculating information matrix...
# Model 1PL (manual constraint)
model_1pl <- '
F1 = 1-37
CONSTRAIN = (1-37, a1)
'
mod_1pl <- mirt(datasampel, model_1pl, itemtype = "2PL", SE = TRUE)
## Iteration: 1, Log-Lik: -2883.820, Max-Change: 0.35974Iteration: 2, Log-Lik: -2850.661, Max-Change: 0.12551Iteration: 3, Log-Lik: -2837.880, Max-Change: 0.09703Iteration: 4, Log-Lik: -2831.309, Max-Change: 0.07668Iteration: 5, Log-Lik: -2827.655, Max-Change: 0.06125Iteration: 6, Log-Lik: -2825.514, Max-Change: 0.04932Iteration: 7, Log-Lik: -2822.137, Max-Change: 0.01299Iteration: 8, Log-Lik: -2822.044, Max-Change: 0.01085Iteration: 9, Log-Lik: -2821.977, Max-Change: 0.00904Iteration: 10, Log-Lik: -2821.800, Max-Change: 0.00389Iteration: 11, Log-Lik: -2821.788, Max-Change: 0.00311Iteration: 12, Log-Lik: -2821.778, Max-Change: 0.00255Iteration: 13, Log-Lik: -2821.732, Max-Change: 0.00223Iteration: 14, Log-Lik: -2821.727, Max-Change: 0.00200Iteration: 15, Log-Lik: -2821.723, Max-Change: 0.00199Iteration: 16, Log-Lik: -2821.706, Max-Change: 0.00692Iteration: 17, Log-Lik: -2821.704, Max-Change: 0.00182Iteration: 18, Log-Lik: -2821.703, Max-Change: 0.00121Iteration: 19, Log-Lik: -2821.696, Max-Change: 0.00101Iteration: 20, Log-Lik: -2821.695, Max-Change: 0.00081Iteration: 21, Log-Lik: -2821.695, Max-Change: 0.00074Iteration: 22, Log-Lik: -2821.692, Max-Change: 0.00052Iteration: 23, Log-Lik: -2821.691, Max-Change: 0.00049Iteration: 24, Log-Lik: -2821.691, Max-Change: 0.00029Iteration: 25, Log-Lik: -2821.691, Max-Change: 0.00027Iteration: 26, Log-Lik: -2821.691, Max-Change: 0.00024Iteration: 27, Log-Lik: -2821.691, Max-Change: 0.00021Iteration: 28, Log-Lik: -2821.691, Max-Change: 0.00029Iteration: 29, Log-Lik: -2821.691, Max-Change: 0.00012Iteration: 30, Log-Lik: -2821.691, Max-Change: 0.00071Iteration: 31, Log-Lik: -2821.691, Max-Change: 0.00026Iteration: 32, Log-Lik: -2821.691, Max-Change: 0.00010Iteration: 33, Log-Lik: -2821.691, Max-Change: 0.00055Iteration: 34, Log-Lik: -2821.690, Max-Change: 0.00024Iteration: 35, Log-Lik: -2821.690, Max-Change: 0.00009
##
## Calculating information matrix...
# Model 2PL
mod_2pl <- mirt(datasampel, 1, itemtype = "2PL", SE = TRUE)
## Iteration: 1, Log-Lik: -2883.820, Max-Change: 0.82253Iteration: 2, Log-Lik: -2779.845, Max-Change: 0.43948Iteration: 3, Log-Lik: -2763.651, Max-Change: 0.28936Iteration: 4, Log-Lik: -2755.711, Max-Change: 0.21147Iteration: 5, Log-Lik: -2751.294, Max-Change: 0.17181Iteration: 6, Log-Lik: -2748.587, Max-Change: 0.14222Iteration: 7, Log-Lik: -2744.085, Max-Change: 0.03571Iteration: 8, Log-Lik: -2743.909, Max-Change: 0.02695Iteration: 9, Log-Lik: -2743.840, Max-Change: 0.02044Iteration: 10, Log-Lik: -2743.717, Max-Change: 0.01599Iteration: 11, Log-Lik: -2743.674, Max-Change: 0.00984Iteration: 12, Log-Lik: -2743.643, Max-Change: 0.00827Iteration: 13, Log-Lik: -2743.545, Max-Change: 0.00529Iteration: 14, Log-Lik: -2743.528, Max-Change: 0.00526Iteration: 15, Log-Lik: -2743.512, Max-Change: 0.00522Iteration: 16, Log-Lik: -2743.432, Max-Change: 0.00593Iteration: 17, Log-Lik: -2743.424, Max-Change: 0.00521Iteration: 18, Log-Lik: -2743.416, Max-Change: 0.00478Iteration: 19, Log-Lik: -2743.378, Max-Change: 0.00349Iteration: 20, Log-Lik: -2743.373, Max-Change: 0.00403Iteration: 21, Log-Lik: -2743.369, Max-Change: 0.00319Iteration: 22, Log-Lik: -2743.357, Max-Change: 0.00423Iteration: 23, Log-Lik: -2743.354, Max-Change: 0.00278Iteration: 24, Log-Lik: -2743.351, Max-Change: 0.00278Iteration: 25, Log-Lik: -2743.336, Max-Change: 0.00214Iteration: 26, Log-Lik: -2743.335, Max-Change: 0.00214Iteration: 27, Log-Lik: -2743.333, Max-Change: 0.00205Iteration: 28, Log-Lik: -2743.326, Max-Change: 0.00183Iteration: 29, Log-Lik: -2743.325, Max-Change: 0.00161Iteration: 30, Log-Lik: -2743.324, Max-Change: 0.00156Iteration: 31, Log-Lik: -2743.321, Max-Change: 0.00125Iteration: 32, Log-Lik: -2743.320, Max-Change: 0.00120Iteration: 33, Log-Lik: -2743.320, Max-Change: 0.00116Iteration: 34, Log-Lik: -2743.318, Max-Change: 0.00090Iteration: 35, Log-Lik: -2743.318, Max-Change: 0.00087Iteration: 36, Log-Lik: -2743.318, Max-Change: 0.00084Iteration: 37, Log-Lik: -2743.317, Max-Change: 0.00065Iteration: 38, Log-Lik: -2743.317, Max-Change: 0.00063Iteration: 39, Log-Lik: -2743.317, Max-Change: 0.00061Iteration: 40, Log-Lik: -2743.316, Max-Change: 0.00032Iteration: 41, Log-Lik: -2743.316, Max-Change: 0.00037Iteration: 42, Log-Lik: -2743.316, Max-Change: 0.00036Iteration: 43, Log-Lik: -2743.316, Max-Change: 0.00029Iteration: 44, Log-Lik: -2743.316, Max-Change: 0.00027Iteration: 45, Log-Lik: -2743.316, Max-Change: 0.00027Iteration: 46, Log-Lik: -2743.316, Max-Change: 0.00026Iteration: 47, Log-Lik: -2743.316, Max-Change: 0.00021Iteration: 48, Log-Lik: -2743.316, Max-Change: 0.00020Iteration: 49, Log-Lik: -2743.316, Max-Change: 0.00018Iteration: 50, Log-Lik: -2743.316, Max-Change: 0.00017Iteration: 51, Log-Lik: -2743.316, Max-Change: 0.00016Iteration: 52, Log-Lik: -2743.316, Max-Change: 0.00017Iteration: 53, Log-Lik: -2743.316, Max-Change: 0.00014Iteration: 54, Log-Lik: -2743.316, Max-Change: 0.00014Iteration: 55, Log-Lik: -2743.316, Max-Change: 0.00012Iteration: 56, Log-Lik: -2743.316, Max-Change: 0.00012Iteration: 57, Log-Lik: -2743.316, Max-Change: 0.00011Iteration: 58, Log-Lik: -2743.315, Max-Change: 0.00012Iteration: 59, Log-Lik: -2743.315, Max-Change: 0.00010
##
## Calculating information matrix...
# Model 3PL
mod_3pl <- mirt(datasampel, 1, itemtype = "3PL", SE = TRUE)
## Iteration: 1, Log-Lik: -2920.204, Max-Change: 1.87712Iteration: 2, Log-Lik: -2756.112, Max-Change: 1.36100Iteration: 3, Log-Lik: -2726.417, Max-Change: 1.11624Iteration: 4, Log-Lik: -2711.942, Max-Change: 0.85281Iteration: 5, Log-Lik: -2703.946, Max-Change: 0.55873Iteration: 6, Log-Lik: -2699.468, Max-Change: 0.45229Iteration: 7, Log-Lik: -2693.920, Max-Change: 0.41679Iteration: 8, Log-Lik: -2692.960, Max-Change: 0.47161Iteration: 9, Log-Lik: -2692.344, Max-Change: 0.64634Iteration: 10, Log-Lik: -2691.357, Max-Change: 0.18062Iteration: 11, Log-Lik: -2691.062, Max-Change: 0.93502Iteration: 12, Log-Lik: -2690.876, Max-Change: 1.24478Iteration: 13, Log-Lik: -2692.339, Max-Change: 0.64506Iteration: 14, Log-Lik: -2690.361, Max-Change: 0.12760Iteration: 15, Log-Lik: -2690.263, Max-Change: 0.28582Iteration: 16, Log-Lik: -2690.124, Max-Change: 0.13390Iteration: 17, Log-Lik: -2690.052, Max-Change: 0.09815Iteration: 18, Log-Lik: -2689.989, Max-Change: 0.10586Iteration: 19, Log-Lik: -2689.853, Max-Change: 0.05304Iteration: 20, Log-Lik: -2689.808, Max-Change: 0.04746Iteration: 21, Log-Lik: -2689.768, Max-Change: 0.05523Iteration: 22, Log-Lik: -2689.628, Max-Change: 0.05650Iteration: 23, Log-Lik: -2689.595, Max-Change: 0.19943Iteration: 24, Log-Lik: -2689.569, Max-Change: 0.03170Iteration: 25, Log-Lik: -2689.562, Max-Change: 0.03219Iteration: 26, Log-Lik: -2689.540, Max-Change: 0.02288Iteration: 27, Log-Lik: -2689.521, Max-Change: 0.02853Iteration: 28, Log-Lik: -2689.459, Max-Change: 0.02237Iteration: 29, Log-Lik: -2689.447, Max-Change: 0.02222Iteration: 30, Log-Lik: -2689.432, Max-Change: 0.15499Iteration: 31, Log-Lik: -2689.416, Max-Change: 0.01569Iteration: 32, Log-Lik: -2689.404, Max-Change: 0.01805Iteration: 33, Log-Lik: -2689.392, Max-Change: 0.01815Iteration: 34, Log-Lik: -2689.334, Max-Change: 0.02360Iteration: 35, Log-Lik: -2689.327, Max-Change: 0.01222Iteration: 36, Log-Lik: -2689.321, Max-Change: 0.02578Iteration: 37, Log-Lik: -2689.313, Max-Change: 0.01340Iteration: 38, Log-Lik: -2689.305, Max-Change: 0.01321Iteration: 39, Log-Lik: -2689.300, Max-Change: 0.01378Iteration: 40, Log-Lik: -2689.288, Max-Change: 0.01661Iteration: 41, Log-Lik: -2689.283, Max-Change: 0.00819Iteration: 42, Log-Lik: -2689.280, Max-Change: 0.01762Iteration: 43, Log-Lik: -2689.274, Max-Change: 0.01283Iteration: 44, Log-Lik: -2689.270, Max-Change: 0.01819Iteration: 45, Log-Lik: -2689.267, Max-Change: 0.01678Iteration: 46, Log-Lik: -2689.255, Max-Change: 0.01501Iteration: 47, Log-Lik: -2689.253, Max-Change: 0.00585Iteration: 48, Log-Lik: -2689.251, Max-Change: 0.01187Iteration: 49, Log-Lik: -2689.248, Max-Change: 0.01184Iteration: 50, Log-Lik: -2689.246, Max-Change: 0.00453Iteration: 51, Log-Lik: -2689.244, Max-Change: 0.00386Iteration: 52, Log-Lik: -2689.243, Max-Change: 0.01263Iteration: 53, Log-Lik: -2689.241, Max-Change: 0.00901Iteration: 54, Log-Lik: -2689.239, Max-Change: 0.00626Iteration: 55, Log-Lik: -2689.238, Max-Change: 0.00745Iteration: 56, Log-Lik: -2689.236, Max-Change: 0.00358Iteration: 57, Log-Lik: -2689.236, Max-Change: 0.00684Iteration: 58, Log-Lik: -2689.234, Max-Change: 0.00240Iteration: 59, Log-Lik: -2689.233, Max-Change: 0.00731Iteration: 60, Log-Lik: -2689.233, Max-Change: 0.00180Iteration: 61, Log-Lik: -2689.232, Max-Change: 0.00282Iteration: 62, Log-Lik: -2689.232, Max-Change: 0.00281Iteration: 63, Log-Lik: -2689.231, Max-Change: 0.00686Iteration: 64, Log-Lik: -2689.230, Max-Change: 0.00179Iteration: 65, Log-Lik: -2689.230, Max-Change: 0.00322Iteration: 66, Log-Lik: -2689.230, Max-Change: 0.00726Iteration: 67, Log-Lik: -2689.229, Max-Change: 0.00156Iteration: 68, Log-Lik: -2689.228, Max-Change: 0.00134Iteration: 69, Log-Lik: -2689.228, Max-Change: 0.00339Iteration: 70, Log-Lik: -2689.228, Max-Change: 0.00741Iteration: 71, Log-Lik: -2689.227, Max-Change: 0.00205Iteration: 72, Log-Lik: -2689.227, Max-Change: 0.00580Iteration: 73, Log-Lik: -2689.226, Max-Change: 0.00131Iteration: 74, Log-Lik: -2689.226, Max-Change: 0.00231Iteration: 75, Log-Lik: -2689.226, Max-Change: 0.00565Iteration: 76, Log-Lik: -2689.225, Max-Change: 0.00113Iteration: 77, Log-Lik: -2689.225, Max-Change: 0.00203Iteration: 78, Log-Lik: -2689.225, Max-Change: 0.00196Iteration: 79, Log-Lik: -2689.225, Max-Change: 0.00288Iteration: 80, Log-Lik: -2689.224, Max-Change: 0.00524Iteration: 81, Log-Lik: -2689.224, Max-Change: 0.00145Iteration: 82, Log-Lik: -2689.224, Max-Change: 0.00169Iteration: 83, Log-Lik: -2689.224, Max-Change: 0.00484Iteration: 84, Log-Lik: -2689.224, Max-Change: 0.00118Iteration: 85, Log-Lik: -2689.224, Max-Change: 0.00134Iteration: 86, Log-Lik: -2689.224, Max-Change: 0.00178Iteration: 87, Log-Lik: -2689.223, Max-Change: 0.00050Iteration: 88, Log-Lik: -2689.223, Max-Change: 0.00037Iteration: 89, Log-Lik: -2689.223, Max-Change: 0.00148Iteration: 90, Log-Lik: -2689.223, Max-Change: 0.00252Iteration: 91, Log-Lik: -2689.223, Max-Change: 0.00158Iteration: 92, Log-Lik: -2689.223, Max-Change: 0.00089Iteration: 93, Log-Lik: -2689.223, Max-Change: 0.00146Iteration: 94, Log-Lik: -2689.223, Max-Change: 0.00066Iteration: 95, Log-Lik: -2689.223, Max-Change: 0.00464Iteration: 96, Log-Lik: -2689.223, Max-Change: 0.00119Iteration: 97, Log-Lik: -2689.223, Max-Change: 0.00116Iteration: 98, Log-Lik: -2689.223, Max-Change: 0.00108Iteration: 99, Log-Lik: -2689.223, Max-Change: 0.00069Iteration: 100, Log-Lik: -2689.223, Max-Change: 0.00048Iteration: 101, Log-Lik: -2689.223, Max-Change: 0.00391Iteration: 102, Log-Lik: -2689.222, Max-Change: 0.00121Iteration: 103, Log-Lik: -2689.222, Max-Change: 0.00104Iteration: 104, Log-Lik: -2689.222, Max-Change: 0.00096Iteration: 105, Log-Lik: -2689.222, Max-Change: 0.00067Iteration: 106, Log-Lik: -2689.222, Max-Change: 0.00047Iteration: 107, Log-Lik: -2689.222, Max-Change: 0.00119Iteration: 108, Log-Lik: -2689.222, Max-Change: 0.00143Iteration: 109, Log-Lik: -2689.222, Max-Change: 0.00204Iteration: 110, Log-Lik: -2689.222, Max-Change: 0.00030Iteration: 111, Log-Lik: -2689.222, Max-Change: 0.00104Iteration: 112, Log-Lik: -2689.222, Max-Change: 0.00023Iteration: 113, Log-Lik: -2689.222, Max-Change: 0.00107Iteration: 114, Log-Lik: -2689.222, Max-Change: 0.00021Iteration: 115, Log-Lik: -2689.222, Max-Change: 0.00104Iteration: 116, Log-Lik: -2689.222, Max-Change: 0.00049Iteration: 117, Log-Lik: -2689.222, Max-Change: 0.00101Iteration: 118, Log-Lik: -2689.222, Max-Change: 0.00036Iteration: 119, Log-Lik: -2689.222, Max-Change: 0.00100Iteration: 120, Log-Lik: -2689.222, Max-Change: 0.00026Iteration: 121, Log-Lik: -2689.222, Max-Change: 0.00098Iteration: 122, Log-Lik: -2689.222, Max-Change: 0.00081Iteration: 123, Log-Lik: -2689.222, Max-Change: 0.00097Iteration: 124, Log-Lik: -2689.222, Max-Change: 0.00126Iteration: 125, Log-Lik: -2689.222, Max-Change: 0.00147Iteration: 126, Log-Lik: -2689.222, Max-Change: 0.00097Iteration: 127, Log-Lik: -2689.222, Max-Change: 0.00080Iteration: 128, Log-Lik: -2689.222, Max-Change: 0.00086Iteration: 129, Log-Lik: -2689.222, Max-Change: 0.00051Iteration: 130, Log-Lik: -2689.222, Max-Change: 0.00033Iteration: 131, Log-Lik: -2689.222, Max-Change: 0.00088Iteration: 132, Log-Lik: -2689.222, Max-Change: 0.00100Iteration: 133, Log-Lik: -2689.222, Max-Change: 0.00112Iteration: 134, Log-Lik: -2689.221, Max-Change: 0.00031Iteration: 135, Log-Lik: -2689.221, Max-Change: 0.00085Iteration: 136, Log-Lik: -2689.221, Max-Change: 0.00036Iteration: 137, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 138, Log-Lik: -2689.221, Max-Change: 0.00030Iteration: 139, Log-Lik: -2689.221, Max-Change: 0.00123Iteration: 140, Log-Lik: -2689.221, Max-Change: 0.00051Iteration: 141, Log-Lik: -2689.221, Max-Change: 0.00095Iteration: 142, Log-Lik: -2689.221, Max-Change: 0.00121Iteration: 143, Log-Lik: -2689.221, Max-Change: 0.00139Iteration: 144, Log-Lik: -2689.221, Max-Change: 0.00095Iteration: 145, Log-Lik: -2689.221, Max-Change: 0.00069Iteration: 146, Log-Lik: -2689.221, Max-Change: 0.00086Iteration: 147, Log-Lik: -2689.221, Max-Change: 0.00043Iteration: 148, Log-Lik: -2689.221, Max-Change: 0.00138Iteration: 149, Log-Lik: -2689.221, Max-Change: 0.00129Iteration: 150, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 151, Log-Lik: -2689.221, Max-Change: 0.00065Iteration: 152, Log-Lik: -2689.221, Max-Change: 0.00086Iteration: 153, Log-Lik: -2689.221, Max-Change: 0.00040Iteration: 154, Log-Lik: -2689.221, Max-Change: 0.00131Iteration: 155, Log-Lik: -2689.221, Max-Change: 0.00122Iteration: 156, Log-Lik: -2689.221, Max-Change: 0.00086Iteration: 157, Log-Lik: -2689.221, Max-Change: 0.00062Iteration: 158, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 159, Log-Lik: -2689.221, Max-Change: 0.00039Iteration: 160, Log-Lik: -2689.221, Max-Change: 0.00125Iteration: 161, Log-Lik: -2689.221, Max-Change: 0.00118Iteration: 162, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 163, Log-Lik: -2689.221, Max-Change: 0.00060Iteration: 164, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 165, Log-Lik: -2689.221, Max-Change: 0.00038Iteration: 166, Log-Lik: -2689.221, Max-Change: 0.00122Iteration: 167, Log-Lik: -2689.221, Max-Change: 0.00114Iteration: 168, Log-Lik: -2689.221, Max-Change: 0.00087Iteration: 169, Log-Lik: -2689.221, Max-Change: 0.00058Iteration: 170, Log-Lik: -2689.221, Max-Change: 0.00088Iteration: 171, Log-Lik: -2689.221, Max-Change: 0.00036Iteration: 172, Log-Lik: -2689.221, Max-Change: 0.00118Iteration: 173, Log-Lik: -2689.221, Max-Change: 0.00111Iteration: 174, Log-Lik: -2689.221, Max-Change: 0.00088Iteration: 175, Log-Lik: -2689.221, Max-Change: 0.00057Iteration: 176, Log-Lik: -2689.221, Max-Change: 0.00440Iteration: 177, Log-Lik: -2689.221, Max-Change: 0.00471Iteration: 178, Log-Lik: -2689.221, Max-Change: 0.00209Iteration: 179, Log-Lik: -2689.221, Max-Change: 0.00080Iteration: 180, Log-Lik: -2689.221, Max-Change: 0.00090Iteration: 181, Log-Lik: -2689.221, Max-Change: 0.00049Iteration: 182, Log-Lik: -2689.221, Max-Change: 0.00448Iteration: 183, Log-Lik: -2689.221, Max-Change: 0.00155Iteration: 184, Log-Lik: -2689.221, Max-Change: 0.00108Iteration: 185, Log-Lik: -2689.221, Max-Change: 0.00090Iteration: 186, Log-Lik: -2689.221, Max-Change: 0.00067Iteration: 187, Log-Lik: -2689.221, Max-Change: 0.00038Iteration: 188, Log-Lik: -2689.221, Max-Change: 0.00444Iteration: 189, Log-Lik: -2689.220, Max-Change: 0.00140Iteration: 190, Log-Lik: -2689.220, Max-Change: 0.00112Iteration: 191, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 192, Log-Lik: -2689.220, Max-Change: 0.00070Iteration: 193, Log-Lik: -2689.220, Max-Change: 0.00040Iteration: 194, Log-Lik: -2689.220, Max-Change: 0.00442Iteration: 195, Log-Lik: -2689.220, Max-Change: 0.00112Iteration: 196, Log-Lik: -2689.220, Max-Change: 0.00112Iteration: 197, Log-Lik: -2689.220, Max-Change: 0.00090Iteration: 198, Log-Lik: -2689.220, Max-Change: 0.00069Iteration: 199, Log-Lik: -2689.220, Max-Change: 0.00043Iteration: 200, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 201, Log-Lik: -2689.220, Max-Change: 0.00027Iteration: 202, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 203, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 204, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 205, Log-Lik: -2689.220, Max-Change: 0.00060Iteration: 206, Log-Lik: -2689.220, Max-Change: 0.00090Iteration: 207, Log-Lik: -2689.220, Max-Change: 0.00038Iteration: 208, Log-Lik: -2689.220, Max-Change: 0.00025Iteration: 209, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 210, Log-Lik: -2689.220, Max-Change: 0.00090Iteration: 211, Log-Lik: -2689.220, Max-Change: 0.00072Iteration: 212, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 213, Log-Lik: -2689.220, Max-Change: 0.00046Iteration: 214, Log-Lik: -2689.220, Max-Change: 0.00029Iteration: 215, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 216, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 217, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 218, Log-Lik: -2689.220, Max-Change: 0.00056Iteration: 219, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 220, Log-Lik: -2689.220, Max-Change: 0.00041Iteration: 221, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 222, Log-Lik: -2689.220, Max-Change: 0.00026Iteration: 223, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 224, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 225, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 226, Log-Lik: -2689.220, Max-Change: 0.00056Iteration: 227, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 228, Log-Lik: -2689.220, Max-Change: 0.00035Iteration: 229, Log-Lik: -2689.220, Max-Change: 0.00022Iteration: 230, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 231, Log-Lik: -2689.220, Max-Change: 0.00089Iteration: 232, Log-Lik: -2689.220, Max-Change: 0.00065Iteration: 233, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 234, Log-Lik: -2689.220, Max-Change: 0.00040Iteration: 235, Log-Lik: -2689.220, Max-Change: 0.00026Iteration: 236, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 237, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 238, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 239, Log-Lik: -2689.220, Max-Change: 0.00051Iteration: 240, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 241, Log-Lik: -2689.220, Max-Change: 0.00038Iteration: 242, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 243, Log-Lik: -2689.220, Max-Change: 0.00024Iteration: 244, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 245, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 246, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 247, Log-Lik: -2689.220, Max-Change: 0.00048Iteration: 248, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 249, Log-Lik: -2689.220, Max-Change: 0.00030Iteration: 250, Log-Lik: -2689.220, Max-Change: 0.00020Iteration: 251, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 252, Log-Lik: -2689.220, Max-Change: 0.00088Iteration: 253, Log-Lik: -2689.220, Max-Change: 0.00054Iteration: 254, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 255, Log-Lik: -2689.220, Max-Change: 0.00034Iteration: 256, Log-Lik: -2689.220, Max-Change: 0.00022Iteration: 257, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 258, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 259, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 260, Log-Lik: -2689.220, Max-Change: 0.00043Iteration: 261, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 262, Log-Lik: -2689.220, Max-Change: 0.00035Iteration: 263, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 264, Log-Lik: -2689.220, Max-Change: 0.00022Iteration: 265, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 266, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 267, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 268, Log-Lik: -2689.220, Max-Change: 0.00041Iteration: 269, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 270, Log-Lik: -2689.220, Max-Change: 0.00025Iteration: 271, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 272, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 273, Log-Lik: -2689.220, Max-Change: 0.00087Iteration: 274, Log-Lik: -2689.220, Max-Change: 0.00046Iteration: 275, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 276, Log-Lik: -2689.220, Max-Change: 0.00028Iteration: 277, Log-Lik: -2689.220, Max-Change: 0.00019Iteration: 278, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 279, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 280, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 281, Log-Lik: -2689.220, Max-Change: 0.00037Iteration: 282, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 283, Log-Lik: -2689.220, Max-Change: 0.00031Iteration: 284, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 285, Log-Lik: -2689.220, Max-Change: 0.00020Iteration: 286, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 287, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 288, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 289, Log-Lik: -2689.220, Max-Change: 0.00035Iteration: 290, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 291, Log-Lik: -2689.220, Max-Change: 0.00022Iteration: 292, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 293, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 294, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 295, Log-Lik: -2689.220, Max-Change: 0.00038Iteration: 296, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 297, Log-Lik: -2689.220, Max-Change: 0.00025Iteration: 298, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 299, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 300, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 301, Log-Lik: -2689.220, Max-Change: 0.00042Iteration: 302, Log-Lik: -2689.220, Max-Change: 0.00086Iteration: 303, Log-Lik: -2689.220, Max-Change: 0.00027Iteration: 304, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 305, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 306, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 307, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 308, Log-Lik: -2689.220, Max-Change: 0.00034Iteration: 309, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 310, Log-Lik: -2689.220, Max-Change: 0.00030Iteration: 311, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 312, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 313, Log-Lik: -2689.220, Max-Change: 0.00016Iteration: 314, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 315, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 316, Log-Lik: -2689.220, Max-Change: 0.00033Iteration: 317, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 318, Log-Lik: -2689.220, Max-Change: 0.00020Iteration: 319, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 320, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 321, Log-Lik: -2689.220, Max-Change: 0.00082Iteration: 322, Log-Lik: -2689.220, Max-Change: 0.00036Iteration: 323, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 324, Log-Lik: -2689.220, Max-Change: 0.00022Iteration: 325, Log-Lik: -2689.220, Max-Change: 0.00016Iteration: 326, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 327, Log-Lik: -2689.220, Max-Change: 0.00085Iteration: 328, Log-Lik: -2689.220, Max-Change: 0.00040Iteration: 329, Log-Lik: -2689.220, Max-Change: 0.00082Iteration: 330, Log-Lik: -2689.220, Max-Change: 0.00025Iteration: 331, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 332, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 333, Log-Lik: -2689.220, Max-Change: 0.00016Iteration: 334, Log-Lik: -2689.220, Max-Change: 0.00082Iteration: 335, Log-Lik: -2689.220, Max-Change: 0.00033Iteration: 336, Log-Lik: -2689.220, Max-Change: 0.00084Iteration: 337, Log-Lik: -2689.220, Max-Change: 0.00027Iteration: 338, Log-Lik: -2689.220, Max-Change: 0.00082Iteration: 339, Log-Lik: -2689.220, Max-Change: 0.00018Iteration: 340, Log-Lik: -2689.220, Max-Change: 0.00017Iteration: 341, Log-Lik: -2689.220, Max-Change: 0.00083Iteration: 342, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 343, Log-Lik: -2689.219, Max-Change: 0.00030Iteration: 344, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 345, Log-Lik: -2689.219, Max-Change: 0.00019Iteration: 346, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 347, Log-Lik: -2689.219, Max-Change: 0.00082Iteration: 348, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 349, Log-Lik: -2689.219, Max-Change: 0.00033Iteration: 350, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 351, Log-Lik: -2689.219, Max-Change: 0.00021Iteration: 352, Log-Lik: -2689.219, Max-Change: 0.00017Iteration: 353, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 354, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 355, Log-Lik: -2689.219, Max-Change: 0.00036Iteration: 356, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 357, Log-Lik: -2689.219, Max-Change: 0.00024Iteration: 358, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 359, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 360, Log-Lik: -2689.219, Max-Change: 0.00017Iteration: 361, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 362, Log-Lik: -2689.219, Max-Change: 0.00030Iteration: 363, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 364, Log-Lik: -2689.219, Max-Change: 0.00026Iteration: 365, Log-Lik: -2689.219, Max-Change: 0.00082Iteration: 366, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 367, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 368, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 369, Log-Lik: -2689.219, Max-Change: 0.00083Iteration: 370, Log-Lik: -2689.219, Max-Change: 0.00029Iteration: 371, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 372, Log-Lik: -2689.219, Max-Change: 0.00019Iteration: 373, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 374, Log-Lik: -2689.219, Max-Change: 0.00082Iteration: 375, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 376, Log-Lik: -2689.219, Max-Change: 0.00031Iteration: 377, Log-Lik: -2689.219, Max-Change: 0.00082Iteration: 378, Log-Lik: -2689.219, Max-Change: 0.00019Iteration: 379, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 380, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 381, Log-Lik: -2689.219, Max-Change: 0.00082Iteration: 382, Log-Lik: -2689.219, Max-Change: 0.00035Iteration: 383, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 384, Log-Lik: -2689.219, Max-Change: 0.00021Iteration: 385, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 386, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 387, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 388, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 389, Log-Lik: -2689.219, Max-Change: 0.00028Iteration: 390, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 391, Log-Lik: -2689.219, Max-Change: 0.00026Iteration: 392, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 393, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 394, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 395, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 396, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 397, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 398, Log-Lik: -2689.219, Max-Change: 0.00024Iteration: 399, Log-Lik: -2689.219, Max-Change: 0.00081Iteration: 400, Log-Lik: -2689.219, Max-Change: 0.00029Iteration: 401, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 402, Log-Lik: -2689.219, Max-Change: 0.00023Iteration: 403, Log-Lik: -2689.219, Max-Change: 0.00019Iteration: 404, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 405, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 406, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 407, Log-Lik: -2689.219, Max-Change: 0.00027Iteration: 408, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 409, Log-Lik: -2689.219, Max-Change: 0.00035Iteration: 410, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 411, Log-Lik: -2689.219, Max-Change: 0.00028Iteration: 412, Log-Lik: -2689.219, Max-Change: 0.00022Iteration: 413, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 414, Log-Lik: -2689.219, Max-Change: 0.00018Iteration: 415, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 416, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 417, Log-Lik: -2689.219, Max-Change: 0.00080Iteration: 418, Log-Lik: -2689.219, Max-Change: 0.00024Iteration: 419, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 420, Log-Lik: -2689.219, Max-Change: 0.00019Iteration: 421, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 422, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 423, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 424, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 425, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 426, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 427, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 428, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 429, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 430, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 431, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 432, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 433, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 434, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 435, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 436, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 437, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 438, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 439, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 440, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 441, Log-Lik: -2689.219, Max-Change: 0.00079Iteration: 442, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 443, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 444, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 445, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 446, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 447, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 448, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 449, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 450, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 451, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 452, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 453, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 454, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 455, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 456, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 457, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 458, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 459, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 460, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 461, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 462, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 463, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 464, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 465, Log-Lik: -2689.219, Max-Change: 0.00078Iteration: 466, Log-Lik: -2689.219, Max-Change: 0.00026Iteration: 467, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 468, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 469, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 470, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 471, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 472, Log-Lik: -2689.219, Max-Change: 0.00026Iteration: 473, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 474, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 475, Log-Lik: -2689.219, Max-Change: 0.00017Iteration: 476, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 477, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 478, Log-Lik: -2689.219, Max-Change: 0.00026Iteration: 479, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 480, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 481, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 482, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 483, Log-Lik: -2689.219, Max-Change: 0.00075Iteration: 484, Log-Lik: -2689.219, Max-Change: 0.00025Iteration: 485, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 486, Log-Lik: -2689.219, Max-Change: 0.00021Iteration: 487, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 488, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 489, Log-Lik: -2689.219, Max-Change: 0.00077Iteration: 490, Log-Lik: -2689.219, Max-Change: 0.00028Iteration: 491, Log-Lik: -2689.219, Max-Change: 0.00075Iteration: 492, Log-Lik: -2689.219, Max-Change: 0.00020Iteration: 493, Log-Lik: -2689.219, Max-Change: 0.00016Iteration: 494, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 495, Log-Lik: -2689.219, Max-Change: 0.00015Iteration: 496, Log-Lik: -2689.219, Max-Change: 0.00075Iteration: 497, Log-Lik: -2689.219, Max-Change: 0.00023Iteration: 498, Log-Lik: -2689.219, Max-Change: 0.00076Iteration: 499, Log-Lik: -2689.219, Max-Change: 0.00028Iteration: 500, Log-Lik: -2689.219, Max-Change: 0.00075
## Warning: EM cycles terminated after 500 iterations.
##
##
## Calculating information matrix...
# Model 4PL
mod_4pl <- mirt(datasampel, 1, itemtype = "4PL", SE = TRUE)
## Iteration: 1, Log-Lik: -3083.659, Max-Change: 4.86216Iteration: 2, Log-Lik: -2762.096, Max-Change: 2.28576Iteration: 3, Log-Lik: -2714.496, Max-Change: 1.33019Iteration: 4, Log-Lik: -2700.297, Max-Change: 1.04665Iteration: 5, Log-Lik: -2693.395, Max-Change: 0.85922Iteration: 6, Log-Lik: -2689.027, Max-Change: 0.70586Iteration: 7, Log-Lik: -2684.520, Max-Change: 0.71040Iteration: 8, Log-Lik: -2682.051, Max-Change: 0.98499Iteration: 9, Log-Lik: -2681.188, Max-Change: 4.92946Iteration: 10, Log-Lik: -2679.694, Max-Change: 0.27113Iteration: 11, Log-Lik: -2679.277, Max-Change: 5.08791Iteration: 12, Log-Lik: -2678.771, Max-Change: 0.71537Iteration: 13, Log-Lik: -2678.725, Max-Change: 4.45057Iteration: 14, Log-Lik: -2678.413, Max-Change: 0.41467Iteration: 15, Log-Lik: -2678.266, Max-Change: 1.37942Iteration: 16, Log-Lik: -2678.106, Max-Change: 0.72648Iteration: 17, Log-Lik: -2678.008, Max-Change: 0.97010Iteration: 18, Log-Lik: -2677.917, Max-Change: 0.69137Iteration: 19, Log-Lik: -2677.816, Max-Change: 1.27667Iteration: 20, Log-Lik: -2677.703, Max-Change: 2.18071Iteration: 21, Log-Lik: -2677.518, Max-Change: 3.12860Iteration: 22, Log-Lik: -2676.951, Max-Change: 2.12239Iteration: 23, Log-Lik: -2676.733, Max-Change: 6.39730Iteration: 24, Log-Lik: -2676.499, Max-Change: 4.17432Iteration: 25, Log-Lik: -2676.366, Max-Change: 0.95612Iteration: 26, Log-Lik: -2676.307, Max-Change: 3.07682Iteration: 27, Log-Lik: -2676.265, Max-Change: 0.65355Iteration: 28, Log-Lik: -2676.258, Max-Change: 1.14222Iteration: 29, Log-Lik: -2676.238, Max-Change: 0.11908Iteration: 30, Log-Lik: -2676.231, Max-Change: 8.49891Iteration: 31, Log-Lik: -2676.176, Max-Change: 0.40728Iteration: 32, Log-Lik: -2676.162, Max-Change: 0.12172Iteration: 33, Log-Lik: -2676.159, Max-Change: 0.01024Iteration: 34, Log-Lik: -2676.159, Max-Change: 0.00968Iteration: 35, Log-Lik: -2676.158, Max-Change: 0.00503Iteration: 36, Log-Lik: -2676.158, Max-Change: 0.00128Iteration: 37, Log-Lik: -2676.158, Max-Change: 0.00060Iteration: 38, Log-Lik: -2676.158, Max-Change: 0.00336Iteration: 39, Log-Lik: -2676.157, Max-Change: 0.00506Iteration: 40, Log-Lik: -2676.157, Max-Change: 0.00254Iteration: 41, Log-Lik: -2676.157, Max-Change: 0.00126Iteration: 42, Log-Lik: -2676.157, Max-Change: 0.00039Iteration: 43, Log-Lik: -2676.157, Max-Change: 0.00032Iteration: 44, Log-Lik: -2676.157, Max-Change: 0.00086Iteration: 45, Log-Lik: -2676.157, Max-Change: 0.00086Iteration: 46, Log-Lik: -2676.157, Max-Change: 0.00038Iteration: 47, Log-Lik: -2676.157, Max-Change: 0.00085Iteration: 48, Log-Lik: -2676.157, Max-Change: 0.00100Iteration: 49, Log-Lik: -2676.157, Max-Change: 0.00046Iteration: 50, Log-Lik: -2676.157, Max-Change: 0.00085Iteration: 51, Log-Lik: -2676.157, Max-Change: 0.00024Iteration: 52, Log-Lik: -2676.157, Max-Change: 0.00102Iteration: 53, Log-Lik: -2676.157, Max-Change: 0.00050Iteration: 54, Log-Lik: -2676.157, Max-Change: 0.00084Iteration: 55, Log-Lik: -2676.157, Max-Change: 0.00054Iteration: 56, Log-Lik: -2676.157, Max-Change: 0.00105Iteration: 57, Log-Lik: -2676.157, Max-Change: 0.00042Iteration: 58, Log-Lik: -2676.157, Max-Change: 0.00034Iteration: 59, Log-Lik: -2676.157, Max-Change: 0.00084Iteration: 60, Log-Lik: -2676.157, Max-Change: 0.00026Iteration: 61, Log-Lik: -2676.157, Max-Change: 0.00109Iteration: 62, Log-Lik: -2676.157, Max-Change: 0.00050Iteration: 63, Log-Lik: -2676.157, Max-Change: 0.00083Iteration: 64, Log-Lik: -2676.157, Max-Change: 0.00067Iteration: 65, Log-Lik: -2676.157, Max-Change: 0.00119Iteration: 66, Log-Lik: -2676.157, Max-Change: 0.00052Iteration: 67, Log-Lik: -2676.157, Max-Change: 0.00041Iteration: 68, Log-Lik: -2676.157, Max-Change: 0.00082Iteration: 69, Log-Lik: -2676.157, Max-Change: 0.00032Iteration: 70, Log-Lik: -2676.157, Max-Change: 0.00025Iteration: 71, Log-Lik: -2676.157, Max-Change: 0.00082Iteration: 72, Log-Lik: -2676.157, Max-Change: 0.00100Iteration: 73, Log-Lik: -2676.157, Max-Change: 0.00041Iteration: 74, Log-Lik: -2676.157, Max-Change: 0.00081Iteration: 75, Log-Lik: -2676.157, Max-Change: 0.00032Iteration: 76, Log-Lik: -2676.157, Max-Change: 0.00026Iteration: 77, Log-Lik: -2676.157, Max-Change: 0.00081Iteration: 78, Log-Lik: -2676.157, Max-Change: 0.00098Iteration: 79, Log-Lik: -2676.157, Max-Change: 0.00044Iteration: 80, Log-Lik: -2676.157, Max-Change: 0.00080Iteration: 81, Log-Lik: -2676.157, Max-Change: 0.00034Iteration: 82, Log-Lik: -2676.157, Max-Change: 0.00027Iteration: 83, Log-Lik: -2676.157, Max-Change: 0.00080Iteration: 84, Log-Lik: -2676.157, Max-Change: 0.00106Iteration: 85, Log-Lik: -2676.157, Max-Change: 0.00046Iteration: 86, Log-Lik: -2676.157, Max-Change: 0.00079Iteration: 87, Log-Lik: -2676.157, Max-Change: 0.00036Iteration: 88, Log-Lik: -2676.157, Max-Change: 0.00029Iteration: 89, Log-Lik: -2676.157, Max-Change: 0.00078Iteration: 90, Log-Lik: -2676.157, Max-Change: 0.00110Iteration: 91, Log-Lik: -2676.157, Max-Change: 0.00049Iteration: 92, Log-Lik: -2676.157, Max-Change: 0.00080Iteration: 93, Log-Lik: -2676.157, Max-Change: 0.00038Iteration: 94, Log-Lik: -2676.157, Max-Change: 0.00030Iteration: 95, Log-Lik: -2676.157, Max-Change: 0.00077Iteration: 96, Log-Lik: -2676.157, Max-Change: 0.00118Iteration: 97, Log-Lik: -2676.157, Max-Change: 0.00051Iteration: 98, Log-Lik: -2676.157, Max-Change: 0.00082Iteration: 99, Log-Lik: -2676.157, Max-Change: 0.00039Iteration: 100, Log-Lik: -2676.157, Max-Change: 0.00032Iteration: 101, Log-Lik: -2676.157, Max-Change: 0.00076Iteration: 102, Log-Lik: -2676.157, Max-Change: 0.00120Iteration: 103, Log-Lik: -2676.157, Max-Change: 0.00053Iteration: 104, Log-Lik: -2676.157, Max-Change: 0.00086Iteration: 105, Log-Lik: -2676.157, Max-Change: 0.00040Iteration: 106, Log-Lik: -2676.157, Max-Change: 0.00032Iteration: 107, Log-Lik: -2676.157, Max-Change: 0.00075Iteration: 108, Log-Lik: -2676.157, Max-Change: 0.00126Iteration: 109, Log-Lik: -2676.157, Max-Change: 0.00053Iteration: 110, Log-Lik: -2676.157, Max-Change: 0.00086Iteration: 111, Log-Lik: -2676.157, Max-Change: 0.00042Iteration: 112, Log-Lik: -2676.157, Max-Change: 0.00033Iteration: 113, Log-Lik: -2676.157, Max-Change: 0.00074Iteration: 114, Log-Lik: -2676.157, Max-Change: 0.00127Iteration: 115, Log-Lik: -2676.157, Max-Change: 0.00055Iteration: 116, Log-Lik: -2676.157, Max-Change: 0.00090Iteration: 117, Log-Lik: -2676.157, Max-Change: 0.00042Iteration: 118, Log-Lik: -2676.157, Max-Change: 0.00033Iteration: 119, Log-Lik: -2676.157, Max-Change: 0.00073Iteration: 120, Log-Lik: -2676.157, Max-Change: 0.00131Iteration: 121, Log-Lik: -2676.157, Max-Change: 0.00055Iteration: 122, Log-Lik: -2676.157, Max-Change: 0.00089Iteration: 123, Log-Lik: -2676.156, Max-Change: 0.00043Iteration: 124, Log-Lik: -2676.156, Max-Change: 0.00034Iteration: 125, Log-Lik: -2676.156, Max-Change: 0.00072Iteration: 126, Log-Lik: -2676.156, Max-Change: 0.00130Iteration: 127, Log-Lik: -2676.156, Max-Change: 0.00056Iteration: 128, Log-Lik: -2676.156, Max-Change: 0.00091Iteration: 129, Log-Lik: -2676.156, Max-Change: 0.00043Iteration: 130, Log-Lik: -2676.156, Max-Change: 0.00034Iteration: 131, Log-Lik: -2676.156, Max-Change: 0.00071Iteration: 132, Log-Lik: -2676.156, Max-Change: 0.00133Iteration: 133, Log-Lik: -2676.156, Max-Change: 0.00055Iteration: 134, Log-Lik: -2676.156, Max-Change: 0.00090Iteration: 135, Log-Lik: -2676.156, Max-Change: 0.00043Iteration: 136, Log-Lik: -2676.156, Max-Change: 0.00034Iteration: 137, Log-Lik: -2676.156, Max-Change: 0.00070Iteration: 138, Log-Lik: -2676.156, Max-Change: 0.00132Iteration: 139, Log-Lik: -2676.156, Max-Change: 0.00056Iteration: 140, Log-Lik: -2676.156, Max-Change: 0.00091Iteration: 141, Log-Lik: -2676.156, Max-Change: 0.00043Iteration: 142, Log-Lik: -2676.156, Max-Change: 0.00034Iteration: 143, Log-Lik: -2676.156, Max-Change: 0.00069Iteration: 144, Log-Lik: -2676.156, Max-Change: 0.00027Iteration: 145, Log-Lik: -2676.156, Max-Change: 0.00108Iteration: 146, Log-Lik: -2676.156, Max-Change: 0.00051Iteration: 147, Log-Lik: -2676.156, Max-Change: 0.00082Iteration: 148, Log-Lik: -2676.156, Max-Change: 0.00065Iteration: 149, Log-Lik: -2676.156, Max-Change: 0.00107Iteration: 150, Log-Lik: -2676.156, Max-Change: 0.00050Iteration: 151, Log-Lik: -2676.156, Max-Change: 0.00040Iteration: 152, Log-Lik: -2676.156, Max-Change: 0.00068Iteration: 153, Log-Lik: -2676.156, Max-Change: 0.00031Iteration: 154, Log-Lik: -2676.156, Max-Change: 0.00125Iteration: 155, Log-Lik: -2676.156, Max-Change: 0.00059Iteration: 156, Log-Lik: -2676.156, Max-Change: 0.00095Iteration: 157, Log-Lik: -2676.156, Max-Change: 0.00080Iteration: 158, Log-Lik: -2676.156, Max-Change: 0.00131Iteration: 159, Log-Lik: -2676.156, Max-Change: 0.00061Iteration: 160, Log-Lik: -2676.156, Max-Change: 0.00049Iteration: 161, Log-Lik: -2676.156, Max-Change: 0.00079Iteration: 162, Log-Lik: -2676.156, Max-Change: 0.00038Iteration: 163, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 164, Log-Lik: -2676.156, Max-Change: 0.00066Iteration: 165, Log-Lik: -2676.156, Max-Change: 0.00116Iteration: 166, Log-Lik: -2676.156, Max-Change: 0.00052Iteration: 167, Log-Lik: -2676.156, Max-Change: 0.00085Iteration: 168, Log-Lik: -2676.156, Max-Change: 0.00040Iteration: 169, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 170, Log-Lik: -2676.156, Max-Change: 0.00065Iteration: 171, Log-Lik: -2676.156, Max-Change: 0.00124Iteration: 172, Log-Lik: -2676.156, Max-Change: 0.00052Iteration: 173, Log-Lik: -2676.156, Max-Change: 0.00085Iteration: 174, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 175, Log-Lik: -2676.156, Max-Change: 0.00033Iteration: 176, Log-Lik: -2676.156, Max-Change: 0.00064Iteration: 177, Log-Lik: -2676.156, Max-Change: 0.00125Iteration: 178, Log-Lik: -2676.156, Max-Change: 0.00053Iteration: 179, Log-Lik: -2676.156, Max-Change: 0.00087Iteration: 180, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 181, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 182, Log-Lik: -2676.156, Max-Change: 0.00063Iteration: 183, Log-Lik: -2676.156, Max-Change: 0.00025Iteration: 184, Log-Lik: -2676.156, Max-Change: 0.00102Iteration: 185, Log-Lik: -2676.156, Max-Change: 0.00048Iteration: 186, Log-Lik: -2676.156, Max-Change: 0.00078Iteration: 187, Log-Lik: -2676.156, Max-Change: 0.00062Iteration: 188, Log-Lik: -2676.156, Max-Change: 0.00102Iteration: 189, Log-Lik: -2676.156, Max-Change: 0.00048Iteration: 190, Log-Lik: -2676.156, Max-Change: 0.00038Iteration: 191, Log-Lik: -2676.156, Max-Change: 0.00062Iteration: 192, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 193, Log-Lik: -2676.156, Max-Change: 0.00119Iteration: 194, Log-Lik: -2676.156, Max-Change: 0.00056Iteration: 195, Log-Lik: -2676.156, Max-Change: 0.00091Iteration: 196, Log-Lik: -2676.156, Max-Change: 0.00076Iteration: 197, Log-Lik: -2676.156, Max-Change: 0.00025Iteration: 198, Log-Lik: -2676.156, Max-Change: 0.00060Iteration: 199, Log-Lik: -2676.156, Max-Change: 0.00028Iteration: 200, Log-Lik: -2676.156, Max-Change: 0.00060Iteration: 201, Log-Lik: -2676.156, Max-Change: 0.00109Iteration: 202, Log-Lik: -2676.156, Max-Change: 0.00048Iteration: 203, Log-Lik: -2676.156, Max-Change: 0.00078Iteration: 204, Log-Lik: -2676.156, Max-Change: 0.00037Iteration: 205, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 206, Log-Lik: -2676.156, Max-Change: 0.00059Iteration: 207, Log-Lik: -2676.156, Max-Change: 0.00115Iteration: 208, Log-Lik: -2676.156, Max-Change: 0.00049Iteration: 209, Log-Lik: -2676.156, Max-Change: 0.00081Iteration: 210, Log-Lik: -2676.156, Max-Change: 0.00038Iteration: 211, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 212, Log-Lik: -2676.156, Max-Change: 0.00059Iteration: 213, Log-Lik: -2676.156, Max-Change: 0.00118Iteration: 214, Log-Lik: -2676.156, Max-Change: 0.00049Iteration: 215, Log-Lik: -2676.156, Max-Change: 0.00080Iteration: 216, Log-Lik: -2676.156, Max-Change: 0.00038Iteration: 217, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 218, Log-Lik: -2676.156, Max-Change: 0.00058Iteration: 219, Log-Lik: -2676.156, Max-Change: 0.00117Iteration: 220, Log-Lik: -2676.156, Max-Change: 0.00049Iteration: 221, Log-Lik: -2676.156, Max-Change: 0.00081Iteration: 222, Log-Lik: -2676.156, Max-Change: 0.00038Iteration: 223, Log-Lik: -2676.156, Max-Change: 0.00030Iteration: 224, Log-Lik: -2676.156, Max-Change: 0.00058Iteration: 225, Log-Lik: -2676.156, Max-Change: 0.00023Iteration: 226, Log-Lik: -2676.156, Max-Change: 0.00095Iteration: 227, Log-Lik: -2676.156, Max-Change: 0.00045Iteration: 228, Log-Lik: -2676.156, Max-Change: 0.00072Iteration: 229, Log-Lik: -2676.156, Max-Change: 0.00058Iteration: 230, Log-Lik: -2676.156, Max-Change: 0.00095Iteration: 231, Log-Lik: -2676.156, Max-Change: 0.00044Iteration: 232, Log-Lik: -2676.156, Max-Change: 0.00035Iteration: 233, Log-Lik: -2676.156, Max-Change: 0.00057Iteration: 234, Log-Lik: -2676.156, Max-Change: 0.00027Iteration: 235, Log-Lik: -2676.156, Max-Change: 0.00111Iteration: 236, Log-Lik: -2676.156, Max-Change: 0.00052Iteration: 237, Log-Lik: -2676.156, Max-Change: 0.00085Iteration: 238, Log-Lik: -2676.156, Max-Change: 0.00071Iteration: 239, Log-Lik: -2676.156, Max-Change: 0.00023Iteration: 240, Log-Lik: -2676.156, Max-Change: 0.00056Iteration: 241, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 242, Log-Lik: -2676.156, Max-Change: 0.00056Iteration: 243, Log-Lik: -2676.156, Max-Change: 0.00102Iteration: 244, Log-Lik: -2676.156, Max-Change: 0.00045Iteration: 245, Log-Lik: -2676.156, Max-Change: 0.00073Iteration: 246, Log-Lik: -2676.156, Max-Change: 0.00035Iteration: 247, Log-Lik: -2676.156, Max-Change: 0.00028Iteration: 248, Log-Lik: -2676.156, Max-Change: 0.00055Iteration: 249, Log-Lik: -2676.156, Max-Change: 0.00106Iteration: 250, Log-Lik: -2676.156, Max-Change: 0.00045Iteration: 251, Log-Lik: -2676.156, Max-Change: 0.00075Iteration: 252, Log-Lik: -2676.156, Max-Change: 0.00035Iteration: 253, Log-Lik: -2676.156, Max-Change: 0.00028Iteration: 254, Log-Lik: -2676.156, Max-Change: 0.00055Iteration: 255, Log-Lik: -2676.156, Max-Change: 0.00109Iteration: 256, Log-Lik: -2676.156, Max-Change: 0.00045Iteration: 257, Log-Lik: -2676.156, Max-Change: 0.00074Iteration: 258, Log-Lik: -2676.156, Max-Change: 0.00035Iteration: 259, Log-Lik: -2676.156, Max-Change: 0.00028Iteration: 260, Log-Lik: -2676.156, Max-Change: 0.00054Iteration: 261, Log-Lik: -2676.156, Max-Change: 0.00108Iteration: 262, Log-Lik: -2676.156, Max-Change: 0.00045Iteration: 263, Log-Lik: -2676.156, Max-Change: 0.00074Iteration: 264, Log-Lik: -2676.156, Max-Change: 0.00035Iteration: 265, Log-Lik: -2676.156, Max-Change: 0.00028Iteration: 266, Log-Lik: -2676.156, Max-Change: 0.00054Iteration: 267, Log-Lik: -2676.156, Max-Change: 0.00022Iteration: 268, Log-Lik: -2676.156, Max-Change: 0.00088Iteration: 269, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 270, Log-Lik: -2676.156, Max-Change: 0.00067Iteration: 271, Log-Lik: -2676.156, Max-Change: 0.00053Iteration: 272, Log-Lik: -2676.156, Max-Change: 0.00087Iteration: 273, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 274, Log-Lik: -2676.156, Max-Change: 0.00033Iteration: 275, Log-Lik: -2676.156, Max-Change: 0.00053Iteration: 276, Log-Lik: -2676.156, Max-Change: 0.00025Iteration: 277, Log-Lik: -2676.156, Max-Change: 0.00102Iteration: 278, Log-Lik: -2676.156, Max-Change: 0.00048Iteration: 279, Log-Lik: -2676.156, Max-Change: 0.00078Iteration: 280, Log-Lik: -2676.156, Max-Change: 0.00065Iteration: 281, Log-Lik: -2676.156, Max-Change: 0.00022Iteration: 282, Log-Lik: -2676.156, Max-Change: 0.00053Iteration: 283, Log-Lik: -2676.156, Max-Change: 0.00024Iteration: 284, Log-Lik: -2676.156, Max-Change: 0.00053Iteration: 285, Log-Lik: -2676.156, Max-Change: 0.00094Iteration: 286, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 287, Log-Lik: -2676.156, Max-Change: 0.00067Iteration: 288, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 289, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 290, Log-Lik: -2676.156, Max-Change: 0.00052Iteration: 291, Log-Lik: -2676.156, Max-Change: 0.00099Iteration: 292, Log-Lik: -2676.156, Max-Change: 0.00042Iteration: 293, Log-Lik: -2676.156, Max-Change: 0.00069Iteration: 294, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 295, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 296, Log-Lik: -2676.156, Max-Change: 0.00052Iteration: 297, Log-Lik: -2676.156, Max-Change: 0.00101Iteration: 298, Log-Lik: -2676.156, Max-Change: 0.00042Iteration: 299, Log-Lik: -2676.156, Max-Change: 0.00069Iteration: 300, Log-Lik: -2676.156, Max-Change: 0.00033Iteration: 301, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 302, Log-Lik: -2676.156, Max-Change: 0.00051Iteration: 303, Log-Lik: -2676.156, Max-Change: 0.00100Iteration: 304, Log-Lik: -2676.156, Max-Change: 0.00042Iteration: 305, Log-Lik: -2676.156, Max-Change: 0.00069Iteration: 306, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 307, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 308, Log-Lik: -2676.156, Max-Change: 0.00051Iteration: 309, Log-Lik: -2676.156, Max-Change: 0.00101Iteration: 310, Log-Lik: -2676.156, Max-Change: 0.00042Iteration: 311, Log-Lik: -2676.156, Max-Change: 0.00068Iteration: 312, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 313, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 314, Log-Lik: -2676.156, Max-Change: 0.00050Iteration: 315, Log-Lik: -2676.156, Max-Change: 0.00099Iteration: 316, Log-Lik: -2676.156, Max-Change: 0.00041Iteration: 317, Log-Lik: -2676.156, Max-Change: 0.00068Iteration: 318, Log-Lik: -2676.156, Max-Change: 0.00032Iteration: 319, Log-Lik: -2676.156, Max-Change: 0.00026Iteration: 320, Log-Lik: -2676.156, Max-Change: 0.00050Iteration: 321, Log-Lik: -2676.156, Max-Change: 0.00020
##
## Calculating information matrix...
# Uji Kecocokan Item (Item Fit) dan Tabel Keputusan
# Rasch
fit_rasch <- itemfit(mod_rasch, fit_stats = "S_X2")
fit_rasch_df <- data.frame(
Butir = fit_rasch$item,
Chi = round(fit_rasch$S_X2, 2),
pvalue = round(fit_rasch$p.S_X2, 2),
Keputusan = ifelse(fit_rasch$p.S_X2 >= 0.05, "Cocok", "Tidak Cocok")
)
print(fit_rasch_df)
## Butir Chi pvalue Keputusan
## 1 b1 16.68 0.34 Cocok
## 2 b2 7.59 0.91 Cocok
## 3 b3 12.72 0.62 Cocok
## 4 b4 27.54 0.04 Tidak Cocok
## 5 b5 14.91 0.46 Cocok
## 6 b6 14.52 0.49 Cocok
## 7 b7 19.48 0.15 Cocok
## 8 b8 24.64 0.05 Cocok
## 9 b9 16.60 0.41 Cocok
## 10 b10 27.74 0.03 Tidak Cocok
## 11 b11 14.61 0.48 Cocok
## 12 b12 21.03 0.14 Cocok
## 13 b13 25.99 0.04 Tidak Cocok
## 14 b14 18.97 0.27 Cocok
## 15 b15 47.14 0.00 Tidak Cocok
## 16 b16 27.18 0.03 Tidak Cocok
## 17 b17 14.68 0.47 Cocok
## 18 b18 16.51 0.42 Cocok
## 19 b19 15.36 0.50 Cocok
## 20 b20 18.47 0.24 Cocok
## 21 b21 26.82 0.03 Tidak Cocok
## 22 b22 15.26 0.36 Cocok
## 23 b23 16.03 0.38 Cocok
## 24 b24 12.87 0.61 Cocok
## 25 b25 10.68 0.64 Cocok
## 26 b26 8.98 0.83 Cocok
## 27 b27 12.48 0.71 Cocok
## 28 b28 16.03 0.31 Cocok
## 29 b29 10.42 0.24 Cocok
## 30 b30 14.83 0.32 Cocok
## 31 b31 16.47 0.29 Cocok
## 32 b32 12.95 0.61 Cocok
## 33 b33 19.53 0.05 Cocok
## 34 b34 62.10 0.00 Tidak Cocok
## 35 b35 44.59 0.00 Tidak Cocok
## 36 b36 13.62 0.55 Cocok
## 37 b37 13.84 0.39 Cocok
cat("Jumlah butir cocok/tidak cocok model Rasch:\n")
## Jumlah butir cocok/tidak cocok model Rasch:
print(table(fit_rasch_df$Keputusan))
##
## Cocok Tidak Cocok
## 29 8
# Model 1PL
fit_1PL <- itemfit(mod_1pl, fit_stats = "S_X2")
fit_1PL_df <- data.frame(
Butir = fit_1PL$item,
Chi = round(fit_1PL$S_X2, 2),
pvalue = round(fit_1PL$p.S_X2, 2),
Keputusan = ifelse(fit_1PL$p.S_X2 >= 0.05, "Cocok", "Tidak Cocok")
)
print(fit_1PL_df)
## Butir Chi pvalue Keputusan
## 1 b1 16.68 0.34 Cocok
## 2 b2 7.59 0.91 Cocok
## 3 b3 12.72 0.62 Cocok
## 4 b4 27.54 0.04 Tidak Cocok
## 5 b5 14.91 0.46 Cocok
## 6 b6 14.52 0.49 Cocok
## 7 b7 19.48 0.15 Cocok
## 8 b8 24.64 0.05 Cocok
## 9 b9 16.60 0.41 Cocok
## 10 b10 27.74 0.03 Tidak Cocok
## 11 b11 14.61 0.48 Cocok
## 12 b12 21.03 0.14 Cocok
## 13 b13 25.99 0.04 Tidak Cocok
## 14 b14 18.97 0.27 Cocok
## 15 b15 47.16 0.00 Tidak Cocok
## 16 b16 27.19 0.03 Tidak Cocok
## 17 b17 14.68 0.47 Cocok
## 18 b18 16.51 0.42 Cocok
## 19 b19 15.36 0.50 Cocok
## 20 b20 18.47 0.24 Cocok
## 21 b21 26.82 0.03 Tidak Cocok
## 22 b22 15.26 0.36 Cocok
## 23 b23 16.03 0.38 Cocok
## 24 b24 12.87 0.61 Cocok
## 25 b25 10.68 0.64 Cocok
## 26 b26 8.98 0.83 Cocok
## 27 b27 12.48 0.71 Cocok
## 28 b28 16.03 0.31 Cocok
## 29 b29 10.43 0.24 Cocok
## 30 b30 14.83 0.32 Cocok
## 31 b31 16.47 0.29 Cocok
## 32 b32 12.95 0.61 Cocok
## 33 b33 19.54 0.05 Cocok
## 34 b34 62.10 0.00 Tidak Cocok
## 35 b35 44.59 0.00 Tidak Cocok
## 36 b36 13.62 0.55 Cocok
## 37 b37 13.84 0.39 Cocok
cat("Jumlah butir cocok/tidak cocok model 1PL:\n")
## Jumlah butir cocok/tidak cocok model 1PL:
print(table(fit_1PL_df$Keputusan))
##
## Cocok Tidak Cocok
## 29 8
# Model 2PL
fit_2PL <- itemfit(mod_2pl, fit_stats = "S_X2")
fit_2PL_df <- data.frame(
Butir = fit_2PL$item,
Chi = round(fit_2PL$S_X2, 2),
pvalue = round(fit_2PL$p.S_X2, 2),
Keputusan = ifelse(fit_2PL$p.S_X2 >= 0.05, "Cocok", "Tidak Cocok")
)
print(fit_2PL_df)
## Butir Chi pvalue Keputusan
## 1 b1 15.37 0.35 Cocok
## 2 b2 7.03 0.93 Cocok
## 3 b3 4.21 0.94 Cocok
## 4 b4 25.13 0.02 Tidak Cocok
## 5 b5 9.36 0.75 Cocok
## 6 b6 13.03 0.45 Cocok
## 7 b7 17.50 0.29 Cocok
## 8 b8 15.41 0.28 Cocok
## 9 b9 15.46 0.42 Cocok
## 10 b10 25.91 0.02 Tidak Cocok
## 11 b11 13.92 0.38 Cocok
## 12 b12 19.70 0.10 Cocok
## 13 b13 26.33 0.05 Tidak Cocok
## 14 b14 13.05 0.22 Cocok
## 15 b15 23.71 0.31 Cocok
## 16 b16 20.45 0.43 Cocok
## 17 b17 14.90 0.67 Cocok
## 18 b18 16.34 0.50 Cocok
## 19 b19 15.12 0.52 Cocok
## 20 b20 12.90 0.74 Cocok
## 21 b21 18.73 0.47 Cocok
## 22 b22 14.01 0.23 Cocok
## 23 b23 16.32 0.50 Cocok
## 24 b24 13.11 0.52 Cocok
## 25 b25 8.19 0.77 Cocok
## 26 b26 8.01 0.89 Cocok
## 27 b27 7.91 0.79 Cocok
## 28 b28 16.71 0.27 Cocok
## 29 b29 7.50 0.59 Cocok
## 30 b30 10.21 0.42 Cocok
## 31 b31 17.66 0.22 Cocok
## 32 b32 10.87 0.54 Cocok
## 33 b33 19.76 0.03 Tidak Cocok
## 34 b34 21.93 0.46 Cocok
## 35 b35 20.40 0.31 Cocok
## 36 b36 11.40 0.65 Cocok
## 37 b37 12.60 0.63 Cocok
cat("Jumlah butir cocok/tidak cocok model 2PL:\n")
## Jumlah butir cocok/tidak cocok model 2PL:
print(table(fit_2PL_df$Keputusan))
##
## Cocok Tidak Cocok
## 33 4
# Model 3PL
fit_3PL <- itemfit(mod_3pl, fit_stats = "S_X2")
fit_3PL_df <- data.frame(
Butir = fit_3PL$item,
Chi = round(fit_3PL$S_X2, 2),
pvalue = round(fit_3PL$p.S_X2, 2),
Keputusan = ifelse(fit_3PL$p.S_X2 >= 0.05, "Cocok", "Tidak Cocok")
)
print(fit_3PL_df)
## Butir Chi pvalue Keputusan
## 1 b1 12.30 0.42 Cocok
## 2 b2 6.98 0.86 Cocok
## 3 b3 3.25 0.92 Cocok
## 4 b4 16.00 0.10 Cocok
## 5 b5 5.52 0.85 Cocok
## 6 b6 13.13 0.44 Cocok
## 7 b7 17.01 0.20 Cocok
## 8 b8 19.48 0.08 Cocok
## 9 b9 15.20 0.23 Cocok
## 10 b10 18.27 0.05 Cocok
## 11 b11 13.92 0.31 Cocok
## 12 b12 17.83 0.12 Cocok
## 13 b13 25.72 0.03 Tidak Cocok
## 14 b14 12.19 0.20 Cocok
## 15 b15 22.33 0.27 Cocok
## 16 b16 15.70 0.55 Cocok
## 17 b17 14.45 0.70 Cocok
## 18 b18 12.02 0.74 Cocok
## 19 b19 15.84 0.46 Cocok
## 20 b20 13.14 0.66 Cocok
## 21 b21 17.60 0.55 Cocok
## 22 b22 13.61 0.33 Cocok
## 23 b23 15.25 0.51 Cocok
## 24 b24 11.00 0.53 Cocok
## 25 b25 8.56 0.66 Cocok
## 26 b26 7.20 0.93 Cocok
## 27 b27 7.22 0.78 Cocok
## 28 b28 18.83 0.09 Cocok
## 29 b29 10.96 0.28 Cocok
## 30 b30 13.96 0.08 Cocok
## 31 b31 13.26 0.35 Cocok
## 32 b32 9.21 0.51 Cocok
## 33 b33 13.23 0.15 Cocok
## 34 b34 20.38 0.50 Cocok
## 35 b35 21.03 0.33 Cocok
## 36 b36 11.89 0.29 Cocok
## 37 b37 12.83 0.54 Cocok
cat("Jumlah butir cocok/tidak cocok model 3PL:\n")
## Jumlah butir cocok/tidak cocok model 3PL:
print(table(fit_3PL_df$Keputusan))
##
## Cocok Tidak Cocok
## 36 1
# Model 4PL
fit_4PL <- itemfit(mod_4pl, fit_stats = "S_X2")
fit_4PL_df <- data.frame(
Butir = fit_4PL$item,
Chi = round(fit_4PL$S_X2, 2),
pvalue = round(fit_4PL$p.S_X2, 2),
Keputusan = ifelse(fit_4PL$p.S_X2 >= 0.05, "Cocok", "Tidak Cocok")
)
print(fit_4PL_df)
## Butir Chi pvalue Keputusan
## 1 b1 12.62 0.32 Cocok
## 2 b2 8.03 0.78 Cocok
## 3 b3 3.85 0.80 Cocok
## 4 b4 16.15 0.10 Cocok
## 5 b5 6.63 0.68 Cocok
## 6 b6 13.98 0.38 Cocok
## 7 b7 17.54 0.18 Cocok
## 8 b8 19.22 0.08 Cocok
## 9 b9 16.25 0.18 Cocok
## 10 b10 17.17 0.05 Tidak Cocok
## 11 b11 15.46 0.16 Cocok
## 12 b12 16.89 0.15 Cocok
## 13 b13 25.89 0.02 Tidak Cocok
## 14 b14 11.15 0.27 Cocok
## 15 b15 25.32 0.09 Cocok
## 16 b16 15.54 0.49 Cocok
## 17 b17 10.58 0.72 Cocok
## 18 b18 12.57 0.56 Cocok
## 19 b19 11.25 0.67 Cocok
## 20 b20 12.11 0.67 Cocok
## 21 b21 16.84 0.60 Cocok
## 22 b22 13.57 0.19 Cocok
## 23 b23 14.90 0.53 Cocok
## 24 b24 11.07 0.44 Cocok
## 25 b25 7.62 0.57 Cocok
## 26 b26 6.71 0.88 Cocok
## 27 b27 7.54 0.67 Cocok
## 28 b28 16.15 0.10 Cocok
## 29 b29 10.50 0.16 Cocok
## 30 b30 12.42 0.19 Cocok
## 31 b31 12.98 0.29 Cocok
## 32 b32 7.93 0.64 Cocok
## 33 b33 4.50 0.72 Cocok
## 34 b34 19.27 0.50 Cocok
## 35 b35 18.39 0.30 Cocok
## 36 b36 13.16 0.21 Cocok
## 37 b37 13.58 0.40 Cocok
cat("Jumlah butir cocok/tidak cocok model 4PL:\n")
## Jumlah butir cocok/tidak cocok model 4PL:
print(table(fit_4PL_df$Keputusan))
##
## Cocok Tidak Cocok
## 35 2
# Plot Empirical untuk Butir ke-9 pada Setiap Model
itemfit(mod_rasch, 'S_X2', empirical.plot = 9)
itemfit(mod_1pl, 'S_X2', empirical.plot = 9)
itemfit(mod_2pl, 'S_X2', empirical.plot = 9)
itemfit(mod_3pl, 'S_X2', empirical.plot = 9)
itemfit(mod_4pl, 'S_X2', empirical.plot = 9)
# Visualisasi Plot Semua Butir
for (coba in 1:nrow(fit_rasch_df)) {
print(itemfit(mod_rasch, 'S_X2', empirical.plot = coba))
}
# Uji Kecocokan Model dengan ANOVA (AIC, BIC)
anova(mod_rasch, mod_1pl)
## AIC SABIC HQ BIC logLik X2 df p
## mod_rasch 5719.360 5713.502 5765.839 5833.764 -2821.68
## mod_1pl 5719.381 5713.522 5765.859 5833.785 -2821.69 -0.02 0 NaN
anova(mod_1pl, mod_2pl)
## AIC SABIC HQ BIC logLik X2 df p
## mod_1pl 5719.381 5713.522 5765.859 5833.785 -2821.690
## mod_2pl 5634.631 5623.222 5725.142 5857.418 -2743.315 156.75 36 0
anova(mod_2pl, mod_3pl)
## AIC SABIC HQ BIC logLik X2 df p
## mod_2pl 5634.631 5623.222 5725.142 5857.418 -2743.315
## mod_3pl 5600.438 5583.324 5736.205 5934.618 -2689.219 108.193 37 0
anova(mod_3pl, mod_4pl)
## AIC SABIC HQ BIC logLik X2 df p
## mod_3pl 5600.438 5583.324 5736.205 5934.618 -2689.219
## mod_4pl 5648.312 5625.494 5829.334 6093.886 -2676.156 26.126 37 0.909
library(mirt)
# Model null (baseline): item independen, satu iterasi
mod_null <- mirt(datasampel, 1, itemtype = 'Rasch', technical = list(NCYCLES = 1), SE = FALSE)
## Iteration: 1, Log-Lik: -2858.002, Max-Change: 0.88713
## Warning: EM cycles terminated after 1 iterations.
logLik_null <- as.numeric(logLik(mod_null))
# Log-likelihood tiap model
logLik_rasch <- as.numeric(logLik(mod_rasch))
logLik_1pl <- as.numeric(logLik(mod_1pl))
logLik_2pl <- as.numeric(logLik(mod_2pl))
logLik_3pl <- as.numeric(logLik(mod_3pl))
logLik_4pl <- as.numeric(logLik(mod_4pl))
# Fungsi NFI
nfi <- function(logLik_null, logLik_model) {
return((logLik_null - logLik_model) / logLik_null)
}
# Hitung NFI masing-masing
nfi_rasch <- nfi(logLik_null, logLik_rasch)
nfi_1pl <- nfi(logLik_null, logLik_1pl)
nfi_2pl <- nfi(logLik_null, logLik_2pl)
nfi_3pl <- nfi(logLik_null, logLik_3pl)
nfi_4pl <- nfi(logLik_null, logLik_4pl)
# Hitung AIC dan BIC via anova
anova_rasch_1pl <- anova(mod_rasch, mod_1pl)
anova_1pl_2pl <- anova(mod_1pl, mod_2pl)
anova_2pl_3pl <- anova(mod_2pl, mod_3pl)
anova_3pl_4pl <- anova(mod_3pl, mod_4pl)
get_anova_ic <- function(anova_res, idx = 1) {
list(AIC = anova_res$AIC[idx], BIC = anova_res$BIC[idx])
}
AIC_rasch <- get_anova_ic(anova_rasch_1pl, 1)$AIC
BIC_rasch <- get_anova_ic(anova_rasch_1pl, 1)$BIC
AIC_1pl <- get_anova_ic(anova_rasch_1pl, 2)$AIC
BIC_1pl <- get_anova_ic(anova_rasch_1pl, 2)$BIC
AIC_2pl <- get_anova_ic(anova_1pl_2pl, 2)$AIC
BIC_2pl <- get_anova_ic(anova_1pl_2pl, 2)$BIC
AIC_3pl <- get_anova_ic(anova_2pl_3pl, 2)$AIC
BIC_3pl <- get_anova_ic(anova_2pl_3pl, 2)$BIC
AIC_4pl <- get_anova_ic(anova_3pl_4pl, 2)$AIC
BIC_4pl <- get_anova_ic(anova_3pl_4pl, 2)$BIC
# Fungsi jumlah cocok/tidak cocok
count_fit <- function(df) {
cocok <- sum(df$Keputusan == "Cocok", na.rm = TRUE)
tidak <- sum(df$Keputusan == "Tidak Cocok", na.rm = TRUE)
na <- sum(is.na(df$Keputusan))
return(c(cocok, tidak, na))
}
fit_counts <- rbind(
count_fit(fit_rasch_df),
count_fit(fit_1PL_df),
count_fit(fit_2PL_df),
count_fit(fit_3PL_df),
count_fit(fit_4PL_df)
)
# Rekapitulasi
n_sampel <- nrow(datasampel)
n_item <- ncol(datasampel)
rekap_model <- data.frame(
Model = c("Rasch", "1PL", "2PL", "3PL", "4PL"),
n = rep(n_sampel, 5),
AIC = c(AIC_rasch, AIC_1pl, AIC_2pl, AIC_3pl, AIC_4pl),
BIC = c(BIC_rasch, BIC_1pl, BIC_2pl, BIC_3pl, BIC_4pl),
Butir_Cocok = fit_counts[,1],
Butir_Tidak_Cocok = fit_counts[,2],
Butir_NA = fit_counts[,3],
Jumlah_Variabel = rep(n_item, 5),
NFI_maks = round(c(nfi_rasch, nfi_1pl, nfi_2pl, nfi_3pl, nfi_4pl), 4)
)
print(rekap_model)
## Model n AIC BIC Butir_Cocok Butir_Tidak_Cocok Butir_NA
## 1 Rasch 150 5719.360 5833.764 29 8 0
## 2 1PL 150 5719.381 5833.785 29 8 0
## 3 2PL 150 5634.631 5857.418 33 4 0
## 4 3PL 150 5600.438 5934.618 36 1 0
## 5 4PL 150 5648.312 6093.886 35 2 0
## Jumlah_Variabel NFI_maks
## 1 37 0.0021
## 2 37 0.0021
## 3 37 0.0298
## 4 37 0.0489
## 5 37 0.0535