সম্ভাবনা ঘনত্ব অন্তর অন্তর সন্ধান করুন

Question

আমার কাছে ভেক্টর আছে

x <- c(1,2,3,4,5,5,5,6,6,6,6,
       7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,
       7,7,7,7,7,7,7,7,8,8,8,8,9,9,9,10)

(আমার প্রকৃত ভেক্টরের দৈর্ঘ্য> 10,000) এবং আমি অন্তরগুলিকে সন্ধান করতে চাই যেখানে 90% ঘনত্ব থাকে। Is quantile(x, probs=c(0.05,0.95), type=5)সবচেয়ে উপযুক্ত বা অন্য কোন উপায় আছে কি?

Comments

আপনার প্রশ্নটি "বিরতি যেখানে ..." সম্পর্কে কিছুটা অস্পষ্ট - একাধিক বিরতি থাকতে পারে। আপনি কি পুরোপুরি অভ্যন্তরীণ 90%, প্রতিটি দিকে প্রতিসাম্যিকভাবে ছাঁটাই করতে আগ্রহী? সর্বোপরি, সর্বনিম্ন থেকে 90% আইল পর্যন্ত, 90% ডেটা ধরা হয়, একইভাবে 10% আইল সর্বাধিক মান পর্যন্ত value

— ইটরেটর

আপনি কি সংক্ষিপ্ত বিরতি, প্রতিসাম্য ব্যবধান (প্রতিটি প্রান্তের সমান সম্ভাবনা) বা অন্য কিছু খুঁজছেন?

— গ্লেন_বি -রিনস্টেট মনিকা

Answer 1

উপরে উল্লিখিত হিসাবে, একটি বিরতি সংজ্ঞায়িত করার বিভিন্ন উপায় রয়েছে যা 90% ঘনত্ব অন্তর্ভুক্ত করে। এখনও যেটি চিহ্নিত করা যায়নি তা হ'ল সর্বোচ্চ [উত্তরোত্তর] ঘনত্ব অন্তর ( উইকিপিডিয়া ), যা সংক্ষিপ্ত ব্যবধান হিসাবে সংজ্ঞায়িত হয় যার জন্য শেষ পয়েন্টগুলির অনুভূতিসংক্রান্ত ঘনত্বের ফাংশন মানগুলির মধ্যে পার্থক্য নামমাত্র সম্ভাবনা "।

library(coda)
HPDinterval(as.mcmc(x), prob=0.9)

Answer 2

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

অন্যান্য প্যাকেজগুলিতে ফাংশন রয়েছে যা বর্ণনামূলক পরিসংখ্যানের আরও সম্পূর্ণ সেট সরবরাহ করে। আমি Hmisc::describeযেটি ব্যবহার করি তা হ'ল তবে describeফাংশন সহ আরও কয়েকটি প্যাকেজ রয়েছে ।

Answer 3

আপনার উপায়টি বোধগম্য মনে হয়, বিশেষত উদাহরণে পৃথক পৃথক ডেটা সহ,

quantile(x,probs=c(0.05,0.95), type=5)
 5% 95% 
2.8 9.0

তবে অন্য উপায়টি হ'ল একটি গণিত ঘনত্ব কার্নেল ব্যবহার করা:

dx <- density(x)
dn <- cumsum(dx$y)/sum(dx$y)
li <- which(dn>=0.05)[1]
ui <- which(dn>=0.95)[1]
dx$x[c(li,ui)]
[1] 2.787912 9.163246

Answer 4

হ্যাঁ. :-)। আপনি stats::densityআরও কার্যকর হতে পারে আউটপুট পেতে পারেন ।