Birkaç güzel, az ya da çok basit seçenekleriniz var. Üniformanız önceden onları basitleştirmeye yardımcı olur.
Seçenek 1: Bağımsızlık örnekleyici. Teklif dağıtımınızı, birim karesi üzerinde tekdüze bir dağılıma eşit olarak ayarlayabilirsiniz; bu, örneklerin siz çağırırken kısıtlı bölgenin dışında kalmamasını sağlar. Potansiyel olumsuz: Posterior birim karenin çok küçük bir bölgesinde yoğunlaşırsa, çok düşük bir kabul oranına sahip olabilirsiniz. OTOH, U (0,1) dağılımından daha hızlı rasgele sayılar üretmek zordur. Potansiyel ters: sizin için daha az iş.
2. Seçenek: Parametrelerinizi sınırlandırılmamış bir şeye dönüştürün, dönüştürülmüş parametreler için önerilerde bulunun, sonra parametreleri olasılık işlevlerinde kullanmak üzere geri dönüştürün. Bu durumda öncekinin dönüştürülmüş parametreler üzerinde olacağına dikkat edin, çünkü bunun için teklifte bulunuyorsunuz, bu yüzden yeni önceliği almak için dönüşümün Jacobian'ıyla uğraşmanız gerekecek. Analiziniz için, elbette, MCMC tarafından üretilen rasgele sayıları orijinal parametrelerine geri dönüştüreceksiniz. Potansiyel olumsuz: sizin için daha fazla başlangıç işi. Yukarı yönlü: teklifleriniz için daha iyi kabul oranı.
Seçenek 3: Birim karesinde bulunan bağımsızlık örnekleyicisinden başka bir teklif dağıtımı oluşturun. Bu, ünitenizi önceden korumanıza izin verir, ancak teklif olasılıklarını hesaplarken daha fazla karmaşıklık pahasına. Bunun bir örneği, parametrelerinizden birinin geçerli değeri olmasına izin vermek , parametrelerle ( n x , n ( 1 - x ) ) bir Beta dağılımı olacaktır . N ne kadar büyük olursa, teklifiniz o anki değerin etrafında o kadar yoğunlaşır. Potansiyel olumsuz: sizin için daha fazla başlangıç işi. Pptential Baş: teklifleriniz için daha iyi kabul oranı - ama yaparsanız nx(nx,n(1−x))nn çok büyükse ve köşeye yaklaşırsanız, çıkmadan önce köşede çok sayıda küçük hamle yapabilirsiniz.
4. Seçenek: Birim karenin dışında kalan tüm teklifleri reddetmeniz yeterlidir (Xian'ın yarı yürekli önerisi). Bunun yalnızca başka bir teklif oluşturmakla aynı şey olmadığını unutmayın; bu durumda teklifi reddedersiniz, yani parametre için bir sonraki değeriniz parametre için geçerli değerle aynıdır. Bu, parametre alanınızın bazı bölgeleri için önceden sıfır olasılıkınız olsaydı ve o bölgeye düşen rastgele bir sayı oluşturduysanız ne olacağıdır. Potansiyel olumsuz: Bir köşeye yaklaşırsanız, düşük bir kabul olasılığınız olabilir ve bir süre takılabilirsiniz. Potansiyel ters: sizin için daha az iş.
Seçenek 5: Birim karede, karşılaştığınız gerçek problemle aynı olan düzlemde genişletilmiş bir sorun oluşturun, her şeyi doğru yapın, ardından MCMC örneklemesinin sonuçlarını işlerken, dışarıdaki tüm örnekleri dışarı atın birim karenin. Potansiyel olumsuz: Bu genişletilmiş problemi oluşturmak çok kolaysa, sizin için daha az iş olabilir. Potansiyel dezavantaj: Markov zinciri bir süre birim karesinin dışında bir yerde dolaşırsa, örneklerinizin çoğunu atacağınız için aslında korkunç kabul olasılıklarına sahip olabilirsiniz.
Şüphesiz başka seçenekler de var, diğer insanların önerilerini görmek isterim!
2 ve 3 arasındaki fark, gerçekte yaptığınız şey için gerçek çıkarımlara rağmen, bir dereceye kadar kavramsaldır. Sadece öneri dağıtım parametresi bazı ayarlama kenara (I R programlama ediyorsam) ve ekstra çaba miktarı R önerisi olasılıkları ne söylemesine izin istesem de muhtemelen, 3 ile giderdim , bakışlar benim için küçük. JAGS veya HATA kullansaydım, elbette, bu tamamen farklı bir mesele olurdu, çünkü bu araçlar kendi önerilerini ele alıyor.n