JMP数据分析论坛

标题: 一个计算单比例置信区间的例子 [打印本页]

---

作者: skyfree    时间: 2013-8-28 10:41
标题: 一个计算单比例置信区间的例子
下面一个计算单比例置信区间的例子:
实现效果：
[attach]568[/attach]

程序如下（参考附件）：数据与程序下载 http://www.kuaipan.cn/file/id_186398518334717989.htm

1. //*********************************************
   
2. // 此程序是用来计算单比例样本的比例估计和对应置信区间。计算基于正态分布，当样本小于30时，计算结果有偏差。
   
3. // 此程序仅用来展示JSL用法，不对算法进行深入研究，如需更精确计算，请参考相关公式。
   
4. // Skyfree
   
5. // 2013.08
   
6. //**********************************************
   
7. 
   
8. Clear Globals();
   
9. 
   
10. :: pConf = Function( {y = ., n = .},
    
11. 
    
12. confVal = Num( conf << get selected );
    
13. z = Normal Quantile( confVal + ((1 - confVal) / 2) );
    
14. perc = (y / n);
    
15. UCL = ((perc + (z ^ 2 / (2 * n))) + z * (Root( ((perc * (1 - perc) / n) + (z ^ 2) / (4 * n ^ 2)), 2 ))) / (1 + z ^ 2 / n);
    
16. LCL = ((perc + (z ^ 2 / (2 * n))) - z * (Root( ((perc * (1 - perc) / n) + (z ^ 2) / (4 * n ^ 2)), 2 ))) / (1 + z ^ 2 / n);
    
17. Show( UCL );
    
18. Show( LCL );
    
19. Show(perc);
    
20. 
    
21. If( Or( y > n, y < 0, n < 0 ),
    
22. UCLanswer << set text( "成功次数不能小于实验总次数/数值不能小于0. 请重试！" );
    
23. LCLanswer << set text( "" );
    
24. UCLanswer << font color( 3 );
    
25. UCLanswer << set font size( 10 );
    
26. ,
    
27. Meananswer << set text( "比例估计: " || Char( Round( perc, 5 ) ) );
    
28. Meananswer << font color( 0 );
    
29. Meananswer << set font size( 10 );
    
30. UCLanswer << set text( "置信区间上限UCL: " || Char( Round( UCL, 5 ) ) );
    
31. UCLanswer << font color( 0 );
    
32. UCLanswer << set font size( 10 );
    
33. LCLanswer << set text( "置信区间下限LCL: " || Char( Round( LCL, 5 ) ) );
    
34. LCLanswer << font color( 0 );
    
35. LCLanswer << set font size( 10 );
    
36. );
    
37. );
    
38. 
    
39. 
    
40. PropCalc = New Window ( "单比例置信区间",
    
41. Panel Box( "比例置信区间计算",
    
42. 
    
43. 
    
44. Lineup Box(
    
45. N Col( 1 ),
    
46. Text Box( "成功次数： " ), yIn = Text Edit Box( 1 ),
    
47. Text Box( "总实验次数： " ), nIn = Text Edit Box( 10 ),
    
48. Text Box( "置信水平： " ), Conf = Combo Box({"0.95", "0.99", "0.90"} ),
    
49. 
    
50. V List Box(
    
51. Border Box(
    
52. Top( 10 ),
    
53. Left( 10 ),
    
54. Bottom( 15 ),
    
55. sides( 143 ),
    
56. V List Box( UCLanswer = Text Box( "置信区间上限UCL: " ),Meananswer = Text Box( "比例估计: " ), LCLanswer = Text Box( "信区间下限LCL: " ),
    
57. Text Box( "注：当样本量小于30，计算结果将不精确 " ),)
    
58. 
    
59. ),
    
60. H List Box(
    
61. Button Box( "确认",
    
62. y = Num( yIn << get text );
    
63. n = Num( nIn << get text );
    
64. ::pConf( y = Num( yIn << get text ), n = Num( nIn << get text ) );
    
65. )
    
66. ),
    
67. Button Box( "退出",
    
68. PropCalc << CloseWindow;
    
69. Clear Globals();
    
70. )
    
71. )
    
72. )
    
73. )
    
74. 
    
75. );

复制代码

---

欢迎光临 JMP数据分析论坛 (http://www.jmpforum.net/)

X3