source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD 2.R')
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error_onerun.R')
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error_onerun.R')
error
# error <- error_table[[1]]
dat <- data.frame(y1 = log10(error[,1]), y2 = log10(error[,2]), y3 = log10(error[,3]), time = 0:20)
dat <- reshape2::melt(data = dat, id.vars = 4, measure.vars = c(1,2,3), value.name = "error", variable.name = "Mode")
g2 <- ggplot2::ggplot(data = dat, aes(x = time, y = error, shape = Mode))+
geom_point(size = 3)+
geom_line()+
theme_bw()+
scale_shape_manual(values = c(0,1,2), labels = c(1,2,3))+
scale_x_continuous(breaks = c(seq(0,20,3),20))+
labs(x = "Iteration", y = TeX("$Log_{10}$(Statistical Error)"))+
theme(axis.text = element_text(size = 14),
axis.title = element_text(size = 14),
title = element_text(size = 14, hjust = 0.5),
legend.text = element_text(size = 17),
legend.title = element_text(size = 17),
legend.key.size = unit(1, "cm"),
legend.position = "bottom")
print(g2)
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error_onerun.R')
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error.R')
error_table[[1]]
error <- error_table
## Plot
dat <- data.frame(y1 = error[,1], y2 = error[,2], y3 = error[,3], time = 0:(NROW(error)-1))
dat <- reshape2::melt(data = dat, id.vars = 4, measure.vars = c(1,2,3), value.name = "error", variable.name = "Mode")
g1 <- ggplot2::ggplot(data = dat, aes(x = time, y = error, shape = Mode))+
geom_point(size = 3)+
geom_line()+
theme_bw()+
scale_shape_manual(values = c(0,1,2), labels = c(1,2,3))+
scale_x_continuous(breaks = c(seq(0,20,3), 20))+
labs(x = "Iteration", y = "Statistical Error")+
theme(axis.text = element_text(size = 14),
axis.title = element_text(size = 14),
title = element_text(size = 14, hjust = 0.5),
legend.text = element_text(size = 17),
legend.title = element_text(size = 17),
legend.key.size = unit(1, "cm"),
legend.position = "bottom")
print(g1)
dat <- data.frame(y1 = log10(error[,1]), y2 = log10(error[,2]), y3 = log10(error[,3]), time = 0:20)
dat <- reshape2::melt(data = dat, id.vars = 4, measure.vars = c(1,2,3), value.name = "error", variable.name = "Mode")
g2 <- ggplot2::ggplot(data = dat, aes(x = time, y = error, shape = Mode))+
geom_point(size = 3)+
geom_line()+
theme_bw()+
scale_shape_manual(values = c(0,1,2), labels = c(1,2,3))+
scale_x_continuous(breaks = c(seq(0,20,3),20))+
labs(x = "Iteration", y = TeX("$Log_{10}$(Statistical Error)"))+
theme(axis.text = element_text(size = 14),
axis.title = element_text(size = 14),
title = element_text(size = 14, hjust = 0.5),
legend.text = element_text(size = 17),
legend.title = element_text(size = 17),
legend.key.size = unit(1, "cm"),
legend.position = "bottom")
print(g2)
error
error_list[[1]]
error_table <- error_list[[1]]
error_table
# for (i in 1:length(a1.list)){
# for (j in 1:length(a2.list)){
for (t in 2:times){
# error_table[[i,j]] <- error_table[[i,j]] + error_list[[t]][[i,j]]
error_table <- error_table + error_list[[t]]
}
# error_table[[i,j]] <- error_table[[i,j]]/times
error_table <- error_table/times
error_table
error <- error_table
log10(error[,1])
dat <- data.frame(y1 = log10(error[,1]), y2 = log10(error[,2]), y3 = log10(error[,3]), time = 0:20)
dat <- reshape2::melt(data = dat, id.vars = 4, measure.vars = c(1,2,3), value.name = "error", variable.name = "Mode")
g2 <- ggplot2::ggplot(data = dat, aes(x = time, y = error, shape = Mode))+
geom_point(size = 3)+
geom_line()+
theme_bw()+
scale_shape_manual(values = c(0,1,2), labels = c(1,2,3))+
scale_x_continuous(breaks = c(seq(0,20,3),20))+
labs(x = "Iteration", y = TeX("$Log_{10}$(Statistical Error)"))+
theme(axis.text = element_text(size = 14),
axis.title = element_text(size = 14),
title = element_text(size = 14, hjust = 0.5),
legend.text = element_text(size = 17),
legend.title = element_text(size = 17),
legend.key.size = unit(1, "cm"),
legend.position = "bottom")
print(g2)
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error_onerun.R')
print(g1)
print(g2)
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/stat_error_onerun.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
output
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
?sqrtm
??sqrtm
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
print(output)
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation3.R')
options(digits = 4)
print(output)
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1_xi1.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1_xi5.R')
sessionInfo()
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
sessionInfo()
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
print(tre_ula)
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
dat.pib
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation2.R')
p
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
s.list.final
debugSource('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/AD.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation1.R')
source('~/OneDrive - Florida State Students/ongoing_projects/GLAA/code_R1/simulation2.R')
