Er zijn verschillende oplossingen. Een is om 'na.rm=TRUE' te specificeren:
### Gemakshalve even de lijst variabelenamen opslaan
varNames <- paste0("ASQ_", c(1, 2, 3, 19, 31, 33, 37, 38);
rowMeans(dat[, varNames], na.rm=TRUE);
Een andere manier is de functie 'validMeans' te gebruiken uit het 'userfriendlyscience' package, dat hetzelfde doet als 'mean.5' uit SPSS:
userfriendlyscience::validMeans(dat[, varNames], requiredValidValues=5);
Hier kan "requiredValidValues" ook een proportie zijn, trouwens.