similar to: codamenu() :Error in coda.options....

Dear R users, I was wondering if anyone knows how to obtain(subset) the first and/or the last record of a subject in a longitudinal setup. Normally in SAS one uses first.variable1 and last.variable1. So my question is that is there an R way of doing this. Regards, -- Luwis Diya, Phd student (Biostatistics), Biostatistical Center, School Of Public Health, Catholic University of Leuven, U.Z. St

Dear R Users, I have a question, how does one print an r function like t.test in R (to get the whole function and not just a summary of the environment)? For example if I type the following: > t.test function (x, ...) UseMethod("t.test") <environment: namespace:stats> > print(t.test) function (x, ...) UseMethod("t.test") <environment: namespace:stats> I

Dear R users I have been encountering an error (in the Winbugs window: Blackbox - index out of range ..... ). So to me it seems R2WinBUGS is not linking up well with the new R version since the same model runs well under r-2.9.2. So is my assumption correct or am I doing something wrong. The data I used and the code are in the help of R2WinBUGS (?bugs). Regards, -- Luwis Diya, Leuven

Dear R users, I would like to fit cross classified or multiple membership logistic models or a 3 level hierarchical logistic model using the Umacs package. Can anyone advise me on how to proceed or better point me to examples of how its done. Regards, -- Luwis Diya, Leuven Biostatistics and Statistical Bioinformatics Centre (L-BioStat), Kapucijnenvoer 35 blok d - bus 7001, 3000 Leuven,

Dear R users: I am a newby to R, under Windows XP. Sorry if my question is too obvious, but I have not found any answer in the archives. I was using the R-Coda package with R version 2.4.1 and everything was going OK. After upgrading to R v.2.5.0, I receive the following error message when I'm trying to run codamenu(): "Error in coda.options(default = TRUE) : cannot change value of

I have been trying to write my own user defined function in Rpart.I imitated the anova splitting rule which is given as an example.In the work I am doing ,I am calculating the concentration index(ci) ,which is in between -1 and +1.So my deviance is given by abs(ci)*(1-abs(ci)).Now when I run rpart incorporating this user defined function i get the following error message: Error in

I am having not a little difficulty in using R-CODA version 0.4 with R-Windows 0.64.2. Basically I want to produce some postscript plots from the updates. It seems the WinBugs output is OK for cut-and-paste operations but don't allow much flexibility for anything else. So far, I can produce the .ind and .out files from WinBugs 1.2 and these seem OK (in particular they are text files not

Dear R users: I am using the package coda (the last verison in CRAN) to analyse the output from a MCMC Bayesian analysis. And I get unconsitented results. I have export the chain using the read.table function and after I have transformed this data frame to an mcmc object using the mcmc function. I am interested in three variables, when I use the function effectiveSize I have these figures:

Dear R users, I am having problems in trying to run R2jags on a Linux Cluster. When I tried to run the model given in the R2jags manual I get the following error error message: Error in FUN(X[[3L]], ...) : object 'J' not found Can any one help me with this problem? Looking to hear from you all. Regards, -- Luwis Diya, Leuven Biostatistics and Statistical Bioinformatics Centre

Dear R users I am trying to load the package lmer4 but it seems the Matrix library is failing to load and giving me the following error message : >library(Matrix) Error in importIntoEnv(impenv, impnames, ns, impvars) : object 'Logic' is not exported by 'namespace:methods' Error: package/namespace load failed for 'Matrix' I have even reinstalled the Matrix

My reply seems not to have made it to either list 24 hours after I sent it, so I'm sending it again. Apologies if you get this twice (or 4 times if you're on both lists, or 6 times if you're John) Martyn -----FW: RE: [R] First (well second) CODA hurdle problem----- Date: Thu, 05 Aug 1999 09:40:01 +0200 (CEST) From: Martyn Plummer <plummer at iarc.fr> To: John Logsdon

The document R-exts contains the following example of using mathematics in R documentation. \deqn{p(x) = {\lambda^x\ \frac{e^{-\lambda}}{x!}} {p(x) = lambda^x exp(-lambda)/x!} There is a syntax error in there, but that's not my point. The problem is that using "R CMD Rd2dvi" I find that putting the alternate forms of the equation on top of each other doesn't work.

HI Everybody Does anyone know of documentation about different ways of obtaining user input in R. I have used readline() but I wondered is there are sophisticated packages that does things like validate answers or generate selection lists. bets regards Christaan [[alternative HTML version deleted]]

Hello I want to use ntlm_auth together with winbindd for RADIUS-authentication of users against the users of a Samba-server. The freeRADIUS-daemon and the Samba-, winbindd- and ntlm_auth binaries are all on the same machine. Samba works fine and the whole setup worked fine with Samba 3.0.22. Actually I have to set up this scenario on a machine with Samba 3.0.24 and it does not work at all :-(.

I have successfully installed coda in R, and I am also able to read output from a Bugs run into coda. When I try to plots from CODA Output Analysis Menu, I do however _get_ problems: CODA Output Analysis Menu ************************* 1:Plots 2:Statistics 3:List/Change Options 4:Return to Main Menu Selection: 1 Error: "missing" illegal use of missing Quitting CODA.... Do you

Hi, I've got a linear mixed model created using lmer: A6mlm <- lmer(Score ~ division + (1|school), data=Age6m) (To those of you to whom this model looks familiar, thanks for your patience with this & my other questions.) Anyway, I was trying this to look at the significance of my fixed effects: A6post <- mcmcsamp(A6mlm, 50000) library(coda) HPDinterval(A6post) ..but I got this

Hi All, I am trying to figure out how to use R-Coda with the output from OpenBugs. I have installed and loaded the packages BRugs and R2WinBUGS. I have successfully run a simple Bayes model in WinBUGS using R2WinBUGS' "bugs" and have used "read.bugs" to build the coda object. I can successfully switch to OpenBugs and run the same model and get the basic summary

The object is a mcmc sample from lmer. I am using R v2.7.1. Please let me know what additional information I can provide, hopefully I am just making a simple mistake. Thanks in advance! > data(ratdrink, package = 'faraway') > rd.er <- lmer(wt ~ weeks*treat + (1 | subject), data = ratdrink) > rd.mc <- mcmcsamp(rd.er, 10000) > library(coda) Loading required package:

Dear all, I just wrote a script to have each of my three mcmc chins running on a different computer core to improve computation speed. To do it I use the function sfLapply from the package snowfall. Before using parallell computation, I did my diagnostic analyses just after the model updating, but now I have to save the coda before. Using the CODA package, it is possible to load the outputs in a

Dear all, I'm not 100% sure if this question is best directed at the r-list, or a mailing list concerned with Bayesian analysis, so please accept my apologies if another audience may be more appropriate. I have been using the rjags package to run Jags models with multiple chains and store the results in a Coda based mcmc list. For instance, having created a jags model and done initial