이 사이트에서 다중 회귀 분석을 하면
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 | 생존분석 [1] | toby | 2015.04.29 |
1253 | RStudio v0.99 Preview: Code Diagnostics [1] | fermat39 | 2015.04.30 |
1252 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
1251 | need finite 'xlim' values [3] | anthleemd | 2015.05.02 |
1250 | 샤이니 앱에서 SQL 명령어 사용 방법 [2] | 김브링이 | 2015.05.06 |
1249 | 표만들기 관련 [2] | 현승협 | 2015.05.07 |
1248 | MS word문서 저장관련 [1] | 아자 | 2015.05.09 |
1247 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1246 | 웹에서 하는 R통계 - R코드 [2] | 김영진 | 2015.06.07 |
1245 | Plot inset 만들기 | cardiomoon | 2015.06.09 |
1244 | 오류 질문 - 표작성시 [1] | 미소천사 | 2015.06.15 |
1243 | 오류 질문 - PDF 변환 안되는 file입니다. [1] | 미소천사 | 2015.06.15 |
1242 | print ztable [1] | 김영진 | 2015.06.17 |
1241 | 표 작성시 질문 [1] | simsim | 2015.06.17 |
1240 | 나무 분석 질문입니다. | 미소천사 | 2015.06.19 |
1239 | Interactive ggplot2 [2] | cardiomoon | 2015.06.22 |
1238 | Rstudio server 서버 오류 관련해서 문의드립니다. [2] | 노상래 | 2015.07.22 |
1237 | Rstudio server 무한 로딩 관련해서 문의 드립니다. [6] | 노상래 | 2015.07.24 |
1236 | http://www.r-statistics.com/ 번역소식 [1] | fermat39 | 2015.07.25 |
1235 | 생존분석 변수선택 [3] | nash | 2015.08.11 |
이 부분은 제가 만든 코드는 아니고 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.
도움이 되시길 바랍니다.^^