Различия
Здесь показаны различия между двумя версиями данной страницы.
Предыдущая версия справа и слева Предыдущая версия Следующая версия | Предыдущая версия | ||
classifydiag:amp_al_feticr [2017/03/27 09:55] admin |
classifydiag:amp_al_feticr [2018/11/12 15:07] (текущий) |
||
---|---|---|---|
Строка 36: | Строка 36: | ||
X=ds$Fe_p3_pfu + ds$Ti_C_pfu + ds$Cr_C_pfu, | X=ds$Fe_p3_pfu + ds$Ti_C_pfu + ds$Cr_C_pfu, | ||
Y=ds$Al_c6_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; | 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; | ||
Строка 47: | Строка 43: | ||
spoints <- ds[maxpoint$P-ds$P < 1,] | spoints <- ds[maxpoint$P-ds$P < 1,] | ||
- | model <- lm(Fe_p3_pfu+Ti_C_pfu+Cr_C_pfu~Al_c6_pfu,spoints) | + | model <- lm(Al_c6_pfu~(Fe_p3_pfu+Ti_C_pfu+Cr_C_pfu),spoints) |
mb <- model$coefficients[1] | mb <- model$coefficients[1] | ||
mk <- model$coefficients[2] | mk <- model$coefficients[2] | ||
- | limits <- c(0,2,0,1); | + | limits <- c(0,2,-1,2); |
- | + | ||
- | labels <- c("Fe3+Ti+Cr", "Al6"); | + | |
scale <- 1.8/(1.8 - (maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu)) | scale <- 1.8/(1.8 - (maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu)) | ||
Строка 59: | Строка 53: | ||
mpX <- maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu | mpX <- maxpoint$Fe_p3_pfu + maxpoint$Ti_C_pfu + maxpoint$Cr_C_pfu | ||
mpY <- maxpoint$Al_c6_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), | lines <- data.frame(X1 = c(mpX, 0, 0), | ||
Y1 = c(mpY, scale*mpY, mb), | Y1 = c(mpY, scale*mpY, mb), | ||
- | X2 = c(1.8, mpX, 2), | + | X2 = c(1.8, mpX, 1.8), |
- | Y2 = c(0, mpY, mk*2+mb)); | + | Y2 = c(0, mpY, mk*1.8+mb)); |
+ | |||
+ | labels <- c("Fe3+Ti+Cr", "Al6"); | ||
ret <- c(); | ret <- c(); |