similar to: pch as single character

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 > useful

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

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

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 "<=" can be enforced) is to rewrite the parameters as: pK1 = exp(theta1) ##

Thanks a lot, Berwin. Unfortunately, pK1 may well be negative and as I understand the literature it may be poorly defined as such, and also seems to be at a boundary, since when lower is set to say rep(-4,3) pK1 is returned as -4 while pK2 and pK3 are undisturbed. Perhaps the point is that pK1 is not carrying any information at the pH around 5. Fair enough, I guess. Only, I believe I need

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

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 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

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

I've been successfully using a dll via dyn.load() with R1.6.0 for Windows, but when I try it under R1.6.1 it manages to crash the program completely. Has there been a change in how R1.6.1 handles dynamic loading? I couldn't spot any such changes in the documentation. This problem occurred on two different machines, and both run the code under R1.6.0 without a problem. Rob Hyndman

Dear friends - below is a subset of a much larger material showing two ways of generating two "lines". The intention 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

> On 14 Aug 2017, at 10:13 , Troels Ring <tring at gvdnet.dk> wrote: > > Dear friends - I hope you will accept a naive question on lm: R version 3.4.1, Windows 10 > > I have 204 "baskets" of three types corresponding to factor F, each of size from 2 to 33 containing measurements, and need to know if the standard deviation on the measurements in each basket,sdd, is

Hello. I have an R program that calls gawk (GNU Awk 3.06 for Windows) from within pipe() to preprocess a large file before it is read into a data frame with read.table(). I've recently upgraded from Win98SE to WinXP, and have also upgraded from R1.5.0 to R1.6.1 over the past month or so. This program worked before the upgrade(s), but now fails. I observe the following sort of behavior with

Dear friends - I'm on winXP, R 2.71 - I have with some help dveloped this multivariate normal model, which gives very plausible results in WinBUGS even without any initial values specified. However, when I then try to run the same model via the bugs function in R2WinBUGS with inits specified as inits=NULL the program stops in a dead end. So I have tried to make inits for the bugs function

Dear friends - I'm on windows, R 2.7.0 I try again asking if anyone can explain why a single pig of 16 makes so wild swings. Warnings are issued, and they are 1: pseudoinverse used at 482.1 2: neighborhood radius 242.1 3: reciprocal condition number 0 4: at 360 5: radius 14400 6: all data on boundary of neighborhood. make span bigger 7: There are other near singularities as well. 14400 8:

dear friends - I have an ordered factor, ID, in a data.frame, labs, and make a table, aa <- with(labs,table(val >150,ID)) - and now the order in the table is according to the alphabetic order in ID, and not as I want it. I have tried adding [order(unique(labs$ID)] but it doesn't help. I have R 2.7.1 on WinXP. Best wishes Troels -- Troels Ring - - Department of nephrology - -

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]

Dear friends, regression towards the mean is interesting in medical circles, and a very recent paper (The American Statistician November 2007;61:302-307 by Krause and Pinheiro) treats it at length. An initial example specifies (p 303): "Consider the following example: we draw 100 samples from a bivariate Normal distribution with X0~N(0,1), X1~N(0,1) and cov(X0,X1)=0.7, We then calculate

Dear friends, I have struggled to do what is likely rather simple but cannot get it working. I have a dataframe with 428 obs. of 5 variables on 26 patients. Two variables are responses, one an index of patients (1:26), one a period-indicator :1 :max 3, and one total time from min 1 to max 36 weeks (max 21 entries). A complete set would be 26*21=546 long but here are only 428 entries. Now,

Dear friends - I hope you will accept a naive question on lm: R version 3.4.1, Windows 10 I have 204 "baskets" of three types corresponding to factor F, each of size from 2 to 33 containing measurements, and need to know if the standard deviation on the measurements in each basket,sdd, is different across types, F. Plotting the observed sdd versus the sizes from 2 to 33, called

Windows 98 R : Copyright 2001, The R Development Core Team Version 1.2.1 (2001-01-15) Dear friends. How comes this works and produce a single prediction: x <- rnorm(15) y <- x + rnorm(15) predict(lm(y ~ x)) new <- data.frame(x = seq(-3, 3, 0.5)) predict(lm(y ~ x), new, se.fit = TRUE) pred.w.plim <- predict(lm(y ~ x), new, interval="confidence") new1 <- data.frame(x=3)