이 사이트에서 다중 회귀 분석을 하면
Relative Importance of Predictor Variables 막대그래프가 나오는데, 각 변수의 weights를 구해서 그래프를 그리는 것으로 추정됩니다.
1. 이그래프를 논문 result에 넣고 싶은데 그래프 도출(통계) 방법을 어떻게 기술해야 할지 고민입니다. (통계 방법이나 R코드가 어떻게 되는지 알고 싶습니다.)
(multiple logistic regression analysis 을 시행한 Nature에 실린 논문:
"We used AICc-based model averaging to assess the relative importance of predictor variables."로 기술하고 있고 그래프도 유사합니다.)
2. 이 그래프에서의 %값 = partial R제곱/total R제곱=각각 variable의 relative weight인가요?
좋은 사이트에서 많이 배워갑니다.
감사합니다^^ 늘 건강하세요 선생님:)
Comment 2
-
cardiomoon
2016.07.13 11:14
-
Ikarus86
2016.07.13 11:25
네 교수님 감사합니다^^
No. | Subject | Author | Date |
---|---|---|---|
1254 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1253 | 초초 초보자 문의 드립니다^^ [2] | 매직하마 | 2017.12.22 |
1252 | mytable() 함수에 대한 질문입니다. [3] | 브리디온 | 2017.11.19 |
1251 | 묻고답하기 게시판을 운영합니다 [1] | cardiomoon | 2015.04.20 |
1250 | ANOVA, Kruskal-Wallis 검정 및 다중비교 [3] | 소오름 | 2016.03.31 |
1249 | ios | bnmy6581 | 2018.12.28 |
1248 | 예제 데이터 열 제목 편집방법 [2] | 신돌석 | 2016.04.25 |
1247 | rstudio 사이트에 연결할 수 없음 문제 [1] | 지원 | 2019.09.10 |
1246 | 로지스틱 회귀분석에 대해 질문 올립니다. [12] | SOL | 2017.09.15 |
1245 | Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 | 김규 | 2017.09.30 |
1244 | rJava 에러 문제 해결에 대해 문의드립니다 [1] | mdhjcho | 2017.06.29 |
1243 | 3D 형태 그래프에 대해 질문드립니다. | 백지백지 | 2019.04.26 |
1242 | 메타분석의 변수지정 [1] | 야스오 | 2021.05.24 |
1241 | Logistic regression [3] | namuhanayeyo | 2017.03.02 |
1240 | 생존 분석 시 Cox-regression vs. Log-rank [2] | begabung | 2017.05.22 |
1239 | 표본설계시 subset 명령어 문의드립니다. [1] | 네잎클로바 | 2020.05.09 |
1238 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
1237 | Survdiff로 log-rank test 시 P-value 소수점 [3] | begabung | 2019.02.15 |
1236 | length 오류질문 [1] | 촙춉 | 2020.05.26 |
1235 | 메타분석-출판편향 [1] | 날아라 | 2019.08.26 |
이 부분은 제가 만든 코드는 아니고 Robert I. Kabacoff의 R in action(2nd edition) p209 - 211 에 있는 내용을 약간 modification 하여 만든 부분입니다. 다음은 이 부분에 대한 설명입니다.
A new method called relative weights shows significant promise. The method closely approximates the average increase in R-square obtained by adding a predictor variable across all possible submodels (Johnson, 2004; Johnson and Lebreton, 2004; LeBreton and Tonidandel, 2008). A function for generating relative weights is provided in the next listing.
relweights <- function(fit,...){ R <- cor(fit$model) nvar <- ncol(R) rxx <- R[2:nvar, 2:nvar]
rxy <- R[2:nvar, 1] svd <- eigen(rxx) evec <- svd$vectors ev <- svd$values delta <- diag(sqrt(ev)) lambda <- evec %*% delta %*% t(evec)
lambdasq <- lambda ^ 2
The code in listing 8.16 is adapted from an SPSS program generously provided by Dr. Johnson. See Johnson (2000, Multivariate Behavioral Research, 35, 1–19) for an explanation of how the relative weights are derived.
도움이 되시길 바랍니다.^^