It is rather hard to try to reproduce this without knowing what and how you created geneExpression (what package did you use? - what is the class structure?) and if you are using R's default boxplot or a package specific function. It would also be helpful to have your sessionInfo( ) - Just on a guess, I would say there might be NA values that makes the ylim infinite - You might try explicitly setting a ylim value.