Çalışmanın veri seti, Tokat ilinde merkez ve
ilçelerinde bulunan kesimhanelerde Haziran 2018 ve
Mayıs 2024 arasında kesilen 9.399 adet dana karkasa
ait fiyatlardan (TL/kg) oluşmuştur
14. Çalışmada
zaman serisi analizinde kullanmak üzere her bir ayın
ortalama dana karkas fiyatı hesaplanarak 72 aylık veri
seti elde edilmiştir.
Verilerin analizinde, zaman serisini düzensiz
dalgalanmaların etkisinden arındırarak yalnızca uzun
vadeli hareketlerin etkisi altındaki seri değerini elde
etmek ve kestirim yapmak amacıyla trend analizi
gerçekleştirilmiştir. Bu çalışma kapsamında zaman serisi
analizinde Box-Jenkins modellerinden ARIMA modeli
kullanılmıştır. Verilerin tanımlayıcı istatistikleri yapılarak
fiyat tahmininde uyum iyiliği yüksek olan modelin
parametreleri tahmin edilerek fiyat kestirimi yapılmıştır.
Zaman serilerinin analizi için SPSS, version 25.0 15
programı kullanılmıştır.
Çalışmada, dana karkas fiyatlarının geleceğe
yönelik tahminleri için ARIMA (1,1,0) modeli
kullanılmıştır. Serilerin otokorelasyon fonksiyon grafiği
incelendiğinde, serinin durağan olmadığı tespit edilmiştir.
Seriyi durağan hale getirmek amacıyla birinci fark
alınmıştır. Dana karkas fiyatlarının tahminleri için uygun
ARIMA (p,d,q) modelini belirlemek üzere, seriye ait
otokorelasyon ve kısmi otokorelasyon fonksiyonları
incelenerek, p ve q seviyeleri belirlenmiş ve
parametrelerin anlamlılığı kontrol edilmiştir.
Modellerde yer alan 16
p: Otoregresif model derecesi,
q: Hareketli ortalama model derecesi,
d: Mevsimsel olmayan fark alma derecesidir.
ARIMA (p, d, q) modelinin ifadesi eşitlik [1]’de
verildiği gibi tanımlanabilir:
Burada:
øP Otoregresif operatör için parametre değerleri,
at: Hata terimi katsayıları, θq: Hareketli ortalama
operatörü için parametre değerleri, Ζt: Orjinal serinin d
dereceden farkı alınmış zaman serisi, Yani,
Wt=Yt - Yt-1 t=1,2,…,t [2]
olup birinci farklar serisi eşitlik [2]’de verildiği gibi
tanımlanır. Burada:
Wt=Birinci farklar serisi,
Yt= Orijinal zaman serisinin tesadüfi değişkenler
kümesidir.
Birinci farklar serisi de durağan değilse; birinci fark
serisinin tekrar farkı alınarak durağanlık kontrolü yapılır.
Bu da eşitlik [3]’de verildiği gibi modellenir.
Zt= Wt Wt-1
t=1,2,...,t
[3]
Fark alma derecesi d = 0 olduğunda (bu orjinal
serinin durağan olması anlamına gelir) ARIMA modeli
AR, MA ya da ARMA modeli haline gelecektir. Bu
özelliğinden dolayı ARIMA modellerinin Box-Jenkins
modellerinin tamamını bünyesinde barındırdığı
söylenebilir.
Box-Jenkins ARIMA modelleri kuruluş aşamasında
dört temel adımı içerir. İlk adımda genel olarak modelin
hangi sınıfta olduğu belirlenir. Genel modelin seçiminde
otokorelasyon ve kısmi otokorelasyon fonksiyonlarına ait
grafiklerden yararlanılır. Her iki fonksiyon göz önünde
bulundurularak, ARIMA modelleri için teorik
fonksiyonların özelliklerinden faydalanılır 17.
İkinci adımda ise verilere uyum sağlayan geçici bir
model kurulur. Bu modelin kurulmasında otokorelasyon ve kısmi otokorelasyon fonksiyonlarından faydalanılır. Modelin belirlenmesinde AR, MA, ARMA, ARIMA veya SARIMA gibi bir modele karar verilir 18.
Üçüncü adımda, geçici modelin parametreleri istatistiksel yöntemler kullanılarak kestirilir ve katsayıların standart hataları göz önünde bulundurularak anlamlılıkları belirlenir. Son aşamada ise karar verilen modelin kestirimi için uygunluğu belirlenir. Bunun için geçici kurulan modelin hatalarının otokorelasyon katsayılarının grafiğinden yararlanılarak otokorelasyon fonksiyonu incelenir. Eğer bu fonksiyon belirli bir yapı gösteriyorsa, hataların rastgele dağılmadığı sonucuna varılır, bu da geçici modelin uygun olmadığını gösterir. Bu durumda ikinci aşamaya dönülerek süreç, uygun bir model bulunana kadar yeni bir geçici modelle tekrarlanır. Bu süreçte uygunluk kontrolünü sağlayan model, artık tahmin veya kestirim yapmak için kullanılabilir 19-21.
Modelin tahmin tutarlılığının değerlendirilmesinde; Kestirim Hatası (Forecast Error), Kestirim Hata Ortalaması (Mean Forecast Error), Mutlak Kestirim Hata Ortalaması (Mean Absolute Deviation), Mutlak Kestirim Hata Ortalaması Yüzdesi (Mean Absolute Percentage Error) ve Hata Kareler Ortalaması (Mean Square Error) parametreleri kullanılmıştır 13.