오늘의 ADP 준비 / 2022-07-31 일(D56)

관리도 공부

22회에서 관리도 문제가 출제되었다고 하여 관리도 관련 유튜브 영상을 몇 개 공부했습니다.

관련 기출문제 : ADP 22회 2. 통계분석 (data: 금속 성분 함유량 데이터)
불량률 관리도에 따른 관리 중심선, 관리 상한선, 하한선 구하기 (각 공식 있음), 관리도 시각화

 

출처: 관리도 선택 https://youtu.be/2i38RRmwpVk

 

관리도 Control chart https://youtu.be/vsTcmfP3u1A

 

 

출처 : SPC 7강 관리도활용 https://youtu.be/50r4xeHkwoY

 

편상관계수

PACF, 상호작용 효과, 조절변수... 이런 개념들과 관련성이 있을 듯 해서 공부를 했습니다.

 

str(mtcars)
mtcars2 <- mtcars[,c("mpg","cyl","hp","wt")]
cor(mtcars2)
library(ggm)
pcor(c(1,3,2,4),cov(mtcars2))
pcor(c("mpg","hp","cyl","wt"),cov(mtcars2))

pcor.test(pcor(c(1,3,2,4),cov(mtcars2)),
          q=2,n=nrow(mtcars2))
library(ppcor)
pcor(mtcars2)
# 귀무가설 : 편상관계수는 0이다
pcor.test(mtcars2["mpg"],mtcars2["hp"],mtcars2[c("cyl","wt")])

출처 : 통계데이터분석 - 상관분석 - 편상관관계(https://youtu.be/LKK9PaSu9zU)

 

 

오타가 있어서 발생한 오류였는데, 원인을 파악하는데 시간이 좀 걸렸습니다.

그리고 아직은 mtcars2["mpg"] 표현도 익숙하지 않습니다.

> pcor.test(mtcars2["mpg"],mtcars2["hp"],mtcars2["cyl"],mtcars2["wt"]) 
Error in match.arg(method) :
   'arg'는 반드시 NULL 또는 문자형 벡터이어야 합니다

> class(mtcars2["mpg"])
[1] "data.frame"
> class(mtcars2$mpg)
[1] "numeric"

 

아직 좀 더 연구를 해볼 내용이 있습니다. 다음에 시간이 되면 살펴보려고 합니다.

• PACF 와 편상관계수 관계를 추가로 확인.
• 편상관계수를 엑셀로 구해보고 생각한 로직이 맞는지 확인.

 

scheffe 사후검정

https://it-freelancer.tistory.com/72