আমি এখানে ডেটাতে একটি বয়েশিয়ান লগিট চালানোর চেষ্টা করছি । আমি ব্যবহার করছি 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()
। তবে আউটপুটটি এখনও বায়েশিয়ান স্পিরিটে থাকা উচিত, তাই না? এখানে মূল্য এবং ভ্যালুগুলি কী কী ? এই ঘনঘনবাদী অনুমান সরঞ্জাম না? তাদের কি এখানে আলাদা ব্যাখ্যা করা হচ্ছে?