রৈখিক মিশ্র প্রভাবগুলির মডেলগুলির জন্য কীভাবে সংগ্রহ / গণনা লিভারেজ এবং কুকের দূরত্ব পাবেন


10

merক্লাস অবজেক্টের ( lme4প্যাকেজের মাধ্যমে প্রাপ্ত ) কীভাবে লিভারেজ এবং কুকের দূরত্ব গণনা (বা এক্সট্রাক্ট) করা যায় তা কি কেউ জানেন ? আমি একটি অবশিষ্টাংশ বিশ্লেষণের জন্য এগুলি পরিকল্পনা করতে চাই।

উত্তর:


13

আপনার আর প্যাকেজটি দেখতে হবে influence.ME। এটা আপনাকে মিশ্র প্রভাব দ্বারা উৎপন্ন মডেলের জন্য প্রভাবশালী তথ্য ব্যবস্থা গনা করতে পারবেন lme4

একটি উদাহরণ মডেল:

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

ফাংশনটি influenceপরবর্তী সমস্ত পদক্ষেপের ভিত্তি:

library(influence.ME)
infl <- influence(model, obs = TRUE)

কুকের দূরত্ব গণনা করুন:

cooks.distance(infl)

প্লট কুকের দূরত্ব:

plot(infl, which = "cook")

এখানে চিত্র বর্ণনা লিখুন


ধন্যবাদ! এটি অবশ্যই সহায়তা করে। কুকের দূরত্ব বনাম লিভারেজ প্লটের লিভারেজ গণনা সম্পর্কে কীভাবে?
রোয়ে অ্যাঞ্জেল

@ রইএঞ্জেল আমি অনুমান করি influence.MEপ্যাকেজটি দিয়ে এটি সম্ভব নয় । দুর্ভাগ্যক্রমে, এই কাজের জন্য আমার কাছে কোনও সমাধান নেই।
সোভেন হোহেনস্টেইন

এটি হওয়া উচিত নয় infl <- influence(model, group = "cyl"), কারণ আপনি এলোমেলো প্রভাব হিসাবে উল্লেখ করেছেন (1|cyl)? আমি জানি না, আমি এটি মোটেও বুঝতে পারি না, আমি কেবল প্রভাব ইনস্টল করেছি ... তবে কখন ব্যবহার করতে হবে obs = TRUEএবং কখন ব্যবহার করতে হবে তা সত্যই আমি জানি না group...
কৌতূহলী

আমি নিম্নলিখিতটি যুক্ত করতে চাই: আপনি যদি কুকের ডি দূরত্বের যে সারি নম্বরটি পেতে চান - প্লটটিতে প্লট না করে একই সংখ্যাটি ঘটতে থাকে, তবে আপনি একটি কাটা দিয়ে রান্নার ডি দূরত্বের সংখ্যাগুলি সম্পর্কে নিম্নলিখিত আর সূত্রটি ব্যবহার করতে পারেন উদাহরণস্বরূপ 0.1 এর মূল্য বন্ধcooksD_data<-as.data.frame(cooks.distance(ft1)) cooksD_data_select<-cooksd[cooksD_data>0.1,drop=FALSE,] cooksD_oultiers<-as.numeric(rownames(cooksD_data_select))]
ইলিয়াস এ্যাসটাস্টিক্স

এখানে প্রস্তাবিতhatvalues() ফাংশনটির চেয়ে এটি কি আরও ভাল ?
কৌতূহলী
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.