eff=c(58.2,52.6,56.2,41.2,65.3,60.8,49.1,42.8,54.1,50.5,51.6,48.4,60.1,58.3,70.9,73.2,39.2,40.7,75.8,71.5,58.2,51.0,48.7,41.4) tuijin=rep(c(1,1,2,2,3,3),4) ranliao=rep(1:4,rep(6,4)) tuijin=as.factor(tuijin) ranliao=as.factor(ranliao) datas=data.frame…