y <- c(22.5 , 23.33 , 20.83 , 24.25 )
w <- c( 3,3,3,2)
# define matrix input
yM <- matrix( 0 , nrow=2 , ncol=4 )
wM <- yM
yM[1,] <- yM[2,] <- y
wM[1,] <- w
wM[2,] <- c(1,3,4, 3 )
# fit rowwise monotone regression
monoreg.rowwise( yM , wM )
# compare results with monoreg function from fdrtool package
library(fdrtool)
fdrtool::monoreg(x=yM[1,] , w=wM[1,])$yf
fdrtool::monoreg(x=yM[2,] , w=wM[2,])$yf
Run the code above in your browser using DataLab