similar to: Cross classified or Multiple membership or Hierarchical (3 level ) logistic models using Umacs

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

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

Dear all, I recently started having some problemd with the coda package. I have also deleted it and then installed the package again. The problem is that I cant access the codamenu. That is ; ********************************************************************************************************************************************* > library(coda) > codamenu() Error in

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

We were promised this package last spring but I can't find it anywhere! Does anyone have any info? Thanks. >From RNews: "Umacs (Universal Markov chain sampler) is an R package (to be released in Spring 2006) that facilitates the construction of the Gibbs sampler and Metropolis algorithm for Bayesian inference ..." -- View this message in context:

Hello all... I am just starting to teach myself Bayesian methods, and am interested in learning how to use UMacs. I've read the documentation, but the single example is a bit over my head at the level I am at right now. I was wondering if anyone has any simple examples they'd like to share. I've successfully done a couple of simple gibbs examples, but have had a hard time

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

I have downloaded the Umacs (Universal Markov chain sampler) and submitted the following sample code from Kerman and Gelman.   s <-Sampler( J=8, sigma.y  =c(15,10,16,11,9,11,10,18),           y  =c(28, 8,-3,7,-1,1,18,12),      theta =Gibbs(theta.update,theta.init),           V =Gibbs(V.update,mu.init),         mu =Gibbs(mu.update,mu.init),         tau =Gibbs(tau.update,tau.init),       

This is a question about how to calculate similarities/distances among items that are classified by hierarchical attributes for the purpose of visualizing the relations among items by means of clustering, MDS, self-organizing maps, and so forth. I have a set of ~260 items that have been classified using two sets of hierarchically-organized codes on the basis of form and content. The data looks

https://bugzilla.mindrot.org/show_bug.cgi?id=1753 Summary: Use -funroll-loops with umac.c Product: Portable OpenSSH Version: -current Platform: Itanium OS/Version: Other Status: NEW Severity: enhancement Priority: P2 Component: Build system AssignedTo: unassigned-bugs at mindrot.org

https://bugzilla.mindrot.org/show_bug.cgi?id=1462 Summary: Unaligned access warnings on IA64 when using umac-64 Classification: Unclassified Product: Portable OpenSSH Version: 5.0p1 Platform: Itanium2 OS/Version: Linux Status: NEW Severity: normal Priority: P2 Component: ssh AssignedTo:

Are there any well developed Classified Ads Script for Rails?

Hi All, I''m not convinced I have my tcrules correctly setup and looking at the counters in the mangle table''s tcpost doesn''t really help much as I can''t tell what is the final match. Is there a way to match packets in iptables based on the classifier? i.e. so I can LOG packets classified with 1:18 for example. I can''t see anything in iptables, except

Dear R-users, I would like to fit a glmm with cross-classified random effects with the function MCMCglmm. Something along the lines: model1<-MCMCglmm(response~pred1, random=~re1+re2, data=data) where re1 and re2 should be crossed random effects. I was wondering whether you could tell me specifying cross-classified random effects in MCMCglmm requires a particular syntax? Are there any

All, I have several years of longitudinal test scores for students (many who switch schools at various points in time). I am using a mixed-effects model with crossed random effects to model student trajectories. The model includes time at level 1 and students crossed with schools at level 2. When I run the model I get the posterior variances on the intercepts and slopes for students and schools,

Dear All, The "coxph" function in the "survival" package allows multiple frailty terms. In all the examples I saw, however, the frailty terms are nested. What will happen if I have non-nested (that is, cross-classified) frailties in the model? Will the model still work? Do I need to take special cares when specifying these models? Thanks! Shige [[alternative HTML

Hello, this is probably a dumb question, but I can't for the life of me figure out why sendmail and courier are classified universally as network-mail-mta whereas dovecot is in places like network-mail-general or some such place. I've used several Linux distros and have noted this keenly (they also normally classify photorec, the HD recovery tool as hardware-photography,so this may just be

Enlightenment succeeded! Thanks, David