Commit fc56b74b authored by de Freitas Costa, Eduardo's avatar de Freitas Costa, Eduardo
Browse files

Thommas comments

parent 19e35987
......@@ -79,7 +79,7 @@ set.seed(13)
C[1]<-mean(flock[[1]][time_flock[[1]]>0])
media[1]<-mean(flock[[1]][time_flock[[1]]>5])
#1 step in slaughterhouse: SCALDING
Nenv1[1]<-( media[1]*aext1*10^(conc1*afec1+log(10)*sig_conc1/2) +
afec1*pfec1[j]*wfec1[j]*(media[1])*10^(conc1+log(10)*sig_conc1/2) )/((benv1[j]+cenv1)-(benv1[j]*cenv1))
......@@ -106,7 +106,7 @@ set.seed(13)
(1-afec3)*wfec3[j]*pfec3[j]*(media[1])*10^(conc3+log(10)*sig_conc3/2)
#5 step in slaughterhouse: CHILLING
#4 step in slaughterhouse: CHILLING
Nenv5[1]<- ( aext5*(Next3[1]) )/((benv5[j]+cenv5)-(benv5[j]*cenv5))
......@@ -124,7 +124,7 @@ set.seed(13)
a[1]<-log10(Next5[1]/car_me) #log10(cfu/g)
prev_slaug[1]<-(media[1]+(1-exp(-Next5[1])))-(media[1]*(1-exp(-Next5[1])))
prev_slaug[1]<-(1-exp(-Next5[1]))
#Sub-module Pre-retail reduction# Evers`s model#
......@@ -308,7 +308,8 @@ set.seed(13)
a[t]<-log10(Next5[t]/car_me) #log10(cfu/g)
prev_slaug[t]<-((1-exp(-Next5[t]))+media[t])-((1-exp(-Next5[t]))*media[t])
prev_slaug[t]<-(1-exp(-Next5[t]))
#Sub-module Pre-retail reduction# Evers`s model#
#################################################
......
......@@ -52,7 +52,7 @@ prev<-list()
for (i in 1:16){
human_prev[[i]]<-read.table(paste(wd,"/","Output","/","output_scenario",i,"/","data1.txt",sep=""),header = TRUE)
human_prev[[i]]<-read.table(paste(here(),"/","Output","/","output_scenario",i,"/","data1.txt",sep=""),header = TRUE)
prev_diff[[i]]<-(human_prev[[i]]$H[[length(human_prev[[i]]$H)]]-
human_prev[[1]]$H[[length(human_prev[[1]]$H)]])
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment