প্রতিস্থাপন ছাড়াই স্বতন্ত্র এলোমেলো নমুনার ছেদগুলির কার্ডিনালটির বিতরণ কী?


10

S উপাদানগুলিতে with সহ কিছু সেট রয়েছে এবং চেয়ে কম বা সমান হিসাবে ধনাত্মক পূর্ণসংখ্যা স্থির হয় ।nNa1,a2,...,amn

এর উপাদানগুলির সমান সম্ভাবনা থাকার কারণে, নমুনাগুলি আলাদাভাবে এবং আলাদাভাবে থেকে প্রতিস্থাপন ছাড়াই আঁকা হয় , যার আকার যথাক্রমে হয়।SmL1,L2,...,LmSa1,a2,...,am

নমুনাগুলিগুলির কার্ডিনালিটিটি, সাধারণভাবে, সমান সমর্থন করে তবে এটি কোন বিতরণ অনুসরণ করে?{0,1,,সর্বনিম্ন{ একটি 1 , একটি 2 ,, একটি মি }}|L1L2 ... Lm|{0,1,...,min{a1,a2,...,am}}


আমি এটি পুনরাবৃত্তভাবে গণনার জন্য একটি রেসিপি সরবরাহ করতে পারি তবে আমি একটি বন্ধ ফর্ম সমাধান সম্পর্কে অবগত নই। এটি কি যথেষ্ট হবে, বা আপনি কি a1,,am এবং n দেওয়া বিতরণ ফাংশনটির একটি স্পষ্ট অভিব্যক্তি চান n?
Bridgeburners

@ ব্রিজবার্নার্স একটি রেসিপিটি দুর্দান্ত হবে, কমপক্ষে এটি এই সমস্যাটি সম্পর্কিত এবং সম্পর্কিত সম্পর্কিত কিছু পদ্ধতি / উপায় সরবরাহ করবে।
llrs

উত্তর:


3

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

আমাদের কাছে

P(|L1L2Lm|=k)=(nk)i=1n(nai)j=0min(a1,,am)k(1)j(nkj)l=1n(njkaljk).

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

j=0min(a1,,am)k(1)j(nj,k,njk)l=1nP(Hyp(n,j+k,al)=j+k).
Hyp(n,j+k,al)alnj+k

শিক্ষাদীক্ষা

সম্মিলিত যুক্তিগুলি ট্র্যাক করার জন্য কিছুটা সহজ করার জন্য আসুন আমরা কিছু স্বীকৃতি পাই (আশা করি)। সামগ্রিকভাবে, আমরা এবং । স্থির বিবেচনা করি। অর্ডার করা টিপলসগুলি L_ সংগ্রহটি বোঝাতে আমরা ব্যবহার করব , যেখানে প্রতিটি , সন্তুষ্ট1 , , এম সি ( আই ) মি ( এল 1 , , এল এম ) এল আইএসSa1,,amC(I)m(L1,,Lm)LiS

  • |Li|=ai ; এবং
  • L1Lm=I

আমরা সমানতার পরিবর্তে প্রয়োজন বাদে আমরা অভিন্ন সংগ্রহের জন্য ব্যবহার । এল 1 এল এমআইC(I)L1LmI

একটি মূল পর্যবেক্ষণ হ'ল গণনা করা তুলনামূলকভাবে সহজ। এর কারণ হল অবস্থা সমতূল্য সবার জন্য , তাই একটা ধারনা বিভিন্ন মধ্যে এই অপসারণ কথাবার্তাও মান। প্রতিটি , প্রয়োজনীয়তা পূরণের of সংখ্যা হ'ল since, যেহেতু আমরা সাইজের এর উপসেট বেছে নিয়ে এই জাতীয় একটি করতেএবং তারপর সঙ্গে unioning । এটা যে অনুসরণ করে এল 1 এল এমআমি এল আইআমিআইআমিআই এল আমি ( | এস | - | আমি |C(I)L1LmILiIiiiLiএলআমিএসআমিএকটিআমি-| আমি| আমি| সি(আই)| =ni=1(|এস|-|আই|)(|S||I|ai|I|)LiSIai|I|I

|C(I)|=i=1n(|S||I|ai|I|).

এখন আমাদের মূল সম্ভাবনাটি follows follows এর মাধ্যমে প্রকাশ করা যেতে পারে : পি ( | এল 1এল 2এল এম | = কে ) = আই : | আমি | = কে | সি ( আই ) |C

P(|L1L2Lm|=k)=I:|I|=k|C(I)|all IS|C(I)|.

আমরা এখনই এখানে দুটি সরলকরণ করতে পারি। প্রথমত, ডিনোনেটরটি দ্বিতীয়ত, একটি অনুগতি যুক্তি দেখায় যেশুধুমাত্র উপর নির্ভর করে cardinality মাধ্যমে। যেহেতু সাবলেট রয়েছে যার কার্ডিনালিটি , এটি অনুসরণ করে যে যেখানে কার্ডিয়ালিটি থাকার একটি নির্বিচারে, স্থির উপসেট| সি(আই)| আমি| আমি| (এন)

|C()|=i=1n(|S|ai)=i=1n(nai).
|C(I)|I|I| এসকেআই:| আমি| =কে| সি(আই)| = ( এন)(nk)Skআমি0এসকে
I:|I|=k|C(I)|=(nk)|C(I0)|,
I0Sk

একটি পদক্ষেপ পিছনে নিয়ে যাওয়া, আমরা এখন সেই সমস্যাটি " showing দেখানোর ক্ষেত্রে কমিয়েছি

|C(I0)|=j=0min(a1,,am)k(1)j(nkj)l=1n(njkaljk).

যাক এর স্বতন্ত্র সাব-সেট নির্বাচন করা ঠিক এক উপাদান যোগ করে গঠিত । তারপরে (এই মাত্র বলছে যে যদি , তারপর রয়েছে তবে এতে কোনও অতিরিক্ত উপাদান থাকে না)) আমরা এখন ম্যাথক্যাল c-গণনা সমস্যাটিকে একটি ম্যাথক্যাল গণনা সমস্যায় রূপান্তর করেছি , যা আমরা কীভাবে পরিচালনা করতে হয় তা আরও জানি। আরও নির্দিষ্টভাবে, আমরা আছে এস আই 0 সি ( আই 0 ) = সি ( আই 0 ) ( এন - কে i = 1 সি ( জে আই ) )এল 1এল এম = আমি 0 এল 1এল এম আমি 0 সেJ1,,JnkSI0

C(I0)=C(I0)(i=1nkC(Ji)).
L1Lm=I0L1LmI0C| সি ( আই 0 ) | = | সি ( আমি 0 ) | - | n - কে i = 1 সি ( জে আই ) | = n l = 1 ( এন - কে)C
|C(I0)|=|C(I0)||i=1nkC(Ji)|=l=1n(nkalk)|i=1nkC(Ji)|.

আমরা উপরের ইউনিয়ন প্রকাশের আকার পরিচালনা করতে অন্তর্ভুক্তি-বর্জন প্রয়োগ করতে পারি। এখানে গুরুত্বপূর্ণ সম্পর্কটি হ'ল যে কোনও অমানবিক , কারণ এটি যদি এর সাথে of এর একটি সংখ্যা থাকে তবে এটিতে তাদের থাকে। আমরা আরও নোট করি যে সেটটি এর আকার। অতএব I{1,,nk}

iIC(Ji)=C(iIJi).
L1LmJiiIJi|I0|+|I|=k+|I|
|i=1nkC(Ji)|=I{1,,nk}(1)|I|1|iIC(Ji)|=j=1nkI:|I|=j(1)j1l=1n(njkaljk)=j=1nk(1)j1(nkj)l=1n(njkaljk).
We (আমরা এখানে মানগুলিকে সীমাবদ্ধ করতে পারি যেহেতু দ্বিপদী সহগের শূন্য হয় যতক্ষণ না সমস্ত , অর্থাৎ ।)jjalkljmin(a1,,am)k

অবশেষে সমীকরণের সমাপ্তি দিয়ে ভাবটি স্থির করে উপরে এবং যোগফলকে একীভূত করে আমরা claimed হিসাবে দাবি করা হয়েছে।|C(I0)|

|C(I0)|=j=0min(a1,,am)k(1)j(nkj)l=1n(njkaljk)

সমস্ত প্রচেষ্টা এবং সমাধানের জন্য +1, তবে এর বেশিরভাগটি (এবং অন্য উত্তর) বুঝতে আমার আমার গণিতগুলিকে পোলিশ করতে হবে। ধন্যবাদ
llrs

4

আমি এটিকে সমাধান করার কোনও বিশ্লেষণ পদ্ধতি সম্পর্কে অবগত নই, তবে ফলাফলটি গণনার জন্য এটি একটি পুনরাবৃত্ত উপায়।

জন্য আপনি বেছে করছি বাইরে উপাদান যার সামনে মনোনীত করা হয়েছে। আপনার দ্বিতীয় ড্রতে এর সাথে ছেদ করে এমন উপাদানগুলি এ choosing বেছে নেওয়ার সম্ভাবনা বিতরণ দ্বারা প্রদত্ত:m=2একটি2এন, একটি1সর্বনিম্ন{একটি1,একটি2}এল1

পি(|এন,একটি1,একটি2)=(একটি1)(এন-একটি1একটি2-)(এনএকটি2)

আমরা ফলাফলটি কল করতেআমরা করতে একই যুক্তি ব্যবহার করতে পারি যেখানে হল তিনটি নমুনার ছেদের কার্ডিনালিয়ালিটি। তারপর,2পি(3=|এন,2,একটি3),3

পি(3=)=Σ=0সর্বনিম্ন(একটি1,একটি2)পি(3=|এন,2=,একটি3)পি(2=|এন,একটি1,একটি2)

প্রত্যেকের জন্য এই খুঁজে । পরবর্তী গণনাটি সংখ্যাগতভাবে কঠিন নয়, কারণ কেবল পূর্ববর্তী গণনার ফলাফল এবং is হাইপারজমেট্রিক বিতরণ।{0,1,2,...,সর্বনিম্ন(একটি1,একটি2,একটি3)}পি(2=|এন,একটি1,একটি2)পি(3=|এন,2=,একটি3)

সাধারণভাবে, আপনি নিম্নলিখিত পুনরাবৃত্ত সূত্রগুলি প্রয়োগ করতে পারেন: for for এবং যা কেবলমাত্রপি(মি)

পি(আমি=)=Σ=0সর্বনিম্ন(একটি1,একটি2,...,একটিআমি-1)পি(আমি=|এন,আমি-1=,একটিআমি)পি(আমি-1=),
i{2,3,,m},P(b1)=δa11,
পি(আমি=|এন,আমি-1=,একটিআমি)=()(এন-একটিআমি-)(এনএকটিআমি),
আমি{2,3,...,মি},
পি(1)=δএকটি11,
1=একটি1

এখানে এটি আর:

hypergeom <- function(k, n, K, N) choose(K, k) * choose(N-K, n-k) / choose(N, n)

#recursive function for getting P(b_i) given P(b_{i-1})
PNext <- function(n, PPrev, ai, upperBound) {
  l <- seq(0, upperBound, by=1)
  newUpperBound <- min(ai, upperBound)
  kVals <- seq(0, newUpperBound, by=1)
  PConditional <- lapply(kVals, function(k) {
    hypergeom(k, ai, l, n)
  })
  PMarginal <- unlist(lapply(PConditional, function(p) sum(p * PPrev) ))
  PMarginal
}

#loop for solving P(b_m)
P <- function(n, A, m) {
  P1 <- c(rep(0, A[1]), 1)
  if (m==1) {
    return(P1)
  } else {
    upperBound <- A[1]
    P <- P1
    for (i in 2:m) {
      P <- PNext(n, P, A[i], upperBound)
      upperBound <- min(A[i], upperBound)
    }
    return(P)
  }
}

#Example
n <- 10
m <- 5
A <- sample(4:8, m, replace=TRUE)
#[1] 6 8 8 8 5

round(P(n, A, m), 4)
#[1] 0.1106 0.3865 0.3716 0.1191 0.0119 0.0003
#These are the probabilities ordered from 0 to 5, which is the minimum of A

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