R语言在不同样本量下的Littles MCAR检验
时间:2022-07-23
本文章向大家介绍R语言在不同样本量下的Littles MCAR检验,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
原文链接:http://tecdat.cn/?p=10134
我进行一个小型仿真,以在不同样本量下测试Little的MCAR检验1。我可以研究线性回归中的异方差。我能够找到一些使用Little's MCAR检验的小样本研究人员的例子,因此我进行了仿真。
library(BaylorEdPsych)library(simglm)library(ggplot2)library(dplyr)library(mice)fixed <- ~1 + age + incomefixed_param <- c(2, 0.3, 1.3)cov_param <- list(dist_fun = c('rnorm', 'rnorm'), var_type = c("single", "single"), opts = list(list(mean = 0, sd = 4), list(mean = 0, sd = 3)))
ggplot(little.mcar.p, aes(x = n, y = p)) + geom_boxplot() + geom_crossbar(aes(ymin = q025, y = q05, ymax = q075), data = summarise( group_by(little.mcar.p, n), q025 = quantile(p, .025, na.rm = TRUE), q05 = quantile(p, .05, na.rm = TRUE), q075 = quantile(p, .075, na.rm = TRUE) )) + geom_hline(yintercept = .05) + scale_y_continuous(breaks = seq(0, 1, .05), limits = c(0, 1)) + labs(x = "Sample size", y = "p-value", title = "Little's MCAR test for data that are MCAR", subtitle = "2000 replications", caption = paste(paste("For the narrow boxes, going from top to bottom, lines", "represent 7.5th, 5th and 2.5th percentiles of p-values."), "Test maintains nominal error rate across wide range of sample sizes.", sep = "n"))
ggplot(little.mcar.p.mar, aes(x = n, y = p)) + geom_boxplot() + geom_crossbar(aes(ymin = q925, y = q95, ymax = q975), data = summarise( group_by(little.mcar.p.mar, n), q925 = quantile(p, .925, na.rm = TRUE), q95 = quantile(p, .95, na.rm = TRUE), q975 = quantile(p, .975, na.rm = TRUE) ), linetype = 2) + geom_hline(yintercept = .05) + scale_y_continuous(breaks = seq(0, 1, .05), limits = c(0, 1)) + labs(x = "Sample size", y = "p-value", title = "Little's MCAR test for data that are MAR", subtitle = "2000 replications", caption = paste(paste("For the dashed boxes, going from top to bottom, lines", "represent 97.5th, 95th and 92.5th percentiles of p-values."), "Test only maintains nominal error rate around sample size of 120.", sep = "n"))
回归接近完美(没有多重共线性)。
- Linux系统内存监控、性能诊断工具vmstat命令详解
- go http 服务器编程(2)
- 利用placeholder属性来添加输入框默认文字提示,提高用户体验
- Linux系统监控、诊断工具之top命令详解
- 【Dev Club分享】iOS黑客技术大揭秘
- Linux终端:用cat命令查看不可见字符
- golang 函数定义及其接口实例
- 分享两种圣诞节雪花特效JS代码(网站下雪效果)
- React 移动 web 极致优化
- golang 高效字符串拼接
- Linux+Nginx/Apache/Tomcat新增SSL证书,开启https访问教程
- golang 使用时间通过md5生成token
- golang中对map操作类
- Nginx在线服务状态下平滑升级或新增模块的详细操作记录
- 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 数组属性和方法
- Python自学成才之路 带有参数的装饰器
- Grafana 将默认的元数据库 sqlite 改为 mysql
- 使用IDEA整合spring4+spring mvc+hibernate
- springmvc中@PathVariable和@RequestParam的区别
- centos7 安装mysql5.6
- bashdb安装及调试shell脚本
- linux awk指令详解
- linux sed指令详解
- Grafana使用zabbix自定义模板
- linux shell之变量的使用规则
- 第11期:压缩表
- grafana-zabbix插件安装和配置zabbix mysql
- grafana使用教程之API key
- Grafana使用教程之安装
- Java基础数据类型之包装类equals和==详解