খয়রাত:
সম্পূর্ণ অনুগ্রহ এমন কাউকে দেওয়া হবে যিনি নীচে অনুমানকারী uses ব্যবহার করে বা উল্লেখ করেছেন এমন কোনও প্রকাশিত কাগজের একটি রেফারেন্স সরবরাহ করে।
প্রেরণা:
এই বিভাগটি সম্ভবত আপনার পক্ষে গুরুত্বপূর্ণ নয় এবং আমি সন্দেহ করি এটি অনুগ্রহ পেতে আপনাকে সহায়তা করবে না, তবে যেহেতু কেউ অনুপ্রেরণা সম্পর্কে জিজ্ঞাসা করেছেন, আমি এখানে কী কাজ করছি তা এখানে's
আমি একটি পরিসংখ্যান গ্রাফ তত্ত্ব সমস্যা নিয়ে কাজ করছি। স্ট্যান্ডার্ড ঘন গ্রাফ সীমিতকরণ বস্তু এই অর্থে একটি প্রতিসম ফাংশন । একটি গ্রাফ স্যাম্পলিং ছেদচিহ্ন স্যাম্পলিং যেমন ভাবা যেতে পারে অভিন্ন মূল্যবোধ ইউনিট ব্যবধান উপর ( জন্য ) একজন প্রান্ত সম্ভাবনা এবং তারপর হয় । ফলে অন্তিক ম্যাট্রিক্স বলা যেতে যাক ।
আমরা চিকিৎসা করতে পারে একটি ঘনত্ব যেমন ত যে \ iint ডব্লিউ> 0 । আমরা যদি অনুমান চ উপর ভিত্তি করে একটি কোনো সীমাবদ্ধতা ছাড়া চ , তাহলে আমরা একটি সামঞ্জস্যপূর্ণ অনুমান পেতে পারে না। ধারাবাহিকভাবে চ এর অনুমান সম্পর্কে একটি আকর্ষণীয় ফলাফল পেয়েছি যখন চ সম্ভাব্য ফাংশনের একটি সীমাবদ্ধ সেট থেকে আসে। এই অনুমানক এবং \ যোগফল থেকে , আমরা ডাব্লু অনুমান করতে পারি ।
দুর্ভাগ্যবশত, পদ্ধতি যে আমি যখন আমরা ঘনত্ব বন্টন থেকে নমুনা দেখায় দৃঢ়তা পাওয়া । পথ নির্মাণ করা হয় প্রয়োজন যে আমি পয়েন্ট একটি গ্রিড নমুনা (যেমন গ্রহণ উল্টোদিকে মূল থেকে স্বপক্ষে )। এই পরিসংখ্যান.এসই প্রশ্নে, আমি কী ঘটতে পারে তার 1 টি মাত্রিক (সরল) সমস্যার জন্য জিজ্ঞাসা করছি যখন আমরা কেবল সরাসরি বন্টন থেকে নমুনা না করে কেবল গ্রিডে নমুনা বার্নোলিসকেই নমুনা দিতে পারি।
গ্রাফ সীমা জন্য রেফারেন্স:
এল। লোভাস এবং বি। ঘন গ্রাফ সিকোয়েন্সের সীমাবদ্ধতা ( আরক্সিভ )।
সি। বর্গস, জে চয়েস, এল। লোভাস, ভি। সোস এবং কে। ভেস্তেরগোম্বি। ঘন গ্রাফের রূপান্তরীয় ক্রম i: উপগ্রাফের ফ্রিকোয়েন্সি, মেট্রিক বৈশিষ্ট্য এবং পরীক্ষার। ( আর্কসিভ )
স্বরলিপি:
সিডিএফ সঙ্গে একটি একটানা বন্টন বিবেচনা এবং PDF যা ব্যবধান ইতিবাচক সমর্থন আছে । ধরুন কোন pointmass আছে, , সর্বত্র differentiable এবং এছাড়াও এর supremum হয় ব্যবধান উপর । যাক গড় যে দৈব চলক বন্টন থেকে নমুনা হয় । আইআইডি ইউনিফর্ম র্যান্ডম ভেরিয়েবল ।
সমস্যা সেট আপ:
প্রায়শই, আমরা দেওয়া যাবে ডিস্ট্রিবিউশনের সাথে র্যান্ডম ভেরিয়েবল হতে চলিত সঙ্গে এবং কাজের গবেষণামূলক বণ্টনের ফাংশনের যেমন যেখানে সূচক ফাংশন। দ্রষ্টব্য যে এই অভিজ্ঞতামূলক বিতরণ নিজেই এলোমেলো (যেখানে স্থির থাকে)।
দুর্ভাগ্যক্রমে, আমি সরাসরি থেকে নমুনা আঁকতে পারছি না । তবে, আমি জানি যে কেবলমাত্র এর ইতিবাচক সমর্থন রয়েছে এবং আমি এলোমেলো ভেরিয়েবলগুলি তৈরি করতে যেখানে সাফল্যের সম্ভাবনা সহ বার্নোলি বন্টন সহ একটি র্যান্ডম ভেরিয়েবল হয় যেখানে উপরে এবং সংজ্ঞায়িত করা হয়েছে। সুতরাং, । এই মানগুলি থেকে আমি অনুমান করার একটি স্পষ্ট উপায় কোথায়
প্রশ্নাবলী:
সবচেয়ে কঠিন থেকে সহজ (যা আমি মনে করি) হওয়া সহজ।
কেউ কি জানেন যে এই (বা অনুরূপ কিছু) এর একটি নাম আছে? আপনি কি এমন কোনও রেফারেন্স দিতে পারেন যেখানে আমি এর কয়েকটি বৈশিষ্ট্য দেখতে পারি?
হিসাবে হয় একটি সামঞ্জস্যপূর্ণ মূল্নির্ধারক (এবং আপনি এটি প্রমাণ করতে পারেন)?
হিসাবে সীমাবদ্ধ বিতরণ কী ?
আদর্শভাবে, আমি নিম্নলিখিতগুলি - এর ফাংশন হিসাবে আবদ্ধ করতে চাই eg যেমন, , তবে আমি জানি না সত্য কী। ঘোরা সম্ভাব্যতা বিগ হে
কিছু ধারণা এবং নোট:
এটি অনেকটা গ্রিড ভিত্তিক স্তরের সাথে গ্রহণযোগ্যতা-প্রত্যাখ্যানের নমুনার মতো দেখায় । দ্রষ্টব্য যে এটি এর কারণ নয় কারণ সেখানে আমরা প্রস্তাবটি প্রত্যাখ্যান করলে আমরা আর একটি নমুনা আঁকিনা।
আমি নিশ্চিত যে এই পক্ষপাতদুষ্ট। আমি মনে করি বিকল্প নিরপেক্ষ, তবে এতে অপ্রীতিকর সম্পত্তি রয়েছে যা ।
আমি প্লাগ-ইন অনুমানকারী হিসাবে ব্যবহার করতে আগ্রহী । আমি এটি দরকারী তথ্য বলে মনে করি না, তবে এটি হতে পারে এমন কোনও কারণ সম্ভবত আপনি জানেন।
আর এর উদাহরণ
আপনি যদি সাথে অভিজ্ঞতা অভিজ্ঞতা ভাগ করতে চান তবে কিছু আর কোড । দুঃখিত কিছু ইনডেন্টেশন ভুল ... আমি কীভাবে এটি ঠিক করব তা দেখছি না।
# sample from a beta distribution with parameters a and b
a <- 4 # make this > 1 to get the mode right
b <- 1.1 # make this > 1 to get the mode right
qD <- function(x){qbeta(x, a, b)} # inverse
dD <- function(x){dbeta(x, a, b)} # density
pD <- function(x){pbeta(x, a, b)} # cdf
mD <- dbeta((a-1)/(a+b-2), a, b) # maximum value sup_z f(z)
# draw samples for the empirical distribution and \tilde{F}
draw <- function(n){ # n is the number of observations
u <- sort(runif(n))
x <- qD(u) # samples for empirical dist
z <- 0 # keep track of how many y_i == 1
# take bernoulli samples at the points s
s <- seq(0,1-1/n,length=n) + runif(n,0,1/n)
p <- dD(s) # density at s
while(z == 0){ # make sure we get at least one y_i == 1
y <- rbinom(rep(1,n), 1, p/mD) # y_i that we sampled
z <- sum(y)
}
result <- list(x=x, y=y, z=z)
return(result)
}
sim <- function(simdat, n, w){
# F hat -- empirical dist at w
fh <- mean(simdat$x < w)
# F tilde
ft <- sum(simdat$y[1:ceiling(n*w)])/simdat$z
# Uncomment this if we want an unbiased estimate.
# This can take on values > 1 which is undesirable for a cdf.
### ft <- sum(simdat$y[1:ceiling(n*w)]) * (mD / n)
return(c(fh, ft))
}
set.seed(1) # for reproducibility
n <- 50 # number observations
w <- 0.5555 # some value to test this at (called t above)
reps <- 1000 # look at this many values of Fhat(w) and Ftilde(w)
# simulate this data
samps <- replicate(reps, sim(draw(n), n, w))
# compare the true value to the empirical means
pD(w) # the truth
apply(samps, 1, mean) # sample mean of (Fhat(w), Ftilde(w))
apply(samps, 1, var) # sample variance of (Fhat(w), Ftilde(w))
apply((samps - pD(w))^2, 1, mean) # variance around truth
# now lets look at what a single realization might look like
dat <- draw(n)
plot(NA, xlim=0:1, ylim=0:1, xlab="t", ylab="empirical cdf",
main="comparing ECDF (red), Ftilde (blue), true CDF (black)")
s <- seq(0,1,length=1000)
lines(s, pD(s), lwd=3) # truth in black
abline(h=0:1)
lines(c(0,rep(dat$x,each=2),Inf),
rep(seq(0,1,length=n+1),each=2),
col="red")
lines(c(0,rep(which(dat$y==1)/n, each=2),1),
rep(seq(0,1,length=dat$z+1),each=2),
col="blue")
সম্পাদনাগুলি:
সম্পাদনা 1 -
আমি @ হুঁশিয়ার মন্তব্যগুলিকে সম্বোধন করতে এটি সম্পাদনা করেছি।
সম্পাদনা 2 -
আমি আর কোড যুক্ত করেছি এবং এটি আরও কিছুটা পরিষ্কার করেছি। আমি পঠনযোগ্যতার জন্য স্বরলিপিটি কিছুটা পরিবর্তন করেছি, তবে এটি মূলত একই। আমি অনুমতি পাওয়ার সাথে সাথে এটিতে একটি অনুগ্রহ দেওয়ার পরিকল্পনা করছি, সুতরাং আপনি আরও স্পষ্টতা চান কিনা দয়া করে আমাকে জানান।
সম্পাদনা 3 -
আমি মনে করি আমি @ কার্ডিনালের মন্তব্যগুলিকে সম্বোধন করেছি। আমি টাইপগুলি মোট ভ্যারিয়েশন স্থির করেছিলাম। আমি একটি অনুগ্রহ যোগ করছি।
সম্পাদনা 4 -
@ কার্ডিনাল এর জন্য একটি "অনুপ্রেরণা" বিভাগ যুক্ত করা হয়েছে।