Трудност при конструирането на доверителен интервал на профила в lmer модел

Когато пасвам на по-мерен модел с моите данни, няма предупредително съобщение. Но когато се опитах да конструирам доверителен интервал от confint, той показва следното предупредително съобщение:

Warning messages:
1: In FUN(X[[i]], ...) : non-monotonic profile
2: In nextpar(mat, cc, i, delta, lowcut, upcut) :
 Last two rows have identical or NA .zeta values: using minstep
3: In FUN(X[[i]], ...) : non-monotonic profile

Има ли някакъв начин да се преодолее предупредителното съобщение, т.е. има ли някакъв аргумент (опция), който може да премахне проблема? Използвах аргумента control=lmerControl(optCtrl=list(maxfun=20000)):

fit <- lmer(Y~X+Z+X:Z+(X||group),data=dat,control=lmerControl(optCtrl=list(maxfun=20000)))

person user81411    schedule 08.08.2015    source източник
comment
Най-вероятно вашият модел не може да бъде идентифициран, така че без данни няма да получите добър отговор. И двойното || обърква ме. Сигурен ли си? Като второстепенна точка, която не е от значение за въпроса, X + Z +X:Z може да се запише като X*Z.   -  person Dieter Menne    schedule 08.08.2015
comment
@DieterMenne да, сигурен съм за ||, тъй като ковариацията между произволно пресичане и произволен наклон е нула.   -  person user81411    schedule 08.08.2015
comment
Помага въпросите ви да са самостоятелни. Мисля, че вероятно има код за симулация/репликация за това в предишен въпрос. Можете ли да направите връзка към него или (за предпочитане) да копирате кода, за да направите този пример възпроизводим?   -  person Ben Bolker    schedule 09.08.2015
comment
@BenBolker, моля, вижте stackoverflow.com /questions/31796299/   -  person user81411    schedule 09.08.2015
comment
@BenBolker и в тази връзка намерих само тези три предупредителни съобщения, когато стартирам кода за комбинацията fit <- replicate(1000,noncoverage(30,20,1,.3,.3,.3,(1/18),0,(1/18))). за друга комбинация, като fit <- replicate(1000,noncoverage(30,5,1,.3,.3,.3,(1/18),0,(1/18))), получавам повече от 50 предупредителни съобщения.   -  person user81411    schedule 09.08.2015