我的数据如下图所示，数据名为hjr，在某一个wave下的ID和number是不变一致，想把不同wave下的ref提取出来放在新的数据集中范围为R325-R1075。并且将wave对应的值作为变量名保存起来。wave范围为【325-1075】。

      
转换成：

本人在R软件中编了些东西（附一），实现了两两的合并，但是由于wave为325-1075共750个，两两合并肯定太累了，
想求一个循环，将数据循环合并出来，请各位大侠们出出主意,谢谢啊！

附一：
> R400<-hjr[hjr$wave==400,];R401<-hjr[hjr$wave==401,];R402<-hjr[hjr$wave==402,];
> names(R400)[4]="R400";names(R401)[4]="R401";names(R402)[4]="R402";
> r1<-merge(R400,R401,by=c("ID"));
> p1<-subset(r1,select = c(ID,R400,R401))
> r2<-merge(p1,R402,by=c("ID"));
> p2<-subset(r2,select= c(ID,R400,R401,R402));
>.........