plot_brain.kdiffnet: plot 3d brain network from kdiffnet result

In JointNets: End-to-End Sparse Gaussian Graphical Model Simulation, Estimation, Visualization, Evaluation and Application

Description

This function plots 3d brain network from kdiffnet result

Usage

## S3 method for class 'kdiffnet'
plot_brain(x, type = "task",
  neighbouroption = "task", subID = NULL, index = NULL,
  hastitle = TRUE, haslegend = TRUE, hasbackground = TRUE, ...)

Arguments

x	output generated from kdiffnet function (kdiffnet class)
type	type of graph. There are two options: "task" (difference graph) "neighbour" (zoom into nodes in the difference graph specified further by parameter "index" (node id)
neighbouroption	not used
subID	not used
index	determines which node(s) to zoom into when parameter "type" is "neighbour". This parameter could either be an integer or vector of integers representing node ids
hastitle	determines whether the graph title is displayed or not (TRUE to display / FALSE to hide)
haslegend	not used
hasbackground	determines whether the reference brain is plotted or not (TRUE to display / FALSE to hide)
...	extra parameters passed to igraph::rglplot()

Details

The function plots brain network using rglplot.igraph

Value

3d (rgl) brain network

Author(s)

Zhaoyang Wang (Author), Zhaoyang Wang (maintainer) [e-mail]

Examples

library(JointNets)
graphics.off()
par(ask=FALSE)
par(mfrow=c(1,1))
data(ABIDE_aal116_timeseries)
data(aal116coordinates)
layout = cbind(aal116coordinates$x.mni + 90,
aal116coordinates$y.mni+126, aal116coordinates$z.mni+72)
result = simulation(p=116, s = 0.001, ss = 0.001, n = c(1,1))$simulatedgraphs
class(result) = "simule"
plot_brain(result, type = "task", neighbouroption = "task",
subID = NULL, index = NULL, layout = layout, hasbackground = FALSE)

result = kdiffnet(ABIDE_aal116_timeseries[[1]], ABIDE_aal116_timeseries[[2]],
W = matrix(1,116,116), g = rep(0,116), epsilon = 0.1, lambda = 0.001)
plot_brain(result, type = "task", neighbouroption = "task",
subID = NULL, index = NULL, layout = layout)

JointNets documentation built on July 30, 2019, 1:02 a.m.