我需要测试一个公式是否是片面的(例如~ a
而不是a~b
)。
现在我正在做这样的事情:
test <- list( ~ a + b, a ~ b + c, b + c ~ a )
isOneSided <- function(form) length(form)==2 && sum(grepl("~",form))==1
> sapply(test,isOneSided)
[1] TRUE FALSE FALSE
有更好的方法吗?我担心有些类型的公式我不知道可以逃避这个测试。
答案 0 :(得分:6)
我会使用terms
函数并提取响应属性:
test <- list( ~ a + b, a ~ b + c, b + c ~ a )
sapply( test , function(x) attr( terms(x) , "response" ) == 0 )
# [1] TRUE FALSE FALSE
由于@Arun指出terms
无法在不知道特殊引用的.
的情况下扩展其中包含data.frame
等特殊项的公式对象。解决此问题的方法是在data.frame
函数调用中包含一个虚拟terms
:
## If we want to expand the '.' in b + c ~ .
test <- list( ~ a + b, a ~ b + c, b + c ~ a , b + c ~ . , . ~ b + c )
sapply( test , function(x) attr( terms(x , data = data.frame(runif(1))) , "response" ) == 0 )
# [1] TRUE FALSE FALSE FALSE FALSE