এক স্তরের এবং অন্যগুলির গড়ের পার্থক্যের জন্য একটি বিপরীতে ম্যাট্রিক্স (আরে) কীভাবে নির্দিষ্ট করবেন?

আমার কাছে একটি রিগ্রেশন মডেল রয়েছে যা দেখতে দেখতে:

Y = β_{0} + β_{1} X_{1} + β_{2} X_{2} + β_{3} X_{3} + β_{12} X_{1} X_{2} + β_{13} X_{1} X_{3} + β_{123} X_{1} X_{2} X_{3}

$Y = \beta_0+\beta_1X_1 + \beta_2X_2 + \beta_3X_3 +\beta_{12}X_1X_2+\beta_{13}X_1X_3+\beta_{123}X_1X_2X_3$

... বা আর সংকেত: y ~ x1 + x2 + x3 + x1:x2 + x1:x3 + x1:x2:x3

ধরা যাক এবং শ্রেণীবদ্ধ ভেরিয়েবল এবং সংখ্যাসূচক। হ'ল এর তিনটি স্তর রয়েছে এবং মান বিপরীতে পরিবর্তনের পরিবর্তে, আমাকে পরীক্ষা করা দরকার: $X_1$ $X_2$ $X_3$ $X_1$ $X_{1a}, X_{1b}, X_{1c}$

স্তরের জন্য পথিমধ্যে কিনা উল্লেখযোগ্যভাবে স্তরের জন্য গড় পথিমধ্যে থেকে পৃথক এবং । $X_{1a}$ $X_{1b}$ $X_{1c}$
কিনা প্রতিক্রিয়া স্তরের মধ্যে উল্লেখযোগ্যভাবে ভিন্ন হয় এবং স্তরের গড় এবং । $X_2$ $X_{1a}$ $X_{1b}$ $X_{1c}$
কিনা ঢাল স্তরের মধ্যে উল্লেখযোগ্যভাবে ভিন্ন হয় এবং স্তরের গড় এবং । $X_3$ $X_{1a}$ $X_{1b}$ $X_{1c}$

এই পোস্টের উপর ভিত্তি করে মনে হচ্ছে যে ম্যাট্রিক্স আমি চাই তা হ'ল ...

 2
-1
-1

তাই আমি করি contrasts(mydata$x1)<-t(ginv(cbind(2,-1,-1)))। অনুমান পরিবর্তন হয় তবে অন্যরাও তাই করে। আমি এই স্তরের দ্বিগুণ মান থেকে এবং গ্রুপের পূর্বাভাসিত মানগুলিকে (যখন এবং তার রেফারেন্স স্তরে থাকে) বিয়োগ করে এর নতুন অনুমান পুনরুত্পাদন করতে পারি । তবে আমি বিশ্বাস করতে পারি না যে আমি আমার বিপরীতে ম্যাট্রিক্সটি সঠিকভাবে নির্দিষ্ট করেছি যদি না আমি একইভাবে অন্যান্য সহগগুলিও অর্জন করতে পারি। $\beta_1$ $beta_1$ $X_1b$ $X_1c$ $X_3=0$ $X_2$ $X_1a$

কোষের মাধ্যম এবং বৈপরীত্যের মধ্যকার সম্পর্কের চারপাশে কীভাবে আমার মাথা মোড়ানো যায় তার জন্য কারও কাছে কোনও পরামর্শ আছে? ধন্যবাদ। এই ধরণের বৈপরীত্যের জন্য কি কোনও আদর্শ নাম রয়েছে?

আহা! গ্লেেন_ব এর উত্তরে পোস্ট করা লিঙ্ক অনুসারে , নীচের লাইনটি হল, আপনি যে কোনও গ্রুপের তুলনা করতে চাইলে আপনি কোনও আর-স্টাইলের বিপরীতে বৈশিষ্ট্যে রূপান্তর করতে চান:

একটি বর্গ ম্যাট্রিক্স তৈরি করুন। সারিগুলি আপনার ফ্যাক্টরের স্তরের প্রতিনিধিত্ব করে এবং কলামগুলি বিপরীতে উপস্থাপন করে। প্রথমটি ব্যতীত, যা মডেলটিকে জানায় যে ইন্টারসেপ্টটি কী উপস্থাপন করবে।
আপনি যদি নিজের ইন্টারসেপটকে গ্র্যান্ড মানে হতে চান তবে প্রথম কলামটি একই অ-শূন্য মানের সাথে পূর্ণ করুন, তাতে কিছু আসে যায় না। আপনি যদি ইন্টারসেপটকে স্তরগুলির অন্যতম হতে চান তবে সেই সারিতে একটি নম্বর দিন এবং বাকীটি শূন্যে পূরণ করুন। আপনি যদি ইন্টারসেপটকে বিভিন্ন স্তরের গড় হিসাবে বোঝাতে চান তবে সেই সারিগুলিতে নম্বর দিন এবং বাকী অংশে শূন্য করুন। আপনি যদি এটি একটি ওজনযুক্ত গড় হতে চান তবে বিভিন্ন সংখ্যা ব্যবহার করুন, অন্যথায় একই সংখ্যাটি ব্যবহার করুন। আপনি এমনকি বিরতি কলামে নেতিবাচক মান রাখতে পারেন এবং এর সম্ভবত খুব কিছু বোঝা যায় তবে এটি অন্যান্য বিপরীতে সম্পূর্ণরূপে পরিবর্তিত হয়, তাই এর জন্য আমার কোন ধারণা নেই
অন্যদের তুলনায় আপনি কোন স্তরের চেয়ে তুলনা করতে চান তা ইতিবাচক এবং নেতিবাচক মান সহ বাকী কলামগুলি পূরণ করুন। আমি ভুলে গিয়েছি কেন শূন্যের যোগফল গুরুত্বপূর্ণ, তবে মানগুলি সমন্বয় করুন যাতে কলামগুলি শূন্যের সমষ্টি হয়।
t()ফাংশনটি ব্যবহার করে ম্যাট্রিক্স স্থানান্তর করুন ।
প্যাকেজ ginv()থেকে ব্যবহার করুন MASSবা solve()ট্রান্সপোজড ম্যাট্রিক্সের বিপরীত পেতে পারেন।
প্রথম কলামটি ফেলে দিন, যেমন mycontrast<-mycontrast[,-1]। আপনার কাছে এখন এপিএক্স পি -১ ম্যাট্রিক্স রয়েছে, তবে আপনি আপনার ইন্টারসেপ্টের জন্য যে তথ্য রেখেছিলেন তা পদক্ষেপ 5 এর সময় পুরো ম্যাট্রিক্সে এনকোড করা হয়েছিল।
যদি আপনি lm()এট আল-এর ডিফল্ট আউটপুটের তুলনায় সংক্ষিপ্ত আউটপুটে লেবেলগুলি পড়তে আরও আনন্দদায়ক চান তবে সেই অনুযায়ী আপনার ম্যাট্রিক্সের কলামগুলির নাম দিন। তবে বিরতি সর্বদা স্বয়ংক্রিয়ভাবে নাম দেওয়া হবে (Intercept)।
আপনার ম্যাট্রিক্সকে প্রশ্নের ফ্যাক্টরের জন্য নতুন বিপরীতে পরিণত করুন contrasts(mydata$myfactor)<-mymatrix
চালান lm()(এবং সম্ভবত অন্যান্য অনেক ফাংশন যে সূত্র ব্যবহার করুন) লোড না করেও মান আর স্বাভাবিক হিসাবে glht, doByঅথবা contrasts।

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

r contrasts

— f1r3br4nd
সূত্র

পরবর্তী সমস্ত ভেরিয়েবলগুলির গড়ের সাথে একটির তুলনা হ'ল (স্কেল বাদে), যা হেলমার্ট কোডিং বা হেলমার্ট বিপরীতে বলে । আপনি যেটি দেবেন সেটি হ'ল প্রথম বৈপরীত্য, অন্যটি of আকারযুক্ত সংস্করণ । $(0, 1, -1)^\top$

আর হেলমার্ট কোডিং কে আর বলে, এটি 'রিভার্স হেলমার্ট' বলে। এগুলি পরিবর্তনশীল ক্রমের পরিবর্তনের সমতুল্য।

— গ্লেন_বি -রাইনস্টেট মনিকা
সূত্র

"পরিবর্তনশীল ক্রম পরিবর্তন পর্যন্ত" সমতুল্য হওয়া উচিত "স্তরের ক্রমের বিপরীত দেওয়া"? সুতরাং এসএএস / এসপিএসএস ধরণের হেলমার্ট (অবশিষ্ট স্তরের গড়ের সাথে তুলনা) পেতে, কোনওটিকে সংশ্লিষ্ট ফ্যাক্টরের স্তরের ক্রমটি বিপরীত করতে হবে, বা প্রতিটি কলামটি পুনরায় () পুনরুদ্ধার করতে হবে (ম্যাট্রিক্সের প্রতিটি সারি) দ্বারা ফিরে contra.helmert?

— টিম

@ টিম "ভেরিয়েবল ক্রমের পরিবর্তন" এর অর্থ " ভেরিয়েবল ক্রমের কোনও পরিবর্তন " বোঝানো নয় , তবে বেশ আক্ষরিক অর্থে

$\:$ একটি পরিবর্তনশীল আদেশের পরিবর্তন করুন (যে, বিশেষ করে এক)।

— গ্লেন_বি -রিনস্টেট মনিকা