같은 병원에서 일하고 있는 타과 선생님과 공동으로 실험한 결과과 나와서 우선은 그래프로 그려 보았는데...
몇개의 그래프를 보여 드리면 다음과 같습니다.
측정한 변수는 많은데 우선 대표적인 그래프 두개만 보여 드립니다, 보시다시피 두개의 약물(HPR/FUT)을 투여한 개(beagle) 8마리를 각 그룹당 4마리씩 할당 후 각종 혈액검사를 1시간 간격으로 5-6회 측정하였습니다.
그래프를 보시면 아시겠지만 기울기의 차이는 크게 없습니다.
여기서 질문을 드리면
1) 두 그룹의 데이터를 보면 산포도(spreadness)의 차이가 있는 것 같아, 변이계수(CV)를 이용해서
그룹간 차이가 있는지를 보려고 합니다. 즉 실험군이 대조군에 비해 분산이 적다. 즉 안정적인 값을 도출한다.라고 결론을 내려고 합니다.
이 방법이 가능할까요?
가능하다면 변이계수의 비교는 어떤 방법으로 하나요? t-test면 자유도가 3?,4?
2) mixed model을 적용해 볼까 합니다. 그룹당 개체수가 4라, 너무 작지 않을까 생각이 됩니다.
이 개체수로도 mixed model 적용이 가능할까요?
첫댓글 아래그래프는 시간을 연속형변수로 회귀분석같이 할수 있지만, 위의 그래프는 전혀 trend가 없어서 그냥 discrete하게 적용해야할듯합니다. 아래경우는 mixed model을 써도 될듯합니다만, 위에것은 시간을 discrete하게 하면 추정해야할 parameter수가 증가하여 될지 않될지는 실제로 fit해보야할듯합니다.
아래경우는 random intercept model을 쓰면 잘 맞을듯합니다. intercept의 분산이 같은지 검정할수 있을지는 모르겠습니다. 아마 intercept가 같은 모형과 다른 모형의 AIC, BIC등을 비교해야할듯합니다.
역시 쉬운것은 아니군요. ㅎㅎ
아 맨 마지막줄 수정...두 그룹의 random intercept "분산"이 같은 모형과 다른 모형의 AIC, BIC를 비교^^
어렵군요. 케이스 사이즈도 4개씩이니...