One way you can create community level variables by cluster is to use the bysort command in Stata. For example we can create a binary variable for education level as follows: recode v106 (0/1=0 "noedu/primary") (2/3=1 "secondary+") (8/9=.), gen(edu)
We can then find the proportion with secondary+ education by cluster as follows: bysort v001: egen educlus=mean(edu)
This will create a new variable, educlus, which would be the proportion of respondents with a secondary+ education level in each cluster. You may also do this by region instead of by cluster. The same procedure can be followed with other variables.
For multilevel models in Stata, you have to install stmixed. Here are a presentation and a journal article, but there are many articles on this type of analysis in the literature:
Presentation: http://www.stata.com/meeting/uk13/abstracts/materials/uk13_c rowther.pdf
Article: http://bioline.org.br/request?ep13024
Another Stata command is mestreg but this is only available with Stata 14.