search for: atot

...ich have two values (files 3,32,33) and these have lines which look like... 01/06/2005,1766, or 15/05/2006,289,114 a1 <- read.csv("file1.csv",header=F) etc... a63 <- read.csv("file63.csv",header=F) names(a1) <- c("mdate","file1.column.description") atot <- merge(a1,a2,all=T) followed by repeatedly doing... atot <- merge(atot, a3,all=T) atot <- merge(atot, a4,all=T) etc... I normally start R with --vanilla. What appears to happen is that atot doubles in size each iteration and just falls over due to lack of memory at about i=17... even...

...he calculations are useful in practice and likely to be right in principle but I cannot see how this occurs: why a calculated value of 7.4 known as numeric is not simply plotted as such. It happened on Windows both 7 and 10 with R version 3.4.1. All best wishes Troels ff <- function(H,SID,ATOT,ka)? H + SID - kw/H - ka*ATOT/(H+ka) ?ka <- 1e-7 ?kw <- 1e-14 ?ATOT <- seq(0,0.3,length=100)*1e-3 ?for (i in 1:length(ATOT))? { ?SID[i] <- uniroot(ff,c(-1,1),tol=.Machine$double.eps,maxiter=100000,ka=ka, ATOT=ATOT[i], H = 10^-7.4)$root} ATOT #confirm pH 0 7.4 H <- c() ?for (...

...d likely to be right in principle but I cannot > see how this occurs: why a calculated value of 7.4 known as numeric is > not simply plotted as such. It happened on Windows both 7 and 10 with > R version 3.4.1. > > All best wishes > > Troels > > ff <- function(H,SID,ATOT,ka)? H + SID - kw/H - ka*ATOT/(H+ka) > ?ka <- 1e-7 > ?kw <- 1e-14 > ?ATOT <- seq(0,0.3,length=100)*1e-3 > > ?for (i in 1:length(ATOT))? { > ?SID[i] <- > uniroot(ff,c(-1,1),tol=.Machine$double.eps,maxiter=100000,ka=ka, > ATOT=ATOT[i], H = 10^-7.4)$root} > ATOT...

Thanks a lot - formatting the ordinate as ylim=c(4,10) before plotting pH also removed the problem, and options(digits=10) confirmed that pH was not all exactly 7.4 - as I knew. Still I wonder just why R chooses to plot(ATOT,pH) as shown with repeated "7.4" instead of some more detailed representation. Thanks a gain and happy New Year! Troels Den 26-12-2017 kl. 01:03 skrev Bert Gunter: > Note that ?all.equal clearly says that it tests for **approximate > equality only** with tolerance "close t...

I hope you will forgive me this simple question on titration. I'm trying to find very small values from the algorithm below, which I believe is correctly formatted, and the constants are also correct. When SID goes over ATOT, fitted vales are much too low compared to the literature. I guess I must be using optimize in a wrong way but cannot find out how to improve it. I'm on windows, R 2.0. Best wishes Troels Ring, Aalborg ATOT <- 0.019 KA <- 3e-7 KW <- 4.4e-14 SID <- 0.01 ff <- function(H,KA,SID...

...le_name_1.csv",header=F) ... a63 <- read.csv("another_silly_filename_63.csv",header=F) and then i'm naming the columns in these like... names(a1)[2] <- "silly column name" ... names(a63)[2] <- "daft column name" then trying to merge()... atot <- merge(a1, a2, all=T) and then using language manipulation to loop atot <- merge(atot, a3, all=T) ... atot <- merge(atot, a63, all=T) etc... followed by more language manipulation for() { rm(a1) } etc... i.e. for (i in 2:63) { atot <- merge(atot, eval(parse(text=paste(&q...

...0 2005-06-10 0.2096477 > a2 mdate value2 1 2005-06-01 0.5347852 2 2005-06-02 0.9322765 3 2005-06-03 0.9106499 4 2005-06-04 0.6810564 5 2005-06-05 0.5871867 6 2005-06-06 0.8123808 7 2005-06-07 0.9675379 8 2005-06-08 0.9470369 9 2005-06-09 0.7493767 10 2005-06-10 0.8864103 > atot <- merge(a1,a2,all=T) However, I find the following results to be quite un-intuitive - are they correct? May I draw your attention to lines 9:12... Should lines 9:11 be there? > atot mdate value value2 1 2005-06-01 NA 0.5347852 2 2005-06-02 0.5287683 0.9322765 3 2...

...Bpy%c1OY< z-Q^bdwMyB|v{IhoHlOG8<!^Mc8OikR>eRyRppr!vR&MLDc3Y2Z*2TJY*WwlD5#?`b zJlz6J59y<bJP$Lm!*cios`;L&@PKL>p|A#2<{s4VQQf>-lD>!Pbq`QJrdLkwAJ+-H z>3+f<dq>Q%*D<NkkH=w4MjKpDRu8VFx;y0c&E*UBKk1L>v^XwnUy=*NbHRA7rq6S) zO5v%6YI30}O&^F#9fR~AH+4I3)w=^%XE%ATotE=kFKLHxI($jdP?Lu|GAmD*mJ8&M zPC at 4z%9flGK7 at y@u!G!@3zFkKIBH6Dzv4X at Si246UyU)n;C5Ib|5}Xq9gSc|*?we? zJI>~-dBY>pyzWeMG}gswYP$I8y3eSguP>Y-QF+uH1az<quim{z(0t)ODq>NXnhwKL zdDCyy7(U_ina_GY`+6v~NHvt?B316lX{jE#d*IfKWkc)E_PdLQ$I!d8caX!2kX<f3 z;iV9sq at U7#=XrgbRcw`nsby>9xw?3+KA...