search for: cblog

...lt;-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 use? Neither a I could find any apckage or function to aply the restricted maximum likelihood est...