如何测试公式是否是片面的?

时间:2013-05-04 16:55:12

标签: r

我需要测试一个公式是否是片面的(例如~ 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

有更好的方法吗?我担心有些类型的公式我不知道可以逃避这个测试。

1 个答案:

答案 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