Nothing
dae.specific <-
function(zpre,rel,zpost,mmat,factorno,componentname,effnames,effcodes,effnandc,comcodes,varcodes,cnames,cnamesie,emat,vmat,icol,iecol,gls,ctable)
# dae.specific() - expectations for a specific component
# zpre, zpost are zi or zm or zc
# returns incremented icol which is next vacant col of W
{
# zaz <- matrix(0,nrows(emat),nrows(emat)
nfactorcodes <- length(effcodes[[factorno]])
k <- 0 # counts estimable components
kie <- 0 # counts inestimable components - pos in comcodes[[fno]][]
if(any(!is.na(match(componentname,ctable$allindenv)))) { # cross-class-cov not estimable
for ( i in 1:nfactorcodes) {
for (j in 1:nfactorcodes) {
if ( i == j) { # within class estimable
k <- k + 1
kie <- kie + 1
specificcomponentname <- paste(varcodes[[factorno]][k],componentname,sep=":")
szpre <- zpre[[effnandc[[factorno]][i]]]
szpost <- zpost[[effnandc[[factorno]][i]]]
zaz <- szpre %*% rel %*% t(szpost)
emat[,icol] <- as.vector(mmat %*% zaz %*% mmat) # one col of W matrix
if(gls) {
vmat[,icol] <- as.vector(zaz) # one col of V matrix
}
cnames[icol] <- specificcomponentname # name for this col
icol <- icol + 1
} # end within class estimable
else { # cross class inestimable - just setup names
kie <- kie + 1
specificcomponentname <- paste(comcodes[[factorno]][kie],componentname,sep=":")
cnamesie[iecol] <- specificcomponentname
iecol <- iecol + 1
} # end cross class inestimable
}
}
} # end cross-class cov not estimable
else { # cross-class-cov estimable
for ( i in 1:nfactorcodes) {
for( j in 1:nfactorcodes) {
k <- k + 1
specificcomponentname <- paste(comcodes[[factorno]][k],componentname,sep=":")
szpre <- zpre[[effnandc[[factorno]][i]]]
szpost <- zpost[[effnandc[[factorno]][j]]]
zaz <- szpre %*% rel %*% t(szpost)
emat[,icol] <- as.vector(mmat %*% zaz %*% mmat) # one col of W matrix
if(gls) {
vmat[,icol] <- as.vector(zaz) # one col of V matrix
}
cnames[icol] <- specificcomponentname # name for this col
icol <- icol + 1
}
}
} # end cross-class cov estimable
daelist <- list(cnames=cnames,cnamesie=cnamesie,emat=emat,vmat=vmat,icol=icol,iecol=iecol)
return(daelist)
}
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.