ฉันจะประมาณช่วงความมั่นใจ 95% โดยใช้การทำโปรไฟล์สำหรับพารามิเตอร์ที่ประเมินได้โดยการเพิ่มฟังก์ชั่นบันทึกความเป็นไปได้โดยใช้ optimize ใน R ได้อย่างไร

ฉันจะประมาณช่วงความมั่นใจ 95% โดยใช้การทำโปรไฟล์สำหรับพารามิเตอร์ที่ประเมินได้โดยการเพิ่มฟังก์ชั่นบันทึกความเป็นไปได้โดยใช้ optimize ใน R ได้อย่างไร

ฉันรู้ว่าฉันสามารถประเมินเมทริกซ์ความแปรปรวนร่วมแบบไม่มีสัญญาณได้โดยการสลับเฮสเซียนแต่ฉันกังวลว่าข้อมูลของฉันไม่ตรงตามสมมติฐานที่จำเป็นสำหรับวิธีการนี้ที่จะถูกต้อง ฉันต้องการประเมินช่วงความมั่นใจโดยใช้วิธีอื่น

วิธีความน่าจะเป็นของโปรไฟล์มีความเหมาะสมตามที่กล่าวไว้ในStryhn และ Christensenและในหนังสือ MASS Venables และ Ripley ของ§8.4, pp. 220-221?

ถ้ามีมีแพ็คเกจใดบ้างที่สามารถช่วยฉันทำสิ่งนี้ใน R หรือไม่? ถ้าไม่เช่นนั้นรหัสหลอกสำหรับวิธีการดังกล่าวจะเป็นอย่างไร

mleฟังก์ชั่นแพคเกจจาก stats4 เป็นกระดาษห่อของoptimซึ่งทำให้มันค่อนข้างง่ายต่อการผลิตการคำนวณรายละเอียดความน่าจะเป็น ดูhelp("profile,mle-method", package = "stats4")ข้อมูลเพิ่มเติม

มีแพ็กเกจ ProfileLikelihoodหากคุณใช้ nlme ส่วนตัวแล้วฉันยังไม่ได้จัดการให้ใช้

การใช้แพ็คเกจ lme4a หรือ lmeEigen มีฟังก์ชั่น profile () ซึ่งมีจุดมุ่งหมายที่จะทำสิ่งที่คุณต้องการอย่างแท้จริง ลองใช้วิธีนี้เพื่อติดตั้งแพ็คเกจเหล่านี้:

install.packages("lme4a",repos="http://lme4.r-forge.r-project.org/repos") 

หรือไปที่เว็บไซต์เพื่อรับไฟล์ zip ในทำนองเดียวกันและน่าเสียดายที่ฉันไม่ได้จัดการให้ใช้ :) บางทีเราควรรอการอัพเดตของ lme4

วิธีการมีรายละเอียดในร่างของหนังสือของดักลาสเบตส์

แก้ไข: เจ๋ง! ฟังก์ชั่น profile () สำหรับรุ่น lmer มีให้บริการในรุ่นล่าสุดของ lme4 ที่จะติดตั้งโดยการพิมพ์:

install.packages("lme4",repos="http://r-forge.r-project.org")