search for: hadfield

...urn, > nitt = iter, > family = "multinomial5", > data = data, We have no way to debug this without the data. Perhaps you should contact the maintainer and in your message attach the data? maintainer('MCMCglmm') [1] "Jarrod Hadfield <j.hadfield at ed.ac.uk>" An equally effective approach would be to post (again with data that reproduces the error) on the R-SIG-mixed-models mailing list since Hadfield is a regular contributor on that list. (To me it suggests not an error since you got output but rather a warning....

...s several minutes--and I didn't investigate whether a process that's hung in this way can be killed.) Sorry I don't have much info about the server's Samba configuration because I'm a mere user. But I'm sure my sysadmin could supply any relevant info if necessary. --- Mark Hadfield m.hadfield@niwa.cri.nz http://katipo.niwa.cri.nz/~hadfield National Institute for Water and Atmospheric Research

...gt; > family = "multinomial5", > > data = data, > > We have no way to debug this without the data. Perhaps you should contact > the maintainer and in your message attach the data? > > maintainer('MCMCglmm') > [1] "Jarrod Hadfield <j.hadfield at ed.ac.uk>" > > > An equally effective approach would be to post (again with data that > reproduces the error) on the R-SIG-mixed-models mailing list since > Hadfield is a regular contributor on that list. (To me it suggests not an > error since you got o...

Hi, Thanks in advance for any help on this question. I'm running multinomial models using the MCMCglmm package. The models have 5 outcome variables (each with count data), and an additional two random effects built into the models. The issue is that when I use the following code, the summary only gives me results for four of the outcome variables. Here is the code for my model: m3.random

...log. I didn't manage to apply the backout patch (didn't try very hard) but I fixed the problem by downgrading to 2.5.1. Lapo, does Cygwin rsync-2.5.2-1 suffer from this bug? Based on your message I should think it would. I strongly advise against upgrading until the bug is fixed. --- Mark Hadfield m.hadfield@niwa.co.nz http://katipo.niwa.co.nz/~hadfield National Institute for Water and Atmospheric Research

...all there in the book. HTH, Simon. Simon Blomberg, PhD Depression & Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ Simon.Blomberg at anu.edu.au +61 (2) 6125 3379 > -----Original Message----- > From: Jarrod Hadfield [mailto:jarrod.hadfield at imperial.ac.uk] > Sent: Wednesday, 23 July 2003 5:52 PM > To: r-help at stat.math.ethz.ch > Subject: [R] animal models and lme > > > > -- > Hi R users, > > I'm trying to fit an animal model using lme and cannot fit the > corelati...

Hi Everyone, I would simply like to coerce a character string into an expression: something like: as.expression(paste(letters[1:3], collapse="+")) but I can't seem to get rid of the quotes. The only way I can get it to work is using as.formula: as.expression(as.formula(paste("~", paste(letters[1:3], collapse="+")))) but this requires the expression to

Dear All, I have spatial data (presence/absence for 4000 squares) on 250 bird species and would like to use a model-based clustering technique to test for species associations. Is there any way of passing a distance/correlation matrix to mclust as with hclust, rather than the actual data? Or alternatively, is there a way of getting mclust to handle binary data? I'd appreciate any

...ed, please accept my apologies. Roger. -- ** R.E.Wolff at BitWizard.nl ** https://www.BitWizard.nl/ ** +31-15-2049110 ** ** Delftechpark 11 2628 XJ Delft, The Netherlands. KVK: 27239233 ** f equals m times a. When your f is steady, and your m is going down your a is going up. -- Chris Hadfield about flying up the space shuttle.

...(G1 = list(V = diag(2), nu = 5), G2 = list(V = diag(2), nu = 5), G3 = list(V = diag(2), nu = 5), G4 = list(V = diag(2), nu = 5), G5 = list(V = diag(2), nu = 5))) This is based on Hadfield's Course Notes, as well as some advice found in this post <https://stackoverflow.com/questions/40617099/mcmcglmm-binomial-model-prior>. It's consistent with how I've specified priors for simpler models (with single outcome variables), but I am obviously missing something that must...

...5), > G2 = list(V = diag(2), nu = 5), > G3 = list(V = diag(2), nu = 5), > G4 = list(V = diag(2), nu = 5), > G5 = list(V = diag(2), nu = 5))) > > This is based on Hadfield's Course Notes, as well as some advice found in this > post > <https://stackoverflow.com/questions/40617099/mcmcglmm-binomial-model-prior>. > It's consistent with how I've specified priors for simpler models (with > single outcome variables), but I am obviously missing...

...ll then be: Va 0 0 Vm or possibly, Va Cov(a,m) Cov(m,a) Vm I am quite new to both mixed effect models and R so would like to know if it is possible to specify specific variance covariance structures and whether non-positive-definite matrices can be used. Many thanks Jarrod Hadfield.

Dear All, I was wondering whether any one knows of a matrix bending function in R that can turn non-positive definite matrices into the nearest positive definite matrix. I was hoping there would be something akin to John Henshall's flbend program (http://agbu.une.edu.au/~kmeyer/pdmatrix.html), which allows the standard errors of the estimated matrix elements to be considered in the

Hello all, In a previous posting regarding glmm.admb it is stated that glmm.admb can handle 2 nested random effects. I can only fit a single random term at the moment, and wondered if anyone could provide me with some information on how to specify a model with 2 (nested or cross-classified) random terms? Thanks, Jarrod.

Hi, I'm trying to compare models, one of which has all parameters fixed using offsets. The log-likelihoods seem reasonble in all cases except the model in which there are no free parameters (model3 in the toy example below). Any help would be appreciated. Cheers, Jarrod x<-rnorm(100) y<-rnorm(100, 1+x) model1<-lm(y~x) logLik(model1) sum(dnorm(y, predict(model1),

Hi Everyone, I'm trying to load some C++ code using dyn.load but I'm getting unresolved symbols associated with some external libraries (CSparse). I gather this is something to do with linking as the the code compiles fine. However, I've passed -L/home/jarrod/My_Programs/SuiteSparse/CSparse/Lib -lcsparse to the complier (g++), either directly using R CMD SHLIB or as

Hi, Does anyone know an easy way of retaining rows in a model.matrix where missing values are present in the predictors. Ideally I'd be able to retain these rows as zeros. Thanks, Jarrod -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

Hi, I was wondering whether there was an easy way to change the attributes of all elements in a data.frame (rather than looping through elements)? Specifically, I would like to set the "dim" attributes to NULL Thanks for any help, Jarrod -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

Hi, I am having trouble getting \alt (or \altenv) to work with Schunk/Sinput and was wondering if anybody had had success? With the slide \begin{frame}[fragile]\frametitle{Basic R} \alt<2>{ <<echo=TRUE>>= 2+2 @ }{ <<echo=TRUE, eval=FALSE>>= 2+2 @ } \end{frame} I get the error message: ! FancyVerb Error: Extraneous input `> 2+2 \end {Sinput} \end

Hi, I am sorry if this is perceived as a C++ question rather than an R question. After uploading an R library to CRAN (MCMCglmm) the C++ code failed to pass the memory checks. The errors come in pairs like: Mismatched free() / delete / delete [] at 0x4A077E6: free (vg_replace_malloc.c:446) by 0x144FA28E: MCMCglmm (MCMCglmm.cc:2184) Address 0x129850c0 is 0 bytes inside a block of size 4