বলুন আমরা সময়ের সাথে সাথে কিছু গণনা উপাত্তের এলোমেলো প্রভাবের মডেলটিতে কাজ করছি এবং আমরা কিছু ট্রেন্ডের জন্য নিয়ন্ত্রণ করতে চাই। সাধারণত, আপনি যেমন কিছু করতে চাই:
lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson")এর জন্য একটি চতুর্ভুজ আকৃতি অন্তর্ভুক্ত করতে t। এই সম্পর্কটিকে মডেল করার জন্য কি কোনও স্বচ্ছ মসৃণ বা স্প্লাইমের মতো আরও কিছু পরিশীলিত স্মুথ টেকনিক ব্যবহার করা সম্ভব?
উত্তর:
আপনি যা দেখান তা যদি lmerএলোমেলো প্রভাবের শর্তের জন্য কোনও সূত্রের জন্য কাজ করে তবে প্রাসঙ্গিক ভিত্তি ফাংশনগুলি সেট আপ করার জন্য আপনার স্প্লিংস প্যাকেজ থেকে ফাংশনগুলি ব্যবহার করতে সক্ষম হবেন ।
require("lme4")
require("splines")
lmer(counts ~ dependent_variable + (bs(t) | ID), family="poisson")
আপনি যা করতে চান তার উপর নির্ভর করে আপনার গ্যাম 4 প্যাকেজ এবং এমজিসিভি প্যাকেজটিও দেখতে হবে। প্রাক্তনটি মূলত উপরের কলটিতে bs()কিছুটা আনুষ্ঠানিকভাবে আনুষ্ঠানিকভাবে চলছে lmer()এবং বিশ্লেষণের অংশ হিসাবে মসৃণতা নির্বাচন সম্পাদনের অনুমতি দেয়। ফাংশন সহ gam()উত্তরোত্তরগুলি এই জাতীয় ফিটিং মডেলগুলিতে কিছুটা নমনীয়তার অনুমতি দেয় (যদি আপনি বুঝতে চান যে আপনি কী করার চেষ্টা করছেন)। দেখে মনে হচ্ছে আপনি এর মধ্যে আলাদা ট্রেন্ড চান ID? আরও স্থির প্রতিক্রিয়া পদ্ধতির কিছু হবে:
gam(counts ~ dependent_variable + ID + s(t, by = ID) , family="poisson")এলোমেলো প্রভাব অন্তর্ভুক্ত করা যেতে পারে gam()ব্যবহার মডেলের s(foo, bs = "re")প্রকার পদ যেখানে fooহবে IDআপনার উদাহরণে। byশব্দ ধারণাটি একটি এলোমেলো প্রভাবের সাথে একত্রিত করার অর্থটি বোধগম্য কিনা সে সম্পর্কে ভাবার বিষয় এবং আমি মন্তব্য করার যোগ্য এমন কিছু নয় not