মিশ্র মডেলগুলির জন্য নমুনা আকারের গণনা


23

আমি ভাবছি মিশ্রিত মডেলগুলিতে নমুনা আকার গণনা করার কোনও পদ্ধতি আছে কিনা? আমি lmerমডেলগুলি ফিট করার জন্য আর এ ব্যবহার করছি (আমার এলোমেলো slালু এবং ইন্টারসেপ্ট রয়েছে)।


3
সিমুলেশন সর্বদা একটি বিকল্প - যেমন কোনও নির্দিষ্ট বিকল্প অনুমান এবং নমুনা আকারের অধীনে ডেটা সিমুলেট করুন এবং আপনি কতবার আগ্রহের নাল অনুমানকে প্রত্যাখ্যান করেন তা দেখতে মডেলটিকে অনেকবার পুনরায় ফিট করে। আমার অভিজ্ঞতা থেকে এটি বেশ (কম্পিউটার) সময় সাশ্রয়ী কারণ এটি প্রতিটি মডেলের ফিট হওয়ার জন্য কমপক্ষে কয়েক সেকেন্ড সময় নেয়।
ম্যাক্রো

উত্তর:


29

longpowerপ্যাকেজ কার্যকরী লিউ এবং লিয়াং (1997) এবং Diggle নমুনা আকার গণনার এট (2002)। ডকুমেন্টেশনের উদাহরণ কোড রয়েছে। এখানে একটি, lmmpower()ফাংশন ব্যবহার করে :

> require(longpower)
> require(lme4)
> fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy) 
> lmmpower(fm1, pct.change = 0.30, t = seq(0,9,1), power = 0.80)

     Power for longitudinal linear model with random slope (Edland, 2009) 

              n = 68.46972
          delta = 3.140186
         sig2.s = 35.07153
         sig2.e = 654.941
      sig.level = 0.05
              t = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9
          power = 0.8
    alternative = two.sided
       delta.CI = 2.231288, 4.049084
           Days = 10.46729
        Days CI = 7.437625, 13.496947
           n.CI = 41.18089, 135.61202

liu.liang.linear.power()যা পরীক্ষা করে দেখুন " লিনিয়ার মিশ্রিত মডেলের জন্য নমুনা আকারের গণনা সম্পাদন করে"

লিউ, জি।, এবং লিয়াং, কেওয়াই (1997)। সম্পর্কযুক্ত পর্যবেক্ষণগুলি সহ অধ্যয়নের জন্য নমুনা আকারের গণনা। বায়োমেট্রিকস, 53 (3), 937-47।

ডিগল পিজে, হেগারটি পিজে, লিয়াং কে, জেগার এসএল। অনুদায়ী তথ্য বিশ্লেষণ। দ্বিতীয় সংস্করণ. অক্সফোর্ড। পরিসংখ্যান বিজ্ঞান সিরিয়াস। 2002

সম্পাদনা করুন: আরেকটি উপায় হ'ল ক্লাস্টারিংয়ের প্রভাবের জন্য "সংশোধন" করা। একটি সাধারণ রৈখিক মডেলটিতে প্রতিটি পর্যবেক্ষণ স্বতন্ত্র থাকে তবে গুচ্ছ পর্যবেক্ষণের উপস্থিতিতে স্বাধীন হয় না যা কম স্বতন্ত্র পর্যবেক্ষণ বলে মনে করা যেতে পারে - কার্যকর নমুনার আকারটি আরও কম। কার্যকারিতা হ্রাস এই নকশা প্রভাব হিসাবে পরিচিত :

DE=1+(m1)ρ
mρDE

3
DEFF=1+(m1)ρxρϵ,
ρxρϵ

আপনি কি আমাকে এই সূত্রটির জন্য উদ্ধৃতি দেওয়ার দিকে নির্দেশ করতে পারেন?
জোশুয়া রোজেনবার্গ

10

সাধারণ 2 নমুনা পরীক্ষার বাইরে যে কোনও কিছুর জন্য আমি নমুনা আকার বা শক্তি অধ্যয়নের জন্য সিমুলেশন ব্যবহার করতে পছন্দ করি। প্রিপেইকেজড রুটিনগুলির সাহায্যে আপনি কখনও কখনও প্রোগ্রামগুলি থেকে যে অনুমানগুলি করছেন তার ফলাফলগুলির মধ্যে বড় পার্থক্য দেখতে পাবেন (এবং এই অনুমানগুলি কী তা আপনি খুঁজে পেতে সক্ষম হবেন না, যদি তারা আপনার গবেষণার পক্ষে যুক্তিযুক্ত হয় তবে ছেড়ে দিন)। সিমুলেশন দিয়ে আপনি সমস্ত অনুমান নিয়ন্ত্রণ করেন।

এখানে একটি উদাহরণের লিঙ্ক:
https://stat.ethz.ch/pipermail/r-sig-mixed-models/2009q1/001790.html


শুধু ভাবছি, এটি কি জিএলএমআর মডেলগুলির জন্য কাজ করে?
চার্লি গ্লেজ

1
@ কার্লোস গ্লেজ, হ্যাঁ, এটি এমন কোনও মডেলটির জন্য কাজ করে যেখানে আপনি ডেটা অনুকরণ এবং বিশ্লেষণ করতে পারেন। আমি GLMER মডেলগুলির জন্য এটি করেছি done
গ্রেগ স্নো

ঠিক আছে, এবং আমি এটিকে যুক্ত করব "অনুমানগুলি নিয়ন্ত্রণ" ছাড়াও, আপনি "কী হতে পারে" প্রশ্নগুলি জিজ্ঞাসা করতে পারেন, এই অনুমানগুলি ভেঙে ফেলতে পারেন এবং দৃ rob়তার কিছু ব্যবহারিক ধারণা নির্ধারণ করতে পারেন, উদাহরণস্বরূপ, অ-স্বাভাবিক এলোমেলো প্রভাবগুলি সত্যই দক্ষতা নষ্ট করে কিনা।
অ্যাডামো
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.