diff --git a/src/sed/magdziarz/components.csv b/sed/magdziarz1997/components.csv
similarity index 100%
rename from src/sed/magdziarz/components.csv
rename to sed/magdziarz1997/components.csv
diff --git a/reference/magdziarz1997.tab b/sed/magdziarz1997/magdziarz1997.tab
similarity index 100%
rename from reference/magdziarz1997.tab
rename to sed/magdziarz1997/magdziarz1997.tab
diff --git a/src/sed/magdziarz/.RData b/src/sed/magdziarz/.RData
deleted file mode 100644
index 9ff2437..0000000
Binary files a/src/sed/magdziarz/.RData and /dev/null differ
diff --git a/src/sed/magdziarz/.Rhistory b/src/sed/magdziarz/.Rhistory
deleted file mode 100644
index cd1272f..0000000
--- a/src/sed/magdziarz/.Rhistory
+++ /dev/null
@@ -1,146 +0,0 @@
-data=read.csv("components.csv",header=T)
-plot(data,log="xy",ylim=c(.0003,.1),xlim=c(.001,1000))
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_sum = data$blackbody_sum_1 + data$comptonization + data$excess
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_error = data$total - data$test_sum
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-warnings()
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-head(data)
-mean(data$total - data$test_error)
-mean((data$total - data$test_error)^2)
-mean(sqrt((data$total - data$test_error)^2))
-data$test_sum = data$blackbody_sum_2 + data$comptonization + data$excess
-mean(sqrt((data$total - data$test_error)^2))
-mean(sqrt((data$total - data$test_sum)^2))
-data$test_sum = data$blackbody_sum_1 + data$comptonization + data$excess
-mean(sqrt((data$total - data$test_sum)^2))
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_error = data$total - data$test_sum
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_sum = data$blackbody_sum_2 + data$comptonization + data$excess
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_error = data$total - data$test_sum
-mean(sqrt((data$total - data$test_sum)^2))
-rms
-?rms
-library(rms)
-data$test_error = data$total - data$test_sum
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_error_sq = (data$total - data$test_sum)^2
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data = read.csv("components.csv",header=t)
-data = read.csv("components.csv",header=T)
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-data$test_sum_1 = data$blackbody_sum_1 + data$comptonization + data$excess
-data$err_1_sq = (data$total - data$test_sum)^2
-data$test_sum_2 = data$blackbody_sum_2 + data$comptonization + data$excess
-data$err_2_sq = (data$total - data$test_sum_2)^2
-data$err_1_sq = (data$total - data$test_sum_1)^2
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-warnings()
-data$err_1_sq = (data$total - data$test_sum)^2
-data$err_1_sq = (data$total - data$test_sum_1)^2
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-attach(data)
-plot(err_1_sq ~ energy..keV.)
-plot(err_2_sq ~ energy..keV.)
-head(err_1_sq)
-head(err_2_sq)
-plot(err_2_sq ~ energy..keV.,log="x")
-plot(err_1_sq ~ energy..keV.,log="x")
-plot(total ~ energy..keV.,log="x")
-contour()
-?read.table
-TOTL1216 = read.table("../../../TOTL1216.tab",header=T)
-head(TOTL1216)
-contour(TOTL1216)
-man contour()
-?contour
-contour(TOTL1216$Hden,TOTL1216$Phi.H.,TOTL1216$Eq_Width.A.)
-?contour
-contour(TOTL1216$Phi.H.,TOTL1216$Hden,TOTL1216$Eq_Width.A.)
-TOTL1216
-contour(TOTL1216$Phi.H.,TOTL1216$Hden,TOTL1216$Eq_Width.A.)
-contour(y = TOTL1216$Phi.H.,x = TOTL1216$Hden,z = TOTL1216$Eq_Width.A.)
-?contour
-?contour
-contour(seq(1:4),seq(5:8),seq(3:6))
-?contour
-?op
-?op
-?opar
-?contour
-filled.contour
-?filled.contour
-filled.contour(seq(1:4),seq(5:8),seq(3:6))
-contour(y = TOTL1216$Phi.H.,x = TOTL1216$Hden,z = TOTL1216$Eq_Width.A.)
-TOTL1216$Hden
-TOTL1216$Phi.H.
-?sort
-TOTL1216[order(TOTL1216$Phi.H.,TOTL1216$Hden)]
-TOTL1216[order(Phi.H.,Hden)]
-TOTL1216[order(TOTL1216$Phi.H.,Hden)]
-TOTL1216[order(TOTL1216$Phi.H.,Hden),]
-TOTL1216[order(TOTL1216$Phi.H.,TOTL1216$Hden),]
-TOTL1216[order(TOTL1216$Hden,TOTL1216$Phi.H.),]
-library(ggplot)
-library(ggplot2)
-library(ggplot2)
-ggplot(Statsrep, aes(x=X, y=Y, z=Zn)) + 
-    geom_density2d()
-ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + 
-ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + geom_density2d()
-library(stats)
-ggplot(Statsrep, aes(x=TOTL1216$Hden, y=TOTL1216$Phi.H., z=TOTL1216$Eq_Width.A.)) + geom_density2d()
-?ggplot
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d()
-TOTL1216$Eq_Width.A.
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_density2d()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.))
-?ggplot
-?interp
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour(alpha=0.5)
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) + geom_contour(alpha=1)
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) +
-  geom_tile() + 
-  coord_equal() +
-  geom_contour(color = "white", alpha = 0.5) + 
-  s
-  theme_bw()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) +
-  geom_tile() + 
-  coord_equal() +
-  geom_contour(color = "white", alpha = 0.5) + 
-  theme_bw()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) +
-  geom_tile("TOTL 1216") + 
-  coord_equal() +
-  geom_contour(color = "white", alpha = 0.5) + 
-  theme_bw()
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) +
-  geom_tile() + 
-  coord_equal() +
-  geom_contour(color = "white", alpha = 0.5) + 
-  theme_bw()
-?ggplot
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-plot(err_1_sq ~ energy..keV.)
-plot(err_1_sq ~ energy..keV.,log="x")
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-plot(err_2_sq ~ energy..keV.,log="x")
-plot(total ~ energy..keV.,log="x")
-mean(sqrt((data$total - data$test_sum)^2))
-mean(sqrt((data$total - data$test_sum_1)/data$total))
-mean((data$total - data$test_sum_1)/data$total)
-mean((data$total - data$test_sum_2)/data$total)
-plot(data,log="xy",ylim=c(.0003,.2),xlim=c(.001,1000))
-ggplot(TOTL1216, aes(x=Hden, y=Phi.H., z=Eq_Width.A.)) +
-  geom_tile() + 
-  coord_equal() +
-  geom_contour(color = "white", alpha = 0.5) + 
-  theme_bw()
-q()