Nothing

```
# <rank.btmm>
# Copyright (C) <2014> <Hsiuying Wang, Yu-Chun Lin>
#
#
# This program is free software; you can redistribute it and/or modify
#
# it under the terms of the GNU General Public License as published by
#
# the Free Software Foundation; either version 2 of the License, or
#
# (at your option) any later version.
#
#
#
# This program is distributed in the hope that it will be useful,
#
# but WITHOUT ANY WARRANTY; without even the implied warranty of
#
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
#
# GNU General Public License for more details.
rank.btmm=function(data){
data=as.matrix(data)
n=dim(data)[1]
c=dim(data)[2]
k=t(apply(data,1, function(x) x==1))
B=matrix(0,c,c)
wm=B;temp=B
for(i in 1:n){
B[k[i,],!k[i,]]=1
wm=wm+B
B=temp
}
nmo=c-1
pi=rep(1,nmo)
change=1
gm=t(wm[,1:nmo])+wm[1:nmo,]
wins=apply(wm[1:nmo,],1,sum)
gind=(gm>0)
z=matrix(0,nrow=nmo,ncol=c)
pisum=z
while(norm(as.matrix(change),"F")>1e-8){
pius=matrix(rep(pi,c),ncol=c)
piust=t(pius[,1:nmo])
piust=cbind(piust,rep(1,nmo))
pisum[gind]=pius[gind]+piust[gind]
z[gind]=gm[gind]/pisum[gind]
newpi=wins/apply(z,1,sum)
change=newpi-pi
pi=newpi
}
result=round(c(pi,1),4)
gamma=result/sum(result)
names(gamma)=c(1:c)
gamma_temp=as.numeric(names(sort(gamma,decreasing=T)))
rank=numeric(c)
for(i in 1:c)
{
rank[gamma_temp[i]]=i
}
result_rank=rbind(gamma,rank)
return(result_rank)
}
```

**Any scripts or data that you put into this service are public.**

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.