Nothing
plot.gexp.fe_lsd <- function(x,
main = NULL,
sub = NULL,
colgrid = 'red',
coltext = 'blue',
ltygrid = 'dotted',
lwdgrid = par('lwd'),
xleftimg = par()$usr[1],
ybottomimg = par()$usr[3],
xrightimg = par()$usr[2],
ytopimg = par()$usr[4],
dynamic = FALSE,
random = TRUE,
...)
{
mfactors <- attr(x$X,
'contrasts')
labelfac <- names(mfactors)[-c(1:2)]
labelrow <- names(mfactors)[1]
labelcol <- names(mfactors)[2]
levelsfac <- lapply(mfactors[-c(1:2)],
rownames)
levelsrow <- rownames(mfactors[[1]])
levelscol <- rownames(mfactors[[2]])
labelinter <- suppressWarnings(do.call('interaction',
levelsfac))
levelsinter <- levels(labelinter)
repp <- dim(x$X)[1]/(length(levelsrow))
if(is.null(main)){
main = 'Factorial Structure: Latin Square Design'
}
if(is.null(sub)){
sub <- paste('Factors: ',
paste(labelfac,
collapse = ', '),
'\n',
paste('Levels: ',
paste(unlist(levelsfac),
collapse = ', '),
sep = ''),
'\n',
paste('Rows: ',
length(levelsrow),
sep = ''),
'\n',
paste('Columns: ',
length(levelscol),
sep = ''))
}
ifelse(random == FALSE,
{
rtreat <- latin(length(levelsinter),
levelss = levelsinter,
nrand = 0)
treat <- c(rtreat)
},
{
rtreat <- latin(length(levelsinter),
levelss = levelsinter)
treat <- c(rtreat)
})
rowsquare <- columsquare <- length(levelsrow)
aux_posxcentro <- 1/rowsquare
aux_posxcentro1 <- aux_posxcentro + ((rowsquare - 1)*2/rowsquare)
posxcentro <- posycentro <- seq(aux_posxcentro,
aux_posxcentro1,
by = 2/rowsquare)
if(!dynamic){
op <- par('xaxs', 'yaxs') # Original par('xaxs', 'yaxs')
par(xaxs = 'i',
yaxs = 'i')
plot(1,
type = 'n',
xlim = c(0, 2),
ylim = c(0, 2),
axes = FALSE,
xlab = '',
ylab = '',
main = main,
sub = sub,
...)
box()
grid(nx = columsquare,
ny = rowsquare,
col = colgrid,
lty = ltygrid,
lwd = lwdgrid)
text(x = rep(posycentro, length(posxcentro)),
y = rep(posxcentro, rep(length(posycentro),
length(posxcentro))),
treat,
col = coltext)
arrows(-0.05,
seq(0,
2,
by = 2/rowsquare),
-0.05,
seq(2/rowsquare,
2,
by = 2/rowsquare),
angle = 90,
xpd = TRUE,
code = 3,
length = 0.06)
arrows(seq(0,
2,
by = 2/rowsquare),
2.05,
seq(2/rowsquare,
2,
by = 2/rowsquare),
2.05,
angle = 90,
xpd = TRUE,
code = 3,
length = 0.06)
text(-0.08,
posxcentro,
levelsrow,
col = colgrid,
xpd = TRUE,
srt = 90)
text(posxcentro,
2.08,
levelscol,
col = colgrid,
xpd = TRUE)
par(op) # Restoring the original par('xaxs', 'yaxs')
} else {
auxin <- tcltk::tk_choose.files()
auxin1 <- gsub('[\\s\\S]*?\\.',
'',
auxin,
perl = TRUE)
auxin2 <- toupper(auxin1)
switch(auxin2,
PNG = {
myimage <- png::readPNG(auxin)
},
JPEG = {
myimage <- jpeg::readJPEG(auxin)
},
JPG = {
myimage <- jpeg::readJPEG(auxin)
})
plot(1,
type = 'n',
xlab = '',
ylab = '',
axes = FALSE,
main = main,
sub = sub,
...)
rasterImage(myimage,
xleft = xleftimg,
ybottom = ybottomimg,
xright = xrightimg,
ytop = ytopimg)
text(x = locator(),
y = NULL,
paste(labelrow,
1:rowsquare),
col = coltext)
text(x = locator(),
y = NULL,
paste(labelcol,
1:columsquare),
col = coltext)
text(x = locator(),
y = NULL,
treat,
col = coltext)
}
}
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.