Nothing
gexp.simple_lsd <- function(x,
...)
{
ifelse(is.null(x$fe),
fe <- list(f1 = rep(1,
3)),
fe <- x$fe)
ifelse(is.null(x$rowe),
rowe <- unlist(fe),
rowe <- x$rowe)
ifelse(is.null(x$cole),
cole <- rowe,
cole <- x$cole)
if(length(fe) != 1){
stop('Use only one factor!')
}
if(x$r != 1){
x$r <- 1
warning('Internaly replicates was set to one (r=1)!')
}
treatments <- makeTreatments(fl = x$fl,
fe = fe,
quali = x$qualiquanti$quali,
quanti = x$qualiquanti$quanti,
posquanti = x$qualiquanti$posquanti)
contrast <- makeContrasts(factors = treatments,
quali = x$qualiquanti$quali,
quanti = x$qualiquanti$quanti,
posquanti = x$qualiquanti$posquanti)
rowcolumn <- list()
n <- length(treatments[[1]])
ifelse(is.null(x$rowl),
{
rowcolumn$Row <- factor(1:dim(as.matrix(rowe))[1])
contrast[["Row"]] <- diag(dim(as.matrix(rowe))[1])
},
{
rowcolumn[[names(x$rowl)]] <-factor(unlist(x$rowl))
contrast[[names(x$rowl)]] <- diag(dim(as.matrix(rowe))[1])
})
ifelse(is.null(x$coll),
{
rowcolumn$Column <- factor(1:dim(as.matrix(cole))[1])
contrast[["Column"]] <- diag(dim(as.matrix(cole))[1])
},
{
rowcolumn[[names(x$coll)]] <- factor(unlist(x$coll))
contrast[[names(x$coll)]] <- diag(dim(as.matrix(cole))[1])
})
if(!is.null(x$contrasts)){
contrast[names(x$contrasts)] <- x$contrasts
contrasts <- contrast
}else{
contrasts <- contrast
}
cformula <- paste('~',
paste(c(names(rowcolumn),
names(treatments)),
collapse = '+'))
sorttreatment <- latin(n,
levelss = treatments[[1]],
nrand = 0)
if(x$qualiquanti$quanti){
sorttreatment <- as.ordered(sorttreatment)
}else{
sorttreatment <- c(sorttreatment)
}
combrowcolumn <- expand.grid(rowcolumn,
KEEP.OUT.ATTRS = FALSE)
dados <- data.frame(combrowcolumn,
T1 = sorttreatment,
stringsAsFactors = TRUE)
names(dados) <- gsub("T1",
names(treatments),
names(dados))
XB <- makeXBeta(cformula,
dados,
mu = x$mu,
fe = fe,
blke = x$blke,
rowe = rowe,
cole = cole,
inte = x$inte,
contrasts = contrasts)
Z <- NULL
if(is.null(x$err)){
e <- mvtnorm::rmvnorm(n = dim(XB$XB)[1],
sigma = diag(length(x$mu)))
}else{
if(!is.matrix(x$err))
stop("This argument must be a matrix n x 1 univariate or n x p multivariate!")
e <- x$err
}
yl <- XB$XB + e
colnames(yl) <- paste('Y',
1:dim(yl)[2],
sep = '')
Y <- round(yl,
x$round)
#J.C.Faria
if(!x$qualiquanti$quali){
dados <- lapply(dados,
function(x) if(is.ordered(factor(x))) as.numeric(as.character(x)) else x)
dados <- as.data.frame(dados)
}
dados <- cbind(dados,
Y)
res <- list(X = XB$X,
Z = Z,
Y = Y,
dfm = dados)
class(res) <- c(paste('gexp',
class(x),
sep = '.'),
'gexp',
'list')
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.