Description Usage Arguments Value Examples
roc
calculate the fpr and tpr for the roc curve
1 | roc(a, a0)
|
a |
p * p estimated graph |
a0 |
p * p true graph |
a list.
tpr |
tpr sequence |
fpr |
fpr sequence |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | p=30
n=120
tmp=runif(p-1,0.5,1)
s=c(0,cumsum(tmp));
s1=matrix(s,p,p)
cov.mat.true=exp(-abs(s1-t(s1)))
prec.mat.true=solve(cov.mat.true);
a=matrix(rnorm(p*n),n,p)
data.sa=a%*%chol(cov.mat.true);
true.graph = outer(1:p,1:p,f<-function(x,y){(abs(x-y)==1)})
methodlist = c("ridge","lasso","sam")
fit = vector(mode="list", length=3)
info = vector(mode="list", length=3)
auc = NULL
plot.new()
for(i in 1:3){
method = methodlist[i]
fit[[i]] = pgraph(data.sa, method = method)
info[[i]] = roc(fit[[i]]$statmat.pearson, true.graph)
auc[i] = sum(-diff(info[[i]]$fpr)*info[[i]]$tpr[-1])
lines(info[[i]]$fpr,info[[i]]$tpr, xlab='FPR',ylab='TPR', type='s', col=i)
cat(method, ': auc=', auc[i],'\n')
}
auc
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.