library(openxlsx) #import data df<-openxlsx::read.xlsx(file.choose()) #Adcock_etal_2020_Pinprick_Data.xlsx table(df$PrickResponseAfter) #8/26=31% of horn buds were fully de-sensitized #6/26=23% fully sensitized #12/26=46% partially de-sensitized