2016-05-24 7 views
0

Ich habe versucht, das Problem zu lösen Lesen anderer Antworten, aber nicht die Lösung bekommen. Ich führe ein lmer Modell:Ich kann Fehler nicht verstehen lmer

MODHET <- lmer(PERC ~ SITE + TREAT + HET + TREAT*HET + (1|PINE), data = PRESU).

Perc ist der Prozentsatz der Prädation. Site ist eine kategoriale Variable, die ich als Blockierungsfaktor verwende. Es ist Site-Identität, wo ich das Experiment durchgeführt habe. TREAT ist kategoriale Variable von 2 Ebenen. HET ist eine kontinuierliche Variable. Die Anzahl der Beobachtungen ist 56 in 7 Standorte unterteilt Vielleicht ist das Problem, wie ich den Zufallsfaktor ausgedrückt. In jedem Standort wählte ich 8 Kiefern unter 15 aus, um das Experiment durchzuführen. Ich habe die Kieferidentität als kategorischen Zufallsfaktor aufgenommen. Zum Beispiel in Standort 1 Kiefern a1 genannt werden, a3, a7 ecc, während in Park 2 b1 genannt werden, b4, b12 ecc ... Die Ausgabe des Modells ist

Error: number of levels of each grouping factor must be < number of observations

ich nicht verstehe, wo der Fehler ist. Könnte es sein, wie ich die Kiefern nannte? Ich habe auch versucht

MODHET <- lmer(PERC ~ SITE + TREAT + HET + TREAT*HET + (1|SITE:PINE), data = PRESU)

aber die Ausgabe ist das gleiche. Ich hoffe, dass ich meine Probleme gut erklärt habe. Ich lese in diesem Forum ähnliche Fragen dazu aber ich bekomme immer noch keine Lösung.

Vielen Dank für Ihre Hilfe

Antwort

0

Verwenden Argument control = lmerControl(check.nobs.vs.nRE = "ignore") in Ihrem lmer-Aufruf, diesen Fehler zu unterdrücken. Ich schätze, das löst das eigentliche Problem nicht. Es scheint mir, dass Ihre Gruppierungsebene keine "Gruppen" enthält, wahrscheinlich ist "SITE" Ihr zufälliger Abschnitt?

Wenn Sie PINIEN verschachtelt als „Subjekte“ innerhalb SITES betrachten, dann würde ich folgende Formel vorschlagen:

MODHET <- lmer(PERC ~ TREAT*HET + (1|SITE), data = PRESU) 

oder

MODHET <- lmer(PERC ~ TREAT*HET + (1 | SITE/PINE), data = PRESU) 

Aber meine Antwort falsch sein, ich bin nicht Sicher, ob ich genug Informationen habe, um zu verstehen, worauf Sie abzielen.

bearbeiten: Sorry, Verschachtelung wurde nicht korrekt angegeben, ich habe es in der obigen Formel behoben. Siehe auch this answer.

+0

Hej Daniel, vielen Dank für Ihre Antwort. Ich werde versuchen, das Ziel des Experiments besser zu erklären. Ich möchte die Prädationsrate an 7 verschiedenen Standorten testen (Site 1,2,3, ecc ...). In jeder Site habe ich zufällig 8 Kiefern ausgewählt. Diese Kiefern wurden unter 15 für ein vorheriges Experiment ausgewählt. Sie heißen also (a1, a3, a7 für Site 1, b2, b6, b7 für SIte2 usw.). Sitte und Pines sind als kategorisch definiert. In jedem Standort wurden die Kiefern durch zwei Behandlungen geteilt (Holzpräsenz oder nicht). Diese mehr Informationen können vielleicht helfen. Da Kiefern "innerhalb" jeder Seite sind, denke ich, dass die zweite funktionieren könnte. Ich werde es versuchen. –

+0

Entschuldigung, ich hatte einen Fehler in meinem Formular und habe es behoben. Wenn Sie 'TREAT * HET' angeben, wird dies außerdem zu' TRET + HET + TREAT: HET' erweitert, so dass Sie 'TRET + HET' nicht hinzufügen müssen. – Daniel

+0

Hej Daniel, danke für die Hilfe. Ich habe das Problem gelöst! –