Различия

Здесь показаны различия между двумя версиями данной страницы.

--- classifydiag:amp_al_feticr [2017/03/27 08:36]
admin
+++ classifydiag:amp_al_feticr [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
 ====== Amp: Testing plot for monomineral P sensor ======
-Plot Amp in Al6 - Fe3+Ti+Cr space.
+Plot Amp in Al6 - Fe3+Ti+Cr space. For [[sensors:p:amp_simakin]].
 === Details ===
@@ Строка 36: / Строка 36: @@
                   X=ds$Fe_p3_pfu + ds$Ti_C_pfu + ds$Cr_C_pfu,
                   Y=ds$Al_c6_pfu);
-marks  <- data.frame(Label=c("Key point\n(1.8,0)"),
-                     X=c(1.8),
-                     Y=c(0.0));
 ds$P <- 16.051*(1.8*ds$Al_c6_pfu/(1.8-ds$Fe_p3_pfu-ds$Ti_C_pfu-ds$Cr_C_pfu))-7.819;
-maxpoint <- ds$P[ds$P == max(ds$P),]
+maxpoint <- ds[ds$P == max(ds$P),]
-limits <- c(0,7,0,3);
+spoints <- ds[maxpoint$P-ds$P < 1,]
-labels <- c("Fe3+Ti+Cr", "Al6");
+model <- lm(Al_c6_pfu~(Fe_p3_pfu+Ti_C_pfu+Cr_C_pfu),spoints)
+mb <- model$coefficients[1]
+mk <- model$coefficients[2]
-lines <- data.frame(X1 = c(maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu),
+limits <- c(0,2,-1,2);
-                    Y1 = c(maxpoint$Al_c6_pfu),
-                    X2 = c(1.8),
+scale <- 1.8/(1.8 - (maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu))
-                    Y2 = c(0));
+mpX <- maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu
+mpY <- maxpoint$Al_c6_pfu
+marks  <- data.frame(Label=c("Key point\n(1.8,0)",
+                             paste0("PmaxLine\n P_max=",format(round(max(ds$P), 1), nsmall = 1) ),
+                             paste0("PmeanLine\n P_mean=",format(round(mean(spoints$P), 1), nsmall = 1))),
+                     X=c(1.8 , 0 , 0),
+                     Y=c(0.0, scale*mpY+0.2, mb-0.2));
+lines <- data.frame(X1 = c(mpX, 0,         0),
+                    Y1 = c(mpY, scale*mpY, mb),
+                    X2 = c(1.8, mpX,       1.8),
+                    Y2 = c(0,   mpY,       mk*1.8+mb));
+labels <- c("Fe3+Ti+Cr", "Al6");
 ret <- c();