메뉴 건너뛰기

웹에서 하는 R 통계

ANCOVA에 대해서 여쭤보고 싶습니다.

2015.09.13 06:36

고든 조회 수:716



 안녕하세요. 

 

 다시 질문있어 교수님께 여쭤봅니다. 


 제가 하려는 것은 IFG2 세 군간의 차이를 변수를 보정해도 차이가 있다는 걸 ANCOVA로 입증하려는 것입니다. 


 > out2=lm(Ep~factor(IFG2)+factor(sex)+age+LVMI+factor(HTN)+bmi+heartrate,data=result)

> anova(out2)

Analysis of Variance Table


Response: Ep

                Df Sum Sq Mean Sq F value    Pr(>F)    

factor(IFG2)     2   3612  1805.9  645.23 < 2.2e-16 ***

factor(sex)      1   2927  2926.8 1045.72 < 2.2e-16 ***

age              1  23029 23029.0 8228.05 < 2.2e-16 ***

LVMI             1   1182  1181.5  422.16 < 2.2e-16 ***

factor(HTN)      1   1554  1554.2  555.32 < 2.2e-16 ***

bmi              1   3855  3854.9 1377.33 < 2.2e-16 ***

heartrate        1   2125  2124.9  759.19 < 2.2e-16 ***

Residuals    21337  59719     2.8                      

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1


 일단 여기까지는 잘 됩니다. 문제는 그 다음 Tukey나 bonferroni로 각 군간의 차이를 밝히는 부분입니다. 일주일 넘게 여러 방법으로 해봤는데 방법이 없습니다. 


 Tukey HSD  의 경우 범주형 변수로 하면 잘되지만 연속 변수를 넣으면 '다음과 같이 요인이 아닌 것들은 무시됩니다' 라는 메세지만 뜹니다. post hoc을 어떻게 할 수 있을까요. 이런 저런 패키지 다 써보고 나름 고민해도 안되서 여쭤 봅니다. 


 긴글 읽어주셔서 감사합니다. 

번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 575
380 rJava 에러 문제 해결에 대해 문의드립니다 [1] mdhjcho 2017.06.29 3444
379 예제 데이터 열 제목 편집방법 [2] 신돌석 2016.04.25 3172
378 생존 분석 시 Cox-regression vs. Log-rank [2] begabung 2017.05.22 2694
377 로지스틱회귀분석 in R [3] wjdEoddl 2018.01.19 2431
376 WebR관련 문의드립니다. [1] neurojang 2018.05.17 2100
375 안녕하세요? mytable에서 나온 결과를 워드 저장하는 방법에 대해 문의드립니다. [2] 박만영 2018.03.12 2018
374 특정 조건을 만족하는 ROW 삭제에 관한 질문입니다! [1] file ksyjudy 2017.07.30 1780
373 로지스틱 회귀분석에 대해 질문 올립니다. [12] SOL 2017.09.15 1661
372 error 문의드립니다. [1] chamchi 2017.07.19 1542
371 web-R에서 사용되는 명령어를 볼 수 있는 방법이 있을까요? [1] 통계장님 2018.03.15 1242
370 r studio server 가 오늘(3/15) 만료된다고 합니다 [2] 스튜던트 2018.03.15 1237
369 메타분석에 관한 문의 드립니다. [1] 봄날 2018.03.08 1214
368 mytable로 만든표를 워드파일로 변환 하는 방법 [3] 정암 2017.06.02 1203
367 R서버 설치 문의드립니다. [4] 화모기 2015.05.13 1176
366 data frame 문의드립니다. [1] hichoi 2018.04.25 1092
365 time ROC 에 대해 질문드립니다. file YH 2018.04.19 954
364 웹에서 하는 ggplot2에서.. [3] suein 2018.03.15 954
363 메타분석 후 결과 forest plot이 위쪽과 아래쪽이 다 출력되지 않는 것에 대한 질문입니다. ^^ [2] file hanclan 2018.08.09 924
362 csv파일 업로드 안되는 문제 [2] file 신돌석 2016.04.25 917
361 결측값 처리와 관련하여 [3] 화모기 2015.10.07 916