23.02.2026
dilek vural
2026-02-25
DERSTE YAPILANLARIN TEKRARI
Öncelikle derste yapılanları en baştan tekrar uyguladım, önemli bilgileri özetledim.
Kayıp veri bulma, atama ve t test denemeleri
Veri setini aktar
TRUSA <- readRDS("TRUSA.RDS")Eksik veri kontrolü
library(naniar)
miss_var_table(TRUSA)Toplam değişkeni oluşturma
SUM <-TRUSA %>% select(starts_with("M0")) %>% rowSums()
TRUSA$SUM <- SUMToplam değişkeninin ülkelere göre, betimleyici analiz sonuçları
TRUSA %>%
group_by(CNT) %>% summarise(n=n(), ort=mean(SUM), sd=sd(SUM), min=min(SUM), max=max(SUM))2.yol
describe.by(TRUSA$SUM, TRUSA$CNT)##
## Descriptive statistics by group
## group: TUR
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 435 13.45 7.57 11 12.74 7.41 2 32 30 0.71 -0.61 0.36
## ------------------------------------------------------------
## group: USA
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 716 17.05 7.53 17 16.92 8.9 1 34 33 0.1 -0.9 0.28# independent sample t test
library(effsize) # to calculate effect size
t_test_result <- t.test(SUM ~ CNT, data = TRUSA, var.equal = TRUE)
print(t_test_result)##
## Two Sample t-test
##
## data: SUM by CNT
## t = -7.8348, df = 1149, p-value = 1.064e-14
## alternative hypothesis: true difference in means between group TUR and group USA is not equal to 0
## 95 percent confidence interval:
## -4.493049 -2.693382
## sample estimates:
## mean in group TUR mean in group USA
## 13.45287 17.04609# calculate cohen's d
cohen_d_result <- cohen.d(TRUSA$SUM[TRUSA$CNT=="TUR"],
TRUSA$SUM[TRUSA$CNT=="USA"],
pooled=TRUE)Veri Setinde Yüzde 5,10,15 değerinde eksik veri oluşturma
library(mvdalab)
TRUSA_5 <- introNAs(TRUSA, percent = 5)
TRUSA_10 <- introNAs(TRUSA, percent = 10)
TRUSA_15 <- introNAs(TRUSA, percent = 15)Oluşturduğumuz veri setlerinde her bir sütununda toplam kaç tane eksik (NA) olduğunu hesaplama:
TRUSA_5%>%is.na() %>% colSums()## IDSTUD IDBOOK M042182 M042081 M042049 M042052 M042076 M042302A
## 60 51 50 65 69 57 59 60
## M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268
## 48 54 56 62 50 60 70 58
## M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D
## 60 56 56 41 58 50 72 61
## M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202
## 59 53 65 55 61 75 59 56
## M062246 M062286 M062325 M062106 M062124 CNT SUM
## 59 52 55 59 50 56 47TRUSA_10%>% is.na() %>% colSums()## IDSTUD IDBOOK M042182 M042081 M042049 M042052 M042076 M042302A
## 104 111 116 119 116 112 114 115
## M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268
## 126 116 123 115 122 108 105 111
## M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D
## 108 123 103 109 113 114 119 110
## M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202
## 118 122 100 126 91 115 114 112
## M062246 M062286 M062325 M062106 M062124 CNT SUM
## 111 149 128 120 125 115 111TRUSA_15%>% is.na() %>% colSums()## IDSTUD IDBOOK M042182 M042081 M042049 M042052 M042076 M042302A
## 161 201 175 173 180 145 166 159
## M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268
## 176 165 173 171 177 157 167 176
## M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D
## 183 194 159 177 178 182 184 180
## M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202
## 184 179 180 161 166 171 176 165
## M062246 M062286 M062325 M062106 M062124 CNT SUM
## 152 157 188 192 153 180 170# Bu testlerin hepsi hata verdi; mcar_test (TRUSA_5); mcar_test (TRUSA_10); mcar_test (TRUSA_15).
# Aşağıdaki kodu çalıştırınca oldu. Tüm sonuçlar istatistiksel olarak anlamlı değil, yani kayıp veri rastgele.
TRUSA_5 %>%
select(starts_with("M0")) %>%
mcar_test()TRUSA_10 %>%
select(starts_with("M0")) %>%
mcar_test()TRUSA_15 %>%
select(starts_with("M0")) %>%
mcar_test()Kayıp veri rastgele ve az sayıda ise liste bazı silme yapılır:
TRUSA_5_lw<- na.omit(TRUSA_5)
TRUSA_10_lw<- na.omit(TRUSA_10)
TRUSA_15_lw<- na.omit(TRUSA_15)Eksik veri olan veri setinde t-test ve etki büyüklüğü hesabı: Eksik veri olan veri setinde(yüzde 15 olanda) yapamadık t test çünkü elimizde t test için yeterli veri kalmadı.
t_test_result_1<- t.test(SUM ~ CNT, data=TRUSA_5_lw, var.equal = TRUE)
print(t_test_result_1)##
## Two Sample t-test
##
## data: SUM by CNT
## t = -2.7177, df = 156, p-value = 0.007318
## alternative hypothesis: true difference in means between group TUR and group USA is not equal to 0
## 95 percent confidence interval:
## -5.9795980 -0.9459339
## sample estimates:
## mean in group TUR mean in group USA
## 13.25000 16.71277# calculate cohen's d
cohen_d_result_1 <- effsize::cohen.d(TRUSA_5_lw$SUM[TRUSA_5_lw$CNT=="TUR"],
TRUSA_5_lw$SUM[TRUSA_5_lw$CNT=="USA"],
pooled=TRUE)
print(cohen_d_result_1)##
## Cohen's d
##
## d estimate: -0.4404269 (small)
## 95 percent confidence interval:
## lower upper
## -0.7642605 -0.1165933Kayıp veri fazlaysa multiple imputation yap:
library(dplyr)
library(mice)
TRUSA_5_im1<- mice(TRUSA_5[,3:37], m=5, maxit=50, method='logreg', seed=500)##
## iter imp variable
## 1 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 1 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 1 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 1 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 1 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 2 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 2 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 2 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 2 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 2 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 3 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 3 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 3 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 3 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 3 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 4 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 4 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 4 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 4 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 4 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 5 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 5 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 5 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 5 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 5 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 6 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 6 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 6 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 6 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 6 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 7 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 7 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 7 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 7 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 7 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 8 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 8 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 8 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 8 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 8 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 9 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 9 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 9 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 9 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 9 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 10 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 10 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 10 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 10 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 10 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 11 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 11 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 11 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 11 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 11 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 12 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 12 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 12 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 12 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 12 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 13 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 13 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 13 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 13 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 13 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 14 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 14 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 14 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 14 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 14 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 15 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 15 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 15 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 15 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 15 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 16 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 16 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 16 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 16 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 16 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 17 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 17 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 17 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 17 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 17 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 18 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 18 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 18 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 18 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 18 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 19 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 19 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 19 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 19 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 19 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 20 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 20 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 20 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 20 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 20 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 21 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 21 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 21 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 21 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 21 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 22 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 22 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 22 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 22 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 22 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 23 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 23 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 23 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 23 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 23 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 24 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 24 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 24 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 24 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 24 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 25 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 25 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 25 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 25 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 25 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 26 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 26 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 26 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 26 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 26 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 27 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 27 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 27 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 27 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 27 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 28 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 28 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 28 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 28 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 28 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 29 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 29 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 29 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 29 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 29 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 30 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 30 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 30 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 30 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 30 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 31 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 31 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 31 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 31 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 31 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 32 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 32 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 32 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 32 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 32 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 33 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 33 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 33 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 33 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 33 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 34 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 34 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 34 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 34 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 34 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 35 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 35 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 35 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 35 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 35 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 36 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 36 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 36 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 36 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 36 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 37 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 37 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 37 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 37 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 37 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 38 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 38 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 38 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 38 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 38 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 39 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 39 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 39 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 39 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 39 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 40 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 40 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 40 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 40 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 40 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 41 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 41 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 41 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 41 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 41 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 42 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 42 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 42 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 42 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 42 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 43 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 43 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 43 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 43 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 43 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 44 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 44 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 44 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 44 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 44 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 45 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 45 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 45 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 45 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 45 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 46 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 46 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 46 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 46 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 46 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 47 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 47 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 47 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 47 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 47 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 48 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 48 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 48 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 48 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 48 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 49 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 49 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 49 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 49 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 49 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 50 1 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 50 2 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 50 3 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 50 4 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124
## 50 5 M042182 M042081 M042049 M042052 M042076 M042302A M042302B M042302C M042100 M042202 M042240 M042093 M042271 M042268 M042159 M042164 M042167 M062208 M062208A M062208B M062208C M062208D M062153 M062111A M062111B M062237 M062314 M062074 M062183 M062202 M062246 M062286 M062325 M062106 M062124completed_data_1<-complete(TRUSA_5_im1,1)
completed_data_2<-complete(TRUSA_5_im1,2)
completed_data_3<-complete(TRUSA_5_im1,3)
completed_data_4<-complete(TRUSA_5_im1,4)
completed_data_5<-complete(TRUSA_5_im1,5)Normallik
Eksik veri doldurduktan sonra 2. aşama uç değer kontrolü yapmak. Tek değişkenlik uç değer kontrolü için Z puan kontrolü yap: -3 ve +3 aralığında olmalı, veri 1000’den fazla ise +4 ve -4 arası kabul edilir.
Eğer bir uç değer örneklem temsili sağlıyorsa tutabilirsin. Uç değeri kullanmak istersen, etkisini azaltmak gerekir. Logaritmik dönüşümü kullanabilirsin ama yorumunu dönüşmüş değişkene göre yapman gerekiyor.
Doğası gereği çarpık dağılan bir değişken olabilir mesela madde kullanımına ilişkin bir madde, onu kabul edebilrisin sınırlılık olarak yazarsın.
Çok değişkenli uç değer kontrolü: İki veya daha fazla değişkenin alışılmadık kombinasyonudur. Mahalanobis değerini hesapladıktan sonra kimin “uç değer” olup kimin olmadığını belirlemek için bir sınır çizmek gerekir. Bunun için Ki-Kare Dağılımı kullanılır.Kritik değeri geçen çok değişkenli uç değerdir.
- aşama normallik varsayımı kontrolüdür. Normalden sapan değişken varsa önemli bir değişken değilse kullanmayabilirsin, eğer önemli bir değişkense dönüşüm yapılabilir ama yorumlama yine buna göre olmalı.Genelde dönüştürmeyi çok yapmıyouz yorumu zor diye.
Normallik kontrolü için: Çarpıklık/Skewness ve Basıklık/Kurtosis değerlerİ, histogram, q,q plot kontrol edilir, doğrusallık kontorolü için scatter plot incelenir.Çarpıklık basıklık:-1,+1 arasında ise normal dağılıma yakın denebilir. -2,+2 aralığı yüksek oran kabul etme.
Çoklu bağlantı sorunu: Birbirleriyle çok yüksek korelasyona sahip değişkenler varsa,sonuçları çarpıtabilir.Bunları temel bileşenler analizi ile bir araya getirip tek değişken haline getirebiliriz.
Kolmogorov-Smirnov (K-S) testi, Shapiro-Wilk testi normallik varsayımını kontrol etmek için artık güncel değil(örneklem büyüklüğünden çok etkileniyor) (Howell, 2012).
library(tidyverse)
library(stevemisc)
library(knitr)
library(haven)
library(summarytools)
library(outliers)
library(ggplot2)
library(plotly)
library(ggpmisc)
library(psych)
library(sur)
library(moments)
library(corrplot)
library(olsrr)library(dplyr)
library(haven) # SPSS dosyalarını R ortamına aktarmak için haven paketini kullanın.
screen <- read_sav("SCREEN.sav")
screen <- expss::drop_var_labs(screen)
head(screen) # Veri setinin ilk birkaç satırını görüntülelibrary(dplyr)
screen <- screen %>%
mutate(INCOME = ifelse(is.na(INCOME), mean(INCOME, na.rm =TRUE),INCOME)) %>% na.omit()
summary(screen)## SUBNO TIMEDRS ATTDRUG ATTHOUSE
## Min. : 1.0 Min. : 0.000 Min. : 5.00 Min. : 2.00
## 1st Qu.:136.8 1st Qu.: 2.000 1st Qu.: 7.00 1st Qu.:21.00
## Median :313.5 Median : 4.000 Median : 8.00 Median :24.00
## Mean :317.3 Mean : 7.914 Mean : 7.69 Mean :23.54
## 3rd Qu.:483.2 3rd Qu.:10.000 3rd Qu.: 9.00 3rd Qu.:27.00
## Max. :758.0 Max. :81.000 Max. :10.00 Max. :35.00
## INCOME EMPLMNT MSTATUS RACE
## Min. : 1.000 Min. :0.000 Min. :1.00 Min. :1.000
## 1st Qu.: 3.000 1st Qu.:0.000 1st Qu.:2.00 1st Qu.:1.000
## Median : 4.000 Median :0.000 Median :2.00 Median :1.000
## Mean : 4.208 Mean :0.472 Mean :1.78 Mean :1.086
## 3rd Qu.: 6.000 3rd Qu.:1.000 3rd Qu.:2.00 3rd Qu.:1.000
## Max. :10.000 Max. :1.000 Max. :2.00 Max. :2.000Kategorik değiken için uç değer kontrolü:
race değişkenin frekans tablosu: Yorum: iki kategorili değişkenlerden race değişkeninin kategorilere dağılımları incelendiğinde kategoriler arasında yaklaşık 10.1:1 (91/9) oranı olduğu görülmektedir.Bir ırk nerdeyse hiç temsil edilmemiş gözüküyor.Normal popülasyonda nasıl temsil onu kontrol et, senin verinde bu temsil sağlanmış mı? Bu oran oldukça yüksektir. Değişken araştırma için önemli değilse çıkarılabilir, aksi halde değişkenle ilgili sonuçlar yorumlanırken bu durum göz önüne alınmalıdır.
table(screen$RACE)##
## 1 2
## 424 40library(summarytools)
freq(screen$RACE,
round.digits=2,report.nas = FALSE,
style = "rmarkdown") ## ### Frequencies
## #### screen$RACE
## **Type:** Numeric
##
## | | Freq | % | % Cum. |
## |----------:|-----:|-------:|-------:|
## | **1** | 424 | 91.38 | 91.38 |
## | **2** | 40 | 8.62 | 100.00 |
## | **Total** | 464 | 100.00 | 100.00 |Mstatus değişkeni yorum:İki kategorili değişkenlerden mstatus değişkeninin kategorilere dağılımları incelendiğinde kategoriler arasında yaklaşık 3.5:1 (78/22) oranı olduğu görülmektedir. Bu oran kabul edilebilir bir orandır.
library(knitr)
freq(screen$MSTATUS,report.nas = FALSE) %>%
kable(format='markdown',
caption="Frekans Tablosu",digits = 2)| Freq | % Valid | % Valid Cum. | % Total | % Total Cum. | |
|---|---|---|---|---|---|
| 1 | 102 | 21.98 | 21.98 | 21.98 | 21.98 |
| 2 | 362 | 78.02 | 100.00 | 78.02 | 100.00 |
| 0 | NA | NA | 0.00 | 100.00 | |
| Total | 464 | 100.00 | 100.00 | 100.00 | 100.00 |
Sürekli değişkenlerde uç değerler
Sürekli değişkenler için tek değişkenli uç değerleri belirlemenin bir yolu, değişkene ait bütün değerlerin ortalama 0, standart sapma 1 olacak şekilde standart değerlere (z puanlarına) dönüştürülmesidir. Tek değişkenli uç değerler çok büyük z puanlarına sahiptirler.
outliers paketinde scores() fonksiyonu ile z değerleri hesaplanabilir.
library(outliers)
z.scores <- screen %>%
select(2:5) %>%
scores(type = "z") %>%
round(2)
head(z.scores)summarytools paketinde descr() fonksiyonu ile z değerlerinin minumum ve maksimum değerleri incelenebilir. Yorum: timedrs değişkeni için z puanlarının maksimum değerin 4.0’ten büyük, atthouse değişkeni z puanlarının içinse minimum değerin 4.0’ten küçük olduğu görülmektedir. Diğer değişkenler için değerler beklenen sınırlar içerisindedir.
summarytools::descr(z.scores,
stats = c("min", "max"),
transpose = TRUE,
headings = FALSE) Z skor için diğer bir yol, interaktif tablo.
library(DT)
DT::datatable(z.scores)Tek değişkenli uç değerleri saptamak için grafiksel yöntemlerden de yararlanılabilir (Örneğin, histogramlar, kutu grafikleri, normal olasılık grafikleri gibi).
# timedrs değişkenin incelenmesi
library(ggplot2)
ggplot(screen, aes(x = TIMEDRS)) +
geom_histogram(bins = 30L, fill = "#0c4c8a")
Grafik üzerinde ortalamanın gösterilmesi
library(ggpmisc)
ggplot(screen, aes(x = TIMEDRS)) + geom_histogram() +
geom_vline(xintercept =7.914, color = "red",
linetype = "dashed") +
annotate("text", label = "Ort = 7.913", x = 10, y = 100, color ="black")
Yoğunluk grafiği
ggplot(screen, aes(x = TIMEDRS)) +
geom_histogram(aes(y=..density..))+
geom_density(alpha=.5, fill="#0c4c8a") +
theme_minimal()
Interaktif Grafik
library(plotly)
plot_ly(x = screen$TIMEDRS, type = "histogram",
histnorm = "probability")Kutu Grafik
ggplot(screen, aes(y = TIMEDRS)) +
geom_boxplot() 
boxplot.stats fonksiyonun out bileşini uç değerleri vermektedir.
out <- boxplot.stats(screen$TIMEDRS)$out
out## [1] 60 23 39 33 38 34 27 30 25 49 60 27 27 52 24 57 52 58 57 43 37 75 29 30 25
## [26] 37 56 29 37 81 27 23Bu uç değerlerin timedrs değişkeninde hangi gözlemlerde olduğu ise aşağıdaki kodlarla belirlenebilir.
out_ind <- which(screen$TIMEDRS %in% c(out))
out_ind## [1] 40 64 67 76 79 96 102 117 150 163 168 170 178 193 203 206 213 249 274
## [20] 278 285 289 309 342 344 362 367 374 388 404 408 443İnteraktif kutu grafiği
plot_ly(y = screen$TIMEDRS, type = 'box') İnteraktif kutu grafiği üzerine uç değerlerin idlerinin işaretlenmesi ise aşağıdaki kodla yapılabilir.
plot_ly(y = screen$TIMEDRS, type = 'box') %>%
layout(title = 'Box Plot',
annotations = list( x = -0.01, y = boxplot.stats(screen$TIMEDRS)$out,
text = paste(out_ind), showarrow = FALSE,
xanchor = "right"))timedrs değişkenin mstatus değişkenine göre incelenmesi
ggplot(screen, aes(x = factor(MSTATUS),
y = TIMEDRS, fill = factor(MSTATUS))) +
geom_boxplot() +
theme_minimal()
Attdrug değişkenin incelenmesi
ggplot(screen) + aes(x = ATTDRUG) +
geom_histogram( bins = 6, fill = "#0c4c8a")+
theme_minimal()
Atthouse değişkenin incelenmesi
ggplot(screen) +
aes(x = ATTHOUSE) +
geom_histogram( bins = 10, fill = "darkgreen") +
theme_minimal() 
Hem z puanı hem grafik incele beraber karar ver.
plot_ly(y = screen$ATTHOUSE, type = 'box') Veri setinde potansiyel tek değişkenli uç değerler tespit edildiğinde, önce uç değerin nedeni araştırılmalıdır. Eğer veri girişinde hata varsa veya kayıp veri kodlanırken hata yapıldıysa düzeltilmelidir.
Bunun dışındaki nedenlerde değişkenin dönüştürülmesinin uygun olup olmayacağına karar verilmelidir.
Dönüşümler hem dağılımların normalliğini geliştirir hem de tek değişkenli uç değerleri dağılımın merkezine çekerler ve etkisini azaltırlar.
Dönüşüme karar verilirse çok değişkenli uç değerler incelenmeden dönüşüm yapılmalıdır. Çünkü çok değişkenli uç değerlerin belirlenmesinde kullanılan istatistikler normal dağılımı gerektirir.
Veri setinde iki değişken – timedrs ve atthouse uç değerlere sahiptir.
timedrs değişkeni için uç değer olarak belirlenen değerlerin beklenen değerlerin üstünde olduğu ancak veri girişinde hata bulunmadığı rapor edilmiş, bu değerlere sahip bireylerin veri setinde kalmasına karar verilmiştir.
atthouse değişkeni için uç değerler olarak belirlenen değerler diğer değerlerden kopuktur. Bu değerlerin evren için beklenen değerler mi olduğuna veya veri girişinde hata yapılıp yapılmadığına karar verilmelidir.
Her iki durumda da veri setinde 260. ve 296. satırda yer alan 2 birey (346 ve 407 subno.lu bireyler) veri setinden çıkarılabilir. 2 bireyin veri setinden çıkarılması sonucu örneklem büyüklüğü 462’ye eşit olacaktır.
screen[c(260,298),]screen2 <- screen[-c(260,298),]Mahalanobis Uzaklığı Çok değişkenli uç değerleri belirlemenin bir yolu Mahalanobis uzaklığını hesaplamaktır. Mahalanobis uzaklığı z puanının çok boyutlu versiyonudur. Bir gözlemin, dağılımın kovaryansı (çok boyutlu varyansı) verildiğinde, dağılımın ağırlık merkezinden (çok boyutlu ortalamasından) uzaklığını ölçer.
Mahalonobis uzaklığı ki-kare dağılımı gösterir (serbestlik derecesi hesaplamada kullanılan değişken sayısına eşittir) ve ki-kare dağılımı kullanılarak değerlendirilebilir. Eğer hesaplanan Mahalonobis uzaklığının gözlenme olasılığı 0.001 veya daha küçükse gözlem uçdeğerdir.
Bu yöntem eşit aralık veya eşit oran düzeyinde ölçülen değişkenler için veya sürekli değişken olarak ele alınan sıralama ölçeğinde ölçülen değişkenler için geçerli olup sınıflama düzeyinde ölçülen değişkenler için geçerli değildir.
# Mahalanobis uzaklığı hesaplama
library(psych)
veri <- screen2[,1:5]
md <- mahalanobis(veri, center = colMeans(veri), cov = cov(veri))
head(md,20)## [1] 3.785517 4.541493 3.501077 7.281365 5.457240 2.896550 5.807898
## [8] 3.879478 4.751166 7.415405 10.602100 5.249121 6.073732 3.271885
## [15] 12.316463 4.440749 4.836160 6.362806 4.126524 10.797545# Mahalanobis uzaklığı kritik değer belirleme
library(psych)
alpha <- .001
cutoff <- (qchisq(p = 1 - alpha, df = ncol(veri)))
cutoff## [1] 20.51501# Mahalonobis uzaklığı değerleri ki-kare ile değerlendirilir (serbestlik derecesi bağımsız değişken sayısına eşittir). Buna göre 20.51501 kritik değerinden büyük olan değerler 0.001 alfa düzeyinde istatistiksel olarak anlamlı olarak değerlendirilir.
ucdegerler <- which(md > cutoff)
veri[ucdegerler, ]veri[40,]data_temiz <- veri[-ucdegerler, ]
# 548, 398, 48, 235, 330, 502, 276, 291 ve 370 subno.lu bireyler için Mahalonobis uzaklık değerleri kritik değerden büyüktür. Bu gözlemler çok değişkenli uç değerler olarak değerlendirilir.Çok değişkenlik normallik sayıltısı
Çok degiskenli normallik sayıltısını test etmek için doğrudan bir test bulunmadığından, genellikle her bir değişken ayrı ayrı test edilir ve eğer her bir değişken normal dağılım gösteriyorsa çok değişkenli normal oldukları varsayılır.
Çok değişkenli normallik her bir değişkenin ve degiskenlerin bütün doğrusal kombinasyonlarının normal dağıldığı sayıltısıdır.
Normalliğin iki bileşeni vardır: Çarpıklık ve basıklık
Bir değişkene ait dağılım normal olduğunda, değişkenin çarpıklık ve basıklık değerleri sıfıra eşittir.
Kural olarak eğer değişkenin çarpıklık ve basıklık değerleri -1.0 ile +1.0 arasındaysa, değişkenin normale oldukça yakın olduğu söylenebilir.
Sayıltının karşılanması durumunda analizin artıkları (hataları-açıklayamadığımız kaynaklar) da normal dağılır.
Çok değişkenli normallik sayıltısı farklı çok değişkenli istatistikler için farklı ele alınır.
library(sur)
attach(screen)
skew(TIMEDRS) # çarpıklık değeri ve çarpıklığın standart hatasını hesaplar## [1] 3.234045skew(data_temiz$TIMEDRS) ## [1] 2.35295se.skew(TIMEDRS) # standart hata## [1] 0.1133494skew.ratio(TIMEDRS) ## [1] 28.53164# çarpıklığın standart hataya bölünmesi -3+3 arasında olmalı, 28 çıktı çok yüksek
# skew(TIMEDRS)/se.skew(TIMEDRS)# jarque.test fonksiyonu veri normal dağılımdan farklılaşmamaktır yokluk hipotezini test etmektedir.(yeni nesil normallik testi)
library(moments)
library(labelled)
jarque.test(remove_labels(TIMEDRS))##
## Jarque-Bera Normality Test
##
## data: remove_labels(TIMEDRS)
## JB = 4034.9, p-value < 2.2e-16
## alternative hypothesis: greaterjarque.test(remove_labels(ATTDRUG))##
## Jarque-Bera Normality Test
##
## data: remove_labels(ATTDRUG)
## JB = 5.0552, p-value = 0.07985
## alternative hypothesis: greaterjarque.test(remove_labels(ATTHOUSE)) # sapma var istatistiksel olarak anlamlı olduğu için ##
## Jarque-Bera Normality Test
##
## data: remove_labels(ATTHOUSE)
## JB = 61.092, p-value = 5.418e-14
## alternative hypothesis: greaterjarque.test(remove_labels(data_temiz$ATTHOUSE))##
## Jarque-Bera Normality Test
##
## data: remove_labels(data_temiz$ATTHOUSE)
## JB = 1.5383, p-value = 0.4634
## alternative hypothesis: greaterset.seed(0)
normal <- rnorm(200)
non_normal<- rexp(200, rate=3)
par(mfrow=c(1,2))
hist(normal, col='steelblue', main='Normal')
hist(non_normal, col='steelblue', main='Non-normal')
par(mfrow=c(1,2))
qqnorm(normal, main='Normal')
qqline(normal)
qqnorm(non_normal, main='Non-normal')
qqline(non_normal)
ggplot(data = screen, aes(sample = ATTHOUSE )) +
geom_qq()+
geom_qq_line( ) 
ÖĞRENME GÜNLÜĞÜ 2
Bugünkü dersimizde geçen haftanın konusunu önce tekrar ettik, ben ilk derse katılamadığım için iyi oldu. Bu haftanın konusu ise, uç değer ve normallik kontrolü oldu.Normallik varsayımını kontrol ederken ezbere bildiğimiz Kolmogorov-Smirnov testinin, güncel literatürde (Howell, 2012) artık yerinin olmadığını öğrenmek ilginçti. Normallik kontrolü için, Q-Q plot, histogramlar ve çarpıklık/basıklık oranlarını incelemek ve grafiklerle bunun kontrolünü sağlamak ve çok değişkenli uç değerleri belirlemek için Mahalanobis uzaklığını hesaplamak aslında bildiğim kısımdı. Fakat jarque.test ile tek değişkenli normallik kontrolü yapmak yeni öğrendiğim bir bilgi oldu. Özetle bu ders, analiz öncesi veri temizliği ve varsayım kontrollerinin ne kadar dikkat gerektiren bir süreç olduğunu bana uygulamalı olarak gösterdi. Analiz öncesi veriyi adeta didik didik ettik..