কোনও বিশ্লেষণাত্মক ফর্ম রাখা যথেষ্ট সহজ হতে পারে যখন কোনও পোস্টারিওর ডিস্ট্রিবিউশন বের করার পদক্ষেপগুলি?


12

এটি কম্পিউটেশনাল সায়েন্সেও জিজ্ঞাসা করা হয়েছিল।

আমি 11 টি নমুনা সহ একটি স্বাবলম্বীকরণের জন্য কিছু সহগের একটি বয়েসিয়ান অনুমান গণনা করার চেষ্টা করছি: যেখানে mean গড় 0 এবং বৈকল্পিক সহ গাউসিয়ান the ভেক্টরের উপর পূর্ব বিতরণ গড় সহ গাউসিয়ান এবং এর সাথে একটি তির্যক কোভারিয়েন্স ম্যাট্রিক্স তির্যক টি এন্ট্রির মধ্যে সমান ।

Yi=μ+αYi1+ϵi
ϵiσe2(μ,α)t(0,0)σp2

স্বাবলম্বন সূত্রের ভিত্তিতে, এর অর্থ হ'ল ডেটা পয়েন্টগুলির বিতরণ ( ) গড় এবং বৈকল্পিক with সহ স্বাভাবিক । সুতরাং, ডেটা পয়েন্টগুলির সকলের ঘনত্ব যৌথভাবে (স্বতন্ত্রতা ধরে নিচ্ছি, যা আমি লিখছি সেই প্রোগ্রামটির জন্য ভাল),Yiμ+αYi1σe2(Y)

p(Y|(μ,α)t)=i=21112πσe2exp(YiμαYi1)22σe2.

বয়েসের উপপাদ্য দ্বারা, আমরা পূর্ববর্তী ঘনত্বের সাথে উপরের ঘনত্বের পণ্যটি নিতে পারি এবং তারপরে আমাদের কেবলমাত্র স্বাভাবিককরণের ধ্রুবক প্রয়োজন। আমার কুণ্ডলীটি হ'ল এটি গাউসীয় বিতরণ হিসাবে কাজ করা উচিত, তাই আমরা এবং সাথে ইন্টিগ্রালগুলি দিয়ে স্পষ্টভাবে গণনা করার চেয়ে শেষে ধীরে ধীরে স্বাভাবিককরণের বিষয়ে চিন্তা করতে পারি ।μα

এই অংশটি নিয়ে আমি সমস্যায় পড়ছি। আমি কীভাবে পূর্বের ঘনত্বের গুণন করতে পারি (যা মাল্টিভারেট হয়) এবং অবিচ্ছিন্ন ডেটা ঘনত্বের এই পণ্যটি? উত্তরোত্তরটি খাঁটি and এবং ঘনত্ব হওয়া দরকার তবে আপনি দেখতে পাচ্ছেন না যে আপনি কীভাবে এই জাতীয় পণ্যটি পেয়ে যাবেন।μα

কোনও পয়েন্টার সত্যই সহায়ক, এমনকি যদি আপনি কেবল আমাকে সঠিক দিকে নির্দেশ করেন এবং তারপরে আমাকে অগোছালো বীজগণিত করা উচিত (যা আমি ইতিমধ্যে বেশ কয়েকবার চেষ্টা করেছি)।

শুরুর দিক হিসাবে, এখানে বয়েসের বিধি থেকে সংখ্যার ফর্মটি রয়েছে:

1(2πσe2)52πσp2exp[12σe2i=211(YiμαYi1)2μ22σp2α22σp2].

সমস্যা দেখেছি কিভাবে যে এই একটি গসিয়ান ঘনত্ব নিচে হ্রাস হয় ।(μ,α)t

যোগ করা হয়েছে

শেষ পর্যন্ত, এটি নিম্নলিখিত সাধারণ সমস্যার দিকে ফোটে। যদি আপনাকে কিছু চতুষ্কোণ প্রকাশ দেওয়া হয় যেমন কীভাবে এটিকে চতুর্ভুজ আকারে রেখেছেন some কিছু 2x2 ম্যাট্রিক্সের জন্য ? এটি সহজ ক্ষেত্রে যথেষ্ট সহজ, তবে গড় অনুমান, এবং পেতে আপনি কোন প্রক্রিয়াটি ব্যবহার করেন ?

Aμ2+Bμα+Cα2+Jμ+Kα+L
(μμ^,αα^)Q(μμ^,αα^)tQμ^α^

দ্রষ্টব্য, আমি ম্যাট্রিক্স সূত্রটি প্রসারণ এবং তারপরে উপরের মতো সহগের সমান করার চেষ্টা করার সরল বিকল্পটি চেষ্টা করেছি। আমার ক্ষেত্রে সমস্যাটি হ'ল ধ্রুব শূন্য এবং তারপরে আমি দুটি অজানাতে তিনটি সমীকরণ পেয়েছি, সুতরাং এটি কেবল সহগের সাথে মিলিয়ে যাওয়ার জন্য নির্ধারিত হয় (এমনকি যদি আমি একটি প্রতিসম চৌম্বক রূপ ম্যাট্রিক্স ধরেও নিই)।L


আমার এই উত্তর [এই প্রশ্ন] ( stats.stackexchange.com/Qestions/22852/… ) সহায়ক হতে পারে। নোট করুন যে আপনার প্রথম পর্যবেক্ষণের জন্য আপনার পূর্বের প্রয়োজন - পুনরাবৃত্তিগুলি সেখানে থামে।
সম্ভাব্যতা ব্লগ

এই ক্ষেত্রে আমার কেন এটি প্রয়োজন তা আমি দেখছি না। আমার পর্যবেক্ষণের পরে শর্তসাপেক্ষে স্বাধীন হওয়ার মতো সময় ব্যবধানগুলিও আমার সাথে আচরণ করার কথা। লক্ষ্য করুন যে যৌথ ঘনত্বের পণ্যটি কেবলমাত্র । আমি মনে করি না যে এখানে আমার অনুক্রমিকভাবে আপডেট হওয়া সূত্রটি পাওয়া যাবে, কেবল উত্তরকালের একটিমাত্র সূত্র । i=2..11p((μ,α)t|Y)
এলী

পূর্ববর্তী এর "মাল্টিভারিয়েট" ডেটা ঘনত্বের মধ্যে "ইউনিভারিটি" এর সাথে বিরোধী নয়, কারণ তারা এর ঘনত্ব । p(α,μ)yi
শি'য়ান

উত্তর:


7

পূর্বের উত্তরের আমার উত্তরে যে ক্লুটি ছিল তা হ'ল আমি কীভাবে প্যারামিটারগুলিকে একীভূত করেছি - কারণ আপনি এখানে ঠিক একই অবিচ্ছেদ্য কাজ করবেন। আপনি প্রশ্নটি ভেরিয়েন্স প্যারামিটারগুলি পরিচিত বলে ধরে নিয়েছেন, তাই তারা ধ্রুবক। আপনাকে কেবলমাত্র উপর নির্ভরতা দেখতে হবে। এটি দেখতে, নোট করুন যে আমরা লিখতে পারি:α,μ

p(μ,α|Y)=p(μ,α)p(Y|μ,α)p(μ,α)p(Y|μ,α)dμdα
=1(2πσe2)52πσp2exp[12σe2i=211(YiμαYi1)2μ22σp2α22σp2]1(2πσe2)52πσp2exp[12σe2i=211(YiμαYi1)2μ22σp2α22σp2]dμdα

লক্ষ্য করুন কিভাবে আমরা প্রথম ফ্যাক্টর টান করতে পারেন আউট ডিনোমিনেটরে ডাবল অবিচ্ছেদ্য এবং এটি সংখ্যার সাথে বাতিল হয়। আমরা স্কোয়ারের বের করতে পারি এবং এটি বাতিল হয়ে যাবে। আমরা যে অবিচ্ছেদ্যের সাথে বাকী রয়েছি তা এখন (স্কোয়ার শব্দটি প্রসারিত করার পরে):1(2πσe2)52πσp2exp[12σe2i=211Yi2]

=exp[10μ2+α2i=110Yi22μi=211Yi2αi=211YiYi1+2μαi=110Yi2σe2μ22σp2α22σp2]exp[10μ2+α2i=110Yi22μi=211Yi2αi=211YiYi1+2μαi=110Yi2σe2μ22σp2α22σp2]dμdα

এখন আমরা সাধারণ পিডিএফ থেকে একটি সাধারণ ফলাফল ব্যবহার করতে পারি।

exp(az2+bzc)dz=πaexp(b24ac)
এই বর্গ সম্পন্ন অনুসরণ এবং তিনি লক্ষ করেন উপর নির্ভর করে না । নোট করুন যে এর অভ্যন্তরীণ অবিচ্ছেদ্য এই ফর্মটির এবং এবং। এই অবিচ্ছেদ্য কাজ করার পরে, আপনি দেখতে পাবেন যে বাকী অবিচ্ছেদ্য overaz2+bzczμa=102σe2+12σp2b=i=211Yiαi=110Yiσe2c=α2i=110Yi22αi=211YiYi12σe2+α22σp2αএই ফর্মটিরও তাই, আপনি আবার এই সূত্রটি ব্যবহার করতে পারেন, একটি আলাদা । তারপরে আপনি pos form আকারে আপনার পোস্টারিয়র লিখতে সক্ষম হবেন যেখানে একটি ম্যাট্রিক্সa,b,c12π|V|12exp[12(μμ^,αα^)V1(μμ^,αα^)T]V2×2

আপনার আরও ক্লু প্রয়োজন হলে আমাকে জানান।

হালনাগাদ

(দ্রষ্টব্য: সঠিক সূত্র, হওয়া উচিত এর পরিবর্তে )10μ2μ2

আপনি যদি আপডেটটিতে লিখেছেন এমন চতুষ্কোণ ফর্মটি পর্যালোচনা করে দেখি তবে আমরা দেখতে পাচ্ছি যে এখানে সহগ রয়েছে ( পোস্টেরিয়রের ক্ষেত্রে অপ্রাসঙ্গিক কারণ আমরা সর্বদা কোনও ধ্রুবক যুক্ত করতে পারি যা ডিনোমিনেটরে বাতিল হয়ে যায়)। আমরা আছে অজানা । সমীকরণ রৈখিকভাবে স্বতন্ত্র হওয়া পর্যন্ত এটি একটি "ভালভাবে উত্থাপিত" সমস্যা। আমরা যদি চতুর্ভুজটিকে প্রসারিত করি আমরা পেয়েছি:5L5μ^,α^,Q11,Q12=Q21,Q22(μμ^,αα^)Q(μμ^,αα^)t

Q11(μμ^)2+Q22(αα^)2+2Q12(μμ^)(αα^)
=Q11μ2+2Q21μα+Q22α2(2Q11μ^+2Q12α^)μ(2Q22α^+2Q12μ^)α+
+Q11μ^2+Q22α^2+2Q12μ^α^

দ্বিতীয় ক্রমের সহগের তুলনা করে আমরা যা আমাদের জানায় যে (বিপরীত) কোভেরিয়েন্স ম্যাট্রিক্স দেখতে কেমন looks এছাড়াও আমাদের কাছে পরিবর্তে for এর জন্য আরও দুটি জটিল সমীকরণ রয়েছে । এগুলি ম্যাট্রিক্স আকারে লিখিত হতে পারে:A=Q11,B=2Q12,C=Q22α^,μ^Q

(2ABB2C)(μ^α^)=(JK)

এইভাবে অনুমানগুলি দেওয়া হয়:

(μ^α^)=(2ABB2C)1(JK)=14ACB2(BK2JCBJ2KA)

দেখানো হচ্ছে যে অনন্য অনুমান নেই । এখন আমাদের রয়েছে: 4ACB2

A=102σe2+12σp2B=i=110Yiσe2C=i=110Yi22σe2+12σp2J=i=211Yiσe2K=i=211YiYi1σe2

মনে রাখবেন যে আমরা যদি জন্য সংজ্ঞায়িত এবং সীমা তবে জন্য অনুমানগুলি সর্বনিম্ন স্কোয়ার দ্বারা দেওয়া হয় অনুমান এবং যেখানে এবং । সুতরাং পূর্ববর্তী অনুমানগুলি ওএলএস অনুমান এবং পূর্ববর্তী প্রাক্কলন মধ্যে ওজনযুক্ত গড় ।Xi=Yi1i=2,,11σp2μ,αα^=i=211(YiY¯)(XiX¯)i=211(XiX¯)2μ^=Y¯α^X¯Y¯=110i=211YiX¯=110i=211Xi=110i=110Yi(0,0)


এটি বিশেষভাবে সহায়ক নয় কারণ আমি সুনির্দিষ্টভাবে উল্লেখ করেছি যে এটি এখানে ডোনমিনেটর নয়। ডিনোমিনিটরটি কেবলমাত্র একটি স্বাভাবিক ধ্রুবক, যা আপনি গাওসিয়ান ফর্মে অঙ্কগুলি হ্রাস করার পরে তা স্পষ্ট হবে। ডিনোমিনেটরে সংবিধাগুলি মূল্যায়নের জন্য কৌশলগুলি গাণিতিকভাবে সত্যিই দুর্দান্ত, তবে আমার আবেদনের প্রয়োজন নেই just আমার কেবলমাত্র সমাধানটির সমাধানের দরকার যা হ'ল সংখ্যার হেরফের করা।
এলী

এই উত্তরটি আপনাকে অংকের এবং ডিনোমিনেটর উভয়ই দেয়। সংখ্যক মধ্যে যথাযথ দ্বিতীয় ডিগ্রি বহুপদী প্রদর্শন করে যা সম্ভাব্যতাবিহীন চাপ দিয়ে সাধারণ চতুর্ভুজ আকারে নিয়ে যায়। (α,μ)
শি'য়ান

@ সিম - সাধারণকরণের ধ্রুবক গণনা করে আপনি প্রয়োজনীয় চতুর্ভুজ ফর্মটি তৈরি করবেন। এতে বর্গটি পূর্ণ করার জন্য প্রয়োজনীয় শর্তাদি থাকবে
সম্ভাব্যতাব্লোগিক

এটি কীভাবে আপনাকে চতুর্ভুজ রূপ দেয় তা আমি বুঝতে পারি না। আপনার পোস্ট করা গাউসিয়ান অবিচ্ছেদ্য পরিচয় ব্যবহার করে ডিনোমিনেটরে দুটি অখণ্ডের কাজ করেছি। শেষ পর্যন্ত, আমি কেবল একটি বিশাল, অগোছালো ধ্রুবক পাই। এই ধ্রুবকটিকে গ্রহণ করার এবং এটিকে 1/2 পাওয়ার ইত্যাদির জন্য নির্ধারকরূপে পরিণত করার কোনও সুস্পষ্ট উপায় বলে মনে হচ্ছে না, উল্লেখ করার জন্য নয় আমি এর কোনটি কীভাবে নতুন গণনা করা যায় তা ব্যাখ্যা করে না গড় ভেক্টর ' .. এটিই মূল প্রশ্নে আমি সাহায্যের জন্য বলছিলাম। (μ^,α^)t
এলী

বিস্তারিত সংযোজনের জন্য ধন্যবাদ। চতুর্ভুজ রূপটি বের করার জন্য বীজগণিত করার চেষ্টা করার সময় আমি কিছু নির্বোধ ত্রুটি করছিলাম। ওএলএস অনুমানের সাথে সম্পর্কিত সম্পর্কে আপনার মন্তব্যগুলি অত্যন্ত আকর্ষণীয় এবং প্রশংসিতও। আমি মনে করি এটি আমার কোডটিকে ত্বরান্বিত করবে কারণ আমি একটি বিশ্লেষণাত্মক ফর্ম থেকে নমুনা আঁকতে সক্ষম হব যা বিল্ট-ইন, অনুকূলিত পদ্ধতিগুলি রয়েছে। আমার আসল পরিকল্পনাটি ছিল এ থেকে নমুনার জন্য মেট্রোপলিস-হেস্টিংস ব্যবহার করা, তবে এটি খুব ধীর ছিল। ধন্যবাদ!
এলী
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.