Nested LMM (Hiyerarşik Karma Model)
v1.0.1
Engineering context — fixed effect (treatment) + nested random effects. urun_dayanim_MPa ~ tedavi + (1|fabrika/uretim_hatti/parti). Difference from VARCOMP: a fixed effect (e.g. treatment/group comparison) is added.
🎯 What is it for?
When interpreting fixed effects in hierarchical data, it also accounts for the variance of the random factors. In Engineering, if the same batches receive different treatments, a classical t-test/ANOVA ignores the clustering → type-I error. Nested LMM computes this correctly.
📌 When is it used?
- Comparison of a fixed effect (treatment, intervention, condition) + hierarchical sampling
- Multi-site clinical research, multi-center field trial
- Adjusting site/group-level variance in RCTs
⚙ Assumptions
- Continuous DV, categorical fixed effect.
- Nested structure: uretim_hatti within fabrika, parti within uretim_hatti.
- Random effects ~ Normal(0, σ²).
- ≥1 observation in each treatment-parti combination.
📊 How to Run It in MerQur
Panel assignments (form fields in the program):
- Dependent Variable:
deger - Replication Variable:
blok - Upper Group:
ust_grup - Inner (Nested) Group:
alt_grup_no - SS Type:
1
🧪 Example Data — Engineering
📂 samples/Ileri_Duzey_v102/muhendislik/108_nested_lmm.xlsx · n = 216 · Treatment: Control/Treatment-A/Treatment-B · Hierarchy: fabrika/uretim_hatti/parti
⬇ Download sample data set (.xlsx)
| fabrika | uretim_hatti | parti | tedavi | urun_dayanim_MPa |
|---|---|---|---|---|
| FABRIKA_01 | URETIM_HATTI_01_01 | PARTI_01_01_01 | Kontrol | 355.218 |
| FABRIKA_01 | URETIM_HATTI_01_01 | PARTI_01_01_01 | Tedavi-A | 376.136 |
| FABRIKA_01 | URETIM_HATTI_01_01 | PARTI_01_01_01 | Tedavi-B | 374.672 |
| FABRIKA_01 | URETIM_HATTI_01_01 | PARTI_01_01_02 | Kontrol | 364.527 |
| FABRIKA_01 | URETIM_HATTI_01_01 | PARTI_01_01_02 | Tedavi-A | 353.466 |
| … | … | … | … | … |
📊 MerQur Output (summary)
==================================================
Fixed effect (tedavi) β̂ SE p
Kontrol (ref) 0.00 – –
Tedavi-A +4.0 1.20 < .001 Tedavi-B +7.0 1.20 < .001 Random components σ² fabrika 324.00 uretim_hatti 144.00 parti 64.00 Residual 100.00
📄 APA 7 Yorumu
⚠ Common Mistakes
- Specifying “crossed” instead of “nested” — the uretim_hatti ids must not repeat across each fabrika.
- Optimizer stalling (lbfgs llf=inf) — MerQur resolves it with a bfgs/cg cascade.
- Treatment × fabrika imbalance — the fixed-effect coefficient may be biased.
📚 Cite MerQur
Örücü, Ö. K. (2026). MerQur: An Integrated Platform for Academic Data Analysis and Reporting [Computer Software] (Version 1.0.0). https://doi.org/10.53463/merqur.2026001
📝 Production Note — The sample data on this page was generated synthetically (fixed SEED=42, generator: samples/Ileri_Duzey_v102/_generate_v102_datasets.py). The page content was prepared with the support of Anthropic Claude, and its academic accuracy was verified by the author.