Ok, I studied the source code of melt.data.frame. With na.rm=T melt operates normally except that it deletes rows from the molten data.frame that have NAs in the value column. NAs in the id.vars are not touched. This could be clearer in the documentation especially as it seems that earlier versions of reshape behaved differently. Best, Titus