Stessa formula con risultati diversi in R ed Excel?


0

Ho calcolato il carico per 3 parametri di qualità dell'acqua in R per determinati giorni per l'intero anno utilizzando questa formula.

enter image description here

Ho calcolato 7 scenari usando la formula precedente. Ogni scenario rappresenta 1 giorno della settimana (lunedì di ogni settimana, martedì di ogni settimana ...) per l'intero anno.

Usando 1 anno dataset (i calcoli dei set di dati e degli excel possono essere scaricati da Qui ), Ho usato (tapply) la funzione in R per calcolare il carico usando la formula precedente

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

Ho calcolato la stessa formula in eccellere prima. Confrontando i risultati ottenuti eccellere e R, sono simili per 2 giorni e diversi per il resto dei giorni. Ti sarei grato se potessi controllare la funzione R e il calcolo excel e potrei aiutarmi a capire perché ho ottenuto risultati diversi usando la stessa formula?

Grazie in anticipo

Risposte:


0

Non so davvero perché vorresti farlo, dato che la formula è per carico annuale. Perché dovrebbe esserci un ciclo settimanale che si adatta esattamente nei giorni feriali? Ma se lo vuoi fare comunque (magari in acque reflue industriali?), Questo è il modo in cui ci proverei:

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

Saluti, Bacca


Molte grazie per il vostro aiuto. In effetti la mia formula per il carico annuale ma il campionamento su base giornaliera è difficile a causa dell'elevato costo del campionamento e dell'analisi. Di conseguenza, utilizziamo campioni settimanali per verificare come i carichi annuali calcolati per (52) campioni differiscano dal carico annuale calcolato da (365) campioni
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.