r 에서 여러개의 모델을 간단한 코딩으로 돌리고 싶어서 질문 드립니다.
예를들어서
변수가 "y1","y2","y3","var1","var2","var3"순서로 있는 데이터 테이블에서
y1=var1+var2+var3
y2=var1+var2+var3
y3=var1+var2+var3
이 세가지 모델을 for문을 이용해서
x=c(1:3)
for(i in 1:length(x)){
model_x[i] = glm(data[,x[i]] ~ var1+var2+var3 ,data =data, family = "binomial")
}
이렇게 해서 3가지모델을 한거번에 해보려고 했는데 for문으로 원래 안되는건지...방법이 잘못된건지 궁금해서 질문드립니다.
또, 설명변수 부분도 데이터 테이블상의 위치나 이런걸 이용해서 할수 있을까요?
Comment 5
-
cardiomoon
2017.01.23 00:21
-
chamchi
2017.01.23 00:54
모델은 그 세가지 모델입니다.
for문에 나오는 x 제가 잘못적은 부분이 있습니다...
x=("y1","y2","y3") 이렇게 3가지 모델의 각 반응변수의 이름을 지정하고
for문을 이용해서 i=1일때 x[1]으로 x에서 y1을 선택하고, 반응변수로 data 에서 y1변수를 선택해서 실행해서 model_y1저장하는 식으로 세번 반복하려고 합니다.
( x[1]=y1 이고, data[,x[1]] 은 data[,"y1"] 이되서 data에서 y1에 해당하는 열과 모든 행을 불러 올 수 있다고 생각했습니다.)
-
cardiomoon
2017.01.23 23:49
문자열을 language로 사용하실 수는 없습니다. 문자열을 language로 사용하시려면 문자열을 parsing하셔야 합니다. 다음을 참조하세요
> fit=lm(mpg~hp*wt,data=mtcars)
> fit
Call:
lm(formula = mpg ~ hp * wt, data = mtcars)
Coefficients:
(Intercept) hp wt hp:wt
49.80842 -0.12010 -8.21662 0.02785
> var<-"mpg"
> dataname="mtcars"
> temp<-paste0("lm(",var,"~hp*wt,data=",dataname,")")
> temp
[1] "lm(mpg~hp*wt,data=mtcars)"
> fit1=eval(parse(text=temp))
> fit1
Call:
lm(formula = mpg ~ hp * wt, data = mtcars)
Coefficients:
(Intercept) hp wt hp:wt
49.80842 -0.12010 -8.21662 0.02785
-
chamchi
2017.01.24 20:57
답변해 주셔서 정말 감사합니다.
알려주신 자료로 문제를 잘 해결 했습니다.
-
cardiomoon
2017.01.25 14:32
힌트를 조금 드렸는데 해결하셨다니 다행입니다^^
No. | Subject | Author | Date |
---|---|---|---|
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 |
1174 | group-based trajectory modeling | swpapa | 2023.03.19 |
세가지 모델은
y1~var1+var2+var3
y2~var1+var2+var3
y3~var1+var2+var3
이런 모델을 말씀하시나요?
아래의 for문은 갑자기 x가 나오는데요