#upload data
grooming <- read.csv ("grooming.csv", header = T, stringsAsFactors = F)

#change headers
colnames(grooming)[colnames(grooming)=="Monkey_ID1"] <- "ID1"
colnames(grooming)[colnames(grooming)=="Monkey_ID2"] <- "ID2"
colnames(grooming)[colnames(grooming)=="Tot_grooming_bouts"] <- "Gr"

#grooming rates
grooming$grooming_rates <- grooming$Gr/grooming$ObsTime

#round grooming rates to two decimal positions
grooming$grooming_rates <-round(grooming$grooming_rates, digits =2)

#convert group from integer to character
grooming$Group <- as.character (grooming$Group)

#change group names
grooming["Group"] [grooming["Group"] == "0"] <- "mg"
grooming["Group"] [grooming["Group"] == "1"] <- "hg"
grooming["Group"] [grooming["Group"] == "2"] <- "rg"

#median grooming rates for each group
summarize_data <-ddply(grooming, .(Group), summarize,              #to summarize the data by mean and sd
                       median = median(grooming_rates),
                       sd = round (sd(grooming_rates), 2))

#only "rg" group
rg <- grooming [which(grooming$Group == "rg"),]
                       
#remove observation time less than 4 hours
rg_nolessfour<- rg [-which(rg$ObsTime < 4),]

#remove "fr"
#if you want to remove a specific ID, you can do it as follows:
rg_nolessfour_nofr <- rg_nolessfour [-c(which (rg_nolessfour$ID1 == "fr")),]