perl 不放回抽样

好吧,上次由于交代不是很清楚,搞了个抽样放回的程序,今天重新弄了一下,把它写成模块的形式:

package Sampling;
use warnings;
use strict;

require Exporter;
our @ISA = qw(Exporter);
our @EXPORT_OK = qw(swor fac); #swor -- sampling without replacement

sub swor{
	my ($sample_num, $sampling_num) = @_;
	my @array = (1..$sample_num);
	my @percom;
	for(my $i = 0; $i < $sampling_num; $i ++)
	{
		my $index = int(rand(@array));
		push @percom, $array[$index];
		my @tmp;
		if($index eq 0)
		{
			@tmp = @array[($index + 1)..$#array];
		}elsif($index eq $#array){
			@tmp = @array[0..($index - 1)];
		}else{
			@tmp = (@array[0..($index - 1)], @array[($index + 1)..$#array]);
		}
		@array = ();
		@array = @tmp;
	}
	return @percom;
}

sub fac{
	my $value = shift;
	if($value eq 1)
	{
		return 1;
	}elsif($value eq 0){
		return 1;
	}else{
		return $value * &fac($value - 1);
	}
}

文件名保存与package相同的名字,后缀为.pm,主程序调用就可以了。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值