我有一个简单的nx1整数数组,我想引导它来评估比例的置信区间。
我找到了IBM SPSS的解决方案,但是我想用matlab做这个分析,你可以在这里找到这个例子: http://publib.boulder.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.cs%2Fbootstrap_telco_frequencies_table.htm
在Matlab中的我有这些数据,它来自名为 c 的数组的列表:
Value | Count | Percent
1 300 2.99%
2 2928 29.16%
3 0 0.00%
4 3244 32.31%
5 0 0.00%
6 2589 25.78%
7 980 9.76%
我尝试使用BOOTFUN作为bootstrap置信区间,如下表达式:
n = histc(c,unique(c))/sum(n);
我的意思是 n 前一个数组的比例。
最后,我使用bootci函数来评估间隔:
ci= bootci(1000,n,c);
我知道我在bootfun的设置中错了,但我不知道如何解决它,因此我想得到你的帮助。我希望问题很清楚。
答案 0 :(得分:0)
你可以看看像这样的东西
>> proportions = @(x) histc(x, unique(x)) / length(x); # NB this is a function
>> x = randsample(1:6, 100, 1); # sample with replacement from 1:6
>> ci = bootci(10000, @proportions, x) # compute 95% confidence intervals
0.0600 0.1100 0.0400 0.1700 0.1300 0.1100
0.1800 0.2700 0.1600 0.3400 0.2900 0.2500
这里第一类的置信区间为(6%,18%),第二类的置信区间为(11%,27%)等。