지난번에 질문을 드렸었는데, 난해한 부분이 생겨서 다시 질문 드립니다.
ANOVA 분석과 함께 Tukey 테스트를 하고 그래프를 그리려고 합니다.
Standard error 에 관한 질문인데요.
첨부한 파일에서 각 genotype 과 density 가 length 에 영향을 주었는지 보려고 합니다.
그런데 density 가 4 인 것은 총 4개에서 구한 평균 length 값이고,
8 인것은 총 8개에서 구한 평균 length 값 입니다. 16은 총 16에서 구했구요.
그렇다면 각각 수가 다른데요.
저는 R에서 group.SEs <- (group.SDs/sqrt(?)) 이렇게 standard error 를 구하는데...
이런 데이터 경우 어떻게 구하나요?
제 생각엔 4의 경우 group.SEs <- (group.SDs/sqrt(4))
8의 경우 group.SEs <- (group.SDs/sqrt(8)) 그리고 16의 경우 group.SEs <- (group.SDs/sqrt(16)) 같은데요.
이렇게 다를경우 R에서 어떻게 구하나요?
자동으로 인식해서 구해주는 함수식이 있나요?
그리고 이곳 웹에서도 첨부한 자료로 ANOVA 분석이랑 Tukey 그리고 그래프 가능한가요?
Comment 7
-
cardiomoon
2016.03.05 22:41
Recommend:1 Comment
-
R까기
2016.03.06 00:44
선생님, 제가 R 스튜디오랑 비교를 해보았는데 ANOVA 분석 값이 다르게 나옵니다.
Degrees of Freedom 이 잘못된 값으로 나오는데...어떻게 해야되는지 모르겠네요.
그리고 SE 값은 각각 density 4 8 16 마다 n 값이 다르므로 제가 질문드린 공식 group.SEs <- (group.SDs/sqrt(?)) 에 적용 할 수가 없네요.
그렇다면 R 스튜디오에서 이런 데이터가 있을경우 SE 값을 어떻게 구하고 그래프에 에러바를 만들 수 있나요?
바쁘신데 질문이 많아서 죄송합니다.
-
R까기
2016.03.06 00:19
자세히 설명주셔서 고맙습니다.
아직 사용해보지 않았는데 정말 간편해 보이네요.
좋은 자료 정말 감사드립니다.
-
cardiomoon
2016.03.06 08:26
ANOVA분석 값이 다르게 나오는 이유는
1) 완전교차되는 설명변수들에 A,B 두개 를 넣은 경우는 A+B+A:B(상호작용)이 같이 계산되어 그렇습니다.
2) 상호작용이 없는 모형을 원하시면 바로 위의 칸에 설명변수를 두개 넣으시면 됩니다. -> 같은 결과를 얻으실 겁니다.
결과확인후 download as pptx버튼을 누르시면 결과를 pptx로 얻으실 수 있습니다.
-
R까기
2016.03.06 10:38
R studio 에서도 * 즉 interaction 을 보고자 했는데 다른 결과가 나왔습니다.
genotype 4개 그리고 density 3개를 처리했으면 DF 값이 각각 (4-1), (3-1) 로 3과 2가 나와야 되는데 3과 1로 나오더라구요.
잘못된듯 합니다만...
-
cardiomoon
2016.03.06 08:27
에러바 있는 그래프는 자동으로 나옵니다.
-
R까기
2016.03.06 10:39
이상하네요. bar 그래프를 찾지 못하겠네요.
제가 지금 외국에 있는데요. 혹시 에러가 나는걸까요?
평균을 구할 변수에 length,그룹으로 나눌 변수에 genotype과 density를 입력하시면 다음의 결과를 얻습니다.
여기서 se 가 표준에러입니다. "전단계로 되돌리기"를 누르신 후 ANOVA는 다표본비교...분산분석에 가신 후 반응변수에 length, 완전교차되는 설명변수에 genotype과 density 를 넣고 분산분석 버튼을 누르면 분산분석과 TUKEY가 같이 진행됩니다. 아직 파워포인트로 다운로드되는 부분은 못 만들었는데 이번주 내에 만들겠습니다. ㅠ,ㅠ