Displaying 1 result from an estimated 1 matches for "vlogse".
Did you mean:
vloge
2007 Nov 30
1
meta-analysis on diagnostic tests (bivariate approach)
...the restricted maximum likelihood estimation.
This is the script (or function) so far....
bivar<-function(TP,FP,FN,TN){
# find a way to sum 0.5 to each null cell if any
Se<-TP/(TP+FN)
Sp<-TN/(TN+FP)
logSe<-log(Se/(1-Se))
logSp<-log(Sp/(1-Sp))
vlogSe<-1/(Se*(1-Se)*(TP+FN))
vlogSp<-1/(Sp*(1-Sp)*(TN+FP))
bsvlogSp<-0
cblog<-0
bsvlogSp<-0
sdata<-rbind(bsvlogSp,cblog,bsvlogSp,vlogSe,vlogSp)
Does anyone know if there is any package or function with this bivariate
random effect model that I cuold try to...