আপনার ম্যাকনামারের পরীক্ষা দরকার ( http://en.wikedia.org/wiki/McNemar%27s_test , http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3346204/ )। নিম্নলিখিত একটি উদাহরণ:
1300 pts এবং 1300 ম্যাচ নিয়ন্ত্রণগুলি অধ্যয়ন করা হয়। ধূমপানের স্থিতি নিম্নরূপে পেশ করা হয়েছে:
Normal
|no |yes|
Cancer|No |1000|40 |
|Yes |200 |60 |
সারণীর প্রতিটি এন্ট্রি একটি CASE- নিয়ন্ত্রণ পেয়ার সম্পর্কিত তথ্য দেখায়: 1000 এর অর্থ 1000 কেস-নিয়ন্ত্রণ জোড়গুলির মধ্যে, না কোনও ধূমপায়ী ছিল। 40 হ'ল কেস-কন্ট্রোল জুটির সংখ্যা যেখানে নিয়ন্ত্রণ ছিল ধূমপায়ী এবং ক্যান্সারের রোগী ছিল না, ইত্যাদি। এই টেবিলটি তৈরি করতে এবং ম্যাকনেমার পরীক্ষা করতে নিম্নলিখিত আর কোড ব্যবহার করা যেতে পারে।
mat = as.table(rbind(c(1000, 40), c( 200, 60) ))
colnames(mat) <- rownames(mat) <- c("Nonsmoker", "Smoker")
names(dimnames(mat)) = c("Cancer", "Normal")
mat
# Normal
# Nonsmoker Smoker
# Cancer
# Nonsmoker 1000 40
# Smoker 200 60
mcnemar.test(mat)
# McNemar's Chi-squared test with continuity correction
#
#data: mat
#McNemar's chi-squared = 105.34, df = 1, p-value < 2.2e-16
ম্যাকনামারের পরীক্ষাটিও বাইনারি ফলাফলের চলকটিতে হস্তক্ষেপের প্রভাব নির্ধারণের জন্য ব্যবহৃত হয়। উপরোক্ত ফলাফলের পূর্বের জুটি উপরে তালিকাভুক্ত এবং পরীক্ষিত।
সম্পাদনা করুন: @ গুং দ্বারা প্রদত্ত উদাহরণস্বরূপ বাড়ানো হচ্ছে, যদি ধূমপানের স্থিতিটি আপনার ডেটাফ্রেম মাইডিএফ-তে তালিকাভুক্ত থাকে:
pairID cancer control
1 1 1
2 1 1
3 1 0
...
ম্যাকনামার্স পরীক্ষা নিম্নলিখিত আর কমান্ড দিয়ে করা যেতে পারে:
> tt = with(mydf, table(cancer, control))
> tt
control
cancer 0 1
0 5 1
1 3 2
> mcnemar.test(tt)
McNemar`s Chi-squared test with continuity correction
data: tt
McNemar`s chi-squared = 0.25, df = 1, p-value = 0.6171