rich club range

[Syam Gadde]

Hi again,

In the manual, you have an example of how to calculate the rich club attrs from a graph:

    g[[1]][[N]] <- rich_club_attrs(g[[1]][[N]], c(rich_core(g[[1]][[N]])$k.r, 12))

where you specify a range of degrees where phi_norm > 1 is significant.  In your example above, the bottom of the range is k.r from the output of the rich_core function, but the top of the range is what seems like an arbitrary 12.  In my case the rich club size is actually greater than 12, so I am wondering if an appropriate upper bound could be calculated programmatically, something like:

  max(rich.dt[threshold==0.2 & Subject.ID=='mySubjectID' & p.fdr < 0.05,]$k)

Do you think that's an appropriate way to calculate this, or should I be doing something more sophisticated?

Thanks,

-syam

[Chris Watson]

Hi Syam,

You can, essentially, do whatever you want here, as long as you can justify it. That specific code example was not meant to state some kind of rule (as you have noticed, it was an arbitrary range) but rather to give a simple example for that function.

Your second example is almost identical to what I have done in the past, so I think it is reasonable.

Chris

--
You received this message because you are subscribed to the Google Groups "brainGraph-help" group.
To unsubscribe from this group and stop receiving emails from it, send an email to brainGraph-he...@googlegroups.com.
To post to this group, send email to brainGr...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/brainGraph-help/d6eb5164-d98c-4d72-af48-bba56b2024e1%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.