def xperiment_sophiexp():
X, Y, head = parser.csv(currentroot + "SophieXP.csv")
plotter = Scatter2D(X,
Y.ravel(),
axlabels=[axtitles["DHI"], axtitles["D13C"]])
add_reference_ellipses(plotter.ax, 2, ["alma"])
# add_reference_ellipses(plotter.ax, 3, ["alma"])
plotter.scatter(label_points=plotter.Y)
plt.show()
def pull_data(source, sep="\t"):
from csxdata.parser import parser
data, labels, header = parser.csv(projectroot + source,
indeps=1,
headers=1,
dehungarize=True,
decimal=True,
sep=sep)
X_C, Y_C, DHI, D13C = data.T
CCode = labels[:, 0]
return Y_C, DHI, D13C, CCode
def load_eu_full():
X, Y, head = parser.csv(currentroot + "EUborMinden.csv", decimal=True)
return X, Y.ravel()
def load_world_nmr():
X, Y, head = parser.csv(currentroot + "WW_NMR.csv", decimal=True)
return X[:, 0], X[:, 1].ravel(), Y.ravel()
def load_eu_nmr():
X, Y, head = parser.csv(currentroot + "EUborNMR.csv", decimal=True)
return X[:, 0], X[:, 1].ravel(), Y.ravel()
def xperiment_czhech_samples():
def plotone(czparam, huparam, pnm):
plt.title(ttl.format(axtitles[pnm]))
ax = plt.gca() # type: plt.Axes
for i, name in enumerate(sorted(cznmr), start=1):
czmin, czmax = czparam[name]
czh = czmax - czmin
czr = Rectangle((i - 0.2, czmin),
0.2,
czh,
facecolor="red",
edgecolor="black")
ax.add_patch(czr)
ax.annotate(f"{czmin:.2f}",
xy=(i, czmax),
xycoords="data",
verticalalignment="bottom",
horizontalalignment="right")
ax.annotate(f"{czmax:.2f}",
xy=(i, czmin - 0.05),
xycoords="data",
verticalalignment="top",
horizontalalignment="right")
humin, humax = huparam[name]
huh = humax - humin
hur = Rectangle((i, humin),
0.2,
huh,
facecolor="green",
edgecolor="black")
ax.add_patch(hur)
ax.annotate(f"{humin:.2f}",
xy=(i, humax),
xycoords="data",
verticalalignment="bottom",
horizontalalignment="left")
ax.annotate(f"{humax:.2f}",
xy=(i, humin - 0.02),
xycoords="data",
verticalalignment="top",
horizontalalignment="left")
ax.autoscale()
ax.set_ylabel(axtitles[pnm])
ax.set_xlim([0, 7])
ax.set_xticklabels([""] + sorted(cznmr))
fm = plt.get_current_fig_manager()
fm.window.showMaximized()
plt.tight_layout()
plt.show()
X, Y, head = parser.csv(currentroot + "cseh.csv", decimal=True)
cznmr = dict(zip(Y.ravel(), X[:, :2]))
czirms = dict(zip(Y.ravel(), X[:, 2:]))
X, Y, head = parser.csv(currentroot + "SZIdata.csv", 4)
Y = Y[:, 3]
hudict = split_by_categories(Y.ravel(), X)
hunmr = {
name: (d[:, 0].min(), d[:, 0].max())
for name, d in hudict.items()
}
huirms = {
name: (d[:, 2].min(), d[:, 2].max())
for name, d in hudict.items()
}
ttl = "Cseh (piros) és magyar (zöld) {} érték tartományainak összevetése"
plotone(cznmr, hunmr, "DHI")
plotone(czirms, huirms, "D13C")