hoxa9 <- "U82759_at" absts <- pm.getabst(hoxa9, "hu6800") gives Warning message: the condition has length > 1 and only the first element will be used in: if is.na(pm)) rval[[i]] <- NA else { Happens when there are multiple Pubmed IDs corresponding to a gene. Thanks, Ritu