From 2eaaff05da0f89294788bd0bed564e1bf257acf2 Mon Sep 17 00:00:00 2001 From: JuanSilva89 <31934980+JuanSilva89@users.noreply.github.com> Date: Mon, 1 Apr 2019 09:57:21 -0400 Subject: [PATCH] Update Script for Probit Analysis v3.R --- Script for Probit Analysis v3.R | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/Script for Probit Analysis v3.R b/Script for Probit Analysis v3.R index 5bf6e02..1dab071 100644 --- a/Script for Probit Analysis v3.R +++ b/Script for Probit Analysis v3.R @@ -174,7 +174,7 @@ original.data output<-matrix(0,nrow = length(dataf$dose),ncol = 9) output<-cbind(dataf$dose,raw.abbott.correction,"",round(dataf$probit,digits = 3),dataf$total,dataf.output$dead,round(new.expected.mortality,digits = 3),"",round(chitt.output,digits = 3),extreme.values) -colnames(output)<-c("Dose","Mort.Corr(%)","Probit","Total","Dead","Dead Expected","X2","Extreme values","","") +colnames(output)<-c("Dose","Mort.Corr(%)","Probit","Total","Dead","Expected","X2","","","") output #This is the mortality corrected after optimization @@ -225,12 +225,13 @@ write.table(format(eq.mat,digits = 3),append = T, exporttab, sep = "\t", col.nam write("", exporttab, append=T) write.table(format(chimatrix,digits = 3),append = T, exporttab, sep = "\t", col.names=F,row.names = F,quote=F) write("", exporttab, append=T) -write(" Abbott", exporttab, append=T) +write(" Abbott Dead ", exporttab, append=T) write.table(output,append = T, exporttab, sep = "\t", col.names=T,row.names = F,quote=F) write("", exporttab, append=T) +write("Extreme values are marked with an asterisk (*)", exporttab, append=T) +write("", exporttab, append=T) write("", exporttab, append=T) write(" Conf. Interv. 95%", exporttab, append=T) write.table(round(lcs,digits = 3),append = T, exporttab, sep = "\t", col.names=T,row.names = F,quote=F) remove(list=c("datac")) -