I have been using MLWin for a while and am used to being able to put out predicted values and their associated standard errors.  I have recently started using r but am struggling to replicate this data.  Any suggestions?

To put my problem into context, I have fitted a three level linear model with a constant and one other independent variable.  I am only modelling the random effects for the constant.  I am interested in the level 3 and level 2 predictions so that I can graph my outcomes for the units at level 3 and level 2.  I know that in R I can get the level 2 predictions for the model using the fixed() command but I don't seem to be able get the associated se, nor can I get the same figure for level 3.  Reading threads on this matter in the R community seems to imply that I am being niaive in my wish but I am confused as to why I can do this in one software package but not in the other.