Вы можете попробовать что-то вроде этого:
#define matrix with no. of columns and rows
m<-matrix(1:6,nrow = 6,ncol = 6)
m_req<-m
for(i in 1:nrow(m_req)){
if(i!=nrow(m_req)){
m_req[i,]<-m_req[i,]+m_req[i+1,]
m_req[,i]<-m_req[,i]+m_req[,i+1]
}
}
req_columns<-seq(1,ncol(m_req),by=2)
m_req<-m_req[req_columns,req_columns]