抽样问题

newfinder

大侠们，我有一个有关bootstrap再抽样的问题想请教一下。

我现在有一组数，@array=(1,2,3..)有20个数，现在想用bootstrap重复1000次进行一个再抽样，产生一组新的20个数。
不知道该怎么写，

或者perl里有没有相关的模块是进行bootstrap再抽样的。

劳烦大家了

xiumu2280

1. Bootstrap（进化树）编辑
   
2. 本词条缺少概述、信息栏、名片图，补充相关内容使词条更完整，还能快速升级，赶紧来编辑吧！
   
3. 即自展值，是用来检验你所计算的进化树分支可信度的。简单地讲就是把序列的位点都重排，重排后的序列再用相同的办法构树，如果原来树的分枝在重排后构的树中也出现了，就给这个分枝打上一分，如果没出现就给0分，这样经过你给定的repetitions次（至少1000次）重排构树打分后，每个分枝就都得出分值，计算机会给你换算成bootstrap值。重排的序列有很多组合，值越小说明分枝的可信度越低，最好根据数据的情况选用不同的构树方法和模型。

复制代码

帮楼主补充下  百度百科


是这个意思？

newfinder

回复 2# xiumu2280

恩，意思是这样子的，但bootstrap是一个抽样方法，叫自助抽样方法。它也可以用于进化树的分析。
   

afukada

回复 1# newfinder


你是要boostrap 1000次

那會產生1000個包含20個element的array

所以你的題目到底是要什麼?

newfinder

回复 4# afukada
是根据原来的20个数，进行1000次重复抽样，从而产生一个经验分布，然后再从这个分布中随机抽取20个数

   

afukada

你的意思是:

假設有 @array=(1,2,3,4,5)

然後做bootstrap 5次

         位置  12345
可以得到 1. 23445
             2. 33412
             3. 23215
             4. 54112
             5. 14412

然後依據第一個位置上各數字的比例(1:1次, 2:2次, 3:1次, 4:0次, 5:1次)來抽在第一個位置上的數字

以此類推

是這樣嗎?

newfinder

回复 6# afukada
对，对，是的，是这样子的

   

afukada

1. @array=(1..20);
   
2. $rtimes=1000;
   
3. 
   
4. map{$list[$_][int(rand(scalar(@array)))]++}(0..$#array) for(0..$rtimes-1);
   
5. @acclist=map{$i=$_;[map{eval(join("+",@{$list[$i]}[0..$_]))}(0..$#{$list[$i]})]}(0..$#array);
   
6. @newarray=map{$array[$_]}map{$i=$_;$cannum=rand($rtimes);@sortlist=sort{$a<=>$b}(@{$acclist[$i]},$cannum);grep{$sortlist[$_]==$cannum}(0..$#sortlist)}(0..$#array);
   

复制代码

請自己好好想一下該怎麼做

newfinder

回复 8# afukada

好的，受教了大侠，以后好好学
   

pitonas

{:2_172:}

1. #!/usr/bin/perl
   
2. 
   
3. use 5.010;
   
4. 
   
5. my @old = ( 14, 30, 26, 26, 34, 18, 19, 20, 29, 21, 29, 24, 24, 25 );
   
6. 
   
7. sub bs(\@$) {
   
8.     my ( $array, $resample ) = @_;
   
9.     my $len = @$array;
   
10.     my @resample = map {
    
11.         my $sum;
    
12.         $sum += $array->[ rand($len) ] for 1 .. $len;
    
13.         int( $sum / $len )    # int or not int
    
14.     } 1 .. $resample;
    
15.     @resample[ map rand($resample), 1 .. $len ];
    
16. }
    
17. 
    
18. my @new = bs @old, 1000;

复制代码

old:
[ 14, 30, 26, 26, 34, 18, 19, 20, 29, 21, 29, 24, 24, 25 ]
new:
[ 24, 25, 27, 24, 24, 23, 23, 24, 25, 23, 26, 25, 24, 22 ]

回复 1# newfinder