Here is a peace of the data and code :

sex <-c(2,2,1,1,1,1,1,1,1,2,1,1,1,1,1,2,2,1,1,2,2)
AGN <-
c("C","C","C","C","C","A","A","C","B","B","C","C","C","C","C","C","B","B","C","C","C")
X <- c(2,2,2,2,1,2,2,2,2,2,1,1,1,1,1,1,1,1,1,2,2)

varqual <- data.frame(sex,AGN,X)

desqual <- summary.formula(varqual$X ~.,  data = subset( varqual, select =
-X), method = "reverse",  overall = T, test = T, long = T, exclude1 = F)

desqual doesn't show the results for sex ==1 as it is redundant.
I also tried long =T wich didn't change anything here.

Bonus question if I may :
This function is a little bit complex for me and  I couldn't figure out how
to get either Yates' continuity correction or Fisher Test instead of
Chisquare. Does it ask a lot of program coding ?

Regards.

David

	[[alternative HTML version deleted]]

