이 사이트에서 다중 회귀 분석을 하면
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 |
---|---|---|---|
1094 | 로지스틱 회귀분석 관련 질문 [2] | alse | 2016.06.18 |
1093 | 나무 분석 결과 문제 [3] | 매일이5월 | 2016.06.19 |
1092 | linear regression에서 회귀모형 적합성 여부 [3] | 매일이5월 | 2016.06.19 |
1091 | 학생강의project2 | cardiomoon | 2016.06.21 |
1090 | ROC 분석중 문제가 생겨서 질문 드립니다 [2] | candid00 | 2016.06.26 |
1089 | moonBook2 install중 error message [2] | hawk1227 | 2016.06.27 |
1088 | 표 만들기에서 에러가 생깁니다 [2] | jaljary | 2016.06.27 |
1087 | 선그래프 질문드립니다. [2] | 도치맘 | 2016.07.02 |
1086 | plot.or그래프 그리는 방법 [4] | 두치님 | 2016.07.04 |
1085 | 교수님 결측치 관련하여 문의드립니다. [1] | tunneled_vision | 2016.07.07 |
1084 | 질문드립니다. [1] | 죵죵죵 | 2016.07.10 |
1083 | ROC에서 여러설명변수 갖는 model간 비교 [5] | alse | 2016.07.11 |
1082 | excel 변환 부탁드립니다. [3] | 죵죵죵 | 2016.07.12 |
» | 다중 선형 회귀 분석에서 "relative importance of predictor variables"에 관해 질문드립니다. [2] | Ikarus86 | 2016.07.13 |
1080 | 시계열 회귀와 회귀 [1] | 반디12 | 2016.07.19 |
1079 | 데이타 선택 [4] | 김현정 | 2016.07.19 |
1078 | 로그변환 후 회귀분석 시 결과물 [1] | 반디12 | 2016.07.19 |
1077 | 데이터 분석 후 예상과 결과가 다를 때, [1] | 반디12 | 2016.07.19 |
1076 | 샘플사이즈 문제 [3] | 반디12 | 2016.07.20 |
1075 | 표만들기에서 이런 에러가 생기는데 뭐가 문제인지요? [2] | 지훈이아빠 | 2016.07.22 |
이 부분은 제가 만든 코드는 아니고 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.
도움이 되시길 바랍니다.^^