search for: troels

The help file provided within the rsept.zip is somewhat outdated compared to the one available on http://www.stat.math.ethz.ch/R/manual/; would it be possible to update the file for a more easy download instead of having either to download individually or going online ? Troels =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=--- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat...

...s the vector in spite of > names(ID) NULL > names(patients) NULL > attr(patients,"names") NULL > attr(ID,"names") NULL One thing, is that inconsistent ? Short of editing my data file, which would not be very difficult in this case, how could I do it ? Sincerely yours Troels Ring, MD -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-reques...

Dear friends - copy paste missed SID <- c() before the first loop - sorry BW Troels Den 25-12-2017 kl. 19:12 skrev Troels Ring: > > Dear friends - merry Christmas and thanks a lot for much help during > the year! > > In the example below I fail to understand how the calculated value pH > is represented in a simple plot - also included. The calculations are &...

Dear friends. When using in windows version the menu function to change directory, what actions are actually performed, i.e. what commands could I write to a file to make a CD ? I couldn't find it mentioned anywhere - undoubtedly an oversigth on my part. Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe...

...a simple plot - also included. The 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...

...) 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 to 1.5 e-8. > > So.. > > > all.equal(z,pH, tol = 1e-15) > [1] "Mean relative difference: 6.732527e-11" &g...

...with centrum in (1,-1) and radius 5 to show my girl that the line y=3*x+1 goes through (1,4) and (-2,-5) of the circle, but on Windows XP, R 2.20 the drawing was not good at all, and the known solutions were not shown in the graph. I guess I got it wrong? Is this use not intended ? Best wishes Troels Ring, MD Aalborg, Denmark symbols(x=1,y=-1,circles=5,inches=FALSE, xlim=c(-10,10),ylim=c(-10,10)) curve(3*x+1,-10,10,1000,add=T) abline(v=c(-2,1)) abline(h=c(-1,-5,4))

Dear friends - I'm using Tinn-R in Windows and find it OK. My son is running MacOS and is just starting using R and I know nothing about MacOS and R but we seem to have difficulties in finding something similar to Tinn for MaC. What are your suggestions? Thanks a lot! Troels Ring, MD Aalborg, Denmark

...tention is to have the colour reflect a variable, pH, but the legend is disturbed. The little part marked "3" above the colour scale is unwelcome. Why did it appear? How could I avoid it? I'm on Windows 7, R version 3.4.1 (2017-06-30) -- "Single Candle" All best wishes Troels Ring Aalborg, Denmark library(ggplot2) DF1 <- structure(list(P = c(0, 0.00222222222222222, 0.00444444444444444, 0.00666666666666667, 0.00888888888888889, 0.0111111111111111, 0.0133333333333333, 0.0155555555555556, 0.0177777777777778, 0.02, 0, 0.00222222222222222, 0.00444444444444444, 0.00666666...

...and wanted to use identify to show some points but was told by identify that the x and y vectors producing a fine graph with 84 points were not equal in length. Below are the Dato for date - and length(Dato) finds 9 but str finds 84 as known. Will identify not work in this context ? Best wishes Troels Ring Aalborg, Denmark > Dato [1] "2000-01-04" "2000-01-07" "2000-01-10" "2000-01-13" "2000-01-17" ....... [81] "2003-04-23" "2003-05-14" "2003-07-30" "2003-08-14" > length(Dato) [1] 9 > str(Dat...

...e or copy the directories mass, nnet and class to the library directory - then execute link.html.help() Now execute library(MASS) and data(petrol) Should work. "Troels Ring" <tring at mail1.stofane To: "R-help" <r-help at stat.math.ethz.ch> t.dk> cc:...

Dear friends. I've made a very simple procedure to make Bland-Altman plots and it works OK except for the inability to take a variate as argument in annotation. I'm sure it is a very simple error, so if you have 5 seconds please tell me. Best wishes Troels bland <- function (x) #accepts two columns { mn <- 0.5*(x[,1]+x[,2]) diff <- x[,1]-x[,2] plot(mn,diff,pch=16,xlab="Mean",ylab="Difference", main = paste("Bland Altman plot",mn)) abline(h=mean(diff),lty=2) abline(h=mean(diff)+c(2,-2)*sqrt(var(diff))) } Troel...

G'day Troels, On Mon, 6 Nov 2023 20:43:10 +0100 Troels Ring <tring at gvdnet.dk> wrote: > Thanks a lot! This was amazing. I'm not sure I see how the conditiion > pK1 < pK2 < pK3 is enforced? One way of enforcing such constraints (well, in finite computer arithemtic only "<=&quo...

...o including all three pK values to be in agreement with the molecular information about HEPES - although even this is contentious. Be as it may, I wonder if not your method might work if only we KNOW that pK1 is either positive OR negative, in which case we have pK1 = -exp(theta1)? Best wishes Troels Den 07-11-2023 kl. 05:10 skrev Berwin A Turlach: > G'day Troels, > > On Mon, 6 Nov 2023 20:43:10 +0100 > Troels Ring <tring at gvdnet.dk> wrote: > >> Thanks a lot! This was amazing. I'm not sure I see how the conditiion >> pK1 < pK2 < pK3 is enforced?...

Hi Troels, Its in the incoming directory: ftp://ftp.ci.tuwien.ac.at/incoming Try this for a nice demo: > library(timeslab) > arma(alpha, beta, x, 2, 4, 3, 0.2,500, 30,0) Bernd On 01-Jan-99 Troels Ring wrote: > Hi Johannes Wuebben and a happy new year. > Where did you actu...

Hi, I'm using windows xp and R 2.8.0 - I wonder what is the command to put in a script that has the same effect as when in a plot you choose menu "History" and "Recording". Best wishes Troels -- Troels Ring - - Department of nephrology - - Aalborg Hospital 9100 Aalborg, Denmark - - +45 99326629 - - tring at gvdnet.dk

...="confidence") while this refuses to take the "new" and predict ? lot <- c(30,20,60,80,40,50,60,30,70,60) hours <- c(73,50,128,170,87,108,135,69,148,132) z1 <- lm(hours~lot) new <- data.frame(x=80) predict(z1,new,interval="confidence",level=90) best wishes Troels -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at sta...

Thanks a lot! This was amazing. I'm not sure I see how the conditiion pK1 < pK2 < pK3 is enforced? - it comes from the derivation via generalized Henderson-Hasselbalch but perhaps it is not really necessary. Anyway, the use of Vectorize did the trick! Best wishes Troels Den 06-11-2023 kl. 19:19 skrev Ivan Krylov: > ? Mon, 6 Nov 2023 17:53:49 +0100 > Troels Ring <tring at gvdnet.dk> ?????: > >> Hence I wonder if I could somehow have non linear regression to find >> the 3 pK values. Below is HEPESFUNC which delivers charge in the >>...

...roduced far from uniform (0,1) p values - but I fail to see how it is wrong. I contacted the authors of the paper but they did not answer. So, please, doesn?t the code below specify a bivariate N(0,1) with covariance 0.7? I get p values = 1 all over - not interesting, but how wrong? Best wishes Troels library(MASS) Sigma <- matrix(c(1,0.7,0.7,1),2,2) Sigma res <- NULL for (i in 1:1000){ ff <-(mvrnorm(n=100, rep(0, 2), Sigma, empirical = TRUE)) res[i] <- t.test(ff[,1],ff[,2],paired=TRUE)$p.value} -- Troels Ring - - Department of nephrology - - Aalborg Hospital 9100 Aalborg, Denma...

...thm 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,ATOT) {H^3 + (KA+SID)*H^2+(KA*(SID-ATOT)-KW)*H- KA*KW} -log10(optimize(ff,c(0,.1),tol=.Machine$double.eps,KA=KA,SID=SID,ATOT=ATOT)$minimum)