x,y,z년도 40대가 각각 2373, 1859, 2521명, 50대가 각각 3065,2837,4044명이라고 가정했을때
x y z
40대 2373, 1859, 2521
50대 3065, 2837,4044
cochrane armitage test 시행할 경우 40대가 시간이 갈수록 유의하게 증가하는 트렌드를 보인다. (50대가 시간이 갈수록 유의하게 감소했다.) (p<0.05)
와 같이 해석해도 무리가 없나요?
또 궁금한 것은,
위의 예제를 백분율로 변경했을 경우
x y z
40대 43.64 39.59 38.40
50대 56.36 60.41 61.60
cochrane armitage test 시행후 p value가 다른데 n수로 트렌드분석할지 %로 트렌드분석할지는 연구자재량인지요
감사합니다.
R코드)
<N수>
test<-matrix(c(2373,1859,2521,3065,2837,4044),byrow=TRUE,nrow=2)
test
Desc(test)
CochranArmitageTest(test)
<%>
age<-matrix(c(43.64,39.59,38.40,56.36,60.41,61.6),byrow=TRUE,nrow=2)
age
Desc(age)
CochranArmitageTest(age)
Comment 1
-
cardiomoon
2018.04.02 23:25
| No. | Subject | Author | Date |
|---|---|---|---|
| 1194 | PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] | 긍정적인마음가짐 | 2023.05.06 |
| 1193 |
첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의
[3]  | 수담이 | 2023.05.04 |
| 1192 | 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] | mkw7654 | 2023.05.02 |
| 1191 |
메타분석 U3 그래프 관련 문의드립니다
 | 화이팅 | 2023.05.02 |
| 1190 |
교수님 안녕하세요~~
| 메시12 | 2023.04.29 |
| 1189 |
연속형 변수인데 범주형으로 인지됩니다.
[2] | shararaku | 2023.04.28 |
| 1188 |
통계마당
[2] | ejo | 2023.04.27 |
| 1187 |
moonBook 오류 관련 재문의드립니다.
[1] | holyhyun | 2023.04.24 |
| 1186 | moonBook 사용시 오류메시지로 문의드립니다. [1] | holyhyun | 2023.04.24 |
| 1185 |
에러 문의
[2] | 화이팅 | 2023.04.19 |
| 1184 |
Cox regression 에 대한 문의입니다.
[2] | saga | 2023.04.18 |
| 1183 |
PSM 문의
[3] | starwars | 2023.04.13 |
| 1182 | ggPS, ggPSMSummary function [6] | 로렌쵸메디치 | 2023.04.12 |
| 1181 | Survival 분석에 대해서 문의 드립니다. | 정훈신 | 2023.04.05 |
| 1180 | 안녕하세요. 분석 툴에 대해 문의드립니다. [1] | 스매직맨 | 2023.04.04 |
| 1179 | 드디어 publish 되었네요. [2] | 정지윤 | 2023.04.04 |
| 1178 | 연구동향분석 시 논문 데이터 정제 관련 [5] | vtlbtl | 2023.03.31 |
| 1177 |
안녕하세요 교수님 에러 문의드립니다
[2] | 소소 | 2023.03.26 |
| 1176 | 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 | 쥬쥬즈주 | 2023.03.22 |
| 1175 |
PSM 매칭 버튼 보이지 않음.
| swpapa | 2023.03.19 |
1) 통계를 돌리기 앞서 중요한 것은 정보의 전처리입니다. 나이는 연속형 변수인데 이것을 왜 인위적인 40대와 50대의 그룹변수로 바꾸어 처리하는지요?
차라리 x년도, y년도, z년도의 나이정보를 그대로 사용하여 평균이 같은지 다른지 보는 것이 낫지 않을까요?
40대, 50대로 보면 10의 자리에서 내림을 하여 41세와 49세를 같게 취급하고 49세와 50세는 다르게 취급하는 것인데 이는 중요한 정보를 일부러 손상시키는 것입니다.
2) prop.trend.test를 시행하는 방법은 다음과 같습니다.
> test<-matrix(c(2373,1859,2521,3065,2837,4044),byrow=TRUE,nrow=2)
> prop.trend.test(test[2,],colSums(test))
Chi-squared Test for Trend in Proportions
data: test[2, ] out of colSums(test) ,
using scores: 1 2 3
X-squared = 32.97, df = 1, p-value = 9.361e-09
여기서 prop.trend.test()의 도움말을 보면 다음과 같습니다.
Test for trend in proportions
Description
Performs chi-squared test for trend in proportions, i.e., a test asymptotically optimal for local alternatives where the log odds vary in proportion with
score. By default,scoreis chosen as the group numbers.Usage
Arguments
xNumber of events
nNumber of trials
scoreGroup score
즉 x는 number of event, n은 number of trials로 전체 숫자를 가지고 하는 것이 맞습니다.