1,389 questions
Advice
0
votes
2
replies
50
views
How to interpret estimates from latent class mixed model (lcmm in R)
I'm learning how to run latent class mixed models using the r lcmm package and am struggling to interpret the estimates.
I've read this paper by Proust-Lima et al as best I can, but am still missing ...
- 3,973
0
votes
0
answers
26
views
What is the weights argument in glmer() of lme4 doing when using gamma response?
I am using a gamma response in glmer() from lme4. Some rows in my data set represents averages from n_i observations, while most observations just are individual measurements. My idea was to ...
- 87
Advice
0
votes
2
replies
52
views
What are 'weights' argument in lme4 glmer() with gamma response doing?
I am using a gamma response in glmer() from lme4. Some rows in my data set represents averages from n_i observations, while most observations just are individual measurements. My idea was to ...
- 87
0
votes
0
answers
85
views
Creating Estimation Plots in R with multiple variables
I have a mixed model that includes three different variables: Treatment, Stimulation, and DRG, as well as a random effect with DRG nested under a group ((1|Group/DRG)). I am looking to assess how each ...
Best practices
1
vote
1
replies
99
views
Calculate sample size in R - Longitudinal repeated measures study using a linear mixed effects model
I'm trying to calculate a sample size for a phychological study where the purpose of the study is to see if patients can maintain a score at baseline during 3 follow up measurements.
The mean baseline ...
- 35
-3
votes
1
answer
42
views
is mmec function of sommer package deprecated in latest version?
I just want to ask whether the function mmec is deprecated in the latest version? because when I am trying to use it, its displaying (could not find function "mmec") message. It is ...
0
votes
1
answer
30
views
`glmm.hp()` error: Error in str2lang(x) : <text>:3:0: unexpected end of input
I am working on a project that needs glmm.hp for a mixed effects logistic regression. I noticed that when I fit models that do not have too many terms, everything runs fine. When I put in all the ...
- 2,736
0
votes
1
answer
38
views
`glmm.hp()` error: Error in seq.default(1, length(varname), 2) : wrong sign in 'by' argument
library(glmm.hp)
library(partR2)
library(lme4)
data(biomass)
mod <- lmer(Biomass~(1 | Population) +
Year + Temperature + Precipitation + SpeciesDiversity, data=biomass)
glmm.hp(mod)...
- 2,736
2
votes
1
answer
164
views
glmmTMB: start a new optimization at the optimum of previous model
Suppose I have the following model:
library(glmmTMB)
Owls <- transform(
Owls,
Nest = reorder(Nest, NegPerChick),
NCalls = SiblingNegotiation,
FT = FoodTreatment
)
fit_zipoisson <- ...
- 5,495
2
votes
1
answer
100
views
Joint Modeling Survival Error “Invalid survival times for this distribution”
I’m trying to fit a joint longitudinal–survival model using the JM package. My longitudinal part is a random-slope model with nlme::lme, and my survival part is a Cox model with survival::coxph. When ...
- 410
1
vote
0
answers
34
views
Simr powerCurve no errors but 0 successes
I am trying to run a power analysis for a new experiment using simr, based on data from a similar previous experiment. The original experiment was a forced choice task (so `answer' is either 0 or 1, ...
- 85
0
votes
0
answers
90
views
Trouble in repeated measurments: lme model fails with "NA/NaN/Inf" error
I'm trying to fit a linear mixed-effects model with lme() from the nlme package to analyze repeated measures. The model runs fine on the full dataset, but after filtering by sex, I get the following ...
- 410
2
votes
1
answer
101
views
Nested random effect using lmerTest package in R
I want to build a nested linear mixed effect model using lmerTest package. The predictors are gender, student, school and neighborhood, and outcome is score. The gender and student will be taken for ...
- 511
0
votes
1
answer
55
views
Sum-coding when variables have three levels: question about interactions
The two variables of my basic lmer() model...
lmer(SCPRE ~ LANGUAGE * ORDER + (1|PID) + (1 + |WORD), data = XP_DATA)
...have three levels each:
LANGUAGE: English, Italian, Spanish
ORDER: A, B, C
...
- 13
4
votes
1
answer
100
views
glmmTMB issue with number of observations and groups
I have a dataset with 125 animals across 3 sites and 100500 observations. Both show up properly when looking at the structure of the data but when I run the model with an updated data frame (I added a ...
- 45