Различия

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

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

Следующая версия
Предыдущая версия
recalc:wtp_at_to_atp [2016/05/17 11:59]
admin создано
recalc:wtp_at_to_atp [2018/11/12 15:08] (текущий)
Строка 38: Строка 38:
   {   {
   wtname <- gsub('​_wtp',​ "",​ wtnamew);   wtname <- gsub('​_wtp',​ "",​ wtnamew);
 +  ​
   econst = subset(elements,​element==wtname) ​   econst = subset(elements,​element==wtname) ​
   mw <- c(mw,​c(1/​econst$weight));​   mw <- c(mw,​c(1/​econst$weight));​
Строка 46: Строка 46:
   data_st1 <- data/100;   data_st1 <- data/100;
   data_st2 <- rowApply('​*',​data_st1,​mw);​   data_st2 <- rowApply('​*',​data_st1,​mw);​
-  st2_sum <- rowSum(data_st2);+  st2_sum <- rowSums(data_st2);
   data_st3 <- (data_st2/​st2_sum)*100;​   data_st3 <- (data_st2/​st2_sum)*100;​
 +
 +  names(data_st3) <- gsub("​_wtp",​ "​_atp",​ names(data_st3));​
   ​   ​
   data_atp <- cbind(data_atp,​ data_st3);   data_atp <- cbind(data_atp,​ data_st3);
Строка 53: Строка 55:
   return(data_atp);​   return(data_atp);​
 } }
 +
 +recalc.wtp_at_to_atp.test <- function () {
 +  td <- data.frame(Name="​test",​ S_wtp=19.87,​ Fe_wtp=34.32,​ Co_wtp=0.17,​
 +                                Ni_wtp=0.16,​ As_wtp=45.29,​ Sb_wtp=0.37);​
 +
 +  td <- rbind(td,​td)
 +
 +  rs <- recalc.wtp_at_to_atp(td);​
 +
 +  expect_equal(object = rs$As_atp[1],​ 32.72, tolerance = 0.01);
 +  expect_equal(object = rs$As_atp[2],​ 32.72, tolerance = 0.01);
 +}
 +
 </​Code>​ </​Code>​
  
 {{tag> order11}} {{tag> order11}}
  • recalc/wtp_at_to_atp.1463486381.txt.gz
  • Последние изменения: 2018/11/12 15:08
  • (внешнее изменение)