qfun <-
function(x,y,coef,q){
x=as.matrix(x)
p1=ncol(x)+1
r=y-coef[1]-x%*%coef[2:p1]
rhoq=sum(r*(q-as.numeric((r<0))))
s=sum(rhoq)
s
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.