similar to: array question on indexing

Dear friends - please give me a hand. I have a dataset of 40 patients in two groups observed on three occasions. I only want to plot a line for each patient in the two groups. I use the ltext function to put the patient number but fail to make lattice understand the numbers as unique since apparently it starts all over with 1:20 for each panel instead of respecting my desires to have 1:20

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

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'm examining the characteristics of two models that both fit the sodium concentration in 16 pigs quite well under treatment or control conditions. The more complicated model is by anova better than the less complicated model. To take it further I have generated replicate data using the independent variables and parameter estimates under the two models. A clinically

Dear friends. I believe this problem has been discussed in various forms now and then, so I hope you will forgive me I ask how to do a truncated model like this, where the observed y is recorded as 10 whenever it is higher or equal to that value x <- rnorm(1000) y <- 10*x + rnorm(1000) y[which(y>10)] <- 10 and recover the "true" model ? Best wishes Troels Troels Ring,

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

> 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

Dear friends - occurring in Windows R2.6.2 I am modeling physical chemistry in collaboration with a friend who has preferred working in Excel. I used uniroot, and find a solution to a two buffer problem in acid-base chemistry which I believe is physiologically sensible. Using "goal seek" in Excel my friend found another plausible root, quite close to zero, and a plot of the function

On 2017-08-14 5:53 AM, peter dalgaard wrote: >> 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

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

Dear friends. I hope you will bear over with a lengthy question. I have a problem of 76 patients on hemodialysis measured repeatedly once a month for nine months for methylmalonic acid and vitamin B12. The vitamin had been administered previously to these patients but was curtailed in the first part and then started again. There was a very clear response in that as vitamin B12 went down

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

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

Dear friends. I wanted to reproduce a sheet of paper useful for plotting the reciprocal of some biochemical variable (creatinine) over time (Bleyer Am J Kidney Dis 34:576-578,1999). Time on X axis with units years and months. On Y axis, 1/creatinine, but with the original creatinine value indicated - to make it easier for people to use it. When using the graph, the linear decrease in 1/crea over

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

> On 14 Aug 2017, at 13:43 , Spencer Graves <spencer.graves at effectivedefense.org> wrote: > > > > On 2017-08-14 5:53 AM, peter dalgaard wrote: >>> 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

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 actually put the timeslab ? > Yours > > Troels Ring, MD > Department of Nephrology

Dear friends - a couple of papers in PNAS (lastly:framework for making better predictions by directly estimating variables' predictivity, Lo et al PNAS 2016; 113:14277-14282) have focused interest on mapping complex traits to multiple loci spread all over the genome. I have been around on the relevant taskview(s) I hope but fail to see that the backward haplotype transmission association