r studio：模拟我的代码1000次，选择p值

这是我的原始代码：

x = rbinom(1000,1,0.5)
z = log(1.3)*x       
pr = 1/(1+exp(-z)) 
y = rbinom(1000,1,pr)     
k=glm(y~x,family="binomial")$coef
t=exp(k)

---

我怎样才能模拟1000次，然后选择一个具有p值的呢？我不会为您这样做，但您可能需要执行以下步骤：

将代码编写为一个函数，返回您感兴趣的值 使用诸如replicate之类的工具多次运行此函数并记录所有答案 使用类似分位数的方法提取您感兴趣的百分比

---

这是tidyverse及其列表列的完美应用程序。请参阅内联注释中的解释

图书馆管理员 图书馆扫帚 为每个模拟创建一个带有id列的TIBLE，并将x包装在列表中 sim% 要生成z、pr、y、k，请使用purrr包中的map和map2在列表列x上循环 `~ ... ` 类似于'function.x{…}` `.x`表示正在映射上使用的变量 mutatez=mapx，~log1.3*.x， pr=mapz，~1/1+exp-.x， y=mappr，~rbinom1000,1,x， k=map2x，y，~glm.y~.x，family=二项式， 使用broom:：tidy以TIBLE的形式获取模型摘要 sum=mapk，扫帚：：整洁%>% 选择id和sum并取消对tibbles的测试 selectid，总和%>% unnestcols=csum%>% 以p<0.05的速度降低截距和每.x一次 过滤条款=拦截 p、 值<0.05 模拟 >A tibble:545x6 >id项估计标准误差统计p值 > >1 3.x 0.301 0.127 2.37 0.0176 >27.x 0.263 0.127 2.06 0.0392 >3.8.x 0.293 0.127 2.31 0.0211 >4 11.x 0.377 0.128 2.96 0.00312 >512.x0.2650.1272.080.0373 >613.x0.3660.1272.880.00403 >7 14.x 0.461 0.128 3.61 0.000305 >8.17.x 0.274 0.127 2.16 0.0309 >918.x0.3940.1273.090.00200 >1019.x0.3710.1272.920.00354 >…还有535行

---

由v0.3.0于2020年5月18日创建的新问题：现在我提取所有p@misayali你忘记接受这个答案了吗？若要将答案标记为已接受，请单击答案旁边的复选标记，将其从灰色切换为已填写。@RyanSYKwan我刚刚完成了！谢谢你的帮助@RyanSYKwan在这里还有一个小问题：我使用broom:：tidy以TIBLE的形式获取模型摘要，但是如果我想打印出在TIBLE函数中创建的原始矩阵，比如只包含x和y，我应该怎么做？首先对z=log1.3*x进行日志转换，然后再将t=expk反向，这似乎是多余的。而我们可以z=1.3*x，直接与coef进行比较。