যেহেতু আমি একজন সফটওয়্যার ইঞ্জিনিয়ার আরও পরিসংখ্যান শেখার চেষ্টা করছি আমি এমনকি শুরু করার আগে আমাকে ক্ষমা করতে হবে, এটি গুরুতর নতুন অঞ্চল ...
আমি পিএমসি শিখছি এবং কিছু সত্যিকারের (সত্যই) সাধারণ উদাহরণ দিয়ে কাজ করছি। আমি কাজ করতে না পারার একটি সমস্যা (এবং এর সাথে সম্পর্কিত কোনও উদাহরণ খুঁজে পাচ্ছে না) দুটি সাধারণ বিতরণ থেকে উত্পন্ন ডেটাতে একটি মডেল ফিট করা।
বলুন আমার 1000 মান রয়েছে; একটি থেকে উত্পন্ন 500 Normal(mean=100, stddev=20)
এবং এ থেকে আরও 500 জেনারেট করা Normal(mean=200, stddev=20)
।
যদি আমি তাদের কাছে কোনও মডেল ফিট করতে চাই, অর্থাৎ পিএমসি ব্যবহার করে দুটি উপায় এবং একক মানক বিচ্যুতি নির্ধারণ করুন। আমি জানি এটা কিছু ...
mean1 = Uniform('mean1', lower=0.0, upper=200.0)
mean2 = Uniform('mean2', lower=0.0, upper=200.0)
precision = Gamma('precision', alpha=0.1, beta=0.1)
data = read_data_from_file_or_whatever()
@deterministic(plot=False)
def mean(m1=mean1, m2=mean2):
# but what goes here?
process = Normal('process', mu=mean, tau=precision, value=data, observed=True)
অর্থাত্, উত্পাদনের প্রক্রিয়াটি সাধারণ, তবে দুটি মানগুলির মধ্যে একটি মি। আমি কেবল জানি না যে কোনও মান আসে m1
বা আসে কিনা এর মধ্যে "সিদ্ধান্ত" কে উপস্থাপন করতে হয় m2
।
সম্ভবত আমি পুরোপুরি এই মডেলিংয়ের জন্য ভুল পদ্ধতি গ্রহণ করছি? কেউ কি আমাকে একটি উদাহরণে নির্দেশ করতে পারেন? আমি বাগ এবং জেএজিএস পড়তে পারি তাই কিছু ঠিক আছে।