এই তিনটি মডেলের প্রত্যেকটির জন্য মডেল সূত্রটি লিখে এটি স্পষ্ট হয়ে উঠতে পারে। যাক ব্যক্তির জন্য পর্যবেক্ষণ হতে আমি সাইটে ঞ প্রতিটি মডেল এবং সংজ্ঞায়িত একটি আমি ঞ , টি আমি ঞ অনুরূপভাবে আপনার মডেল ভেরিয়েবল উল্লেখ করতে।YijijAij,Tij
glmer(counts ~ A + T, data=data, family="Poisson")
মডেল হয়
log(E(Yij))=β0+β1Aij+β2Tij
যা কেবলমাত্র একটি সাধারণ পিসন রিগ্রেশন মডেল।
glmer(counts ~ (A + T|Site), data=data, family="Poisson")
মডেল হয়
log(E(Yij))=α0+ηj0+ηj1Aij+ηj2Tij
যেখানে এলোমেলো প্রভাব যা সাইট জে থেকে পৃথক পৃথক প্রতিটি পর্যবেক্ষণ দ্বারা ভাগ করা হয় । আপনার নির্ধারিত মডেলটিতে এই এলোমেলো প্রভাবগুলি নিখরচায় সম্পর্কযুক্ত (যেমন Σ এর উপর কোনও বিধিনিষেধ তৈরি হয় না) অনুমোদিত। স্বাধীনতা চাপিয়ে দেওয়ার জন্য আপনাকে এগুলি বিভিন্ন বন্ধনীগুলির মধ্যে রাখতে হবে , যেমন এটি করবে। এই মডেলটি ধরে নেয় যে লগ ( E ( Y i j।)ηj=(ηj0,ηj1,ηj2)∼N(0,Σ)jΣ(A-1|Site) + (T-1|Site) + (1|Site)
হয় α 0 সব সাইটগুলির জন্য কিন্তু প্রতিটি সাইটের একটি র্যান্ডম অফসেট (হয়েছে η ঞ 0 ) ও একটি র্যান্ডম রৈখিক উভয় সঙ্গে সম্পর্ক একজন আমি ঞ , টি আমি ঞ ।log(E(Yij))α0ηj0Aij,Tij
glmer(counts ~ A + T + (T|Site), data=data, family="Poisson")
মডেল হয়
log(E(Yij))=(θ0+γj0)+θ1Aij+(θ2+γj1)Tij
তাই এখন সঙ্গে কিছু "গড়" সম্পর্ক রয়েছে একটি আমি ঞ , টি আমি ঞ , স্থায়ী প্রভাব কর্তৃক প্রদত্ত θ 0 , θ 1 , θ 2 কিন্তু সেই সম্পর্ক প্রতিটি সাইটের এবং যারা পার্থক্য জন্য ভিন্ন র্যান্ডম প্রভাব, দ্বারা ধরা রয়েছে γ ঞ 0 , γ ঞ 1 , γ ঞ 2log(E(Yij))Aij,Tijθ0,θ1,θ2γj0,γj1,γj2। অর্থাৎ, বেসলাইনটি এলোমেলোভাবে স্থানান্তরিত হয় এবং দুটি ভেরিয়েবলের theালু এলোমেলোভাবে স্থানান্তরিত হয় এবং একই সাইট থেকে প্রত্যেকে একই র্যান্ডম শিফট ভাগ করে দেয়।
টি কি? এটি কি এলোমেলো প্রভাব? একটি নির্দিষ্ট প্রভাব? দুটি জায়গায় টি রেখে আসলে কী সম্পাদিত হচ্ছে?
আপনার সহকারীদের মধ্যে একটি। এটি কোনও এলোমেলো প্রভাব নয় -এটি একটি এলোমেলো প্রভাব। উপরের মডেলটিতে- γ j 1 দ্বারা প্রদত্ত এলোমেলো প্রভাবের উপর নির্ভর করে টি এর একটি স্থির প্রভাব রয়েছে। এই র্যান্ডম এফেক্টটি অন্তর্ভুক্ত করে যা সম্পাদিত হয় তা হ'ল টি এবং লগ ( E ( Y i j ) ) এর মধ্যে সম্পর্কের সাইটগুলির মধ্যে বৈচিত্র্যের অনুমতি দেওয়া।TSite
TSite
γj1Tlog(E(Yij))
মডেল সূত্রের এলোমেলো প্রভাব বিভাগে কেবল কখন কিছু উপস্থিত হওয়া উচিত?
এটি আবেদনের প্রসঙ্গে কী বোঝায় তা বিষয় of
ইন্টারসেপ্ট সম্পর্কিত - অনেক কারণের জন্য আপনার সেখানে স্থির বিরতি রাখা উচিত (দেখুন, উদাহরণস্বরূপ, এখানে ); Re: এলোমেলো পথিমধ্যে, , এই প্রাথমিকভাবে একই সাইট এ প্রণীত পর্যবেক্ষণ মধ্যে পারস্পরিক প্রবৃত্ত কাজ করে। যদি এই জাতীয় পারস্পরিক সম্পর্কের উপস্থিতি বোঝায় না, তবে এলোমেলো প্রভাবটি বাদ দেওয়া উচিত।γj0
এলোমেলো slালু সম্পর্কে, কেবলমাত্র এলোমেলো slালু এবং কোনও স্থির opালু সহ একটি মডেল একটি বিশ্বাসকে প্রতিফলিত করে যে, প্রতিটি সাইটের জন্য প্রতিটি সাইটের জন্য এবং আপনার কোভেরিয়ের মধ্যে কিছুটা সম্পর্ক রয়েছে , তবে আপনি যদি সেগুলি গড় করেন তবে সমস্ত সাইটের উপর প্রভাব, তারপরে কোনও সম্পর্ক নেই। উদাহরণস্বরূপ, যদি আপনার টিতে একটি এলোমেলো ope াল থাকে তবে কোনও স্থির opeাল না থাকে, তবে এটির মতো সময়টি বলতে গেলে গড়ের কোনও প্রভাব হয় না (উদাহরণস্বরূপ তথ্যগুলিতে কোনও ধর্মনিরপেক্ষ প্রবণতা নেই) তবে প্রতিটি সময়ের সাথে সাথে এলোমেলো পথে চলেছে, যার অর্থ হতে পারে। আবার এটি আবেদনের উপর নির্ভর করে।log(E(Yij))TSite
নোট করুন যে এটি ঘটছে কিনা তা দেখতে আপনি এলোমেলো প্রভাব সহ এবং মডেলটিকে ফিট করতে পারেন - আপনার স্থির মডেলের কোনও প্রভাব দেখা উচিত নয় তবে পরবর্তী মডেলটিতে উল্লেখযোগ্য এলোমেলো প্রভাব দেখা উচিত। আমি আপনাকে সাবধান করে দিতে পারি যে এই জাতীয় সিদ্ধান্তগুলি মডেল নির্বাচনের পরিবর্তে অ্যাপ্লিকেশনটির বোঝার ভিত্তিতে প্রায়শই ভাল করা হয়।