সম্পাদনা 2: আমি প্রাথমিকভাবে ভেবেছিলাম যে আমার একটি ফ্যাক্টরের উপর বারবার ব্যবস্থা নিয়ে একটি দ্বি-গুণক আনোভা চালানো দরকার, তবে আমি এখন মনে করি একটি লিনিয়ার মিশ্র-প্রভাব মডেল আমার ডেটার জন্য আরও ভাল কাজ করবে। আমার মনে হয় আমি কী হতে হবে তা প্রায় জানি, তবে এখনও কয়েকটি পয়েন্টে আমি বিভ্রান্ত।
বিশ্লেষণের জন্য আমার যে পরীক্ষাগুলি প্রয়োজন তা এ জাতীয় চেহারা:
- বিষয়গুলি বেশ কয়েকটি চিকিত্সা গ্রুপকে দেওয়া হয়েছিল
- প্রতিটি বিষয়ের পরিমাপ একাধিক দিনে নেওয়া হয়েছিল
- তাই:
- বিষয় চিকিত্সার মধ্যে নেস্ট করা হয়
- চিকিত্সা দিনের সাথে পার হয়
(প্রতিটি বিষয় কেবল একটি চিকিত্সার জন্য বরাদ্দ করা হয়, এবং প্রতিটি বিষয়ে প্রতিটি বিষয়ে পরিমাপ নেওয়া হয়)
আমার ডেটাসেটে নিম্নলিখিত তথ্য রয়েছে:
- বিষয় = ব্লকিং ফ্যাক্টর (এলোমেলো ফ্যাক্টর)
- দিন = বিষয় বা পুনরাবৃত্তি ব্যবস্থা ফ্যাক্টরের মধ্যে (স্থির ফ্যাক্টর)
- চিকিত্সা = বিষয় ফ্যাক্টরের মধ্যে (স্থির ফ্যাক্টর)
- ওবস = পরিমাপক (নির্ভরশীল) পরিবর্তনশীল
আপডেট ঠিক আছে, তাই আমি গিয়ে একজন পরিসংখ্যানবিদদের সাথে কথা বললাম, তবে সে একজন এসএএস ব্যবহারকারী। তিনি মনে করেন যে মডেলটি হওয়া উচিত:
চিকিত্সা + দিন + বিষয় (চিকিত্সা) + দিন * বিষয় (চিকিত্সা)
স্পষ্টতই তাঁর স্বরলিপিটি আর সিনট্যাক্স থেকে পৃথক, তবে এই মডেলটির জন্য অ্যাকাউন্টিং করার কথা রয়েছে:
- চিকিত্সা (স্থির)
- দিন (স্থির)
- চিকিত্সা * দিনের মিথস্ক্রিয়া
- বিষয় চিকিত্সার মধ্যে নেস্টেড (এলোমেলো)
- "চিকিত্সার মধ্যে সাবজেক্ট" (এলোমেলো) দিয়ে দিনটি পার হয়ে গেল
সুতরাং, এটি কি সঠিক বাক্য গঠন?
m4 <- lmer(Obs~Treatment*Day + (1+Treatment/Subject) + (1+Day*Treatment/Subject), mydata)আমি বিশেষত উদ্বিগ্ন যে দিনটি "চিকিত্সার মধ্যে সাবজেক্টের" অংশটি সঠিকভাবে পার হয়ে গেছে কিনা। এসএএসের সাথে পরিচিত কেউ, বা আত্মবিশ্বাসী যে তাঁর মডেলটিতে কী চলছে তা তারা বুঝতে পেরেছেন, আর সিনট্যাক্সে আমার দুঃখজনক প্রচেষ্টা মেলে কিনা সে বিষয়ে মন্তব্য করতে সক্ষম?
এখানে একটি মডেল তৈরি এবং সিনট্যাক্স লেখার আমার আগের প্রচেষ্টাগুলি (উত্তর এবং মন্তব্যে আলোচিত):
m1 <- lmer(Obs ~ Treatment * Day + (1 | Subject), mydata)বিষয়টিকে চিকিত্সার মধ্যে বাসা বেঁধে দেওয়া হয়েছে কীভাবে আমি তা মোকাবিলা করব? এর থেকে কীভাবে m1আলাদা হয়:
m2 <- lmer(Obs ~ Treatment * Day + (Treatment|Subject), mydata)
m3 <- lmer(Obs ~ Treatment * Day + (Treatment:Subject), mydata)
এবং হয় m2এবং m3সমতুল্য (এবং যদি না, কেন)?
এছাড়াও, আমি যদি সম্পর্কের কাঠামোটি (যেমন correlation = corAR1) নির্দিষ্ট করতে চান তবে lme4 এর পরিবর্তে nlme ব্যবহার করা দরকার ? পুনরাবৃত্তি ব্যবস্থা অনুসারে , একটি ফ্যাক্টরের উপর বারবার ব্যবস্থা নিয়ে পুনরাবৃত্তি-বিশ্লেষণের জন্য, সমবায় কাঠামো (একই বিষয়ের পরিমাপের মধ্যে পারস্পরিক সম্পর্কের প্রকৃতি) গুরুত্বপূর্ণ।
যখন আমি বারবার ব্যবস্থা আনোভা করার চেষ্টা করছিলাম তখন আমি টাইপ II এসএস ব্যবহার করার সিদ্ধান্ত নিয়েছি; এটি কি এখনও প্রাসঙ্গিক, এবং যদি তাই হয় তবে আমি কীভাবে তা উল্লেখ করব?
ডেটা দেখতে কেমন লাগে তার উদাহরণ এখানে:
mydata <- data.frame(
Subject = c(13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 29, 30, 31, 32, 33,
34, 35, 36, 37, 38, 39, 40, 62, 63, 64, 65, 13, 14, 15, 16, 17, 18,
19, 20, 21, 22, 23, 24, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39,
40, 62, 63, 64, 65, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 62, 63, 64, 65),
Day = c(rep(c("Day1", "Day3", "Day6"), each=28)),
Treatment = c(rep(c("B", "A", "C", "B", "C", "A", "A", "B", "A", "C", "B", "C",
"A", "A", "B", "A", "C", "B", "C", "A", "A"), each = 4)),
Obs = c(6.472687, 7.017110, 6.200715, 6.613928, 6.829968, 7.387583, 7.367293,
8.018853, 7.527408, 6.746739, 7.296910, 6.983360, 6.816621, 6.571689,
5.911261, 6.954988, 7.624122, 7.669865, 7.676225, 7.263593, 7.704737,
7.328716, 7.295610, 5.964180, 6.880814, 6.926342, 6.926342, 7.562293,
6.677607, 7.023526, 6.441864, 7.020875, 7.478931, 7.495336, 7.427709,
7.633020, 7.382091, 7.359731, 7.285889, 7.496863, 6.632403, 6.171196,
6.306012, 7.253833, 7.594852, 6.915225, 7.220147, 7.298227, 7.573612,
7.366550, 7.560513, 7.289078, 7.287802, 7.155336, 7.394452, 7.465383,
6.976048, 7.222966, 6.584153, 7.013223, 7.569905, 7.459185, 7.504068,
7.801867, 7.598728, 7.475841, 7.511873, 7.518384, 6.618589, 5.854754,
6.125749, 6.962720, 7.540600, 7.379861, 7.344189, 7.362815, 7.805802,
7.764172, 7.789844, 7.616437, NA, NA, NA, NA))