Dear all,
I?m specifying the fisher.exact test for use with summary.formula as follows:
u<-function(a,b){
j<-fisher.test(a)
p<-list(P=j$p.value,stat=NA,df=NA,testname=j$method,statname="")
return(p)
}
However I?m also required to specify stat & df. However this doesnt apply to the fisher test. I?ve tried specifying them as NA and "" without success-throws either a blank or an error msg trying to round a non-numeric value respectively.
reproducible example:
ex<-pbc
summary(trt~sex+ascites,data=ex,test=T,method="reverse")
summary(trt~sex+ascites,data=ex,test=T,method="reverse",catTest=u)
The closest I get is
u<-function(a,b){
j<-fisher.test(a)
p<-list(P=j$p.value,stat=1,df=1,testname=j$method,statname="")
return(p)
}
However then I manually have to edit the output. Is there a smart way of doing this?
//M