ช่วงความเชื่อมั่นสำหรับผลิตภัณฑ์ของสองพารามิเตอร์


11

ให้เราสมมติเรามีสองพารามิเตอร์และP_2นอกจากนี้เรายังมีตัวประมาณค่าความน่าจะเป็นสูงสุดสองตัวและและสองช่วงความมั่นใจสำหรับพารามิเตอร์เหล่านี้ มีวิธีสร้างช่วงความมั่นใจสำหรับหรือไม่p1p2p1^p2^p1p2

คำตอบ:


13

คุณสามารถใช้วิธีเดลต้าในการคำนวณผิดพลาดมาตรฐานของ{2}} วิธีการเดลต้าระบุว่าการประมาณความแปรปรวนของฟังก์ชันมอบให้โดย: การประมาณความคาดหวังของในอีกทางหนึ่งจะได้รับจาก: ดังนั้นความคาดหวังจึงเป็นเพียงฟังก์ชัน ฟังก์ชั่นของคุณคือ{2} ความคาดหวังของ น่าจะเป็น:p1^p2^g(t)

Var(g(t))i=1kgi(θ)2Var(ti)+2i>jgi(θ)gj(θ)Cov(ti,tj)
g(t)
E(g(t))g(θ)
g(t)g(p1,p2)=p1p2g(p1,p2)=p1p2p1p2{2} สำหรับความแปรปรวนเราต้องการอนุพันธ์บางส่วนของ : g(p1,p2)
p1g(p1p2)=p2p2g(p1p2)=p1

การใช้ฟังก์ชั่นสำหรับความแปรปรวนข้างต้นเราได้:

Var(p1^p2^)=p2^2Var(p1^)+p1^2Var(p2^)+2p1^p2^Cov(p1^,p2^)
ข้อผิดพลาดมาตรฐานก็จะเป็นราก sqare ของการแสดงออกดังกล่าวข้างต้น เมื่อคุณได้รับข้อผิดพลาดมาตรฐานมันจะส่งตรงไปยังการคำนวณช่วงความมั่นใจ 95% สำหรับ :p1^p2^p1^p2^±1.96SE^(p1^p2^)

ในการคำนวณข้อผิดพลาดมาตรฐานของคุณต้องมีความแปรปรวนของและซึ่งคุณสามารถรับได้ โดยVariance-covariance matrixซึ่งจะเป็นเมทริกซ์ 2x2 ในกรณีของคุณเพราะคุณมีค่าประมาณสองค่า องค์ประกอบในแนวทแยงในเมทริกซ์ความแปรปรวนร่วม - ความแปรปรวนร่วมคือความแปรปรวนของและในขณะที่องค์ประกอบนอกแนวทแยงมุมเป็นความแปรปรวนร่วมของและ (เมทริกซ์มีความสมมาตร) @gung กล่าวถึงในความคิดเห็นเมทริกซ์แปรปรวนร่วมสามารถแยกได้โดยโปรแกรมทางสถิติส่วนใหญ่ บางครั้งอัลกอริทึมการประมาณให้p1^p2^p1^p2^ Σ ^ p 1 ^ p 2 ^ p 1 ^ p 2 Σp1^p2^p1^p2^Hessian matrix (ฉันจะไม่ลงรายละเอียดเกี่ยวกับที่นี่) และเมทริกซ์ความแปรปรวนร่วมแปรปรวนสามารถประเมินได้โดยค่าผกผันของ Hessian เชิงลบ (แต่ถ้าคุณเพิ่มความน่าจะเป็นบันทึก!; ดูโพสต์นี้ ) ศึกษาเอกสารของซอฟต์แวร์เชิงสถิติและ / หรือเว็บของคุณเกี่ยวกับวิธีแยก Hessian และวิธีคำนวณค่าผกผันของเมทริกซ์

หรือคุณสามารถรับผลต่างของและจากช่วงความเชื่อมั่นในวิธีต่อไปนี้ (ใช้ได้กับ 95% -CI):{ต่ำกว่าขีด จำกัด สำหรับ -CI ข้อผิดพลาดมาตรฐานโดยประมาณคือ:โดยที่คือ quantile ของการแจกแจงแบบปกติมาตรฐาน (สำหรับ , ) จากนั้นp1^p2^SE(p1^)=(upper limitlower limit)/3.92100(1α)%SE(p1^)=(upper limitlower limit)/(2z1α/2)z1α/2(1α/2)α=0.05z0.9751.96Var(p1^)=SE(p1^)2. เช่นเดียวกับที่เป็นจริงสำหรับความแปรปรวนของ{2}} เราจำเป็นต้องมีความแปรปรวนร่วมของและด้วย (ดูย่อหน้าด้านบน) หากและเป็นอิสระความแปรปรวนร่วมเป็นศูนย์และเราสามารถวางคำศัพท์ได้p2^p1^p2^p1^p2^

บทความนี้อาจให้ข้อมูลเพิ่มเติม


4
+1 ความแปรปรวนของพารามิเตอร์และความแปรปรวนร่วมของพวกเขาสามารถพบได้โดยการตรวจสอบเมทริกซ์ความแปรปรวนร่วม - ความแปรปรวนร่วมของซึ่งซอฟต์แวร์ทางสถิติส่วนใหญ่สามารถให้ เช่นใน R คืออะไร? vcov ; และในเอสเอจะถูกเพิ่มเป็นตัวเลือกคำสั่งรูปแบบใน PROC REG βcovb
gung - Reinstate Monica

1
@gung ในเรื่องของ pedantry มันอาจคุ้มค่าที่จะชี้ให้เห็น (เพราะฉันรู้ว่ามันทำให้บางคนสับสน) ว่ามันเป็นเมทริกซ์ความแปรปรวนร่วม - ความแปรปรวนร่วมของมากกว่า (และอันที่จริงมันไม่ได้เป็นอย่างนั้นจริงๆ เนื่องจากค่าเบี่ยงเบนมาตรฐานจะต้องมีการประมาณจากกลุ่มตัวอย่างดังนั้นจึงเป็นเมทริกซ์ความแปรปรวนร่วม - ความแปรปรวนร่วมโดยประมาณ .. )บีตาβ^β
Silverfish

3
@Silverfish ลงโทษอย่างถูกต้อง ครั้งต่อไปฉันจะพูดว่า "เมทริกซ์ความแปรปรวนร่วม - ความแปรปรวนร่วมของ " β^
gung - Reinstate Monica

1
คุณสามารถลองสร้างฟังก์ชั่นความเป็นไปได้ของโพรไฟล์! และสร้างช่วงความมั่นใจจากนั้น
kjetil b halvorsen

ไม่เนื่องจากเป็นพารามิเตอร์ใช่หรือไม่ var(p1)=0
user0

1

ฉันพบสมการที่แตกต่างกันสำหรับการคำนวณความแปรปรวนของผลิตภัณฑ์

หาก x และ y กระจายอย่างอิสระความแปรปรวนของผลิตภัณฑ์ค่อนข้างตรงไปตรงมา: V (x * y) = V (y) * E (x) ^ 2 + V (x) * E (y) ^ 2 + V ( x) * V (y) ผลลัพธ์เหล่านี้ยังสรุปกรณีที่เกี่ยวข้องกับตัวแปรสามตัวหรือมากกว่า (Goodman 1960) ที่มา: การควบคุมยาฆ่าแมลง (1980), ภาคผนวก F

Coolserdash: องค์ประกอบสุดท้าย V (x) * V (y) หายไปในสมการของคุณ หนังสืออ้างอิง (สารควบคุมกำจัดศัตรูพืช) ผิดหรือเปล่า?

นอกจากนี้สมการทั้งสองอาจไม่สมบูรณ์แบบ " ... เราแสดงให้เห็นว่าการแจกแจงผลิตภัณฑ์ของตัวแปรอิสระสามตัวนั้นไม่ปกติ " ( แหล่งที่มา ) ฉันคาดหวังว่าจะมีความเบ้เป็นบวกแม้จะอยู่ในผลิตภัณฑ์ของตัวแปรที่กระจายตัวสองแบบ


0
  1. ความยาวของ CI / 2 / 1.96 = se คือข้อผิดพลาดมาตรฐานของ A หรือ B
  2. se ^ 2 = var คือความแปรปรวนของการประมาณ A หรือ B
  3. ใช้ A หรือ B โดยประมาณเป็นวิธี A หรือ B เช่น E (A) หรือ E (B)
  4. ติดตามหน้านี้http://falkenblog.blogspot.se/2008/07/formula-for-varxy.htmlเพื่อรับ var (A * B), เช่น var (C)
  5. สแควร์รูทของ var (C) คือ se ของ C
  6. (C - 1.96 * se (C), C + 1.96 * se (C)) คือ 95% CI ของ C

โปรดทราบว่าหาก A และ B ของคุณมีความสัมพันธ์กันคุณต้องพิจารณาความแปรปรวนร่วมของพวกเขาเช่นกัน

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.