স্ট্যানে সংজ্ঞায়িত প্রিয়ার ছাড়া প্যারামিটার


15

আমি কেবল স্ট্যান এবং ব্যবহার করতে শিখতে শুরু করেছি rstan। জেজিএস / বিজিজি কীভাবে কাজ করেছে সে সম্পর্কে আমি সবসময় বিভ্রান্ত না হয়ে আমি ভেবেছিলাম মডেলটির প্রতিটি প্যারামিটারটি আঁকানোর জন্য আপনাকে সর্বদা কোনও প্রকারের পূর্ব বিতরণটি সংজ্ঞায়িত করতে হবে। এটি স্ট্যান্ডে এটির ডকুমেন্টেশনের উপর ভিত্তি করে আপনাকে এটি করতে হবে না বলে মনে হচ্ছে। এখানে একটি নমুনা মডেল যে তারা দিতে এখানে

data {
  int<lower=0> J; // number of schools 
  real y[J]; // estimated treatment effects
  real<lower=0> sigma[J]; // s.e. of effect estimates 
} 
parameters {
  real theta[J]; 
  real mu; 
  real<lower=0> tau; 
} 
model {
  theta ~ normal(mu, tau); 
  y ~ normal(theta, sigma);
} 

উভয়ই muবা tauপ্রাইরির সংজ্ঞা দেওয়া হয়নি। আমার কিছু জেএজিএস মডেল স্ট্যানে রূপান্তর করার সময়, আমি খুঁজে পেয়েছি যে আমি যদি অপরিবর্তিত প্রিয়ারগুলির সাথে অনেকগুলি বা সর্বাধিক পরামিতি রেখে যাই তবে তারা কাজ করে।

সমস্যাটি হ'ল আমি যখন সংজ্ঞায়িত প্রিয়ার ছাড়া পরামিতিগুলি পাই তখন স্ট্যান কী করছে তা আমি বুঝতে পারি না। এটি কি অভিন্ন বিতরণের মতো কিছুতে ডিফল্ট হয়? এটি কি এইচএমসির বিশেষ বৈশিষ্ট্যগুলির মধ্যে এটি যে প্রতিটি প্যারামিটারের জন্য পূর্বনির্ধারিত প্রয়োজন হয় না?

উত্তর:


18

স্ট্যান রেফারেন্স ম্যানুয়াল থেকে (এর পূর্ববর্তী সংস্করণ) থেকে:

পূর্বের উল্লেখ না করা পূর্ববর্তী ইউনিফর্ম নির্দিষ্ট করার সমতুল্য।

প্যারামিটারটি সীমাবদ্ধ থাকলে একটি ইউনিফর্ম পূর্ববর্তী [...]

স্ট্যান প্রোগ্রামগুলিতে অনুচিত প্রিয়ারদেরও অনুমতি দেওয়া হয়; তারা নমুনা বিবৃতি ছাড়াই নিয়ন্ত্রণহীন পরামিতি থেকে উত্থিত হয়। কিছু ক্ষেত্রে, একটি অনুপযুক্ত পূর্বে সঠিক পোস্টারিয়র হতে পারে তবে প্যারামিটার (গুলি) বা উপাত্তগুলিতে সীমাবদ্ধতা বা উপাত্তের প্রাপ্যতা নিশ্চিত করার বিষয়টি গ্যারান্টি দিয়ে দেওয়া ব্যবহারকারীদের উপর নির্ভর করে।

( 1.0.1 সংস্করণে C3 বিভাগটিও দেখুন )।

স্ট্যানে এটি ঠিক আছে তবে BUGS এ নয় এর অন্তর্নিহিত কারণটি এই কারণের সাথে হতে পারে যে BUGS- এ আপনার মডেল "প্রোগ্রাম" একটি প্রথাগত গ্রাফিক্যাল মডেল নির্দিষ্ট করছে, যখন স্টানে আপনি যৌথ সম্ভাবনা গণনা করার জন্য একটি সামান্য ফাংশন লিখছেন ঘনত্ব ফাংশন। সমস্ত ভেরিয়েবলের জন্য একটি যথাযথ পূর্ব নির্ধারণ না করা গ্রাফিকাল মডেলগুলির দুর্দান্ত ফর্মাল বৈশিষ্ট্যগুলি স্ক্রু করতে পারে।

তবে হ্যামিলটোনিয়ান এমসির জন্য আপনাকে কেবলমাত্র (সংখ্যায়) যৌথ ঘনত্বের কার্য গণনা করতে হবে। পূর্বে একটি ফ্ল্যাট (এমনকি অনুপযুক্ত) কেবল ঘনত্বের জন্য একটি স্থির মেয়াদ অবদান রাখে, এবং যতক্ষণ না উত্তরাধিকার যথাযথ (সীমাবদ্ধ মোট সম্ভাবনার ভর) - এটি যে কোনও যুক্তিসঙ্গত সম্ভাবনা কার্যক্রমে থাকবে - এটি এইচএমসিতে সম্পূর্ণ উপেক্ষা করা যেতে পারে পরিকল্পনা.


8

থেকে স্ট্যান রেফারেন্স v1.0.2 (PG 6, পাদটীকা 1)

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

উভয় muএবং sigmaঅনুপযুক্ত ইউনিফর্ম প্রিয়ার আছে।

ফণা অধীনে, muএবং sigmaঅন্যরকম চিকিত্সা করা হয়। sigmaএকটি নিম্ন সীমা সঙ্গে সংজ্ঞায়িত করা হয়; স্ট্যান নমুনাগুলি log(sigma)(রূপান্তরটির জন্য জ্যাকবীয় সামঞ্জস্য সহ)। রূপান্তর সম্পর্কিত আরও তথ্যের জন্য, অধ্যায় 27 (পৃষ্ঠা 153) দেখুন।


ঠিক তাই আমি এ সম্পর্কে পরিষ্কার, যদি লগ (সিগমা) স্তরের স্ট্যান নমুনাগুলি, ফ্ল্যাট পূর্ববর্তীটি এখনও সিগমার উপরে এবং লগ (সিগমা) এর উপরে নয়?
রাসমুস বুথ

1
হ্যাঁ, ইউনিফর্ম আগে এখনও শেষ sigma, না log(sigma)। স্ট্যান প্যারামিটারটি রূপান্তর করে এবং জ্যাকবীয়িয়ান ব্যবহার করে ভেরিয়েবল সামঞ্জস্যের সঠিক পরিবর্তন প্রয়োগ করে এটি সম্পাদন করে।
21 '13 এ সিক্লিক
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.