r - R 或 Stata 中分类变量的自举比例

Question

我需要帮助在 R 或 Stata 软件中进行引导。我想计算说“是”和“不”的人的比例，例如政策的有效性

在Stata我有这个代码

bs "summarize y1" "r(mean)", reps(200) size(770)

r(mean)估计比例的值应该是多少？

另外，我在 R 中有这个代码：

test <- function (q13){
    test13 <- table(q13)
    rel_freq <- test13/sum(test13)
    return(rel_freq)
      }

results <- boot(data=q13, statistic=test,
                R=200)

如何更正代码？我收到错误

统计错误（数据，原始，...）：未使用的参数（原始）

Answer 1

在Stataproportion中，如果一个变量有两个以上的类别，您可以使用：

//样本数据

sysuse auto, clear
keep if (headroom==2.0 |headroom==2.5)
gen prop=.
replace prop=0 if headroom==2.0
replace prop=1 if headroom==2.5

//说0是是，1是不是

set seed 123
bootstrap _b, reps(100):proportion prop

根据@Nick 更新：对于二进制变量，以下内容就足够了

bootstrap r(mean), reps(100): summarize prop, meanonly

..................................................... ..................................................... ..................................................... .....................

在R中，您可以使用bootpackage 和 mtcars 数据执行以下操作：

library(boot)
set.seed(123)
x<-mtcars$vs
myprop<-function(x,i){
sum(x[i]==0)/length(x)
}

bootprop <- boot(x,myprop,100)