আমি এমন একটি বিতরণ সন্ধান করছি যেখানে সম্ভাব্যতা ঘনত্বটি গড় থেকে কিছুটা দূরে পরে বা আমার নিজের কথায় একটি "মালভূমি-আকৃতির বিতরণ" পরে দ্রুত হ্রাস পায়।
গাউসিয়ান এবং ইউনিফর্মের মধ্যে কিছু।
আমি এমন একটি বিতরণ সন্ধান করছি যেখানে সম্ভাব্যতা ঘনত্বটি গড় থেকে কিছুটা দূরে পরে বা আমার নিজের কথায় একটি "মালভূমি-আকৃতির বিতরণ" পরে দ্রুত হ্রাস পায়।
গাউসিয়ান এবং ইউনিফর্মের মধ্যে কিছু।
উত্তর:
আপনি জেনারালাইজড নরমাল (সংস্করণ 1) , সাববোটিন বিতরণ , বা ক্ষতিকারক শক্তি বিতরণের নামে পরিচিত বিতরণটির সন্ধান করতে পারেন । এটা তোলে অবস্থান দ্বারা parametrized হয় , স্কেল σ এবং আকৃতি β পিডিএফ সঙ্গে
হিসাবে আপনি লক্ষ্য করতে পারেন জন্য এটা বর্ণনার অনুরূপ এবং এগোয় Laplace বিতরণের সঙ্গে β = 2 এটা স্বাভাবিক-র দিকে এগোয়, এবং যখন β = ∞ অভিন্ন বন্টন করা হয়।
আপনি যদি এমন সফ্টওয়্যার সন্ধান করেন যা এটি কার্যকর করেছে, আপনি normalp
আর (মাইনো এবং রুগিগেরি, 2005) এর জন্য লাইব্রেরিটি পরীক্ষা করতে পারেন । এই প্যাকেজটি সম্পর্কে যা সুন্দর তা হ'ল অন্যান্য বিষয়গুলির মধ্যে এটি সাধারণভাবে বিতরণ করা ত্রুটিগুলি, অর্থাৎ আদর্শকে হ্রাস করে রিগ্রেশন প্রয়োগ করে ।
মাইনো, এএম, এবং রুগিগেরি, এম (2005)। তাত্পর্যপূর্ণ শক্তি বিতরণের জন্য একটি সফ্টওয়্যার সরঞ্জাম: নরমালপ প্যাকেজ। পরিসংখ্যান সফটওয়্যার জার্নাল, 12 (4), 1-24।
@ স্ট্রংব্যাডের মন্তব্যটি সত্যই একটি ভাল পরামর্শ। আপনি যদি প্যারামিটারগুলি সঠিকভাবে বেছে নেন তবে একটি অভিন্ন আরভি এবং গাউসিয়ান আরভি এর যোগফল আপনাকে ঠিক যা খুঁজছেন তা দিতে পারে। এবং এটি আসলে একটি যুক্তিসঙ্গত সুন্দর বদ্ধ ফর্ম সমাধান আছে।
এই ভেরিয়েবলের পিডিএফ এক্সপ্রেশন দ্বারা দেওয়া হয়:
অসীম সংখ্যক "মালভূমি-আকারের" বিতরণ রয়েছে।
আপনি কি "গাউসিয়ান এবং ইউনিফর্মের মধ্যে" চেয়ে আরও নির্দিষ্ট কিছু করার পরে? এটা কিছুটা অস্পষ্ট।
এখানে একটি সহজ একটি: আপনি সর্বদা একটি ইউনিফর্মের প্রতিটি প্রান্তে একটি অর্ধ-স্বাভাবিক আটকে রাখতে পারেন:
আপনি সাধারণের স্কেলের তুলনায় ইউনিফর্মের "প্রস্থকে" নিয়ন্ত্রণ করতে পারেন যাতে আপনার আরও বিস্তৃত বা সংকীর্ণ মালভূমি থাকতে পারে, পুরো বন্টন বিতরণ করা যায়, যার মধ্যে গাউসিয়ান এবং ইউনিফর্ম সীমিত ক্ষেত্রে রয়েছে।
ঘনত্বটি হ'ল:
আমরা সম্ভবত এই ঘনত্বটিকে "গাউস-লেজযুক্ত ইউনিফর্ম" বলি।
আমার "শয়তানের টাওয়ার" বিতরণটি এখানে দেখুন [1]:
"স্লিপ-পোষাক" বিতরণ আরও আকর্ষণীয়।
আপনি যে আকারে চান তা নিয়ে বিতরণগুলি তৈরি করা সহজ।
[1]: ওয়েস্টফল, পিএইচ (2014)
"কুর্তোসিস হিসাবে পিকনেস, 1905 - 2014. আরআইপি"
আমি। তাত্ক্ষণিকবাজার। 68 (3): 191–195। doi: 10.1080 / 00031305.2014.917055
সর্বজনীন অ্যাক্সেস পিডিএফ: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4321753/pdf/nihms-599845.pdf
অন্য এক ( সম্পাদনা : আমি এখন এটা সরলীকৃত। EDIT2 : আমি এটা এমনকি আরও সরলীকৃত, যদিও এখন ছবি সত্যিই ঠিক এই সমীকরণ প্রতিফলিত না):
এখানে কিছু নমুনা কোড দেওয়া হল:
f = function(x, a, b, alpha){
y = log((cosh(2*alpha*pi*a)+cosh(2*alpha*pi*x))/(cosh(2*alpha*pi*b)+cosh(2*alpha*pi*x)))
y = y/pi/alpha/6
return(y)
}
f
আমাদের বিতরণ হয়। এর ক্রম জন্য এটি প্লট করা যাকx
plot(0, type = "n", xlim = c(-5,5), ylim = c(0,0.4))
x = seq(-100,100,length.out = 10001L)
for(i in 1:10){
y = f(x = x, a = 2, b = 1, alpha = seq(0.1,2, length.out = 10L)[i]); print(paste("integral =", round(sum(0.02*y), 3L)))
lines(x, y, type = "l", col = rainbow(10, alpha = 0.5)[i], lwd = 4)
}
legend("topright", paste("alpha =", round(seq(0.1,2, length.out = 10L), 3L)), col = rainbow(10), lwd = 4)
কনসোল আউটপুট:
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = 1"
#[1] "integral = NaN" #I suspect underflow, inspecting the plots don't show divergence at all
#[1] "integral = NaN"
#[1] "integral = NaN"
এবং চক্রান্ত:
আপনি যথাক্রমে opeালের শুরু এবং শেষের দিকে যথাক্রমে পরিবর্তন করতে পারেন a
এবং b
এরপরে আরও সাধারণীকরণ প্রয়োজন হবে, এবং আমি এটি গণনা করিনি (এজন্য আমি ব্যবহার করছি a = 2
এবং b = 1
চক্রান্তে)।
আপনি যদি একটি কেন্দ্রীয় মালভূমি এবং ত্রিভুজ বিতরণের পক্ষগুলির সাথে খুব সাধারণ কিছু সন্ধান করছেন তবে আপনি উদাহরণস্বরূপ এন ত্রিভুজ বিতরণ, এন মালভূমি এবং উত্থানের মধ্যে পছন্দসই অনুপাতের উপর নির্ভর করে একত্রিত করতে পারেন। ত্রিভুজ কেন, কারণ বেশিরভাগ ভাষায় তাদের নমুনা ফাংশন ইতিমধ্যে বিদ্যমান। আপনি এগুলির মধ্যে যেকোন একটি থেকে এলোমেলোভাবে সাজান।
আর এ যে দেবে:
library(triangle)
rplateau = function(n=1){
replicate(n, switch(sample(1:3, 1), rtriangle(1, 0, 2), rtriangle(1, 1, 3), rtriangle(1, 2, 4)))
}
hist(rplateau(1E5), breaks=200)
এখানে একটি চমত্কার এক: দুটি লজিস্টিক ফাংশন এর পণ্য।
(1/B) * 1/(1+exp(A*(x-B))) * 1/(1+exp(-A*(x+B)))
এতে টুকরোড়া না হওয়ার সুবিধা রয়েছে।
বি প্রস্থ সামঞ্জস্য করে এবং এ ড্রপ অফের খাড়াতা সামঞ্জস্য করে। নীচে এ = 2 সহ বি = 1: 6 দেখানো হয়েছে। দ্রষ্টব্য: এটি কীভাবে সঠিকভাবে করা যায় তা নির্ধারণ করার জন্য আমি সময় নিইনি।