R/D.m.R

#calculates the Dm-penalty matrix
D.m <- function(penden.env) {
 assign("Dm",kronecker(t(get("L",penden.env)),diag(1,get("p",penden.env)))%*%((kronecker(diag(1,(get("K",penden.env)-get("m",penden.env))),crossprod(get("Z",penden.env))/get("n",penden.env)))%*%kronecker(get("L",penden.env),diag(1,get("p",penden.env)))),penden.env)
}

Try the pendensity package in your browser

Any scripts or data that you put into this service are public.

pendensity documentation built on May 2, 2019, 3:58 a.m.