이 사이트에서 다중 회귀 분석을 하면
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 |
---|---|---|---|
214 | 웹에서 하는 R통계 주성분 분석 function error [1] | 김재성 | 2016.08.18 |
213 | 이전에 생존분석 PPT다운 안되어 글쓴이 입니다. [2] | 따거81 | 2016.08.17 |
212 | R 관련 질문입니다. [7] | kenny | 2016.08.16 |
211 | 생존분석 결과다운로드 [1] | 파랑파랑하늘 | 2016.08.16 |
210 | 인공지능 IBM Watson의 API를 연결할 수 있는 cognizeR [1] | Ben | 2016.08.16 |
209 | 무엇이 문제일까요? [1] | 김현정 | 2016.08.09 |
208 | 교수님 질문드립니다. [1] | emailhere | 2016.08.05 |
207 | 생존분석에서 생존표... [2] | 파랑파랑하늘 | 2016.08.02 |
206 | 생존분석후 파일 다운 받기가 잘 안되요... [3] | 따거81 | 2016.08.02 |
205 | 메타분석에서 zero cell에 대하여 질문드립니다. [1] | mwg11 | 2016.08.02 |
204 | 비교 통계 - matched t-test 질문 [1] | 근면성실 | 2016.08.02 |
203 | 교수님 propensity score관련하여 질문드립니다. [1] | tunneled_vision | 2016.07.26 |
202 | 누락된 자료의 처리 중 " 다중대입법" 에 대한 질문 [1] | kaiross | 2016.07.26 |
201 | 안녕하세요 propensity score관련하여 문의드립니다. [1] | tunneled_vision | 2016.07.22 |
200 | 표만들기에서 이런 에러가 생기는데 뭐가 문제인지요? [2] | 지훈이아빠 | 2016.07.22 |
199 | 샘플사이즈 문제 [3] | 반디12 | 2016.07.20 |
198 | 데이터 분석 후 예상과 결과가 다를 때, [1] | 반디12 | 2016.07.19 |
197 | 로그변환 후 회귀분석 시 결과물 [1] | 반디12 | 2016.07.19 |
196 | 데이타 선택 [4] | 김현정 | 2016.07.19 |
195 | 시계열 회귀와 회귀 [1] | 반디12 | 2016.07.19 |
이 부분은 제가 만든 코드는 아니고 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.
도움이 되시길 바랍니다.^^