Prezados,
Eu tenho n Empresas diferentes em suas caracteristicas e cada empresa desta tem os mesmos cargos
Selecionei dados de diversos integrantes de cada cargo de cada Emresa, e foram anotados variaveis comportamentais que variam de 0 a 10, a qual para facilitar chamarei de x1,x2,x3, estas variaveis foram coletados mais de uma vez para cada pessoa.
O Objetivo é prever a probabilidade de ocorrencia de um fato, sendo que cada pessoa foi anotado tambem se este fato ocorreu ou não.
Utilizando
formMod1= fato~ x1+x2+x3+(1 | empresa/ cargo)
Mod1 <- glmer(
formMod1
, data =
dadosord
, family = binomial,
control = glmerControl(optimizer="bobyqa"))
Este modelo eu consigo fazer uma "predição" sem "problemas"
n11 <- data.frame(empresa=factor("M1", levels = levels(dadosord$empresa),ordered=FALSE),
cargo=factor("P1", levels = levels(dadosord$cargo),ordered=FALSE),
x1=1, x2=7,x3=7)
predict(
Mod1 , n11, type="response")
Mas eu fiquei preocupado pois tendo em vista que tenho mais de uma observação das caracteristicas por individuo, se eu nao teria que ou coloca-lo tambem como aleatorio
ou de alguma maneira analisar esta estrtura de possivel correlação
Pensando no primeiro caso fiz
formMod2= fato~ x1+x2+x3+(1 | empresa/ cargo) + (1 | ID)
Mod2 <- glmer( formMod2 , data = dadosord , family = binomial,
control = glmerControl(optimizer="bobyqa"))
Só que neste caso ao tentar prever a probabilidade do fato, a função me pede quem é o id proposto o que foge ao interesse pois devo estudar para uma pessoa aleatoria que tenha as mesmas
caracteristicas x1, x2, x3.
E a segunda forma, se for o caso, eu não sei como faria
Agradeço a orientação dos amigos