를 사용하여 시뮬레이션()함께 lmer()

Question

나는 두 개의 다중 로지스틱 회귀를 사용하는 동 예측,하지만에 두 개의 서로 다른 응답:

fruitMLM <- lmer(InsuffFruit ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))  
fuelMLM <- lmer(Pollution ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))

Income 은 이산으로 값 c(-2,-1,0,1,2), HDI 은 지속적인 0 과 1 사이 Country 범주 및 반응은 모두 1/0.

을 줄거리 신뢰 밴드 실행 시뮬레이션을 사용하여 시뮬레이션()함수에서 팔 패키지:

sim(fruitMLM,100)  
sim(fuelMLM,100)

첫 번째 중 하나를 계산합니다.두 번째 하나는 다음 오류를 반환합니다:

Error in mvnorm(n.sims, bhat[j,], V.beta) :  
  'Sigma' is not positive definite

나는 실제로 이렇게는 8 개의 다른 응답합니다.여섯 그들이 잘 작동하고 그 두 가지 반환되는 이 오류가 있습니다.

누군가는 방법을 알고 이를 해결하려?

Answer 1

할 수 있는 경우에 당신이 볼 수있는 오류 메시지가 올바른지,즉,확인"에 의해 손으로는"공분산 행렬에서 실패하는 모형은 실제로 긍정적인 반 정?나는 확실하지 않다,하지만 내가 상상하는 몇 가지 검사를 수행할 수 있습니다-아마 여기에 뭔가가 있다: http://en.wikipedia.org/wiki/Positive-definite_matrix