data(manu_abs)
### Compositional tables approach
### analysis of the relative structure
result <- tabCoordWrapper(manu_abs, obs.ID='country',row.factor = 'output',
col.factor = 'isic', value='value', test = TRUE)
result$Bootstrap
### Classical approach
### generalized linear mixed effect model
if (FALSE) {
library(lme4)
m <- glmer(value~output*as.factor(isic)+(1|country),
data=manu_abs,family=poisson)
summary(m)
}
Run the code above in your browser using DataLab