বায়সিয়ান গ্ল্যামের


13

আমি এখানে ডেটাতে একটি বয়েশিয়ান লগিট চালানোর চেষ্টা করছি । আমি ব্যবহার করছি bayesglm()armআর দ্য কোডিং প্যাকেজ সহজবোধ্য যথেষ্ট:

df = read.csv("http://dl.dropbox.com/u/1791181/bayesglm.csv", header=T)
library(arm)
model = bayesglm(PASS ~ SEX + HIGH, family=binomial(link="logit"), data=df)

summary(model) নিম্নলিখিত আউটপুট দেয়:

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  0.10381    0.10240   1.014    0.311    
SEXMale      0.02408    0.09363   0.257    0.797    
HIGH        -0.27503    0.03562  -7.721 1.15e-14 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1   1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 2658.2  on 1999  degrees of freedom
Residual deviance: 2594.3  on 2000  degrees of freedom
AIC: 2600.3

এই মাধ্যমে আমাকে চলুন দয়া করে। আমি বুঝতে পারি যে এই কোডটি একটি খুব দুর্বল প্রাক ব্যবহার করে (যেহেতু আমি পূর্বের উপায়গুলি উল্লেখ করছি না) সুতরাং আউটপুটটি ব্যবহারিকভাবে একই হতে চলেছে যদি আমি glm()পরিবর্তে ব্যবহার করি তবে bayesglm()। তবে আউটপুটটি এখনও বায়েশিয়ান স্পিরিটে থাকা উচিত, তাই না? এখানে মূল্য এবং ভ্যালুগুলি কী কী ? এই ঘনঘনবাদী অনুমান সরঞ্জাম না? তাদের কি এখানে আলাদা ব্যাখ্যা করা হচ্ছে?pz


এটি একটি মন্তব্য এবং একটি উত্তর নয়, তবে এটি আমার কাছে কিছুটা অর্থ বোধ করবে। আপনি প্রাক্কলন বিতরণ সর্বাধিকতর যে মানগুলির জন্য সম্ভবত এটি অনুমান পান। এটিও সম্ভব হতে পারে যে তারা কেবল উত্তরোত্তর মাধ্যম? আপনি যদি পারেন তা পরীক্ষা করে নেওয়ার যোগ্য। তবে সঠিক বিশদটি কী তা নয়, একবার যদি আপনার কিছু অনুমান হয় তবে আপনি সেগুলি যথারীতি প্রাক্কলন / ધોરણের মাধ্যমে পরীক্ষা করতে পারেন। ত্রুটি -> জেড-স্কোর পদ্ধতি যা কাজ করে যদি পোস্টারিয়রটি কোনও স্বাভাবিকের কাছে পর্যাপ্ত থাকে (এটি কিছু পরিস্থিতিতে থাকে যা সাধারণত ধারণ করে তবে এটি একটি সাধারণ দিকে যায়)।
এরিক

এরিক ... আপনি সঠিক: সহগগুলি অবশ্যই উত্তরীয় ঘনত্বের মাধ্যম। আমার প্রশ্নটি p- এবং z-মান সম্পর্কে। তারা এখানে কি উপস্থাপন করে?
ব্যবহারকারী3671

ঠিক আছে. আপনার যদি এমন ঘনত্ব থাকে যা প্রায় সাধারণত বিতরণ করা হয় তবে আপনি জেড-স্কোর = গড় / বন্টনটির স্ট্যান্ডার্ড বিচ্যুতি গ্রহণের মাধ্যমে এটির মানটি সাধারণ বন্টনের সাথে তুলনা করে পরীক্ষা করতে পারবেন mean তারপরে আপনি দেখতে পাচ্ছেন যে মান বা তার চেয়ে বড় মান স্ট্যান্ডার্ড সাধারণ বিতরণ -> পি-মানের আওতায় আসবে। বিস্তারিত জানার জন্য উইকিপিডিয়ায় জেড-স্কোর দেখুন।
এরিক

হ্যাঁ ঠিক. তবে কেন এটি কোনও বায়েশিয়ান সেটিংয়ে বিরক্ত করবেন? বায়েশিয়ান অনুমান অনুসারে, বিন্দু অনুমানটি এলোমেলো পরামিতি সম্পর্কে আমার সেরা অনুমান, সুতরাং এটি পরীক্ষা করার দরকার নেই। সর্বাধিকত, আমি একটি "বিশ্বাসযোগ্য অন্তর্বর্তী" অন্তর্ভুক্ত করতে পারি যা একটি ঘন ঘনবাদী "আত্মবিশ্বাসের ব্যবধান" এর সমতুল্য তবে যার পরিসংখ্যানগত ব্যাখ্যাটি সম্পূর্ণ আলাদা different সংক্ষিপ্তসার () আউটপুট সম্পর্কে এটি বিভ্রান্তিকর অংশ। আত্মা বায়েশিয়ান, কিন্তু আউটপুট ঘন ঘন?
ব্যবহারকারী3671

একটি বিষয় হ'ল আপনি যে অনুমানটি পেয়েছেন তা অন্যরকম হবে, যেহেতু আপনি পূর্ব ব্যবহার করেছিলেন used আপনি যখন কোনও বায়সিয়ান উপায়ে কিছু দেখতে পেলেন যে কোনও কিছুর প্রভাব রয়েছে তা দেখানোর জন্য যদি পয়েন্টটি অনুমান করা হয় "সেরা অনুমান" তবে আপনি এটি দেখানোর চেষ্টা করবেন যে বিশ্বাসযোগ্য অন্তরটি শূন্য থাকে না। যখন আপনি একই গড় এবং এসডি (অ্যাসিপোটোটিকভাবে সঠিক) দিয়ে স্বাভাবিকের মধ্যবর্তী অংশটি আনুমানিক করেন তখন (1-পি / 2) বিশ্বাসযোগ্যতা ব্যবধানটি শূন্যের সমন্বিত বৃহত্তম প্রতিসম বিশ্বাসযোগ্যতা অন্তর, তাই আপনার উত্তরটি মূলত একই। P উপরে উল্লিখিত পি-মান।
এরিক

উত্তর:


16

দুর্দান্ত প্রশ্ন! যদিও বায়সিয়ান পি-মান রয়েছে এবং আর্ম প্যাকেজটির অন্যতম লেখক একজন অ্যাডভোকেট, আপনি নিজের আউটপুটে যা দেখছেন তা কোনও বায়সিয়ান পি-মান নয়। এর ক্লাস চেক করুনmodel

class(model)
"bayesglm" "glm"      "lm" 

এবং আপনি দেখতে পারেন যে ক্লাস বয়েসগ্ল্যামের উত্তরাধিকার সূত্রে প্রাপ্ত হয়েছে। তদ্ব্যতীত, আর্ম প্যাকেজ পরীক্ষা একটি বেইসগ্ল্যাম বস্তুর জন্য কোনও নির্দিষ্ট সংক্ষিপ্ত পদ্ধতি দেখায় না। সুতরাং যখন আপনি করবেন

summary(model)

আপনি আসলে করছেন

summary.glm(model)

এবং ফলাফলগুলির ঘনঘনবাদী ব্যাখ্যা পাওয়া। আপনি যদি আরও বায়েশীয় দৃষ্টিভঙ্গি চান তবে বাহুতে কাজটি হবেdisplay()


+1 দুর্দান্ত উত্তর! এটি আর এর সাথে সমস্যা, এমন অনেক উচ্চ বুদ্ধিমান পরিসংখ্যানবিদ আছেন যারা ভয়াবহ কোড লেখেন যা এই ধরণের ল্যান্ডমাইনগুলি প্রায় পড়ে থাকে।
বোগদানোভিস্ট

এটি কোনও তদারকি না করে ডিজাইনারদের অংশে ইচ্ছাকৃত পছন্দ বলে মনে হচ্ছে ।
এটিরেটু -

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

2
অনেক ধন্যবাদ, আতিরেতু। এটি আরও একটি প্রশ্ন উত্থাপন করে। প্রদর্শন () এবং সংক্ষিপ্তসার () এর মধ্যে পার্থক্য কী? আমার কাছে মনে হয় প্রাক্তন থেকে আউটপুটটি কেবলমাত্র পরবর্তী থেকে আউটপুট, কম দুটি কলাম, এবং 2 টি সংখ্যায় বৃত্তাকার। এটি উপরে প্রদর্শিত হবে, জেলম্যানের পোস্ট থেকে আপনি উপরে লিঙ্ক করেছেন।
user3671

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