Log In Sign Up
Models:
Daniel Galbraith
DanielG/
voomDDA
0
Downloads: 1
[0aeb11]: / GetDnb.R

Download this file

14 lines (10 with data), 286 Bytes 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
GetDnb <- function(ns, x, y,beta){



    uniq <- sort(unique(y))

    ds <- matrix(1, nrow=length(uniq), ncol=ncol(x))

    for(k in 1:length(uniq)){

      a <- colSums(x[y==uniq[k],])+beta

      b <- colSums(ns[y==uniq[k],])+beta

      ds[k,] <- a/b

    }

    return(ds)



}