선생님, repeated measures anova를 링크로 알려주신 바 대로 했는데요... 시간에 따른 그래프만 하나 그려주고 p값 등은 나오지 않네요.
error값으로 null value가 있다고 빨간 글씨로 뜨구요. 빈 데이타가 하나라도 있으면 안되는지요?
동시에 여러개 반응 변수를 넣으니 하나로 연결되서 그려지네요..
ID는 No
반응 변수는 HR, SBP, BNP 등이며...
그룹간 변수는 T group이며
그룹내 변수는 phase입니다.
제가 뭘 잘못 했을까요? 화일 첨부드립니다.
그리고, 이런 데이타에서도 multivariate analysis도 가능한지요? 시간에 따른 차이가 두 T group에서 차이가 있는지 (다른 변수의 영향을 배제하고) 분석가능할지요?
제가 통계 초보라... 질문만 많이 드린 것 같아 죄송합니다.
선생님, 미리 감사드립니다. 은인이십니다.
Comment 2
-
cardiomoon
2018.01.23 14:05
일단 BNP는 T.group간의 차이가 있어 보입니다. 다음을 참조하세요.잘못하신 것은 아니고 누락치가 있는 경우 통계 가설 검정을 할때 에러가 나네요. 프로그램을 수정해야겠습니다. 시간이 조금 걸릴 것 같습니다.> data=read.csv("./Downloads/transfusion_data.csv",stringsAsFactors = FALSE)> out<-aov(HR~T.group*Phase.no+Error(No/(Phase.no)),data=data)> summary(out)Error: NoDf Sum Sq Mean SqT.group 1 532.8 532.8Error: No:Phase.noDf Sum Sq Mean SqT.group 1 34.27 34.27Error: WithinDf Sum Sq Mean Sq F value Pr(>F)T.group 1 479 478.7 2.012 0.159Phase.no 1 85 85.2 0.358 0.551T.group:Phase.no 1 276 276.2 1.161 0.284Residuals 105 24977 237.9> out1<-aov(sBP~T.group*Phase.no+Error(No/(Phase.no)),data=data)> summary(out1)Error: NoDf Sum Sq Mean SqT.group 1 75.44 75.44Error: No:Phase.noDf Sum Sq Mean SqT.group 1 5.747 5.747Error: WithinDf Sum Sq Mean Sq F value Pr(>F)T.group 1 59 58.56 0.357 0.552Phase.no 1 91 90.51 0.551 0.459T.group:Phase.no 1 22 21.71 0.132 0.717Residuals 110 18059 164.18> out2<-aov(BNP~T.group*Phase.no+Error(No/(Phase.no)),data=data)> summary(out2)Error: NoDf Sum Sq Mean SqT.group 1 213.2 213.2Error: No:Phase.noDf Sum Sq Mean SqT.group 1 903.6 903.6Error: WithinDf Sum Sq Mean Sq F value Pr(>F)T.group 1 16173 16173 3.969 0.049 *Phase.no 1 9 9 0.002 0.963T.group:Phase.no 1 0 0 0.000 0.999Residuals 103 419716 4075---Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 -
merr
2018.01.23 16:06
감사합니다.
