Nothing
BETA.U.saery.indep <-
function(X, ydi, D, md, sigma2edi, sigmau1, sigmau2) {
p <- ncol(X)
a <- list(1:md[1])
mdcum <- cumsum(md)
for(d in 2:D)
a[[d]] <- (mdcum[d-1]+1):mdcum[d]
yd <- Xd <- Vd.inv <- list()
Q.inv <- matrix(0, nrow=p, ncol=p)
XVy <- 0
for(d in 1:D) {
yd[[d]] <- ydi[a[[d]]]
Xd[[d]] <- X[a[[d]],]
uno <- rep(1,md[d])
V1d <- uno%*%t(uno)
V2d <- diag(uno)
Vd <- (sigmau1*V1d + sigmau2*V2d + diag(sigma2edi[a[[d]]]))
Vd.inv[[d]] <- solve(Vd)
Q.inv <- Q.inv + t(Xd[[d]])%*%Vd.inv[[d]]%*%Xd[[d]]
XVy <- XVy + t(Xd[[d]])%*%Vd.inv[[d]]%*%yd[[d]]
}
Q <- solve(Q.inv)
beta <- Q%*%XVy
u1 <- u2 <- list()
for(d in 1:D){
u1[[d]] <- sigmau1*( t(uno)%*%Vd.inv[[d]]%*%(yd[[d]]-Xd[[d]]%*%beta) )
u2[[d]] <- sigmau2*( Vd.inv[[d]]%*%(yd[[d]]-Xd[[d]]%*%beta) )
}
u1 <- as.matrix(unlist(u1))
u2 <- as.matrix(unlist(u2))
return(list(beta,u1,u2))
}
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.