Search Archives

Hello. I'm estimating some home ranges using the kernelUD function in adehabitatHR. For each individual, I want to estimate the degree of overlap between the home range (kernel 95%) and a SpatialPolygonsDataFrame which corresponds to a protected area. Basically...

Hi, I have a very simple Cox regression model in which I need to include a nested random effect: individual nested in treatment. I know how to pass a single random effect - e.g. frailty(id)- but how can I...

Hi. I have a script in which I export some plots to a folder in my PC. The plots are exported with the function ggsave (in ggplot2). Whenever I re-run the code and export updated plots, the name of...

Hi, I have been using the *sim *function in the arm package to simulate glm models in R and then extract random effects. mymodel =lmer(y~x1+x2+x3+(1|ID),data=mydataset) mysim=sim(mymodel,n.sims=1000) simulated...

Hi. I was wondering how data has to be structured in the database to be able to run a mixed-model with the MCMCglmm funcion using censored family distributions, like "cengaussian" or "cenpoisson". I know my response variable should have...

Hello. I ran a PCA analysis on a dataset with 5 variables and retained two components. I rotated them and now I want to predict the scores in a new data set for which the original variables are available. I...

Hello. I ran a PCA analysis on a dataset with 5 variables and retained two components. I rotated them and now I want to predict the scores in a new data set for which the original variables are available. I...

Hi, For a number of individuals, I have measured several behavioral traits in the wild. Those traits (e.g. home range) can be estimated on different temporal scales, for example daily, weekly or monthly. I want to estimate repeatability of...

Dear all, For a number of individuals, I have measured several behavioral traits in the wild. Those traits (e.g. home range) can be estimated on different temporal scales, for example daily, weekly or monthly. I want to estimate repeatability...

Dear list. I?m fitting some logistic bam models, and for some of them I get the following warning: Warning message: In bgam.fitd(G, mf, gp, scale, nobs.extra = 0, rho = rho, coef = coef, : fitted probabilities numerically 0 or...

Dear list, I?m running a gam model (package mgcv) with a binary response variable (y), and two continuous explanatory variables (x and z), plus their interaction (x:z). I, therefore, obtain four coefficients from my model (intercept, slope of...

Dear all. I'm running a gamm model. From previous tests in lme, I know that a correlation structure is needed, in this case I'm trying to use a corAR1 (the one I used in my tests in lme...

Dear list I have a conceptual question. I run a MCMCglmm model (model 1) and extracted the posterior distribution of the random effects. This is a matrix of 75 (number of individuals) x 5000 (number of iteractions of the model...

This is what I was looking for, thanks Roman! David 2016-01-20 13:44 GMT+01:00 Roman Lu?trik <roman.lustrik at gmail.com>: > If I understand your question, you are after functionality of `rgeos` > package. You can...

Hi all I?m running a MCMCglmm mixed model in which the random effect is the interaction between the individual identity and the season of the year (categorical variable, four levels). The response variable is a behavioral trait (dvm). Basically...

Dear list, I'm running a model like this using lme (nlme): model=lme(res~t1+t2+poly(month,3)+location,random=~1|ID,data=dataset,method="REML") where "res" is the response variable, "t1", "t2" and "month" are explanatory...

Hi all, I'r trying to run a multivariate MCMCglmm with 3 traits. I was suggested to use the long format since I have unequal number of replicates per trait. Trait 1 and 2 were replicated twice, trait 3 was...

Hi all, I'm running this model with this prior: prior1=list(R=list(V=diag(4),nu=1.002),G=list(G1=list(V=diag(4),nu=1.002))) mcmc1=MCMCglmm(cbind(t1,t2,t3,t4)~(trait-1)+trait...

Dear list, I am investigating the effect of the interaction between two continuous variables (A and B) on a behavioral trait. I have repeated measures from 64 individuals. The number of measures per individual varies a lot with a minimum...

Hello, I'm fitting some random regression models to investigate variation over time of a response trait. The time variable is "month", and by fitting a random regression model I want to investigate variation in plasticity across individuals, i.e...