স্প্লিংসের সাথে মিশ্রিত প্রভাবগুলির মডেল

Question

আমি একটি অ্যাপ্লিকেশনটিতে একটি স্প্লাইন টার্মের সাথে একটি মিশ্র প্রভাবগুলির মডেলটি ফিট করছি যেখানে সময়ের সাথে প্রবণতা কার্ভি-লিনিয়ার হিসাবে পরিচিত। তবে, আমি যা মূল্যায়ন করতে চাই তা হ'ল লাইনারিটি থেকে স্বতন্ত্র বিচ্যুতির কারণে কার্ভি-লিনিয়ার প্রবণতাটি ঘটে কিনা, বা এটি গ্রুপ স্তরের এমন একটি প্রভাব যা কোনও গ্রুপ স্তরের ফিটকে কার্ভি-লিনিয়ার প্রদর্শিত করে। আমি জেএম প্যাকেজ থেকে একটি ডেটাसेट বোরিং একটি পুনরুত্পাদনযোগ্য উদাহরণ দিতে।

library(nlme)
library(JM)
data(pbc2)

fitLME1 <- lme(log(serBilir) ~ ns(year, 2), random = ~ year | id, data = pbc2)
fitLME2 <- lme(log(serBilir) ~ year, random = ~ ns(year, 2) | id, data = pbc2)

মূলত আমি জানতে চাই যে এর মধ্যে কোনটি আমার ডেটার সাথে আরও ভাল ফিট করে। তবে তুলনা anovaআমাকে একটি অশুভ সতর্কতা দেয়:

        Model df      AIC      BIC    logLik   Test  L.Ratio p-value
fitLME1     1  7 3063.364 3102.364 -1524.682                        
fitLME2     2  9 2882.324 2932.472 -1432.162 1 vs 2 185.0399  <.0001
Warning message:
In anova.lme(fitLME1, fitLME2) :
  fitted objects with different fixed effects. REML comparisons are not meaningful.

এখন আমি সচেতন যে সর্বাধিক সম্ভাব্য পদ্ধতির মাধ্যমে এই ধরণের তুলনা করতে অসুবিধা রয়েছে - তবে বিকল্পটি কী?

Comments

পদ্ধতি = "এমএল" বা সঠিক বাক্যাংশ যাই হোক না কেন পুনরায় ফিট করার চেষ্টা করুন এবং সতর্কতাটি চলে গেছে কিনা দেখুন।

— mdewey

Answer 1

যেমনটি মডিউই বলেছে তারপরে আরএএমএল অনুমান পদ্ধতি ছাড়াই মডেলটিকে রিফাইট করুন। সতর্কবার্তাটি যেমন বলেছে, আপনার যখন বিভিন্ন স্থির প্রভাবের কাঠামো থাকে তখন তুলনাগুলি অর্থবহ হয় না।

পরবর্তী সমস্যাগুলি হ'ল মডেলগুলি বাসা বাঁধে না তাই এফ-টেস্টটি সম্ভবত বোঝায় না। আপনি তথ্যের মানদণ্ডটি দেখতে পারেন। উভয় পক্ষই fitLME2।