Podając model skrzyżowanych efektów mieszanych, próbuję uwzględnić interakcje. Jednak pojawia się następujący komunikat o błędzie:
Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) :
nlminb problem, convergence error code = 1
message = iteration limit reached without convergence (10)
Model ma następujące cechy: 1. 3 typy dysz (efekt stały) 2. 5 operatorów, każdy z 3 powtarzającymi się pomiarami przepływu paliwa z 3 typów dysz.
Poproszono mnie o uwzględnienie w modelu interakcji między typem dyszy a operatorem. Oto mój kod dla modelu:
flow.lme <- lme(rate ~ nozzle, error= nozzle|operator, data=Flow)
Dlaczego dostaję ten komunikat o błędzie?
mixed-model
lme4-nlme
Tal Bashan
źródło
źródło
operator|nozzle
losu?Odpowiedzi:
Nie słyszałem o tej
error
argumentacjilme
i nie widzę jej w dokumentacji. Jesteś pewien, że to nie literówka? Ale aby odpowiedzieć na zadane pytanie:Próbować
?lmeControl
Ustawianie
maxIter
,msMaxIter
,niterEM
i / lubmsMaxEval
argumenty do wartości wyższych niż domyślne może rozwiązać ten problem. Przechwyć dane wyjściowelmeControl
do obiektu, a następnie przekaż ten obiekt docontrol
argumentulme
.Lub...
Nowe domyślne ustawienia optymalizatora
lme
są niestabilne. Połowę czasu tego rodzaju problemy są dla mnie rozwiązywane, kiedy zmieniam je z powrotem na stary optymalizator. Można to zrobić poprzez ustawienieopt
argumentu dolmeControl
celu'optim'
.Tak więc, składając to w całość:
źródło
lmeControl
jest to funkcja znlme
pakietuPo pierwsze, jest to model ANOVA, a nie model mieszany.
Po drugie, wydaje mi się, że twój model nie został zidentyfikowany. W formie równania masz
Ostatni termin ma 15 oddzielnych wartości dla 15 twoich obserwacji. Nie ma już żadnych stopni swobody, aby uzyskać inne warunki w modelu. Włączenie interakcji było kiepską radą. Będziesz musiał je porzucić; nawet włączenie ich jako skrzyżowanych efektów nie pomoże, ponieważ będą one idealnie współliniowe z ustalonymi efektami i nie będą możliwe do oszacowania. Maksymalne prawdopodobieństwo lub model REML z 15 obserwacjami nie ma sensu; asymptotyczne wyniki teorii największego prawdopodobieństwa po prostu nie działają: to Ferrari, którym próbujesz jechać na zaoranym polu.
źródło
random=~nozzle|operator
zamiastrandom=~1|operator
.this-does-not-have-to-be-a-mixed-model
tagu. W około dwóch trzecichmixed-models
pytania, które widzę, powiedzenie czegoś takiego jest częścią mojej odpowiedzi.