Filter löschen
Filter löschen

Which formula fits better for this Linear Mixed-Effects Model?

2 Ansichten (letzte 30 Tage)
Miquel Bosch Bruguera
Miquel Bosch Bruguera am 10 Mär. 2020
Beantwortet: Peng Li am 5 Aug. 2020
I am currently analyzing a dataset that contains a list of flight simulator tests performed by different pilots. I want to analyze if a certain flight parameter (i.e. amount of input errors during flight, lateral deviation to ideal path, etc.) are affected by the categorical variables in the following list:
  1. Experimental Campaign: the pilots flew the same flights but in different places (and environmental conditions, such as lack of oxygen, isolation, etc). 5 different campaigns were done. A performance difference is likely to appear depending on the campaigns.
  2. Group: in each campaign, the pilots were divided in two groups: Frequend and infrequent flyers. A performance difference is expected between both groups.
  3. Session: During each campaign, the same amount of flights were performed, each month for FF pilots, or every three months for IF Pilots. In total, 10 sessions were made. A variation of performance might happen throughout the experiment, also affected by Group and Campaign.
  4. Flight Scenario: three different flight scenarios were flown, which required different skill levels. The performance is also expected to vary between type of scenario.
Additionally, an extra list of categorical variables could be considered (Gender, Age, Background, etc.).
Could you please tell me which LME Model formula would you better implement in order to understand the dataset presented? And if you wish, how would you better plot the results of such an analysis?
Thank you!!
  1 Kommentar
Paul am 24 Jul. 2020
Bearbeitet: Paul am 24 Jul. 2020
From what you describe it seems you only have fixed effects, no random effects. Is that correct? Random effects are those factors which can vary between datapoints, and which can have an effect on your data, but about which you are not interested, for instance the ID of the pilot (one pilot might have produced multiple datapoints because he participated in multiple tests) or time of day when the measurement happened (morning/afternoon/evening) or something like that.
If you have no random effects, your formula could be this:
'FlightParameter ~ ExperimentalCampaign*Group*Session*FlightScenario'
In this case you are comparing interactions between all four fixed effects, which would produce a rather complicated model, but LMEs should be able to handle this fairly easily as long as you have enough datapoints.
How you can best plot this depends entirely on what the goal of the visualization is. What is the goal of your visualization?

Melden Sie sich an, um zu kommentieren.

Antworten (1)

Peng Li
Peng Li am 5 Aug. 2020
To understand your dataset, I still believe that descriptive statics can give you the broad general picture. In terms of statistical analysis, the hypothesis always comes first and then models, not the other way around. So get back to your question, I don't really recommend you do the linear regressions, whatever types, just blindly before you have any hypotheses in mind. And for sure that the software will throw out results whatsoever, even you design has problems. The results will then make no sense.
It might be okay to explore around. But the four-way interaction is way too complicated to be interpretable! It seems that in your question, the same pilot may have completed multiple sessions. That will be where random effects are needed.
To showcase a simpler scenario, for example you'd like to test whether frequent and infrequent pilots (group factor) perform differently in different places, and you'd like to control for demographic varations, you may want to apply this lme model: outcome ~ group * place + age + sex + background + (1|pilot)
the (1|pilot) part in the formula is to take the within-pilot correlation into consideration (random effect).




Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by