mixed-models

specifying multiple separate random effects in nlme

筅森魡賤 提交于 2021-01-24 07:17:28
问题 I am analysing some whale tourism data and am trying to construct linear mixed effect models in the nlme package to see if any of my explanatory variables affect encounter time between whales and tourists. (I am also open to running this model in lme4 .) My variables are: mins : encounter time (response variable) Id : individual whale ID (random effect) Vessel : vessel Id (random effect) Sex : sex of the animal Length : length of the animal Year Month (nested within Year ). So my random

Singularity in backsolve at level 0, block 1 in LME model

China☆狼群 提交于 2020-12-05 11:02:33
问题 dput for data, copy from https://pastebin.com/1f7VuBkx (too large to include here) data.frame': 972 obs. of 7 variables: $ data_mTBS : num 20.3 22.7 0 47.8 58.7 ... $ data_tooth: num 1 1 1 1 1 1 1 1 1 1 ... $ Adhesive : Factor w/ 4 levels "C-SE2","C-UBq",..: 2 2 2 2 2 2 2 2 2 2 ... $ Approach : Factor w/ 2 levels "ER","SE": 1 1 1 1 1 1 1 1 1 1 ... $ Aging : Factor w/ 2 levels "1w","6m": 1 1 1 1 1 1 2 2 2 2 ... $ data_name : Factor w/ 40 levels "C-SE2-1","C-SE2-10",..: 11 11 11 11 11 11 11 11

suppressing some of the plots from plot(allEffects(model)) output

拥有回忆 提交于 2020-07-10 07:47:08
问题 I'm not sure this is possible, but I'm wondering if you can remove some of the not-as-interesting plots from the plot(allEffects(model)) output. I assume this is possible by editing some of the background mechanics that R uses to make the array of graphs. I've already used multiline to reduce the number of graphs from 20, but I'm only interested in 2 of 10 of my remaining plots because it most accurately represents one of my interaction effects from my model (generalized linear mixed model).