R语言多臂试验 - 我们应该考虑多重性吗?
时间:2022-07-22
本文章向大家介绍R语言多臂试验 - 我们应该考虑多重性吗?,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
原文链接:http://tecdat.cn/?p=6408
我们应该担心多重性吗?
事实证明(或许不足为奇)在多臂试验的设定中,某种多样性调整是否合适的问题一直备受争议 。一种观点认为,如果不同的假设代表不同的研究问题,那么不允许进行多重比较是合理的。
模拟研究
为了进行实证研究,我根据Parmar 等人的假设进行了一项非常小的模拟研究(在R中)。
运行模拟
mu < - c(0,0)
tausq < - 0
在顶部,这将所有真正的治疗效果设置为零。这使我们能够检查类型1错误。
贝叶斯方法
遵循安德鲁·格尔曼与之前相关的论文,似乎处理前面问题的方法是贝叶斯分析。在这种情况下,每个治疗效果的后验平均值将朝着估计的总体平均治疗效果缩小,取决于真实治疗效果的估计方差和每个估计效果的精确度。
用于仿真
如果有人发现错误,请在评论中告诉我。
nSims < - 1000
#specify真实效果的均值和方差
mu < - c(-1,-1)
tausq < - 2
#mu < - c(0,0)
#tausq < - 0
specify真正治疗之间的相关性rho
rho < - 0
trueCov < - c(tausq,rho * tausq,rho * tausq,tausq),nrow = 2)
由于共同的控制臂,误差相关性为0.5
errorRho < - 0.5
errorVariance < - 1
errorCov < matrix(c(errorVariance,errorRho * errorVariance,errorRho * errorVariance,errorVariance),nrow = 2)
sigResult < - array(0,dim = c(nSims,2))
trueEffects < array(0,dim = c(nSims,2))
estEffects < - array(0,dim = c(nSims,2))
for(i in 1:nSims){
#generate真正的治疗效果
trueEffects [i,] < mvrnorm(n = 1,mu = mu,Sigma = trueCov)
estEffects [i,] < - trueEffects [i,] + mvrnorm(n = 1,mu = c(0,0),Sigma = errorCov)
testStat < - estEffects [i,] /(errorVariance ^ 0.5)
p_value < - 2 * pnorm(abs(testStat),lower.tail = FALSE)
sigResult [i,] < 1 *(p_value <0.05)
}
proportions of trials中找到每个治疗组
colMeans(sigResult)
beneficial effect与对照相比 #is假设两个治疗组相同(对照组)
bestTrt < - array(0,dim = c(nSims,1))
bestEstEff < - array(0,dim = c(nSims,1))
bestEstMinusTrue < - array(0,dim = c(nSims,1))
for(i in 1:nSims){
bestTrt [i] < - which.min(estEffects [i,])
#find对应的真实效果
bestTrtTrueEff [i] < - trueEffects [i,bestTrt [i]]
bestEstEff [i] < - estEffects [i,bestTrt [i]]
ciCov [i] < - 1 *(((bestEstEff [i] -1.96 * errorVariance ^ 0.5)<bestTrtTrueEff [i])&((bestEstEff [i] + 1.96 * errorVariance ^ 0.5)> bestTrtTrueEff [i]))
}
- 《企鹅电竞weex实践—— iOS SDK的小九九》
- Linux操作系统DNS解析(nameserver)监控脚本
- 不给“爸爸”添麻烦 - iTOP iOS 动态库改造
- 移动SEO分享:php自动提交复合型Sitemap到百度搜索
- 《Android外部存储》
- Android JNI出坑指南
- 《iPhone X ARKit Face Tracking》
- 结合标签广告,定制一个QQ邮箱订阅
- SecureCRT全局发送相同命令,快速抓取服务器信息的方法
- [不定期更新]简单的shell脚本练习实例
- 超简单的MySQL主从复制配置步骤
- 解决Centos下vsftp无法上传文件的问题,附vsftp配置详解
- 为iFrame添加动态载入效果,提高用户体验
- 分享超炫的表白页面和爱的纪念日源码
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 如何处理WordPress网站404状态死链
- (三)RecyclerView简单滑动场景分析
- AWS EC2实例开启IPv6方法 阅读模式
- centos7宝塔环境安装VeryNginx waf系统实现高级防火墙功能
- 腾讯云中关于授权子用户QCloudResourceFullAccess权限后使用api接口创建购买cvm没有支付权限的解决办法
- centos7中systemctl系统服务和资源限制故障
- centos7 安装zabbix并监控windows服务器
- 宝塔面板插件-shadowsocks可视化管理
- Android面试知识笔记:那些年面试官常问的知识点...(附详细解析)
- 聊聊dubbo-go的failbackCluster
- Docker环境准备
- js实现HTML页面时钟动态显示
- WordPress标题的分隔符-被转义为–的解决办法
- WordPress网站实现下雪特效,云体验下雪的赶脚
- 小知识:后台执行Oracle创建索引免受会话中断影响