R中的Sample函数

喜欢ヅ旅行 2022-06-13 03:24 254阅读 0赞

今天介绍一些运算函数，它们的使用很简单，没有什么难度，但是也会用的着。

在医学统计学或者流行病学里的现场调查、样本选择经常会提到一个词：随机抽样。随机抽样是为了保证各比较组之间均衡性的一个很重要的方法。那么今天介绍的第一个函数就是用于抽样的函数sample：

> x=1:10

> sample(x=x)

\[1\] 3 5 9 6 10 7 2 1 8 4

第一行代码表示给x向量x赋值，第二行代码表示对x向量进行随机抽样。结果输出为每次抽样抽得的结果，可以看出该抽样为无放回抽样\-最多抽n次，n为向量中元素的个数。

如果想指定在该向量中抽取元素的个数，需要加一个参数size：

> x=1:1000

> sample(x=x,size=20)

\[1\] 66 891 606 924 871 374 879 573 284 305 914 792 398 497 721 897 324 437

\[19\] 901 33

这是在1~1000的正整数中抽样，其中size指定抽样的次数，抽了20次，结果如上所示。

这些都是无放回抽样。所谓无放回抽样，也就是说某个元素一旦被选择，该总体中就不会再有该元素。如果是有放回抽样，则需添加一个参数repalce=T：

> x=1:10

> sample(x=x,size=5,replace=T)

\[1\] 4 7 2 4 8

“replace”就是重复的意思。即可以重复对元素进行抽样，也就是所谓的有放回抽样。我们看上面的结果，元素4在5次随机抽样的过程中被抽取了两次。

如果我们输入代码的位置与某个函数中参数的位置一一对应的话，我们可以不写该函数的参数，如：

> x=1:10

> sample(x,20,T)

\[1\] 1 2 2 1 5 5 5 9 9 5 2 9 8 3 4 8 8 8 1 1

在上述代码中我们省略了参数x，size和repalce，但是仍然可以运算并且表示对x向量有放回随机抽取20次。我们之所以尽量在每次编写代码时带上参数是因为我觉得这个习惯比较好，而且看起来清楚明白。另外，省略参数的前提是你非常熟悉某个函数参数的位置，否则一旦没有“对位”,那么结果肯定是错误的。而且很多函数有较多的参数，想记住它们的位置是困难的。而如果带上参数，那么即使位置不对应，也一样可以运算：

> x=1:10

> sample(size=20,replace=T,x=x)

\[1\] 4 9 2 6 4 5 4 7 10 5 2 2 3 4 2 4 6 8 7 8

这种优点显而易见，不仅清楚，而且无需对应。另外我们也可以看出，有放回抽样的话size可以无穷大，而无放回抽样size的大小就取决于总体的容量了。

对于掷骰子，投硬币（这可能是介绍抽样必介绍的内容），都属于有放回抽样。

这里要说明，对于sample函数，参数x可以是数值，也可以是字符，实际上参数x代表任意一个向量：

> a=c("A","B")

> sample(x=a,size=10,replace=T)

\[1\] "B" "A" "A" "A" "B" "A" "A" "B" "A" "A"

上述代码可以理解为掷硬币，抛了10次，其中正面（A）与反面（B）出现的次数是可以重复的。

上述抽样过程，每一个元素被抽取的概率相等，称为随机抽样。

有时候我们的抽取元素的概率未必相等（如常见的二项分布概率问题），此时我们需要添加一个参数prob，也就是“probability”（概率）的缩写。假设一名医生给患者做某手术成功的概率是80%，那么现在他给20例病人做手术，可能有哪几次是成功的呢？代码如下：

> x=c("S","F")

> sample(x,size=20,replace=T,prob=c(0.8,0.2))

\[1\] "F" "S" "S" "S" "S" "S" "S" "S" "S" "S" "S" "S" "F" "S" "S" "F" "S" "S"

\[19\] "F" "S"

其中“S”代表成功，“F”代表失败。

> x=c(1,3,5,7)

> sample(x,size=20,replace=T,prob=c(0.1,0.2,0.3,0.9))

\[1\] 3 5 7 3 7 3 7 5 3 7 7 7 1 5 7 5 7 7 3 7

这些代码告诉我们，对每一个元素都可以给定一个概率，且每个概率是独立的，即在参数prob中，不一定所有元素的概率加起来等于1，它只代表某元素被抽取的概率而已。

对于sample函数，其参数x可以是R中任意一个对象（如上述对字符的抽样）。其中还有一个功能相同的函数sample.int，“int”即“intger”的缩写，也就是“整数”。它的参数n必须是正整数：

> x=-10.5:7.5

> sample(x=x,size=3);sample.int(n=x,size=3)

\[1\] -5.5 -7.5 0.5

Error in sample.int(x, size = 3) : invalid first argument

第一行代码生成了\-10.5到7.5的等差数列，结果输出的第一行是sample的结果；第二行是sample.int的结果，提示错误：“第一个自变量无效”，因为它不是正整数。其余的用法与ample是一样的。

发表评论取消回复

表情：

评论列表（有 0 条评论，254人围观）

还没有评论，来说两句吧...

相关阅读

相关 r语言中paste函数_R中的paste（）函数-简要指南

r语言中paste函数 Using the paste() function in R will be straight and simple. In this tutori

左手的ㄟ右手/ 2023年02月28日 06:31/ 0 赞/ 118 阅读

相关【善用工具】googletest中的/samples

> 善用工具 > [http://blog.csdn.net/aganlengzi/article/details/64921496][http_blog.csdn.net

太过爱你忘了你带给我的痛/ 2022年09月30日 13:48/ 0 赞/ 210 阅读

相关 R语言中predict()函数

目录 predict()的用法函数形式参数介绍 predict.glm()的用法函数形式参数介绍输出介绍实例 predict.lm()的用法参

曾经终败给现在/ 2022年09月07日 06:24/ 1 赞/ 1494 阅读

相关 R语言中的curve函数

文章目录 curve函数简介实例绘制y=x的直线绘制 y = 2x + 3 的直线 y= 2x + 3 幂函

女爷i/ 2022年08月30日 01:57/ 0 赞/ 325 阅读

相关 A simple explanation of rejection sampling in R

(This article was first published on [theoretical ecology » Submitted to R-bloggers][th

悠悠/ 2022年08月05日 12:20/ 0 赞/ 165 阅读

相关 R中的Sample函数

今天介绍一些运算函数，它们的使用很简单，没有什么难度，但是也会用的着。在医学统计学或者流行病学里的现场调查、样本选择经常会提到一个词：随机抽样。随机抽样是为了保证各比较组之

喜欢ヅ旅行/ 2022年06月13日 03:24/ 0 赞/ 255 阅读

相关 maven 中settings.xml的配置(sample)

<?xml version="1.0" encoding="UTF-8"?> <settings xmlns="http://maven.apache.org/

谁践踏了优雅/ 2022年04月12日 12:54/ 0 赞/ 183 阅读

相关 print_r() 函数

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2xlbmd5

我不是女神ヾ/ 2022年02月16日 07:27/ 0 赞/ 225 阅读

相关 keras中sample_weight的使用

百度了好久，没有找到与sample\_weight相关的博客，于是自己摸索一下。 sample\_weight是keras中的fit的参数，中文文档介绍如下： ![在这

我会带着你远行/ 2021年11月01日 04:48/ 0 赞/ 532 阅读

相关 R语言中的quantile()函数

在R语言中取百分位比用quantile()函数，下面举几个简单的示例： 1、求某个百分位比 > data <- c(1,2,3,4,5,6,7,8,9,1

绝地灬酷狼/ 2021年09月15日 15:36/ 0 赞/ 1175 阅读