esteem<-c(3, 4, 2, 5, 3, 4, NA, 2, 2, 3, 2, 1, 2, NA, 2, 2, 1, 2,
1, 0, 0, 0)
pesteem<-c("A", "A", "A", "A", "A", "A",
NA, "B", "B", "B", "B", "B",
"B", NA, "B1", "B1", "B1", "B1", "B1", "B1", "B1", "B1")
medabove(esteem, pesteem,"A","B1")
Run the code above in your browser using DataLab