2017-08-07 3 views
1

Ich arbeite mit einem gemischten Modell, um die Auswirkungen von Variablen zu sehen. Der Code, den ich verwenden ist:Warum gibt SAS keine p-Werte für die fixierten Effekte im gemischten Modell?

proc mixed data=pb2; 
class treat_a treat_b hoknr_ day; 
model conc=treat_a|treat_b hoknr_/outp=residuals1 residual; 
repeated day/subject=hoknr_(treat_a treat_b)type=vc; 
run; 

Das Ergebnis hat keine p-Werte für treat_a, treat_b oder treat_a | treat_b aber es hat für hoknr_. Ich habe die wiederholte Aussage ausgeschlossen, ich habe das Modell vereinfacht, ich habe die Klasse gewechselt, aber ich habe immer noch keine p-Werte für alle meine fixen Effekte. Ich habe dieses Modell schon früher benutzt und es funktionierte, jetzt passte ich es zu diesem Datensatz an und es funktioniert nicht vollständig.

Edit1 Die Tabelle des Typs 3 Tests festen Effekte zeigen wie folgt aus: Type 3 Table. Die Behandlungen könnten nicht schätzbar sein (treat_a ist ja oder nein, ebenfalls für treat_b). Ich habe Ja/Nein zu 0 oder 1 geändert, die Typ 3-Tabelle nicht geändert. Ich habe schon früher mit Behandlungen gearbeitet, die in Worten ausgedrückt wurden, was nicht zu einem Tisch wie diesem führte.

Edit2 Wenn der Modellanweisung eine Lösung hinzugefügt wird, ist dies das Ergebnis: Solution for Fixed effects.

Was ist falsch an diesem Modell, dass es keine P-Werte für alle festen Effekte zeigt?

Antwort

2

Möglicherweise müssen Sie mit der Option htype= den gewünschten Testtyp angeben. Es klingt wie einer dieser Prozesse, bei denen jemand die Funktion zunächst nicht programmiert hat, und es war ein nachträglicher Einfall spät in der Entwicklung (nicht unähnlich der showpvalues Option in proc glmselect; bis heute denke ich, das ist die seltsamste Option in einem Regressionsprozess) .

https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect015.htm#statug.mixed.mixedmodelhtype

Typ 3 Tests festen Effekte

Sie die HTYPE = Option in der MODEL-Anweisung Tabellen erhalten von Typ 1 (sequenziell) Tests und Typ 2 (bereinigtes können) testet zusätzlich zu oder anstelle der Tabelle von Typ 3 (Teil-) Tests.

Die ODS-Tabellennamen lauten "Tests1" für die Tests des Typs 1, "Tests2" für die Tests des Typs 2 und "Tests3" für die Tests des Typs 3.

Oder es könnte sein, dass einige Ihrer festen Effekte nicht schätzbar sind.

+0

Vielen Dank für Ihre Antwort. Ich habe der Frage basierend auf Ihrer Antwort eine Bearbeitung hinzugefügt. –

+0

Ist Ihr Modell nicht vollständig? d. h. gibt es mehr Design-Spalten als Freiheitsgrade für Ihre fixen Effekte? Versuchen Sie, die 'Lösung'-Option in Ihrer Modellanweisung zu verwenden und zeigen Sie uns, was Sie bekommen. –

+0

Eine zweite Änderung bezüglich Ihrer zweiten Antwort wurde hinzugefügt. –

Verwandte Themen