Lme [বন্ধ] এ একাধিক (পৃথক) এলোমেলো প্রভাব উল্লেখ করা হচ্ছে


25

আমি আর প্যাকেজ nlme এবং lme4 এ কাজ করছিলাম , একাধিক এলোমেলো প্রভাব সহ মডেলগুলি নির্দিষ্ট করার চেষ্টা করছিলাম। আমি খুঁজে পেয়েছি, যে কেবল এনএলএমই বৈকল্পিকের ভিন্ন ভিন্ন কাঠামো নির্দিষ্ট করতে দেয়। অতএব, আমি একটি মডেল পেয়েছি, যেখানে তাপমাত্রা (Y) সময় (ঘন্টা) উপর নির্ভর করে, বিরতি তারিখ এবং বছর অনুসারে পরিবর্তিত হয়, এবং বৈচিত্রটিও বছরের পর বছর পরিবর্তিত হয়:

fit1 <- lme(Y ~ time, random=~1|year/date, data=X, weights=varIdent(form=~1|year))

যাইহোক, যদি আমাকে আরও একটি এলোমেলো শব্দ যুক্ত করতে হয় (সময় অনুসারে সময় পরিবর্তিত হয়), এবং এটির মতো মডেলটি নির্দিষ্ট করুন:

fit2 <- lme(Y ~ time, random=list(~1|year, ~time-1|date,  ~1|date), data=X, 
            weights=varIdent(form=~1|year))

এলোমেলো প্রভাব একে অপরের মধ্যে নেস্টেট হয়ে যায়: বছরের তারিখ; এবং তারপরে তারিখ এবং বছরে।

আমি চেষ্টাও করেছি

one  <- rep(1, length(Y))
fit3 <- lme(Y ~ time, random=list(one=pdBlocked(list(pdSymm(~1|year/date), 
            pdSymm(~time-1|year)))), data=X, weights=varIdent(form=~1|year))

তবে এটি একটি ত্রুটি দেয়:

Error in pdConstruct.pdBlocked(object, form = form, nam = nam, data = data,  :
  cannot have duplicated column names in a "pdMat" object

আমি বুঝতে পারি যে ইতিমধ্যে একই ধরণের সমস্যা সম্পর্কিত অনেক প্রশ্ন রয়েছে তবে আমি সত্যিই আমার মামলার উত্তর পাইনি। আপনি কি মডেলটির সঠিক স্পেসিফিকেশন সম্পর্কে আমাকে সহায়তা করতে পারেন?

উত্তর:


30

অনেক লড়াইয়ের পরে আমি আমার সমস্যার সমাধান খুঁজে পেয়েছি, যা আমি এখানে এখানে পোস্ট করছি কারও কারওর মতো প্রশ্ন থাকলে:

fit <- lme(Y ~ time, random=list(year=~1, date=~time), data=X, weights=varIdent(form=~1|year))
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.