simulation/Revision1/analysis_simulation_greedy_conditions_poskeyed.R
In susanne-frick/MFCblockInfo: compute block information for multidimensional forced-choice questionnaires
####--------------- analysis simulation testinfo ---------------####
library(psych)
devtools::load_all()
devtools::load_all("../DataAnalysisSimulation/")
####------------------ read in and check results ---------------####
#because 1/400 replications did not converge, results were saved separately
#and combined with combine_results_optimal_testdesign.R
res <- readRDS("simulation/results_simulation_greedy_conditions_poskeyed_part1.rds")
mean.frame(dvs=c("rec","MAB","RMSE"), ivs=c("algorithm","target"), results=res) #check
nrow(res)/(5*6*2*3) #traits x algorithm x target x blocksize
#some are missing
head(res)
anyNA(res)
#some NAs: MIP did not converge
table(is.na(res$MAB))/5
res[is.na(res$MAB),]
#5 for blocksize 2
#all conditions included?
table(res$blocksize, useNA="always")
table(res$keying)
table(res$intercepts)
table(res$loads)
table(res$target)
table(res$ntraits)
#30 rows for the weighted target and blocksize 2 are missing
table(res$rep)
which(table(res$rep)<180) #only one
#identifier for which ones are missing
id.res <- paste(res$blocksize, res$target, res$constraints, res$rep, sep="_")
# saveRDS(id.res, file="simulation/missing_rows_greedy_conditions_poskeyed.rds")
#add missing rows which were run by now
res.miss <- readRDS("simulation/results_simulation_greedy_conditions_poskeyed_missingrows.rds")
res <- rbind(res, res.miss)
tail(res.miss)
tail(res)
#check that no runs were duplicated
id.res <- paste(res$blocksize, res$target, res$constraints, res$rep, res$trait, res$algorithm, sep="_")
any(duplicated(id.res))
####--------------- data preparation ----------------------####
#remove variables that were not varied
res <- res[, ! colnames(res) %in% c("keying","intercepts","loads","length")]
#also remove ntraits because they are computationally too intensive
res <- res[, ! colnames(res) %in% c("ntraits")]
head(res)
#remove rows with convergence issues
res <- res[rowSums(is.na(res))==0,]
#add true reliability
res$rel <- res$rec^2
#add Fisher Z of r(true,est)
res$fisherz.r <- fisherz(res$rec)
#move column target to the beginning
res <- res[,c("rep","trait","blocksize","algorithm","target","rec","fisherz.r","rel","RMSE","MAB")]
head(res)
tail(res)
#make algorithm a factor
str(res)
head(res$algorithm)
res$algorithm <- factor(res$algorithm, levels=c("greedy-a","greedy-d","opt-t","r2","loads","random"))
head(res$algorithm)
####------------------ descriptives in text -------------####
# algorithm vs. random
round(mean(res$MAB[! res$algorithm %in% c("random")]), 2)
round(mean(res$MAB[res$algorithm %in% c("random")]), 2)
# Greedy + R2 vs. MIP T + loads
round(mean(res$MAB[res$algorithm %in% c("greedy-a","greedy-d","r2")]), 2)
round(mean(res$MAB[res$algorithm %in% c("opt","loads")]), 2)
#var: greedy vs. others
# round(sd(res$MAB[res$algorithm %in% c("greedy-a","greedy-d")]), 2)
# round(sd(res$MAB[res$algorithm %in% c("opt","loads","r2")]), 2)
#no real difference in this version
# Equal vs. Weighted
round(mean(res$MAB[res$target %in% c("equal")]), 2)
round(mean(res$MAB[res$target %in% c("weighted")]), 2)
round(mean(res$rec[res$target %in% c("equal")]), 2)
round(mean(res$rec[res$target %in% c("weighted")]), 2)
#blocksize
round(mean(res$MAB[res$blocksize %in% 2]), 2)
round(mean(res$MAB[res$blocksize %in% 3]), 2)
round(mean(res$MAB[res$blocksize %in% 4]), 2)
####------------------ absolute recovery ----------------####
means.paper <- mean.frame(dvs=c("rec","rel","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res, na.rm=T, rn=2)
means.paper$algorithm <- recode.df(means.paper$algorithm,
c("greedy-a","greedy-d","opt-t","r2","loads","random"),
c("Greedy A","Greedy D","MIP T","Block $R^2$","Mean Loadings","Random"))
means.paper$target <- recode.df(means.paper$target, c("equal","weighted"), c("Equal","Weighted"))
#re-order columns
factors <- 1:3
nlevels <- apply(means.paper[,factors], 2, function(cl) length(unique(cl)))
#only every first occurence of factor level
for(cl in factors) {
if(cl > 1) {
means.paper[,cl] <- as.character(means.paper[,cl])
means.paper[-c(seq(1, nrow(means.paper), prod(nlevels[1:(cl-1)]))), cl] <- ""
}
}
means.paper[,factors] <- means.paper[,rev(factors)]
#add SDs
means.paper <- add.brackets.SDs(means.paper)
means.paper
header <- list()
header$pos <- list(-1, nrow(means.paper))
header$command <- c("\\hline \n Algorithm & Target & Algorithm & \\multicolumn{2}{c}{$r(\\theta,\\hat{\\theta})$} & \\multicolumn{2}{c}{$r(\\theta,\\hat{\\theta})^2$} & \\multicolumn{2}{c}{MAB} & \\multicolumn{2}{c}{RMSE} \\\\",
"\\hline \n \\multicolumn{10}{l}{\\small \\textit{Note.} MAB = Mean Absolute Bias, RMSE = Root Mean Squared Error, A = A-optimality,} \\\\
\\multicolumn{10}{l}{\\small D = D-optimality, MIP T = Mixed Integer Programming with T-optimality. Standard} \\\\
\\multicolumn{10}{l}{\\small deviations are given in parentheses.}\n")
print(xtable::xtable(means.paper, digits=2,
caption="Mean trait recovery by algorithm in the simulation on test extension with all positively keyed items",
label="tb:means_ext"), include.colnames = F, include.rownames=F, hline.after=c(0,6,12,18,24,30),
sanitize.rownames.function=function(x){x}, sanitize.colnames.function = function(x){x},
sanitize.text.function = function(x){x},
NA.string = "", table.placement = "htp", add.to.row = header,
caption.placement = "top", latex.environments = NULL,
file="../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/textable_means_extension_poskeyed.tex")
####------------------ differences between algorithms -----------####
contrasts(res$algorithm)
contrasts(res$algorithm) <- matrix(c(rep(1,5),-5, rep(1,4),-4,0, rep(1,3),-3,rep(0,2), rep(1,2),-2,rep(0,3), 1,-1,rep(0,4)), 6, 5,
dimnames=list(c("greedy-a","greedy-d","opt-t","r2","loads","random"),
c("algovsrandom","infovsloadings","optgreedyvsr2","greedyvsopt","avsd")))
contrasts(res$target)
contrasts(res$target) <- matrix(c(1,-1), 2, 1,
dimnames=(list(c("weighted","equal"),
c("weightedvsequal"))))
head(res$blocksize)
res$blocksize <- as.factor(res$blocksize)
head(res$blocksize)
contrasts(res$blocksize)
contrasts(res$blocksize) <- matrix(c(2,-1,-1,0,1,-1), 3, 2,
dimnames=list(c("2","3","4"),
c("2vs34", "3vs4")))
lm.algo.main <- calc.lms.main(dvs=c("fisherz.r","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res)
var.expl(lm.algo.main)
lm.algo <- calc.lms(dvs=c("fisherz.r","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res)
var.expl.algo <- var.expl(lm.algo)
rm.0rows(var.expl.algo)
#format for latex
var.paper <- var.expl.algo
var.paper <- rm.0rows(var.paper)
rownames(var.paper) <- c("Algorithm vs. Random", "Info vs. Mean Loadings","Weighted vs. Equal","2 vs. 3 and 4","3 vs. 4",
"Info vs. Mean Loadings $\\times$ 2 vs. 3 and 4",
"Opt and Greedy vs. $R^2$ $\\times$ 2 vs. 3 and 4",
"Greedy-A vs. Greedy-D $\\times$ 2 vs. 3 and 4",
"Residuals")
var.paper
header <- list()
header$pos <- list(-1, nrow(var.paper))
header$command <- c("\\hline \n Factor & $r(\\theta, \\hat{\\theta})$ & MAB & RMSE\\\\",
"\\hline \n \\multicolumn{4}{l}{\\small \\textit{Note.} MAB = Mean Absolute Bias, RMSE = Root Mean Squared Error. } \\\\
\\multicolumn{4}{l}{\\small $r(\\theta, \\hat{\\theta})$ was Fisher \\textit{Z} transformed.}\n")
print(xtable::xtable(var.paper, digits=0,
caption="Variance in trait recovery explained in \\% by algorithm, target and blocksize in the simulation on
test extension with all positively keyed items",
label="tb:var_sim_ext"), include.colnames = F, include.rownames=T, hline.after=c(0, nrow(var.paper)-1),
sanitize.rownames.function=function(x){x}, sanitize.colnames.function = function(x){x},
sanitize.text.function = function(x){x},
NA.string = "", table.placement = "htp", add.to.row = header,
caption.placement = "top", latex.environments = NULL,
file="../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/textable_var_extension_poskeyed.tex")
####---------------- plots --------------####
library(ggplot2)
library(gridExtra)
library(colorspace)
res.equal3 <- res[res$target=="equal" & res$blocksize==3,]
res.equal2 <- res[res$target=="equal" & res$blocksize==2,]
plot.MAB3 <- ggplot(data=res.equal3, aes(y=MAB, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="MAB", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.RMSE3 <- ggplot(data=res.equal3, aes(y=RMSE, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="RMSE", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.rec3 <- ggplot(data=res.equal3, aes(y=rec, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y=expression(r(theta,hat(theta))), x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
ggsave("../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/plot_greedy_poskeyed_B3.pdf",
grid.arrange(plot.MAB3, plot.RMSE3, plot.rec3, nrow=2, ncol=2),
width=15, height=12, units="in")
###----------- block size 2
plot.MAB2 <- ggplot(data=res.equal2, aes(y=MAB, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="MAB", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.RMSE2 <- ggplot(data=res.equal2, aes(y=RMSE, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="RMSE", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.rec2 <- ggplot(data=res.equal2, aes(y=rec, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y=expression(r(theta,hat(theta))), x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
ggsave("../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/plot_greedy_poskeyed_B2.pdf",
grid.arrange(plot.MAB2, plot.RMSE2, plot.rec2, nrow=2, ncol=2),
width=15, height=12, units="in")
susanne-frick/MFCblockInfo documentation built on Oct. 20, 2024, 8:26 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
####--------------- analysis simulation testinfo ---------------####
library(psych)
devtools::load_all()
devtools::load_all("../DataAnalysisSimulation/")
####------------------ read in and check results ---------------####
#because 1/400 replications did not converge, results were saved separately
#and combined with combine_results_optimal_testdesign.R
res <- readRDS("simulation/results_simulation_greedy_conditions_poskeyed_part1.rds")
mean.frame(dvs=c("rec","MAB","RMSE"), ivs=c("algorithm","target"), results=res) #check
nrow(res)/(5*6*2*3) #traits x algorithm x target x blocksize
#some are missing
head(res)
anyNA(res)
#some NAs: MIP did not converge
table(is.na(res$MAB))/5
res[is.na(res$MAB),]
#5 for blocksize 2
#all conditions included?
table(res$blocksize, useNA="always")
table(res$keying)
table(res$intercepts)
table(res$loads)
table(res$target)
table(res$ntraits)
#30 rows for the weighted target and blocksize 2 are missing
table(res$rep)
which(table(res$rep)<180) #only one
#identifier for which ones are missing
id.res <- paste(res$blocksize, res$target, res$constraints, res$rep, sep="_")
# saveRDS(id.res, file="simulation/missing_rows_greedy_conditions_poskeyed.rds")
#add missing rows which were run by now
res.miss <- readRDS("simulation/results_simulation_greedy_conditions_poskeyed_missingrows.rds")
res <- rbind(res, res.miss)
tail(res.miss)
tail(res)
#check that no runs were duplicated
id.res <- paste(res$blocksize, res$target, res$constraints, res$rep, res$trait, res$algorithm, sep="_")
any(duplicated(id.res))
####--------------- data preparation ----------------------####
#remove variables that were not varied
res <- res[, ! colnames(res) %in% c("keying","intercepts","loads","length")]
#also remove ntraits because they are computationally too intensive
res <- res[, ! colnames(res) %in% c("ntraits")]
head(res)
#remove rows with convergence issues
res <- res[rowSums(is.na(res))==0,]
#add true reliability
res$rel <- res$rec^2
#add Fisher Z of r(true,est)
res$fisherz.r <- fisherz(res$rec)
#move column target to the beginning
res <- res[,c("rep","trait","blocksize","algorithm","target","rec","fisherz.r","rel","RMSE","MAB")]
head(res)
tail(res)
#make algorithm a factor
str(res)
head(res$algorithm)
res$algorithm <- factor(res$algorithm, levels=c("greedy-a","greedy-d","opt-t","r2","loads","random"))
head(res$algorithm)
####------------------ descriptives in text -------------####
# algorithm vs. random
round(mean(res$MAB[! res$algorithm %in% c("random")]), 2)
round(mean(res$MAB[res$algorithm %in% c("random")]), 2)
# Greedy + R2 vs. MIP T + loads
round(mean(res$MAB[res$algorithm %in% c("greedy-a","greedy-d","r2")]), 2)
round(mean(res$MAB[res$algorithm %in% c("opt","loads")]), 2)
#var: greedy vs. others
# round(sd(res$MAB[res$algorithm %in% c("greedy-a","greedy-d")]), 2)
# round(sd(res$MAB[res$algorithm %in% c("opt","loads","r2")]), 2)
#no real difference in this version
# Equal vs. Weighted
round(mean(res$MAB[res$target %in% c("equal")]), 2)
round(mean(res$MAB[res$target %in% c("weighted")]), 2)
round(mean(res$rec[res$target %in% c("equal")]), 2)
round(mean(res$rec[res$target %in% c("weighted")]), 2)
#blocksize
round(mean(res$MAB[res$blocksize %in% 2]), 2)
round(mean(res$MAB[res$blocksize %in% 3]), 2)
round(mean(res$MAB[res$blocksize %in% 4]), 2)
####------------------ absolute recovery ----------------####
means.paper <- mean.frame(dvs=c("rec","rel","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res, na.rm=T, rn=2)
means.paper$algorithm <- recode.df(means.paper$algorithm,
c("greedy-a","greedy-d","opt-t","r2","loads","random"),
c("Greedy A","Greedy D","MIP T","Block $R^2$","Mean Loadings","Random"))
means.paper$target <- recode.df(means.paper$target, c("equal","weighted"), c("Equal","Weighted"))
#re-order columns
factors <- 1:3
nlevels <- apply(means.paper[,factors], 2, function(cl) length(unique(cl)))
#only every first occurence of factor level
for(cl in factors) {
if(cl > 1) {
means.paper[,cl] <- as.character(means.paper[,cl])
means.paper[-c(seq(1, nrow(means.paper), prod(nlevels[1:(cl-1)]))), cl] <- ""
}
}
means.paper[,factors] <- means.paper[,rev(factors)]
#add SDs
means.paper <- add.brackets.SDs(means.paper)
means.paper
header <- list()
header$pos <- list(-1, nrow(means.paper))
header$command <- c("\\hline \n Algorithm & Target & Algorithm & \\multicolumn{2}{c}{$r(\\theta,\\hat{\\theta})$} & \\multicolumn{2}{c}{$r(\\theta,\\hat{\\theta})^2$} & \\multicolumn{2}{c}{MAB} & \\multicolumn{2}{c}{RMSE} \\\\",
"\\hline \n \\multicolumn{10}{l}{\\small \\textit{Note.} MAB = Mean Absolute Bias, RMSE = Root Mean Squared Error, A = A-optimality,} \\\\
\\multicolumn{10}{l}{\\small D = D-optimality, MIP T = Mixed Integer Programming with T-optimality. Standard} \\\\
\\multicolumn{10}{l}{\\small deviations are given in parentheses.}\n")
print(xtable::xtable(means.paper, digits=2,
caption="Mean trait recovery by algorithm in the simulation on test extension with all positively keyed items",
label="tb:means_ext"), include.colnames = F, include.rownames=F, hline.after=c(0,6,12,18,24,30),
sanitize.rownames.function=function(x){x}, sanitize.colnames.function = function(x){x},
sanitize.text.function = function(x){x},
NA.string = "", table.placement = "htp", add.to.row = header,
caption.placement = "top", latex.environments = NULL,
file="../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/textable_means_extension_poskeyed.tex")
####------------------ differences between algorithms -----------####
contrasts(res$algorithm)
contrasts(res$algorithm) <- matrix(c(rep(1,5),-5, rep(1,4),-4,0, rep(1,3),-3,rep(0,2), rep(1,2),-2,rep(0,3), 1,-1,rep(0,4)), 6, 5,
dimnames=list(c("greedy-a","greedy-d","opt-t","r2","loads","random"),
c("algovsrandom","infovsloadings","optgreedyvsr2","greedyvsopt","avsd")))
contrasts(res$target)
contrasts(res$target) <- matrix(c(1,-1), 2, 1,
dimnames=(list(c("weighted","equal"),
c("weightedvsequal"))))
head(res$blocksize)
res$blocksize <- as.factor(res$blocksize)
head(res$blocksize)
contrasts(res$blocksize)
contrasts(res$blocksize) <- matrix(c(2,-1,-1,0,1,-1), 3, 2,
dimnames=list(c("2","3","4"),
c("2vs34", "3vs4")))
lm.algo.main <- calc.lms.main(dvs=c("fisherz.r","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res)
var.expl(lm.algo.main)
lm.algo <- calc.lms(dvs=c("fisherz.r","MAB","RMSE"), ivs=c("algorithm","target","blocksize"), results=res)
var.expl.algo <- var.expl(lm.algo)
rm.0rows(var.expl.algo)
#format for latex
var.paper <- var.expl.algo
var.paper <- rm.0rows(var.paper)
rownames(var.paper) <- c("Algorithm vs. Random", "Info vs. Mean Loadings","Weighted vs. Equal","2 vs. 3 and 4","3 vs. 4",
"Info vs. Mean Loadings $\\times$ 2 vs. 3 and 4",
"Opt and Greedy vs. $R^2$ $\\times$ 2 vs. 3 and 4",
"Greedy-A vs. Greedy-D $\\times$ 2 vs. 3 and 4",
"Residuals")
var.paper
header <- list()
header$pos <- list(-1, nrow(var.paper))
header$command <- c("\\hline \n Factor & $r(\\theta, \\hat{\\theta})$ & MAB & RMSE\\\\",
"\\hline \n \\multicolumn{4}{l}{\\small \\textit{Note.} MAB = Mean Absolute Bias, RMSE = Root Mean Squared Error. } \\\\
\\multicolumn{4}{l}{\\small $r(\\theta, \\hat{\\theta})$ was Fisher \\textit{Z} transformed.}\n")
print(xtable::xtable(var.paper, digits=0,
caption="Variance in trait recovery explained in \\% by algorithm, target and blocksize in the simulation on
test extension with all positively keyed items",
label="tb:var_sim_ext"), include.colnames = F, include.rownames=T, hline.after=c(0, nrow(var.paper)-1),
sanitize.rownames.function=function(x){x}, sanitize.colnames.function = function(x){x},
sanitize.text.function = function(x){x},
NA.string = "", table.placement = "htp", add.to.row = header,
caption.placement = "top", latex.environments = NULL,
file="../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/textable_var_extension_poskeyed.tex")
####---------------- plots --------------####
library(ggplot2)
library(gridExtra)
library(colorspace)
res.equal3 <- res[res$target=="equal" & res$blocksize==3,]
res.equal2 <- res[res$target=="equal" & res$blocksize==2,]
plot.MAB3 <- ggplot(data=res.equal3, aes(y=MAB, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="MAB", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.RMSE3 <- ggplot(data=res.equal3, aes(y=RMSE, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="RMSE", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.rec3 <- ggplot(data=res.equal3, aes(y=rec, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y=expression(r(theta,hat(theta))), x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
ggsave("../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/plot_greedy_poskeyed_B3.pdf",
grid.arrange(plot.MAB3, plot.RMSE3, plot.rec3, nrow=2, ncol=2),
width=15, height=12, units="in")
###----------- block size 2
plot.MAB2 <- ggplot(data=res.equal2, aes(y=MAB, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="MAB", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.RMSE2 <- ggplot(data=res.equal2, aes(y=RMSE, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y="RMSE", x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
plot.rec2 <- ggplot(data=res.equal2, aes(y=rec, x=algorithm, fill=algorithm)) +
geom_violin(show.legend=FALSE) +
labs(y=expression(r(theta,hat(theta))), x="Algorithm") +
scale_x_discrete(labels = c('greedy-a' = "Greedy A",
'greedy-d' = "Greedy D",
'opt-t' = "MIP T",
'r2' = expression(plain(Block)~R^2),
"loads" = "Loadings",
"random" = "Random")) +
scale_fill_manual(values=qualitative_hcl(6)) +
theme(axis.text=element_text(size=16),
axis.title=element_text(size=16))
ggsave("../../Projekte/MFC_blocks/paper/Revision1_Psychometrika/SOM/plot_greedy_poskeyed_B2.pdf",
grid.arrange(plot.MAB2, plot.RMSE2, plot.rec2, nrow=2, ncol=2),
width=15, height=12, units="in")
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.