Nothing
MakeModel <-
function(RespDist="gaussian",BinomialDen=NULL, DataMain, MeanModel,DispersionModel=NULL) {
mc <- match.call()
formulaMean<-MeanModel[3][[1]]
fr <- HGLMFrames(mc, formulaMean,contrasts=NULL)
namesX <- names(fr$fixef)
namesY <- names(fr$mf)[1]
y <- matrix(fr$Y, length(fr$Y), 1)
x <- fr$X
n<-nrow(x)
p<-ncol(x)
random_mean<-findbars(formulaMean)
if (!is.null(random_mean)) {
FL <- HGLMFactorList(formulaMean, fr, 0L, 0L)
namesRE <- FL$namesRE
z <- FL$Design
nrand <- length(z)
q <- rep(0, nrand)
for (i in 1:nrand) q[i] <- dim(z[[i]])[2]
z<-zz<-z[[1]]
} else {
z <- NULL
nrand <- 1
q <- rep(0, nrand)
for (i in 1:nrand) q[i] <- 0
}
res<-list(y,x,z,namesX,namesY,n,p,q)
return(res)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.