R语言实现Xbar-R控制图
Xbar-R控制图在质量管理中主要用于对计量数据进行检测,以达到控制对象质量的目的。
虽然用Excel可以轻松实现控制图的操作,不过作为R软件初学者,我试着用仅有的一点R语言知识进行了控制图的绘制操作。虽然画出来了但是代码实在是低效,以后再慢慢改进吧。
> data <-
read.csv("kzt.csv") #打开数据
> D3 <-
0; D4 <- 2.114; A2 <-0.577
#三个系数,下文会用到
1.###计算各样本的极差
> r1 <-
c()
> for(i in 1:20) r1[i]
<- max(data[i,])-min(data[i,])
2.###计算极差上下界,并画出R控制图
> UCL_R
<- D4*mean(r1)
> CL_R <-
mean(r1)
> LCL_R
<- D3*mean(r1)
> plot(r1,
type="o",ylim=c(0,40),main="Range")
> abline(h=UCL_R,
lty='dashed')
> abline(h=CL_R)
> abline(h=LCL_R,
lty='dashed')
存在有出界值。利用which()函数确定出界位置(虽然能直接看出是样本7)
>
which(r1>UCL_R)
[1] 7
3.###把样本7从数据中删掉然后对样本进行重新编号
> data <-
data[-7,]
> rownames(data)
<- 1:nrow(data)
4.###重新计算各样本的极差
> r11 <-
c()
> for(i in 1:19) r11[i]
<- max(data[i,])-min(data[i,])
5.###重新计算极差上下界,并画出R控制图(代码如2,将r1换成r11即可)
此时R图判稳。接下来作Xbar图。
6.###计算各样本的均值
> m1 <-
c()
> for(i in 1:19) m1[i]
<- apply(data[i,],1,mean)
7.###计算均值上下界,并画出Xbar控制图(套路跟画极值控制图差不多)
> UCL_M
<- mean(m1) A2*mean(r11)
> CL_M <-
mean(m1)
> LCL_M
<- mean(m1)-A2*mean(r11)
>
plot(m1,type="o",ylim=c(60,90),main="Mean")
> abline(h=UCL_M,
lty='dashed')
> abline(h=LCL_M,
lty='dashed')
> abline(h=CL_M)
有出界值,找出出界值
> which(m1
[1] 13
8.###把样本13从数据中删掉然后对样本进行重新编号
> data <-
data[-13,]
> rownames(data)
<- 1:nrow(data)
9.###重新计算各样本均值、极差和均值上下界,并画出R控制图和Xbar控制图
此时Xbar与R图都判稳,生产过程的均值与变异度都处于稳态。延长统计过程状态下的Xbar-R图的控制限,即可进入控制用控制图阶段,实现对过程的日常控制。
最新文章
- AOP详解
- 【leetcode】Simplify Path
- 关于(object sender, EventArgs e)
- Ajax&;Java
- Android BLE 蓝牙低功耗教程,中央BluetoothGatt和周边BluetoothGattServer的实现
- Linux内核内存管理子系统分析【转】
- nginx+tomcat集群配置(3)---获取真实客户端IP
- 10个鲜为人知的WordPress函数
- CentOS 配置本地yum源
- 【Python千问 2】Python核心编程(第二版)-- 欢迎来到Python世界
- android之PackageManager简单介绍
- Codility上的问题 (17) PI 2012
- kafkaspout以及kafkabolt的最简实例
- Android可以换行的布局
- [js高手之路]html5 canvas动画教程 - 边界判断与小球粒子模拟喷泉,散弹效果
- (2017浙江省赛E)Seven Segment Display
- 10分钟理解JS引擎的执行机制
- df
- spring boot 开发非web应用程序
- TCP滑动窗口协议
热门文章
- 移动端网站如何开发(电脑端网站到手机端网站我们需要在html代码中添加哪个meta标签)
- Python学习笔记之常用函数及说明
- css背景属性整理
- Eureka自我保护机制、健康检查的作用、actuator模块监控
- ConnectionString连接字符串-密码丢失的解决方法
- Luogu P1092 虫食算(枚举+剪枝)
- Leetcode532.K-diff Pairs in an Array数组中的K-diff数对
- Promise的源码实现(符合Promise/A+规范)
- js作用域的销毁、不立即销毁、不销毁
- Ubuntu 16.04 配置 L2tp 客户端