আমি আমার দ্বিমুখী আনোয়ার পরে আর-এর সাথে টুকিএইচএসডি পোস্ট-হক পরীক্ষা করতে পছন্দ করব, তাৎপর্যপূর্ণ পার্থক্যের ভিত্তিতে দলবদ্ধভাবে সাজানো জোড়যুক্ত একটি টেবিল পেয়েছি। (শব্দটি সম্পর্কে দুঃখিত, আমি এখনও পরিসংখ্যান সহ নতুন।)
আমি এই জাতীয় কিছু পেতে চাই:
সুতরাং, তারা বা অক্ষর সহ গ্রুপবদ্ধ।
কোন ধারণা? আমি প্যাকেজ HSD.test()থেকে ফাংশনটি পরীক্ষা করেছি agricolae, তবে মনে হচ্ছে এটি দ্বি-মুখী সারণীগুলি পরিচালনা করে না।
উত্তর:
agricolae::HSD.testফাংশন ঠিক করে, কিন্তু আপনি তা জেনে রাখা প্রয়োজন যে করতে হবে আপনি একটি মিথষ্ক্রিয়া মেয়াদে আগ্রহী । স্টাটা ডেটাসেটের সাথে এখানে একটি উদাহরণ দেওয়া হয়েছে:
library(foreign)
yield <- read.dta("http://www.stata-press.com/data/r12/yield.dta")
tx <- with(yield, interaction(fertilizer, irrigation))
amod <- aov(yield ~ tx, data=yield)
library(agricolae)
HSD.test(amod, "tx", group=TRUE)
এটি নীচে প্রদর্শিত ফলাফল দেয়:
Groups, Treatments and means
a 2.1 51.17547
ab 4.1 50.7529
abc 3.1 47.36229
bcd 1.1 45.81229
cd 5.1 44.55313
de 4.0 41.81757
ef 2.0 38.79482
ef 1.0 36.91257
f 3.0 36.34383
f 5.0 35.69507
তারা নীচের কমান্ডগুলির সাথে আমরা কী অর্জন করব তা মিলছে:
. webuse yield
. regress yield fertilizer##irrigation
. pwcompare fertilizer#irrigation, group mcompare(tukey)
-------------------------------------------------------
| Tukey
| Margin Std. Err. Groups
----------------------+--------------------------------
fertilizer#irrigation |
1 0 | 36.91257 1.116571 AB
1 1 | 45.81229 1.116571 CDE
2 0 | 38.79482 1.116571 AB
2 1 | 51.17547 1.116571 F
3 0 | 36.34383 1.116571 A
3 1 | 47.36229 1.116571 DEF
4 0 | 41.81757 1.116571 BC
4 1 | 50.7529 1.116571 EF
5 0 | 35.69507 1.116571 A
5 1 | 44.55313 1.116571 CD
-------------------------------------------------------
Note: Margins sharing a letter in the group label are
not significantly different at the 5% level.
Multcomp প্যাকেজ এছাড়াও সিম্বলিক ঠাহর উপলব্ধ করা হয় ( 'কম্প্যাক্ট চিঠি প্রদর্শন', দেখতে কম্প্যাক্ট পত্র প্রদর্শনের জন্য আলগোরিদিম: তুলনা এবং মূল্যায়ন আরো বিস্তারিত জানার জন্য) উল্লেখযোগ্য pairwise তুলনা, যদিও এটি তাদের একটি ট্যাবুলার বিন্যাসে উপস্থাপন করা হয় না। যাইহোক, এটির একটি প্লট করার পদ্ধতি রয়েছে যা বক্সপ্লটগুলি ব্যবহার করে সুবিধামত ফলাফল প্রদর্শন করতে দেয়। উপস্থাপনা ক্রমের পাশাপাশি (বিকল্প decreasing=) পরিবর্তন করা যেতে পারে , এবং এতে একাধিক তুলনার জন্য আরও অনেক বিকল্প রয়েছে। রয়েছে multcompView প্যাকেজ, যার সুসংবাদ বৈশিষ্ট্য প্রসারিত করে।
এখানে বিশ্লেষণ করা একই উদাহরণ glht:
library(multcomp)
tuk <- glht(amod, linfct = mcp(tx = "Tukey"))
summary(tuk) # standard display
tuk.cld <- cld(tuk) # letter-based display
opar <- par(mai=c(1,1,1.5,1))
plot(tuk.cld)
par(opar)
একই স্তরের ভাগ করে নেওয়ার চিকিত্সা নির্বাচিত স্তরে (ডিফল্ট, 5%) উল্লেখযোগ্যভাবে আলাদা নয়।
ঘটনাচক্রে, একটি নতুন প্রকল্প রয়েছে যা বর্তমানে আর-ফোর্জে হোস্ট করা হয়েছে, যা আশাব্যঞ্জক দেখাচ্ছে: ফ্যাক্টরপ্লট । এটিতে লাইন এবং বর্ণ-ভিত্তিক প্রদর্শনগুলি পাশাপাশি সমস্ত জোড়াযুক্ত তুলনার একটি ম্যাট্রিক্স ওভারভিউ (একটি স্তর প্লটের মাধ্যমে) অন্তর্ভুক্ত রয়েছে। একটি কার্যকারী কাগজ এখানে পাওয়া যাবে: ফ্যাক্টরপ্লট: জিএলএমগুলিতে সাধারণ বৈপরীত্যগুলির উপস্থাপনা উন্নত করা
insert_absorb()যা থেকে চিকিত্সার জুড়ি বের করার চেষ্টা করে। আপনি সম্ভবত আপনার মিথস্ক্রিয়া শব্দটির কোডিং স্তরের জন্য যে বিভাজকটি ব্যবহার করেছেন তা পরিবর্তন করার চেষ্টা করতে পারেন? একটি কাজের উদাহরণ না থাকলে, কী হয়েছে তা বলা শক্ত।
console=TRUEমধ্যে HSD.test()অর্ডার টেবিল পেতে, কেস কেউ এই চেষ্টা করে এবং কোন রেজাল্ট হবে। সম্ভবত একটি আপডেট agricolae।
TukeyHSDহেল্প ফাইল অনুসারে একটি ফাংশন রয়েছে যা কভারেজের নির্দিষ্ট পরিবার ভিত্তিক সম্ভাবনার সাথে একটি ফ্যাক্টরের মাত্রার মাধ্যমগুলির মধ্যে পার্থক্যের উপর আত্মবিশ্বাসের ব্যবস্থার একটি সেট গণনা করে। বিরতিগুলি স্টুডেনটেজড রেঞ্জের পরিসংখ্যান, টুকির "সত্যবাদী গুরুত্বপূর্ণ পার্থক্য" পদ্ধতির উপর ভিত্তি করে। এটি কি আপনি চান কি না?
http://stat.ethz.ch/R-manual/R-patched/library/stats/html/TukeyHSD.html