#data(TexasData)
#nu <- 3
#zeta <- nu/(nu-1)
#heter <- list(sigma2y=VarMeanIng/zeta, sigma2xi=cbind(VarPopLabor16,
# VarP_Desemp)/zeta)
#model <- fmem(MeanIng/100 ~ PopLabor16 + P_Desemp | log(DensPobla) +
# P_PbFem + P_PNegra + bsp(P_PbHisp), data=TexasData, family="Slash",
# eta=nu, burn.in=10000, post.sam.s=10000, heter=heter,thin=10)
#summary(model)
Run the code above in your browser using DataLab