추가 할 수식과 물건이 모두 인수 인 적절한 수식 추가

Question

나는 여기에서 상당한 양의 읽기를 수행했으며 일반적으로 formula objects을 문자열로 조작하지 않아야한다는 것을 알았지 만 나는 그렇지 않습니다. 매우 안전한 방법으로이 작업을 수행하는 방법을 발견

tf <- function(formula = NULL, data = NULL, groups = NULL, ...) { 
# Arguments are unquoted and in the typical form for lm etc 
# Do some plotting with lattice using formula & groups (works, not shown) 
# Append 'groups' to 'formula': 
# Change y ~ x as passed in argument 'formula' to 
# y ~ x * gr where gr is the argument 'groups' with 
# scoping so it will be understood by aov 
new_formula <- y ~ x * gr 
# Now do some anova (could do if formula were right) 
model <- aov(formula = new_formula, data = data) 
# And print the aov table on the plot (can do) 
print(summary(model)) # this will do for testing 
} 

아마도 내가 온 가장 가까운 reformulate을 사용했지만 그것은 단지 RHS에 +하지 * 제공합니다.

p <- tf(carat ~ color, groups = clarity, data = diamonds) 

및 캐럿 ~ 색상 * 명확성을 위해 AOV 결과를 가지고 :이 같은 기능을 사용하고 싶습니다. 미리 감사드립니다. 여기에 솔루션

는 무슨 일이 일어나고 있는지 보여줍니다 아론의 코멘트 @ 기반으로 작동 버전 :

tf <- function(formula = NULL, data = NULL, groups = NULL, ...) { 
print(deparse(substitute(groups))) 
f <- paste(".~.*", deparse(substitute(groups))) 
new_formula <- update.formula(formula, f) 
print(new_formula) 
model <- aov(formula = new_formula, data = data) 
print(summary(model)) 
} 

Answer 1

내가 update.formula이 문제를 해결할 수 있다고 생각,하지만 난에 문제를 했어 함수 호출 내에서 업데이트. 아래에 코딩 했으므로 작동하지만 변수 이름이 아니라 그룹에 열을 전달하고 있습니다. 그런 다음 해당 열을 함수 데이터 세트에 추가 한 다음 업데이트 작업을 수행합니다.

두 번째 방정식에서 원하는 내용을 정확히 수행하는지는 알 수 없지만 update.formula에 대한 도움말 파일을 살펴보고 약간 혼란 스럽습니다.

http://stat.ethz.ch/R-manual/R-devel/library/stats/html/update.formula.html

tf <- function(formula,groups,d){ 
    d$groups=groups 
    newForm = update(formula,~.*groups) 
    mod = lm(newForm,data=d) 
} 

dat = data.frame(carat=rnorm(10,0,1),color=rnorm(10,0,1),color2=rnorm(10,0,1),clarity=rnorm(10,0,1)) 
m = tf(carat~color,dat$clarity,d=dat) 
m2 = tf(carat~color+color2,dat$clarity,d=dat) 

tf2 <- function(formula, group, d) { 
    f <- paste(".~.*", deparse(substitute(group))) 
    newForm <- update.formula(formula, f) 
    lm(newForm, data=d) 
} 
mA = tf2(carat~color,clarity,d=dat) 
m2A = tf2(carat~color+color2,clarity,d=dat) 

편집 : @Aaron 지적 , 그것은 내 문제를 해결 deparse 및 substitute이다 : 나는 코드 예제에 더 나은 옵션으로 tf2을 추가 한 어떻게 두 작업을 볼 수 있도록 .

Answer 2

범위 지정 및 함수 내에서 함수 호출에 문제가있을 때 사용하는 한 가지 기술은 매개 변수를 문자열로 전달한 다음 해당 문자열에서 함수 내에서 호출을 생성하는 것입니다. 여기에 그 모습이 보입니다.

tf <- function(formula, data, groups) { 
    f <- paste(".~.*", groups) 
    m <- eval(call("aov", update.formula(as.formula(formula), f), data = as.name(data))) 
    summary(m) 
} 

tf("mpg~vs", "mtcars", "am") 

또 다른 예 : https://stackoverflow.com/a/7668846/210673에 대한 대답을 참조하십시오. 내가 업데이 트 '와 같은 문제를했을 수 있습니다 생각이 @slammaster보고에 대한 https://stackoverflow.com/a/14858661/210673