[R] MCMCglmm and iteration behaviour (new attempt)

2016-02-16 Thread Rémi Lesmerises
6.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émi

[R] MCMCglmm and iteration

2016-02-16 Thread Rémi Lesmerises
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

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

2016-02-04 Thread Rémi Lesmerises
| 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}

[R] Error with adehabitatHR and kernelbb

2013-05-17 Thread Rémi Lesmerises
- 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 Lesmerises, biol. M.Sc., Candidat

Re: [R] Error with adehabitatHR and kernelbb

2013-05-17 Thread Rémi Lesmerises
.%m.%d %H:%M:%S)) Rémi Lesmerises, biol. M.Sc., PH.D candidate Université du Québec à Rimouski De : David Winsemius dwinsem...@comcast.net À : Rémi Lesmerises remilesmeri...@yahoo.ca Cc : r-help@r-project.org r-help@r-project.org Envoyé le : vendredi 17 mai 2013

Re: [R] Error with adehabitatHR and kernelbb

2013-05-17 Thread Rémi Lesmerises
: 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 dwinsem...@comcast.net À : Rémi

[R] RDA permutest envfit

2013-04-24 Thread Rémi Lesmerises
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

Re: [R] RDA permutest envfit

2013-04-24 Thread Rémi Lesmerises
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

[R] Mancova with R

2013-04-17 Thread Rémi Lesmerises
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-project.org mailing list https

Re: [R] Mancova with R

2013-04-17 Thread Rémi Lesmerises
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

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

2013-04-16 Thread Rémi Lesmerises
to improve {coxme} (r-forge.r-project.org/scm/viewvc.php/pkg/R/dredge.R?view=logroot=mumin), but I did not understand all it works and I'm not sure it fixes my problems... Is there someone that can help me with that? Rémi Lesmerises, biol. M.Sc., Candidat Ph.D. en Biologie Université du Québec à