Nothing
lat.fitres <- function(object) {
fit <- object
## estraggo parametri
(MI <- inspect(fit,"est")$nu) # means
(LA <- inspect(fit,"est")$lambda) # factor loadings
(BE <- inspect(fit,"est")$beta) # structural coefficients
## factor scores
(ETA <- predict(fit))
## dati originali
(X2 <- fit@Data@X[[1]])
colnames(X2) <- fit@Data@ov.names[[1]]
## valori attesi
(hatY <- matrix(as.vector(MI),nrow=nrow(ETA),ncol=nrow(LA),byrow=TRUE) + ETA %*% t(LA))
colnames(hatY) <- paste("hat.",colnames(hatY),sep="")
if (!is.null(BE)) {
(hatETA <- ETA %*% t(BE))
etavar <- colnames(hatETA)[!(apply(hatETA,2,sum)==0)]
hatETA <- hatETA[,etavar]
colnames(hatETA) <- paste("hat.",colnames(hatETA),sep="")
}
## stima dei residui
resY <- X2-hatY
colnames(resY) <- paste("e.",colnames(resY),sep="")
if (!is.null(BE)) {
resETA <- ETA[,etavar] - hatETA
colnames(resETA) <- paste("e.",colnames(resETA),sep="")
}
X2 <- data.frame(X2)
if (!is.null(BE)) {
X2 <- cbind(X2,hatY,resY,hatETA,resETA)
} else {
X2 <- cbind(X2,hatY,resY)
}
return(X2)
}
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.