আর 4 ব্যবহার করে 4 টি মিশ্রিত এএনওওএ-এর মধ্যে আর-বিষয়গুলির সাথে আর ব্যবহার করে বি টি বিপরীতে কীভাবে ব্যবহার করবেন?

Question

এখানে আর এর প্রাথমিক ব্যবহারকারীর বারবার ব্যবস্থা আনোভা নিয়ে লড়াই করা AN

আমার একটি ডেটাসেট রয়েছে যা 4 টি মাত্রার ('গ্রুপ' নামে পরিচিত একটি একক ভেরিয়েবলের মধ্যে কোডেড) এবং চারটি স্তরযুক্ত বিষয় ফ্যাক্টরের মধ্যে একটি থাকে (চারটি পৃথক ভেরিয়েবল 'ডিভি 1', 'ডিভি 2', 'ডিভি 3 ',' ডিভি 4 ')।

আমার নিম্নলিখিত উদ্দেশ্যগুলি রয়েছে:

1. আনোভা সামগ্রিক পুনরাবৃত্তি ব্যবস্থা চালান।
2. কাস্টম বিপরীতে ব্যবহার করে গোষ্ঠীগুলির তুলনা করুন (এসপিএসএস-এ এলএম্যাট্রিক্স কমান্ড হিসাবে)।
3. কাস্টম বিপরীতে ব্যবহার করে বিভিন্ন স্তরের ডিভিয়ের তুলনা করুন (এসপিএস-তে একটি এমএম্যাট্রিক্স কমান্ড হিসাবে)।
4. 2) এবং 3) এর সংমিশ্রণটি একই সাথে করুন তাই আমি কেবলমাত্র কয়েকটি গোষ্ঠীর সাথে-বিষয় ফ্যাক্টরের নির্দিষ্ট স্তরে তুলনা করছি।
5. বিপরীতে একটি সেট চালান যা শূন্যের সমান হয় না।

আমি জানি যে আমি অনেক সমস্যা ছাড়াই এসপিএসে এটি করতে পারি, তবে আর কীভাবে এটি করা যায় সে সম্পর্কে আমি একটি স্পষ্ট ধারণা পেতে পারি না seen এখন পর্যন্ত দেখা গেছে কীভাবে এটি এক পদ্ধতিতে বা আর সম্পর্কিত সম্পর্কিত পদ্ধতির সেট-এর মধ্যে কাজ করতে পারে

Answer 1

একটি সমাধানের স্কেচ (অন্যটির জন্য নীচে দেখুন):

1. এসপিএসএসের (যেমন, প্রতি সারি প্রতি একটি বিষয়) বিস্তৃত ফর্ম্যাটের পরিবর্তে ডেটা দীর্ঘ বিন্যাসে (অর্থাত্ প্রতি সারি মূল্য অনুযায়ী) হওয়া দরকার, পুনরায় আকার প্যাকেজ দেখুন বা ?reshape। এর মধ্যে অন্তর্ভুক্ত রয়েছে যে সাবজেক্ট আইডেন্টিফায়ার (অর্থাত্ বিষয় আইডি) নির্দেশক একটি পরিবর্তনশীল হওয়া দরকার।
2. সমস্ত উপাদান (বিষয় শনাক্তকারী সহ) অবশ্যই শ্রেণীর ফ্যাক্টরের হতে হবে ( strএটি পরীক্ষা করতে আপনার ডেটা ফ্রেমে চালিত করুন)। আপনি যদি এটি না করেন তবে আপনার ফলাফলগুলি ভুল হবে।
3. আপনি যদি ধরণের III- র স্কোয়ারগুলি পেতে চান তবে কোডিংয়ের ক্ষেত্রে ডিফল্ট বিপরীতে সেট করুন:
   options(contrasts=c("contr.sum","contr.poly"))
4. সঙ্গে আকাঙ্ক্ষিত মডেল উল্লেখ lmeথেকে nlmeপ্যাকেজ (ইনস্টল এবং মাধ্যমে প্যাকেজ beforehands লোড install.packages("nlme")এবং library(nlme)) একটি compund প্রতিসম পারস্পরিক সম্পর্ক গঠন ব্যবহার করে। এই প্রশ্নের উত্তর গৃহীত হয়েছে এবং উত্তরটি বিশেষত আমার মন্তব্য দেখুন আপনার ক্ষেত্রে এটির মতো কিছু হতে পারে (যদি আপনি নমুনা ডেটা সরবরাহ করেন, যা দৃ strongly়ভাবে সুপারিশ করা হয় তবে আপনি সঠিক কোডটি গ্রহণ করতে পারবেন):
   my.anova <- lme(dv ~ group*within, data = your.df, random = ~1|id, correlation = corCompSymm(form = ~1|id))
5. anovaআনোভা টেবিলটি পাওয়ার জন্য জেনেরিক ফাংশনটি ব্যবহার করুন (দেখুন ?anova.lme):
   anova(my.anova)
   ধরণের তৃতীয় স্তরের স্কোয়ারগুলি পেতে anovaআর্গুমেন্ট typeসেট করে কমান্ডটি ব্যবহার করুন (কোডিংয়ের ক্ষেত্রে বিপরীতগুলি সেট করা থাকলে "marginal"এটি কেবল কাজ করে, পয়েন্ট 3 দেখুন):
   anova(my.anova, type = "marginal")
6. টাইপযুক্ত লাগানো বস্তুটি lmeএখন বিভিন্ন ফাংশনগুলিকে বৈপরীত্য সম্পাদন করতে দেয় allows সর্বাধিক নমনীয় সমাধান (তবে বরং অচল একটি) Lএকটি কলটিতে anova.lmeআবার যুক্তি (আবার দেখুন ?anova.lme)।
   অন্যান্য সমাধানগুলির জন্য lmeএকটি যুক্তিযুক্ত জিনিস হিসাবে একটি উপযুক্ত জিনিসও প্রয়োজন :
   এছাড়াও gmodels প্যাকেজ estimableথেকে ফাংশনটি খুব নমনীয় । এই প্যাকেজটি ফাংশনটিও সরবরাহ করে। Multcomp প্যাকেজ ব্যবহার করতে পারবেন আলফা-ত্রুটি সমন্বয় ব্যবহার বৈপরীত্য (কিন্তু আপনি শুধুমাত্র আপনার কারণের একটি ব্যবহার বৈপরীত্য সম্পাদন করতে পারবেন), ফাংশন। একটি নতুন এবং প্রতিশ্রুতিবদ্ধ পদ্ধতির বিপরীতে প্যাকেজ, যদিও, এখনও পর্যন্ত এটি সমস্ত সম্ভাব্য বৈসাদৃশ্যকে কমাতে দেখায় না।fit.contrasts
glht
   

---

একটি বিকল্প সমাধান সমন্বয় মাধ্যমে মান ANOVA ব্যবহার করা afexএবং lsmeansরূপরেখা যেমন afex-vignette ।

Comments

(+1) দুর্দান্ত এবং খুব শিক্ষামূলক প্রতিক্রিয়া। ব্লগ পোস্টের জন্য অপেক্ষা করছি ...

— chl

?anova.lmeআমার জন্য কাজ করে না এবং methods(anova)এটি একটি অ-দৃশ্যমান ফাংশন হিসাবে তালিকাবদ্ধ করে।

— জন

@ জন আপনি কি nlmeআগে লোড করেছেন ? যদি না হয়, চালান library(nlme), তবে এটি কাজ করা উচিত। এটি এখনও কাজ না করে থাকলে install.packages("nlme")প্রথমে।

— হেনরিক

আহ ... আমি আপাতদৃষ্টিতে ছিলাম না, কেবল lme4

— জন

@ হেন্রিক, আমি ৫ ম পদক্ষেপে যেতে পারি, তবে আমি বিষয়গুলির মধ্যে নির্দিষ্ট বিষয়গুলির মধ্যে নির্দিষ্ট বিষয়গুলির মধ্যে তুলনামূলক তুলনা করতে সক্ষম নই। আমি যে অনলাইন অনলাইন ডকুমেন্টেশন দেখেছি তার সবগুলি বিষয়-বিষয়গুলির বিভিন্ন স্তরের পার্থক্য পরীক্ষার সাথে সম্পর্কিত। আপনার কোন উদাহরণ কোড আছে? আপনি যদি এটি সেখানে কভার করতে চলেছেন তবে আমি কেবল ব্লগ পোস্টের জন্য অপেক্ষা করতে পারি।

— জলস্তর