andreas/daDoctoR
1
0
Fork 0
mirror of https://github.com/agdamsbo/daDoctoR.git synced 2024-11-21 19:30:22 +01:00
Code Issues Projects Releases Wiki Activity

u

Browse Source
This commit is contained in:
agdamsbo 2018-10-11 16:57:29 +02:00
parent 142274192b
commit f25d7ae34f
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
R/strobe_pred.R
View File

@ -84,20 +84,20 @@ strobe_pred<-function(meas,adj,data,dec=2){
ns<-names(dat2)[i] ns<-names(dat2)[i]
for (r in 1:length(levels(vec))){ for (r in 1:length(levels(vec))){
vr<-levels(vec)[r] vr<-levels(vec)[r]
dr<-vec[vec==vr] dr<-vec[vec==vr&!is.na(vec)]
n<-as.numeric(length(dr)) n<-as.numeric(length(dr))
nall<-as.numeric(nrow(dat[!is.na(dat2[,c(ns)]),])) nall<-as.numeric(nrow(dat[!is.na(dat2[,c(ns)]),]))
nl<-paste0(ns,levels(vec)[r]) nl<-paste0(ns,levels(vec)[r])
pro<-round(n/nall*100,0) pro<-round(n/nall*100,0)
rt<-paste0(n," (",pro,"%)") rt<-paste0(n," (",pro,"%)")
nq<-rbind(nq,cbind(nl,n)) nq<-rbind(nq,cbind(nl,rt))
} }
} }
if (!is.factor(dat2[,i])){ if (!is.factor(dat2[,i])){
num<-dat2[,i] num<-dat2[,i]
ns<-names(dat2)[i] nl<-names(dat2)[i]
nall<-as.numeric(nrow(dat[!is.na(dat2[,c(ns)]),])) rt<-as.numeric(nrow(dat[!is.na(dat2[,c(nl)]),]))
nq<-rbind(nq,cbind(ns,nall)) nq<-rbind(nq,cbind(nl,rt))
} }
} }
@ -114,7 +114,7 @@ strobe_pred<-function(meas,adj,data,dec=2){
res<-cbind(aor_ci,apv) res<-cbind(aor_ci,apv)
rest<-data.frame(names=row.names(res),res,stringsAsFactors = F) rest<-data.frame(names=row.names(res),res,stringsAsFactors = F)
numb<-data.frame(names=nq[,c("nl")],N=nq[,c("n")],stringsAsFactors = F) numb<-data.frame(names=nq[,c("nl")],N=nq[,c("rt")],stringsAsFactors = F)
namt<-data.frame(names=rnames,stringsAsFactors = F) namt<-data.frame(names=rnames,stringsAsFactors = F)
coll<-left_join(left_join(namt,numb,by="names"),rest,by="names") coll<-left_join(left_join(namt,numb,by="names"),rest,by="names")