Улучшить визуализацию большого сетевого графика в пакете qgraph
Я столкнулся с пакетом qgraph, и я так благодарен, что сделал, потому что у меня есть много красивых и информативных сюжетов. В настоящее время я работаю с более широкой базой данных из 755 метаболитов, измеренных у 15 видов бактерий в 2 временных точках, и у меня огромный график, который я не могу прикрепить здесь, поскольку у меня нет 10 пунктов репутации.
Я хотел бы видеть, как метаболиты коррелируют между собой, но я хотел бы видеть только существенные корреляции. Поэтому я использовал код, который Sacha Epskam любезно предоставил на своем сайте:
oralnet<-read.csv2('net24h.csv',dec='.')##file with the 755 metabolites
oralnet2<-oralnet[,2:756]##remove the first column with name of bacteria
namesnet<-read.csv2('netID.csv',dec='.') #file with the names of the bacteria
nameschars <- as.character(t(namesnet)[2:nrow(t(namesnet)),])
uniquechars <- unique(as.character(t(namesnet)[2:nrow(t(namesnet)),]))
nameslist<-list("Carbon" =which(nameschars==uniquechars[1]),
"Nitrogen"=which(nameschars==uniquechars[2]),
"PS"=which(nameschars==uniquechars[3]),
"Nutrient"=which(nameschars==uniquechars[4]),
"Peptide"=which(nameschars==uniquechars[5]))
rownames(oralnet2)<-c("Aa","Fn","Pg","Pi","Tf","Smut","Ssob","An","Csput","Sgord",
"Avisc","Ssal","Ssang","Vparv","Smitis")
qgraph(cor(oralnet2, use="pairwise"),minimum=0.25,groups=nameslist,
legend=TRUE,borders=FALSE,label.cex=0.4, labels=names(oralnet2),
layout="spring",line = 2.5,legend.cex=0.5,label.scale=FALSE,
graph="sig2",alpha=cc(0.0001,0.001))
Как я уже упоминал, существует много корреляций, и я хотел бы удалить некоторые узлы. Я знаю, что qgraph не позволяет вам удалять узлы, и поэтому мне было интересно, есть ли способ улучшить визуализацию графика. Я буду очень признателен за любое предложение! Еще раз, я благодарю вас за ваше время и внимание!
С наилучшими пожеланиями,
Эмма