as.data.frame(table(df$var1,df$var2,df$var3))

table(inwoner$relhhh,inwoner$identificatiehuis,inwoner$IDNR)

subset(table(inwoner$relhhh,inwoner$identificatiehuis,inwoner$IDNR),table(inwoner$relhhh,inwoner$identificatiehuis,inwoner$IDNR)$Freq!='0')

droplevels(table(inwoner$relhhh,inwoner$identificatiehuis,inwoner$IDNR))

relhhh identificatiehuis      IDNR   Freq
0               a1            1234     1
1               a1            1234     0
0               a1            1235     1
1               a1            1235     5
0               a2            1234     7
1               a2            1234     0
0               a2            1235     0
1               a2            1235     3

relhhh identificatiehuis      IDNR   Freq
0               a1            1234     1
0               a1            1235     1
1               a1            1235     5
0               a2            1234     7
1               a2            1235     3