if (FALSE) {
data(Orthodont)
attach(Orthodont)
sex = c()
sex[Sex=="Male"]=0
sex[Sex=="Female"]=1
y = distance #response
x = cbind(1,sex,age) #design matrix for fixed effects
z = cbind(1,age) #design matrix for random effects
#A median regression
median_reg = QRLMM(y,x,z,Subject,MaxIter = 500)
}
Run the code above in your browser using DataLab