如何解决如何使用for循环获取多个样本?
设置为:“假设您有10,000人口,而6,000人口是民主党人。我们从该假设人口中模拟调查抽样。首先,在R中生成人口向量。然后,创建500个随机样本大小为n = 1,000的样本。“
到目前为止,我的代码是:
pop<-c(rep("Democrat",6000),rep("Republican",4000))
nTrials <- 500
n <- 1000
results <- rep(NA,nTrials)
for(i in 1:nTrials)
{
sampled <- sample(x=pop,size=n,replace=FALSE)
results[i]<- sampledpop<-c(rep(1,rep(0,4000))
nTrials <- 500
n <- 1000
results <- matrix(data=NA,ncol = nTrials,nrow = n)
Y<-matrix(data=NA,ncol=nTrials,nrow=1)
for(i in 1:nTrials)
{
sampled <- sample(x=pop,replace=TRUE)
results[,i]<- sampled
Y[,i]<- sum(results[,i])
}
我认为这段代码有效,但是我担心如何判断矩阵是否正确填充。
解决方法
您可以使用查看功能轻松查看已保存的对象。 Link to View Function in R。
我们还可以在R代码中添加一行代码,直到按键后才停止执行。 Stack exchange thread covering this
将两者放在一起,然后我们可以将两行放入循环中,其中一行向我们显示最终输出的当前版本,另一行将循环暂停直到我们继续。这将使我们逐步探索循环的行为。以您的一个循环为例:
for(i in 1:nTrials)
{
sampled <- sample(x=pop,size=n,replace=TRUE)
results[,i]<- sampled
Y[,i]<- sum(results[,i])
View(Y)
readline(prompt="Press [enter] to continue")
}
请记住,这将继续进行指定次数的试验。
您可以限制试验的次数,但是您不能确保获得相同的结果,因此我们可以在代码中加一个中断。 link to info about the break statement。假设我们对事情的发展感到满意,这可以让我们及早中断for循环。为使中断真正发挥作用,让它与一些用户输入配对,您可以选择是否继续。 link for collecting user input in r
因此将所有这些内容进行梳理,我们得到的结果如下:
for(i in 1:nTrials)
{
sampled <- sample(x=pop,i])
View(Y,)
interrupt = readline(prompt="Enter 1 for next loop,0 to exit: ")
if (interrupt == 0) {break}
}
对于我来说,到目前为止,您的代码对我来说还算不错。
,尝试一下
replicate(500L,sample(c("Democrat","Republican"),1000L,replace = TRUE,prob = c(0.6,0.4)),simplify = FALSE)
或者这个
pop <- c(rep("Democrat",6000L),rep("Republican",4000L))
replicate(500L,sample(pop,1000L),simplify = FALSE)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。