Se as variáveis dicotômicas tiverem valores 0 ou 1, pode tentar da seguinte maneira: dados <- data.frame(x1=rbinom(100,1,.5),x2=rnorm(100)) sapply(dados, function(x) all(x%in%c(0,1))) x1 x2 TRUE FALSE