ข้อผิดพลาด LME () - ถึงขีด จำกัด การวนซ้ำ

ในการระบุรูปแบบเอฟเฟกต์แบบผสมข้ามฉันกำลังพยายามรวมการโต้ตอบ อย่างไรก็ตามฉันได้รับข้อความแสดงข้อผิดพลาดต่อไปนี้:

Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : 
nlminb problem, convergence error code = 1
message = iteration limit reached without convergence (10)

แบบจำลองมีดังต่อไปนี้: 1. หัวฉีด 3 แบบ (เอฟเฟกต์คงที่) 2. ผู้ปฏิบัติงาน 5 คนแต่ละคนมีการวัดซ้ำ 3 ครั้งสำหรับการไหลของเชื้อเพลิงจากหัวฉีดทั้ง 3 แบบ

ฉันถูกขอให้รวมการทำงานร่วมกันระหว่างประเภทของหัวฉีดและตัวดำเนินการในโมเดล นี่คือรหัสของฉันสำหรับรูปแบบ:

flow.lme <- lme(rate ~ nozzle, error= nozzle|operator, data=Flow)

ทำไมฉันถึงได้รับข้อความแสดงข้อผิดพลาดนี้?

mixed-model lme4-nlme

— ทัลบาชาน
แหล่งที่มา

คุณไม่ต้องการoperator|nozzleสุ่ม

— Olivia Grigg

ไม่โอเปอเรเตอร์คือเอฟเฟกต์แบบสุ่ม

— f1r3br4nd

คุณสามารถใช้> crtl = lmeControl (opt = 'optim', optimMethod = "SANN")

— AliReza Afshari Safavi

คุณควรตอบกลับของ @ f1r3br4nd เป็นคำตอบ

— JetLag

@AliRezaAfshariSafavi ประโยชน์ของการใช้ "SANN" กับ BFGS เริ่มต้นคืออะไร

— gcamargo

คำตอบ:

ฉันไม่เคยได้ยินerrorเรื่องโต้แย้งlmeและฉันไม่เห็นมันในเอกสารประกอบ คุณแน่ใจหรือว่าไม่ใช่การพิมพ์ผิด แต่เพื่อตอบคำถามที่คุณถาม:

ลอง ?lmeControl

การตั้งค่าmaxIter, msMaxIter, niterEMและ / หรือmsMaxEvalข้อโต้แย้งที่จะสูงกว่าค่าเริ่มต้นอาจแก้ไขปัญหานี้ จับออกจากlmeControlวัตถุแล้วผ่านวัตถุที่กับข้อโต้แย้งของcontrollme

หรือ...

เครื่องมือเพิ่มประสิทธิภาพเริ่มต้นใหม่lmeใช้เป็นขุย ครึ่งเวลาที่ปัญหาเหล่านี้ได้รับการแก้ไขสำหรับฉันเมื่อฉันเปลี่ยนกลับไปเป็นเครื่องมือเพิ่มประสิทธิภาพแบบเก่า คุณทำเช่นนี้โดยการตั้งค่าoptอาร์กิวเมนต์สำหรับการที่จะlmeControl'optim'

ดังนั้นการรวมเข้าด้วยกัน:

ctrl <- lmeControl(opt='optim');
flow.lme <- lme(rate ~ nozzle, error= nozzle|operator, control=ctrl, data=Flow);

— f1r3br4nd
แหล่งที่มา

ในบางกรณีมันอาจคุ้มค่าที่จะรู้ว่าlmeControlมันเป็นฟังก์ชั่นจากnlmeแพ็คเกจ

— Qaswed

อย่างแรกนี่คือโมเดล ANOVA ไม่ใช่โมเดลผสม

ประการที่สองดูเหมือนว่าฉันไม่ได้ระบุรุ่นของคุณ ในรูปแบบสมการคุณมี

{คำตอบ}_{ผม J} = β_{1} {ประเภทหัวฉีด}_{1 ผม J} + β_{2} {ประเภทหัวฉีด}_{2 ผม J} + β_{3} {ประเภทหัวฉีด}_{3 ผม J} + {ผู้ประกอบการ}_{ผม} + {หัวฉีดภายในผู้ประกอบการ}_{ผม J}

$\mbox{response}_{ij} = \beta_1 \mbox{nozzle type}_{1ij} + \beta_2 \mbox{nozzle type}_{2ij} + \beta_3 \mbox{nozzle type}_{3ij} + \mbox{operator}_i + \mbox{nozzle within operator}_{ij}$

คำสุดท้ายมี 15 ค่าแยกต่างหากสำหรับการสังเกต 15 ครั้งที่คุณมี ไม่มีองศาอิสระเหลือสำหรับรับเงื่อนไขอื่น ๆ ในโมเดล รวมถึงการโต้ตอบก็เป็นคำแนะนำที่ไม่ดี คุณจะต้องวางพวกเขาใด ๆ แม้การรวมไว้เป็นเอฟเฟ็กต์แบบกากบาทก็ไม่ได้ช่วยอะไรเพราะมันจะมีสีที่สมบูรณ์แบบกับเอฟเฟกต์คงที่และจะไม่สามารถประเมินได้ ความน่าจะเป็นสูงสุดหรือโมเดล REML ที่มีการสำรวจ 15 ครั้งนั้นไม่สมเหตุสมผล ผลลัพธ์แบบอะซิมโทติคของทฤษฎีความน่าจะเป็นสูงสุดจะไม่ทำงาน: นี่คือเฟอร์รารี่ที่คุณพยายามจะขับบนสนามที่ไถ

— StasK
แหล่งที่มา

หากมีทั้งเอฟเฟกต์แบบสุ่มและแบบคงที่ในโมเดลจากนั้นนิยามเป็นโมเดลเอฟเฟกต์ผสม ไม่ว่าคุณจะเรียกมันว่า ANOVA หรือการถดถอยเป็นปัญหาที่แยกต่างหากและเรียงลำดับคำถามเชิงความหมาย แม้ว่าฉันจะสับสนเล็กน้อย แต่สิ่งที่ OP หมายถึงโดยการโต้ตอบ เท่าที่ผมสามารถบอกได้ว่าเขาทำอยู่แล้วว่าโดยใช้แทนrandom=~nozzle|operator random=~1|operator

— f1r3br4nd

วรรณกรรมบางเล่มอ้างถึงเอฟเฟกต์แบบซ้อนหลายระดับเป็นการโต้ตอบระหว่างการทำรังในระดับต่างๆ ฉันคิดว่าฉันเคยเห็นสิ่งนี้ใน Pinheiro & Bates ฉันเห็นด้วยว่าการบอกสิ่งนี้อย่างถูกต้องเป็นเรื่องของความหมาย แต่ฉันแค่คิดว่าจะแนะนำthis-does-not-have-to-be-a-mixed-modelแท็ก เมื่อประมาณสองในสามของmixed-modelsคำถามที่ฉันได้เห็นการพูดอะไรบางอย่างกับเอฟเฟกต์นั้นเป็นส่วนหนึ่งของคำตอบของฉัน

— StasK

ตลกฉันใช้เวลาไปกับการบอกคนอื่นว่าพวกเขาไม่ได้ใช้โมเดลผสมกัน ที่จริงฉันอยากจะผิดเพราะมันจะทำให้ชีวิตของฉันง่ายขึ้นบ้าง สิ่งที่คุณจะบอก OP กฎของหัวแม่มือคือการพิจารณาเมื่อมีความต้องการรูปแบบผสม?

— f1r3br4nd

โอ้คุณเป็นคนเลวแล้ว อันนี้มีตัวทำนายหมวดหมู่เดียวดังนั้นมันจึงเป็นแบบจำลอง ANOVA สำหรับฉันตามที่ฉันพูดไว้ก่อนหน้านี้ หากคุณมีข้อมูลในระดับที่แตกต่างกัน (เช่นรัฐ \ โรงเรียน \ นักเรียนที่มีข้อมูลเกี่ยวกับรัฐโรงเรียนและนักเรียน) นั่นจะฟังดูเป็นแบบผสมกับฉันมากกว่า โดยทั่วไปถ้าคุณสามารถทำเช่นนี้เป็นผลบวกกำลังสองนั่นก็คือ ANOVA หากคุณสามารถทำสิ่งนี้เป็นรูปแบบการถดถอยนั่นคือรูปแบบการถดถอย ถ้าการทำโอกาสสูงสุด / REML นั้นหลีกเลี่ยงไม่ได้อย่างแน่นอน (เพราะมันอยู่ในกรณีตอบกลับแบบไบนารี่) นั่นคือโมเดลผสมที่เหมาะสมสำหรับฉัน

— StasK