Различия

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

Ссылка на это сравнение

Предыдущая версия справа и слева Предыдущая версия
Следующая версия
Предыдущая версия
classifydiag:amp_al_feticr [2017/03/27 10:02]
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,ds)+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,3,-5,5);+limits <- c(0,2,-1,2);
  
 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))
Строка 57: Строка 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,       1), +                    X2 = c(1.8, mpX,       1.8), 
-                    Y2 = c(0,   ​mpY, ​      ​mk*1+mb));​+                    Y2 = c(0,   ​mpY, ​      mk*1.8+mb));
  
-labels <- c(paste0("​Fe3+Ti+Cr"​,​toString(mb),​toString(mk)), "​Al6"​);​+labels <- c("​Fe3+Ti+Cr",​ "​Al6"​);​
  
 ret <- c(); ret <- c();
  • classifydiag/amp_al_feticr.1490608953.txt.gz
  • Последние изменения: 2018/11/12 15:08
  • (внешнее изменение)