Та же формула с разными результатами в R и Excel?

Question

Используя эту формулу, я рассчитал нагрузку для 3 параметров качества воды в R для конкретных дней на весь год.

Я рассчитал 7 сценариев, используя предыдущую формулу. Каждый сценарий представляет 1 день недели (понедельник каждой недели, вторник каждой недели ....) на весь год.

Используя набор данных за 1 год (отсюда можно загрузить наборы данных и Excel), я использовал (tapply) функцию в R для расчета нагрузки по предыдущей формуле

FWNH4Week=tapply( 1:length( NH4 ), as.factor( Weekday ), function( x ) ((sum(Flow[x]*NH4[x])/sum(Flow[x]))*sum(Flow)))

Я рассчитал ту же формулу в Excel раньше. Сравнивая результаты, полученные от Excel и R, они схожи в течение 2 дней и отличаются в остальное время. Я был бы признателен, если бы вы могли проверить функцию R и вычислить Excel и могли бы помочь мне понять, почему я получил разные результаты, используя одну и ту же формулу?

заранее спасибо

Answer 1

Я действительно не знаю, почему вы хотели бы сделать это, так как формула для годовой нагрузки. Зачем нужен еженедельный цикл, который точно соответствует рабочим дням? Но если вы хотите сделать это в любом случае (может быть, в условиях промышленных сточных вод?), Вот как я бы это сделал:

weekday <- rep(1:7, 25)
set.seed(007)
Flow <- cumsum(rnorm(25*7)) + 20
NH4 <- Flow + rnorm(25*7)
plot(Flow, type="l")
lines(NH4, col=2)
# anual sum:
sum(Flow)*  sum(Flow*NH4) / sum(Flow) # 228 660
# per weekday
sapply( unique(weekday), function(i)
  {
  wd <- which(weekday==i)
  sum(Flow[wd])*  sum(Flow[wd]*NH4[wd]) / sum(Flow[wd])
  } )
# each ca 32 000 , a seventh of annual load

С уважением, Берри