Hi,
It must be an easy question but how to boxplot a subset of data:
data = read.table("my_data.txt", header = T)
boxplot(data$var1[data$loc == "nice"]~data$loc_type[data$loc == "nice"])
#in this case, i want to display only the boxplot loc == "nice"
#doesn't display the boxplot of only loc == "nice". It also displays loc ==
"mice"
Thanks
Carol
loc type loc_type var1 var2
nice SA nice_SA 2.96E+03 3.33E+03
nice SA nice_SA 4.07E+03 4.44E+03
nice SA nice_SA 3.33E+03 4.81E+03
nice Am nice_Am 3.70E+02 0.00E+00
nice Am nice_Am 7.41E+02 0.00E+00
nice Am nice_Am 1.48E+03 0.00E+00
nice Ap nice_Ap 3.70E+02 0.00E+00
nice Ap nice_Ap 3.70E+02 0.00E+00
nice Ap nice_Ap 3.70E+02 0.00E+00
nice Am_Ap nice_Am_Ap 0.00E+00 0.00E+00
nice Am_Ap nice_Am_Ap 0.00E+00 0.00E+00
nice Am_Ap nice_Am_Ap 0.00E+00 0.00E+00
mice SA mice_SA 9.26E+05 6.67E+04
mice SA mice_SA 8.52E+05 7.78E+04
mice SA mice_SA 9.63E+05 8.15E+04
mice Am mice_Am 5.93E+03 5.56E+03
mice Am mice_Am 5.56E+03 7.04E+03
mice Am mice_Am 7.78E+03 7.41E+03
mice Ap mice_Ap 8.89E+03 4.07E+03
mice Ap mice_Ap 9.63E+03 3.33E+03
mice Ap mice_Ap 1.07E+04 2.59E+03
mice Am_Ap mice_Am_Ap 4.07E+03 1.11E+03
mice Am_Ap mice_Am_Ap 4.81E+03 1.11E+03
mice Am_Ap mice_Am_Ap 5.56E+03 3.70E+02
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
