search

SachaEpskamp / qgraph

Developmental version of qgraph
GNU General Public License v2.0
68 stars 21 forks source link

centralityPlot with confidence intervals #40

Closed ariannechr closed 4 years ago

ariannechr commented 4 years ago

Hey Sacha,

Just trying to reproduce your Figure 2 in the Supplementary Materials (i.e. 95% quantile regions plotted in centralityPlot (Epskamp, Borsboom, Fried, 2018)

Bootstrap 95% CIs on centrality indices: boot3 <- bootnet(Data, 100, default = "pcor", type = "parametric")

Plot the 95% quantile regions: plot(boot3, statistics = c("strength","closeness","betweenness"),CIstyle = "quantiles")

Instead, I appear to be getting the node strength plot, like figure 6 (below) from the (Epskamp, Borsboom, Fried, 2018) tutorial paper

Arianne

diegogaldinof commented 4 years ago

Hi, Arianne!

I guess you need to specify the centrality measures in argument statistics.

Try this: boot3 <- bootnet(Data, 100, default = "pcor", type = "parametric", statistics = c("edge", "strength", "betweenness", "closeness"))

I hope to be contributing with your analysis.

Cheers!

ariannechr commented 4 years ago

thanks!

On Wed, Jun 17, 2020 at 3:49 PM Diego Galdino França < notifications@github.com> wrote:

Hi, Arianne!

I guess you need to specify the centrality measures in argument statistics .

Try this: boot3 <- bootnet(Data, 100, default = "pcor", type = "parametric", statistics = c("edge", "strength", "betweenness", "closeness"))

I hope to be contributing with your analysis.

Cheers!

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/SachaEpskamp/qgraph/issues/40#issuecomment-645665937, or unsubscribe https://github.com/notifications/unsubscribe-auth/AMOLD6CVJVTRG6FEBKP6ZN3RXFB5ZANCNFSM4KY3EFYQ .

alexp872 commented 10 months ago

Hello!

This fix doesn't work for me - I still get the node strength plot (as in Fig 6, referenced above), but I want the 95% quantile regions plotted in centralityPlot

boot1 <- bootnet(results, nCores = 8, nBoots = 1000, type = 'nonparametric', statistics= c("edge", "strength","closeness","betweenness"))

plot(boot1, statistics = c("strength"),CIstyle = "quantiles")

Any help much appreciated!

SachaEpskamp commented 10 months ago

This is deliberatly not made possible as bootstrapping CIs around centrality metrics is not possible.