আমি একটি নকশা করা পরীক্ষা থেকে চোখের ট্র্যাকিং ডেটা বিশ্লেষণ করছি। আমার ডেটার একটি সরলিকৃত সংস্করণ দেখতে এইরকম দেখাচ্ছে (আপনি এখানে dp () ডেটা পেতে পারেন ),
head(lookDATA)
participant fixationImage fixationCount
1 9 Automobile 81
2 9 Bird 63
3 9 Chair 82
4 9 Dog 64
5 9 Face 90
6 9 Plant 75
যেখানে অংশগ্রহণকারী প্রতিটি বিষয়ের জন্য একটি অনন্য সনাক্তকারী, সেখানে ফিক্সেশনইমেজ তারা কোন চিত্র বিভাগে স্থির করেছেন এবং ফিক্সেশনকাউন্টটি সেই চিত্রের বিভাগে তারা কতবার নির্ধারণ করেছেন।
আমি lme4 প্যাকেজ থেকে গ্ল্যামার () ব্যবহার করে ডেটাতে একটি পিসন মডেল ফিট করি ।
model<-glmer(fixationCount ~ fixationImage + (1|participant), family = poisson, data = lookDATA)
আমি ফ্যাক্সার স্তরের মধ্যে পার্থক্য পরীক্ষা করতে lsmeans প্যাকেজ থেকে lsmeans () ব্যবহার করেছি ,
cld(lsmeans(model,"fixationImage"))
যা নিম্নলিখিত আউটপুট সরবরাহ করে:
fixationImage lsmean SE df asymp.LCL asymp.UCL .group
Chair 3.786022 0.05764923 NA 3.673018 3.899026 1
Bird 3.866201 0.05750641 NA 3.753476 3.978925 2
Dog 3.868768 0.05751010 NA 3.756037 3.981500 2
Body 3.883644 0.06040952 NA 3.765230 4.002059 23
Plant 3.893327 0.05746744 NA 3.780679 4.005975 23
Automobile 3.901939 0.05745528 NA 3.789315 4.014563 23
Face 3.946848 0.05832549 NA 3.832519 4.061178 3
Lsmeans vignette ব্যবহার সম্পর্কে আমার (সম্ভবত সীমাবদ্ধ) বোধ অনুসারে lsmean কলামটি মডেলের দ্বারা পূর্বাভাস দেওয়া কোনও নির্দিষ্ট বিভাগের গড় গড় সংখ্যা উপস্থাপন করবে।
যাইহোক, এই মানগুলি এই সংখ্যাগুলির সহজ বর্ণনামূলক পরিসংখ্যান থেকে অস্বস্তিকর মনে হয়,
summaryBy(fixationCount ~ fixationImage, data = lookDATA)
fixationImage fixationCount.mean
1 Automobile 55.18750
2 Bird 53.25000
3 Body 57.12821
4 Chair 50.39450
5 Dog 53.82883
6 Face 56.76389
7 Plant 54.71429
সুপারিশ করে যে lsmeans এখানে প্রতিনিধিত্ব করে আমি সঠিকভাবে বুঝতে পারছি না, বা সম্ভবত আমি মডেলটি ভুল বানিয়েছি।
কোন সহায়তা ব্যাপকভাবে প্রশংসা হবে।
$<-.data.frame
(*tmp*
, "sep", মান = ",") এ ত্রুটি হয়েছে : প্রতিস্থাপনের 1 টি সারি রয়েছে, ডেটা রয়েছে 0. রেকর্ডের জন্য আমি আর সংস্করণটি 3.1.2 (2014-10-31) 'কুমড়ো হেলমেট' এবং lsmeans সংস্করণ 2.17 ব্যবহার করছি। তবুও, আপনি আমার প্রশ্নের উত্তর দিয়েছেন এবং আমি নিজেই আউটপুটটি রূপান্তর করব। আবার ধন্যবাদ!