Улучшить визуализацию большого сетевого графика в пакете 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 не позволяет вам удалять узлы, и поэтому мне было интересно, есть ли способ улучшить визуализацию графика. Я буду очень признателен за любое предложение! Еще раз, я благодарю вас за ваше время и внимание!

С наилучшими пожеланиями,

Эмма

0 ответов

Другие вопросы по тегам