###################################################################################################################
## Different plots for SensMixed package
###################################################################################################################
### plot results
plotSensMixed <- function(resSensMixed, mult = FALSE, dprime = FALSE, sep = FALSE, cex = 2,
interact.symbol = ":", isFixed = TRUE, isRand = TRUE, isScaling = TRUE)
{
dens <- function(x)
{
if(x < 0.01)
return(dens=500)
if(x < 0.05)
return(dens=100)
if(x < 0.1)
return(dens=50)
return(dens=10)
}
col.bars.F <- function(x)
{
gr <- gray.colors(3)
if(x < 0.05)
return(gr[1])
if(x < 0.1)
return(gr[2])
return(gr[3])
}
col.bars.Chi <- function(x)
{
gr <- gray.colors(2)
if(x<0.05)
return(gr[1])
return(gr[2])
}
Chi <- sqrt(resSensMixed$random$Chi)
pvalueChi <- resSensMixed$random$pvalueChi
Fval <- sqrt(resSensMixed$fixed$Fval)
pvalueF <- resSensMixed$fixed$pvalueF
if("scaling" %in% names(resSensMixed)){
FScaling <- sqrt(resSensMixed$scaling$FScaling)
pScaling <- resSensMixed$scaling$pvalueScaling
}
if(mult == FALSE){
if(isRand)
return(.plotSensMixed(Chi, pvalueChi, title = expression(paste("Barplot for ",
sqrt(chi^2))),
mult=FALSE, sep = FALSE, interact.symbol = interact.symbol,
cex = cex, ylab = expression(sqrt(chi^2))))
if(isFixed){
if(dprime){
if(!"dprimeav" %in% names(resSensMixed$fixed))
stop("Averaged d primes are not available")
return(.plotSensMixed(resSensMixed$fixed$dprime, pvalueF,
title = expression(paste("Barplot for averaged d-primes")),
mult = FALSE, sep = FALSE,
interact.symbol = interact.symbol,
cex = cex, ylab = expression(paste(tilde(d), "-primes"))))
}
else
return(.plotSensMixed(Fval, pvalueF, title = expression(paste("Barplot for ",
sqrt(F), " values")),
mult=FALSE, sep = FALSE, interact.symbol = interact.symbol,
cex = cex, ylab = expression(sqrt(F))))
}
}
else{
if(isRand)
return(.plotSensMixed(Chi, pvalueChi, mult = TRUE, sep = sep,
interact.symbol = interact.symbol, cex = cex, ylab = expression(sqrt(chi^2))))
if(isFixed){
if(dprime){
if(!"dprimeav" %in% names(resSensMixed$fixed))
stop("Averaged d primes are not available")
return(.plotSensMixed(resSensMixed$fixed$dprime, pvalueF, mult = TRUE, sep = sep,
interact.symbol = interact.symbol, cex = cex,
ylab = expression(paste(tilde(d), "-primes"))))
}
else
return(.plotSensMixed(Fval, pvalueF, mult = TRUE, sep = sep,
interact.symbol = interact.symbol, cex = cex, ylab = expression(sqrt(F))))
}
}
if(("scaling" %in% names(resSensMixed)) && isScaling)
return(.plotSensMixed(FScaling, pScaling, title = expression(paste("Barplot for ",
sqrt(F), " values")),
mult=FALSE, sep = FALSE, interact.symbol = interact.symbol,
cex = cex, ylab = expression(sqrt(F))))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.