Markov zincirlerini depolama
\(m\), 3 günlük uyku yoksunluğundan sonra tepki süresindeki ortalama değişim olsun. Önceki bir egzersizde, \(m\)'nin ardıl modelinden yaklaşık 10.000 çekilişlik bir örnek elde ettin. Ortaya çıkan mcmc.list nesnesini sleep_sim olarak kaydettin ve bu nesne çalışma alanında yüklü:
sleep_sim <- coda.samples(model = sleep_jags, variable.names = c("m", "s"), n.iter = 10000)
Aslında, sleep_sim içindeki \(m\) değerlerinin örneklemi, dağılımı ardıla yakınsayan bağımlı bir Markov zinciridir. sleep_sim'in içeriğini inceleyecek ve analiz üzerinde daha ince kontrol sağlamak için içeriği bir veri çerçevesinde saklayacaksın.
Bu egzersiz, kursun bir parçasıdır
RJAGS ile Bayesçi Modelleme
Egzersiz talimatları
sleep_simliste nesnesininhead()çıktısına bak.İlk
sleep_simliste öğesi \(m\) ve \(s\) zincirlerini içerir. Bunlarısleep_chainsadlı bir veri çerçevesinde sakla. Her zincir elemanı için karşılık gelen yineleme numarasını kaydeden1:10000değerlerinden oluşan biriterdeğişkeni ekle.sleep_chains'in ilk 6 satırına bak.
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# Check out the head of sleep_sim
# Store the chains in a data frame
sleep_chains <- data.frame(___, iter = ___)
# Check out the head of sleep_chains