আমি দ্বিপদী প্রতিক্রিয়া ভেরিয়েবল এবং একটি স্পষ্টিকর ভবিষ্যদ্বাণী নিয়ে একটি গ্ল্যাম চালাচ্ছি। এলোমেলো প্রভাব ডেটা সংগ্রহের জন্য ব্যবহৃত নেস্টেড ডিজাইন দ্বারা দেওয়া হয়। ডেটা দেখতে এমন দেখাচ্ছে:
m.gen1$treatment
[1] sucrose control protein control no_injection .....
Levels: no_injection control sucrose protein
m.gen1$emergence
[1] 1 0 0 1 0 1 1 1 1 1 1 0 0....
> m.gen1$nest
[1] 1 1 1 2 2 3 3 3 3 4 4 4 .....
Levels: 1 2 3 4 5 6 8 10 11 13 15 16 17 18 20 22 24আমি চালিত প্রথম মডেলটি দেখতে এমন দেখাচ্ছে
m.glmm.em.<-glmer(emergence~treatment + (1|nest),family=binomial,data=m.gen1)
আমি দেখতে দুটি সতর্কতা পেয়েছি যা দেখতে দেখতে:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.0240654 (tol = 0.001, component 4)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?মডেল সংক্ষিপ্তসারটি দেখায় যে চিকিত্সাগুলির একটিতে অস্বাভাবিকভাবে বড় স্ট্যান্ডার্ড ত্রুটি রয়েছে, যা আপনি এখানে দেখতে পারেন:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 2.565 1.038 2.472 0.0134 *
treatmentcontrol -1.718 1.246 -1.378 0.1681
treatmentsucrose 16.863 2048.000 0.008 0.9934
treatmentprotein -1.718 1.246 -1.378 0.1681 আমি গ্লোমার কন্ট্রোল এবং অন্যান্য প্যাকেজগুলির ফাংশন থেকে বিভিন্ন অপটিমাইজার চেষ্টা করেছি এবং আমি একই রকম আউটপুট পাই। আমি এলোমেলো প্রভাবটিকে উপেক্ষা করে গ্ল্যাম ব্যবহার করে মডেলটি চালিয়েছি, এবং সমস্যাটি অব্যাহত রয়েছে। ডেটা অন্বেষণের সময় আমি বুঝতে পারি যে উচ্চমানের সাথে চিকিত্সা করা হয়েছে। ত্রুটির প্রতিক্রিয়া পরিবর্তনশীলটিতে কেবল সাফল্য রয়েছে। এই সমস্যাটি হতে পারে কিনা তা খতিয়ে দেখার জন্য আমি সেই চিকিত্সার জন্য "ব্যর্থতা" সহ একটি জাল ডেটা পয়েন্ট যুক্ত করেছি এবং মডেলটি সুচারুভাবে চলে এবং যুক্তিসঙ্গত মানক ত্রুটি দেয়। আপনি এখানে দেখতে পারেন:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.4090 1.6712 2.040 0.0414 *
treatmentcontrol -1.8405 1.4290 -1.288 0.1978
treatmentsucrose -0.2582 1.6263 -0.159 0.8738
treatmentprotein -2.6530 1.5904 -1.668 0.0953 .আমি ভাবছিলাম যে আমার অনুজ্ঞানটি যদি এই চিকিত্সাটির জন্য কোনও ভাল অনুমানের প্রতিরোধকারী ব্যর্থতার অভাব সম্পর্কে সঠিক হয় এবং আমি কীভাবে এই সমস্যাটি নিয়ে কাজ করতে পারি।
আগাম ধন্যবাদ!