Похоже, вы отлично начали.Пара небольших синтаксических вещей;Я думаю, что это работает так, как вы хотите.
initial<- c(11.4,9.6,10.1,8.5,10.3,10.6,11.8,9.8,10.9,10.3,10.2,11.4,9.2,10.6,10.8,8.2,9.1,8.7,9.7,10.8,10.9,10.6,10.1,12.3,8.8,10.4,10.9,10.4,11.6,10.9)
final<-c(138.3,104.0,96.4,89.0,88.0,103.8,147.3,97.1,172.6,146.3,99.0,122.3,103.0,117.8,121.5,93.0,9.3,8.8,8.8,10.1,9.6,8.6,10.4,12.4,9.3,9.5,8.4,8.7,12.5,9.1)
group<-rep(c("treatment","control"), times=c(16,14))
group
selenium_data<-data.frame(initial,final,group)
selenium_data
boxplot(formula=initial~group, data=selenium_data)
