similar to: Rycas help

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

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, 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 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. I have trobpe with path. Below I have tried to paste a path and file name but it will not work > path <- "D:/rw0641/own/own/procedures/GFR_TIME" > path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good

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

Dear friends. I just downloaded tseries dated 27.10.99 (from Aalborg, Denmark) and it seems it lacked something, library ts was not available ? What to do ? > library(tseries) Loading required package: ts Warning: There is no package called `ts' > > Best wishes Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Dear friends - I am trying to understand fda and working with the 2009 book from Springer. I run the scripts directly from a new installed fda library and include sessionInfo() below. This is from the script fdarm-ch09.R - I seem to have got not fda_2.2.8 but 2.2.7 - is that the problem? Where do I get .8 as I just reinstalled? Best wishes Troels Ring, Nephrology 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

Dear friends. I would like to link an object file made in C. Dyn.load rather wants to link a whole library and I only wanted to link this small object file. I followed instructions given for S - by Clive Loader in his software for changepoint handling and it is his object file I made from http://cm.bell-labs.com/cm/ms/departments/sia/clive/smcp.shar Any dirty tricks or easy oversights ? Yours

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. In the documentation for PAR it is specified that "las" allowing turning axis labels away from default=parallel to axis is not implemented. I need to make graphs with horizontal Y axis labels - how to ? I'm using win98. Thank you in advance. Sincerely yours Troels Ring, MD Department of Nephrology Aalborg Hospital Aalborg Denmark

Dear friends. I downloaded the most recent of everything. Now the directory structure when installing VR5_3pl037.zip in library seem to have been disturbed - at least on my old machine. ....... R : Copyright 1999, The R Development Core Team Version 0.64.1 (May 8, 1999) ......... > library(VR) Warning: Package `VR' contains no R code > data(petrol) Warning: Data set `petrol' not

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

Dear friends. I tried the session below with 10 MB in both vsize and nsize but didn't get the example work. Is this a problem in LME or in me or both or somewhere else or undefined ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type