Fitting mixed models in R (with lme4)

Learning Objectives:
* Understand lmer syntax (fixed, random, cluster)
* Understand how to interpret fixed effect parameters
Here's the dataset I'm using: quantpsych.net/data/jedi.csv
Here's my video that shows you how to identify your cluster variable: • How to identify your c...
Here's a video about how to determine whether an effect is fixed or random: • How to decide whether ...
Link about EDA versus CDA: • Ethics in Statistics P...
My Multivariate playlist: • Multivariate Statistics
And here's a paper I wrote about my eight step approach to data analysis: psyarxiv.com/r8g7c/
Undergraduate curriculum playlist (GLM-based approach): kzread.info?list...
Graduate curriculum playlist (also GLM-based approach): kzread.info?list...
Exonerating EDA paper: psyarxiv.com/5vfq6/
Download JASP (and visual modeling module): www.jasp-stat.org

Пікірлер: 39

  • @QuantPsych
    @QuantPsych Жыл бұрын

    Yes, I know my head is in the way of the output. Sorry! But you should still get an idea of how to do these things in R.

  • @igorcastro6776

    @igorcastro6776

    Жыл бұрын

    When will you upload the next video? I'm desperately waiting for it!

  • @vinitalec
    @vinitalec29 күн бұрын

    Your videos are excellent!Thanks for helping me understand this subject.

  • @martinabautista
    @martinabautista5 ай бұрын

    You're funny. I feel glad to came across your content!! I'm using lme4 for my dissertation project

  • @msimister4261
    @msimister4261 Жыл бұрын

    This was extremely helpful... Thank you so much for this!

  • @farmz0r
    @farmz0r Жыл бұрын

    07:40 unfortunately your head is blocking the summary table you are talking about :d 5 referred to fixed effect intercept estimate 3.261 referred to random effects residual Variance 1.8 referred to random effects residual SD cba to do it for the next table atm

  • @mykiawiggins3318
    @mykiawiggins33186 ай бұрын

    Not a fistbump...😳You punched my face!🤣😂 Thank goodness for your videos you're so frrreakin FUNNIEEEEE!

  • @chrislloyd5415
    @chrislloyd5415 Жыл бұрын

    Interested in your view about the following question. Your third model allows slope(darkness,anger) to depend on jedi_id. You could also just add an interaction term darkness*factor(jedi_id) right? Assuming this model is identifiable (and I think it is because anger is varying over time within jedi_id) what is the value of the random effects model which makes the arbitrary assumption that the jedi-specific slopes follow a normal distribution?

  • @EW-to9sr
    @EW-to9sr Жыл бұрын

    Hello, thanks for uploading these tutorial videos. I'm a uni student trying to understand statistics in R and I found your channel is extremely helpful, thanks! Besides, I'd like to know is it reasonable to consider a factor as fixed and random at the same time? Looking forward to seeing you how explain the mixed model in following videos!

  • @QuantPsych

    @QuantPsych

    Жыл бұрын

    Yes. All random effects also have fixed effects. You cannot have a random effect without a fixed effect.

  • @zimmejoc
    @zimmejoc Жыл бұрын

    In the CSV file, the variable is named dojo_id not jedi_id. I always tell my students, "KNOW THY DATA" and that applies to me too. Bad zim for blindly typing what was in the video instead of verifying with the data first :)

  • @QuantPsych

    @QuantPsych

    Жыл бұрын

    Good catch!

  • @EmmaCalikanzaros
    @EmmaCalikanzaros5 ай бұрын

    Thanks, I found your video super helpful. Could you explain what it means when the error "singular model" apprears?

  • @nachete34
    @nachete34 Жыл бұрын

    As always, love ur videos, specially if they involve R. However, while it did not prevent me from following you, your face camera occluded the results window..just being picky :) Different topic...do u consider uploading some machine learning tutorial for dummies?

  • @Lello991
    @Lello991 Жыл бұрын

    Hi prof Fife, I have a question for you: let's assume you want to add age_started as a predictor of darkness. Of course age_started doesn't vary within each level of the jedi_id cluster variable (i.e., each jedi has just a unique value of age_started). Could you add it with no problems? Is it possible to use age_started as fixed effect + jedi_id as random effect together? Or would you encounter some separation / convergence issues? I'm struggling to understand this point but I feel kind of lost. Thank you for your enlightening videos!

  • @elena.s.v.
    @elena.s.v.5 ай бұрын

    Hello! Thank you for your videos, they help a lot. Would it be possible to get more information about the warnings one gets when fitting a random effect model, called "full" in your videos (e.g. singularity, did not converge, etc.) and how to solve them? I used the model comparison to verify whether I should fit a variable as a fixed effect only or as a fixed/random effect and some of my "full" models do not converge or "are singular". I also get warnings when trying to fit my final models (linear mixed model and generalized mixed models (poisson and negative binomial). Thank you for your help!

  • @ALI_B
    @ALI_B9 ай бұрын

    awesome video. I have a question : I have two clustering variables for a finance dataset where data is nested in banks and scores are time series reported quarterly. How to include the variable date in the script for the fixed model : fixed_y = lmer(y ~ 1 + x1 + x2 + x4 + x4 + (1 | Bank))

  • @QuantPsych

    @QuantPsych

    4 ай бұрын

    Good question. I think the proper notation is "...(1|Bank:Time)"

  • @marcellberto2538
    @marcellberto2538 Жыл бұрын

    Head's in the way bro! 😜 Still luvs ya videos though, as always thanks for sharing 😊

  • @nikashomeil2369
    @nikashomeil23693 ай бұрын

    Thank you so much for your wonderful video, I have a question, I have a model in which my fixed effect has 3 levels but when I ask for the summary one of the levels is not reported in the fixed effect, do you have any idea why is that?

  • @trueperson22
    @trueperson22 Жыл бұрын

    Thank you.. If you can explain the difference between the model with ~1 and the last one, I 'd really appreciate that..

  • @QuantPsych

    @QuantPsych

    Жыл бұрын

    ~1 just means to fit an intercept. If it's omitted, R will fit an intercept anyway. The first model (baseline) has to have it because there are no predictors (i.e., fitting this will throw an error: lmer(y~ ( | id), data=d)). For the remaining models, it's redundant, but I put it there so you can easily track what I've added to the model.

  • @jorgemmmmteixeira
    @jorgemmmmteixeira Жыл бұрын

    In there a scenario where makes sense to have random slope and intercept for jedi_id? If so, how would you code it? thx

  • @QuantPsych

    @QuantPsych

    Жыл бұрын

    I'm not sure I understand your question...cluster variables don't have random slopes and intercepts. Variables do.

  • @ropflpfopfl2555
    @ropflpfopfl2555 Жыл бұрын

    Hey, i really like your video and your way of explaining things. Way clearer and simpler than other channels out here :D it is a great example on how mixed models work. However, i may have an understanding problem when it comes to the data format. At 3:40 you say that in their first year, when they are 5, they already killed someone. Does the dataset offers any time-specific variable? because the 5 stands for the point in time they started training, which should be always the same when looking at one jedi. (here jedi_id = jedi_1). Did I miss something? I'm only a poliscience student who stumbled upon that one :)

  • @QuantPsych

    @QuantPsych

    4 ай бұрын

    Ah, good point. I was getting really worried that five year olds were murdering people, but you're right. They might have waited until they were seven :)

  • @mind_palace
    @mind_palace4 ай бұрын

    Repeated measures....to make sense of it, I'd say it could be that every measure per jedi_id, could be the therapist measuring their anger level for every session?

  • @QuantPsych

    @QuantPsych

    4 ай бұрын

    That works :)

  • @ladynoluck
    @ladynoluck Жыл бұрын

    Definitely not watching this as a SAS to R transplant for my dissertation 😂

  • @taranaferdous2858
    @taranaferdous28582 ай бұрын

    Hi! Thank you for this clear explanation. I followed the same for my data (data is in the same format that you have shown). My models worked fine up to fixed part. The moment I added the random part (1 + *** | ID), I got this error : number of observations (=100) What am I missing here or doing wrong?

  • @IbrahimKwakuDuah

    @IbrahimKwakuDuah

    2 ай бұрын

    It is not 1+*, I guess the * would mean everything but you can't do everything

  • @hanswurst4728
    @hanswurst4728 Жыл бұрын

    Cool, but if we could actually see the output while you're interpreting it and not your face, that would be tremendously helpful. Informative nonetheless 👍.

  • @ceciliacocucci8288
    @ceciliacocucci8288 Жыл бұрын

    Dear Dustin I really enjoy your videos and love flexplot, but you always fit linear models. I'm a neonatologist, and almost 90% of my outcomes are binary o dichotomous. Is it possible to give usa video about non linear mixed models, aka logistic, ordinal etc.. ? I guess you'll use the nlme function but I'd really love to have your explanations. Thanks!!!

  • @QuantPsych

    @QuantPsych

    Жыл бұрын

    I'm pretty sure I made a generalized mixed model video. I think it was poisson mixed models, but it might have been logistic. Just take my logistic regression videos, combine those with my mixed models videos, and you'll get it.

  • @dle3528
    @dle35289 ай бұрын

    Your head is in front of the output 😕

  • @QuantPsych

    @QuantPsych

    4 ай бұрын

    I can't help having a fat head :)

  • @chrislloyd5415
    @chrislloyd5415 Жыл бұрын

    Why is it a 1? You reallj don't know? C'mon man! (Just call me Joe Biden). It obviously refers to a column of 1's in the design matrix.

  • @diptarshis
    @diptarshis3 ай бұрын

    Could you please stop with the clowning and get on with the explanations. Your attempt at deadpan humor sucks. You're way better at teaching, do that please !

  • @QuantPsych

    @QuantPsych

    3 ай бұрын

    No, it's so deeply ingrained in my personality. I do me, you do you. And it seems you're outvoted. Everyone else seems to like it :)