বিসিএ পদ্ধতি ব্যবহার করে আত্মবিশ্বাসের ব্যবধান গণনা করার সময় আর বুট প্যাকেজ থেকে ত্রুটিটি "আনুমানিক সমন্বয়" এ "এনএ হয় কেন?

Question

আমার কাছে সংখ্যার ভেক্টর রয়েছে যা আমি এখানে (... / কোড / মাইডাটা.আরডিটা) ডুপুট ব্যবহার করে আপলোড করেছি।

আমি বিসিএ সিআই পেতে চাই তাই আমি এই কোডটি লিখেছি:

my.mean <- function(dat, idx){
 return (mean(dat[idx], na.rm = TRUE))
} 

boot.out<-boot(data=my.data, statistic = my.mean, R=1000)

তবে আমি যখন নিম্নলিখিতগুলি চালনা করি আমি তা পাই:

> boot.ci(boot.out)
Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o,  : 
estimated adjustment 'a' is NA
In addition: Warning message:
In boot.ci(boot.out) : bootstrap variances needed for studentized intervals

এই ত্রুটিটি কী কারণে ঘটছে তা বুঝতে কেউ আমাকে সহায়তা করতে পারে? আপনার সাহায্যের জন্য ধন্যবাদ!

Comments

আমি আর = 1000 এর পরিবর্তে আর = 10000 দিয়ে চেষ্টা করেছি (এবং টাইপ = 'বিসিএ' এর পরিবর্তে = 'সমস্ত' টাইপ করে যা আপনি যা চান তা তাই) এবং এটি কার্যকর হয়) তাই আমার উপসংহার কি এই ত্রুটি হয় যার ফলে যে হবে যথেষ্ট বুটস্ট্র্যাপ প্রতিলিপিকৃত নয় ... কিন্তু আমি আরো বলতে পারবে না ...

— ocram

Answer 1

যেমন আপনি নিজের ত্রুটি বার্তা, boot.ciকলগুলি দেখতে পাচ্ছেন bca.ci। যেহেতু boot.outবস্তুটি সরবরাহ করে না L, তথ্যের জন্য আপনি যে পরিসংখ্যানগুলি গণনা করছেন তার পরীক্ষামূলক প্রভাবের মানগুলি ফাংশনটি bca.ciব্যবহার করে সেগুলি গণনা করার চেষ্টা করে empinfএবং তারপরে (মাইকেল যেমন বলেছেন) এটি ত্বরণ ধ্রুবক গণনা করতে তাদের ব্যবহার করে:

L <- empinf(boot.out, index = index, t = t.o, ...)
a <- sum(L^3)/(6 * sum(L^2)^1.5)

তবে অল্প সংখ্যক প্রতিলিপি সহ, empinfকখনও কখনও ব্যর্থ হয় এবং NAমানগুলির একটি ভেক্টর প্রদান করে । ফলাফলটি হ'ল আপনার কোনও মান নেই L, aগণনা করা যায় না এবং আপনি নিজের ত্রুটি পান। যেমন ওক্রাম বলেছেন, বুস্ট্র্যাপের প্রতিরূপের সংখ্যা বাড়ানো এটি ঠিক করবে। এমনকি R2000 এ দ্বিগুণ করারও সম্ভবত এটি করা উচিত।

Comments

+1 @ কিরান ধন্যবাদ! আমি অনুমান করি যে আমার অনুমানের কিছুটা বৈধতা ছিল। জড়িত আর প্রোগ্রামগুলি সম্পর্কে আপনার জ্ঞান থেকে আপনার সঠিক উত্তর রয়েছে। সুতরাং আমার উত্তর আর প্রয়োজন নেই। কেউ হয় আমাকে অনুমান করা পছন্দ করেনি বা ভেবেছিলেন আমার উত্তরটি একটি মন্তব্য হওয়া উচিত। এই সমস্ত কারণে এটি এখন মুছে ফেলা আমার পক্ষে বোধগম্য।

— মাইকেল আর চেরনিক

যারা সমস্যায় ভুগছেন তাদের ক্ষেত্রে: stat.ethz.ch/pipermail/r-help/2011-Feb February/269006.html অনুসারে পুনরাবৃত্তির সংখ্যা অবশ্যই ডাটা সারিগুলির সংখ্যার চেয়ে বড় হতে হবে। আমি হালকাভাবে ডুবেছি যে এতগুলি পুনরাবৃত্তি প্রয়োজনীয় ...

— বার্নিনলিয়ো