rm(list=ls())
Sys.setenv(LANG = "en")
library("devtools")
source("../R/hatmatrix.R")
#install_version("netmeta",version="0.9-6")
# indata = read.csv("griselda_wide.csv",header=TRUE,sep=",")
# C = getContributionMatrix(indata,type="netwide_binary",model="random",sm="OR")
#
# indata = read.csv2("diabetes_new.csv",header=TRUE,sep=";")
# C = getContributionMatrix(indata,type="netwide_binary",model="random",sm="OR")
indata = read.csv("Macfayeden.csv",header=TRUE,sep=";")
C = getHatMatrix(indata,type="netwide_binary",model="random",sm="OR")
# Direct effects
md = -C$hatMatrix$Pairwise[,1]
newnmd = lapply(names(md), function(n){gsub(" vs ",":",n)})
names(md) = newnmd
md = md[order(names(md))]
#network effects
mn = C$hatMatrix$NMA[,1]
#hatmatrix
HM = C$hatMatrix$H2bu
nmn = HM %*% md
# test that hmatrix works
diff = all(lapply((mn - nmn), function(n){n == 0}))
print("hatmatrix gives correct network measures")
print(diff)
#Check if Contribution matrix * direct effectes is 100
CM = C$contributionMatrix
ch2 = (CM * sign(HM)) %*% md
test2 = all(lapply((ch2), function(n){n == 100}))
print("Check if Contribution matrix * direct effectes is 100")
print(test2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.