[R] MCMCglmm and iteration behaviour (new attempt)

.9136.761 13.020 333.7 <4e-04 *** small_raod -16.945 -20.694 -13.492 194.9 <4e-04 *** I'm then wandering if it is because more iteration produce better estimates and then a model that had a better fit with the data. Anyone can help me

[R] MCMCglmm and iteration

Hi everyone, I'm running a bayesian regression using the package MCMCglmm (Hadfield 2010) and to reach a normal posterior distribution of estimates, I increased the number of iteration as well as the burnin threshold. However, it had unexpected outcomes. Although it improved posterior distributi

[R] Conditional (paired) design for binomial regression in MCMCglmm

| This is my first attempt asking question in this forum and I do it because thorough research on the web didn't give me any answer.I am trying to accomodate a conditional regression in a Bayesian generalized linear mixed model using Monte Carlo Markov Chain. I am using the package {MCMCglmm} i

Re: [R] Error with adehabitatHR and kernelbb

: >coordinates(fece)<-~UTMeasting+UTMnorthin >loc@data$time <- as.POSIXct(strptime(as.character(loc@data$time),"%Y.%m.%d >%H:%M:%S")) Rémi Lesmerises, biol. M.Sc., PH.D candidate Université du Québec à Rimouski De : David Winsemius À :

Re: [R] Error with adehabitatHR and kernelbb

ce@data$time),"%Y.%m.%d >%H:%M:%S")) Rémi Lesmerises, biol. M.Sc., PH.D candidate Université du Québec à Rimouski De : David Winsemius À : Rémi Lesmerises Cc : "r-help@r-project.org" Envoyé le : vendredi 17 mai 2013 13h53 Objet :

[R] Error with adehabitatHR and kernelbb

me? ## commands BBtraj <- list() for (i in 1:(nrow(loc@data)-1)) { BBtraj[[i]] <- kernelbb(as.ltraj(loc@coords[i:(i+1),], date=loc@data$time[i:(i+1)], id = as.character(loc@data$id[i:(i+1)]), typeII = TRUE), sig1=as.numeric(as.character(loc@data$sig1[i])), sig2= 5, grid = 1000) }   Rémi

Re: [R] RDA permutest envfit

Thank you! A last question: Is it still the same explanation if I remove the condition "first=TRUE" (and then testing for all axis) and permutest gives the same result?   Rémi Lesmerises, biol. M.Sc., Candidat Ph.D. en Biologie Université du Québec à Rimouski remilesmeri...@yahoo.ca On W

[R] RDA permutest envfit

Dear all, I did a RDA and when I looked to the signification of the test with permutest, the output was non-significant. But when I used the envfit function, some of the vectors are significant. All the test's conditions are respected. What it means? Is it an error in the script? Commands an

Re: [R] Mancova with R

Dear John, Thanks for your comments! But when I tried your suggestion, the output was as the following:  Response Dist_arbre :             Df     Sum Sq    Mean Sq F value    Pr(>F)     Poids        1 0.00010398 0.00010398  6.2910 0.0364733 *   Age          1 0.520

[R] Mancova with R

o stay in R! Someone have any idea? Thanks in advance for your help!   Rémi Lesmerises, biol. M.Sc., Candidat Ph.D. en Biologie Université du Québec à Rimouski remilesmeri...@yahoo.ca [[alternative HTML version deleted]] __ R-help@r-proje

[R] Model ranking (AICc, BIC, QIC) with coxme regression

ion of interaction between two random variables (R3|Individual), and it's why I have to use coxme. I found a new update in R-forge to improve {coxme} (r-forge.r-project.org/scm/viewvc.php/pkg/R/dredge.R?view=log&root=mumin), but I did not understand all it works and I'm not s