how to present the statistical results on the brain template

Viewing 30 reply threads
  • Author
    Posts
    • #2049
      Anand Joshi
      Moderator

        The visualization is described here:
        http://brainsuite.org/wp-content/uploads/2017/06/BrainSuite_Workshop_2017_BSSR_Shantanu_Joshi.pdf
        Please let me know how it goes.

      • #2050
        glf19780301
        Participant

          Thank you.We’re dealing with the surface area of the brain, how do we deal with that? I didn’t see any relevant content in the document you gave me.

        • #2051
          Anand Joshi
          Moderator

            Are you running ROI wise analysis? You can check slides 43-45.

          • #2052
            glf19780301
            Participant

              Thank you.

              “Volumetric analysis – sMRI” and “volumetric analysis – dMRI” how do I calculate?

            • #2053
              Anand Joshi
              Moderator

                There are slides in the link on diffusion analysis (FA). sMRI analysis (T1) can be done using tensor based morphometry (TBM). Please check slides on those.

              • #2054
                glf19780301
                Participant

                  I want to compare the surface and area of each area of the whole brain. How do I do that?I didn’t see it in “slides 43-45.”
                  Thank you.

                • #2055
                  shjoshi
                  Moderator

                    Thanks Anand for pointing to the right slides.
                    As Anand replied,
                    volumetric MRI analysis can be found starting at slide 33.
                    volumetric diffusion analysis (analysis of diffusion parameters such as FA, MD, RD) can be found from slide 38 onwards.

                    Have you installed bssr and are you able to load the processed data in R?

                    Shantanu

                  • #2056
                    glf19780301
                    Participant

                      I ran into a problem during the installation, so I followed the instructions to run the command
                      “install.packages(c(‘devtools’, ‘htmlwidgets’, ‘iterators’, ‘magrittr’, ‘gtable’, ‘R6’))
                      devtools::install_url(‘http://brainsuite.org/wp-content/uploads/2018/05/bssr_0.2.1.tar.gz’)”

                      show:

                      > install.packages(c(‘devtools’, ‘htmlwidgets’, ‘iterators’, ‘magrittr’, ‘gtable’, ‘R6’))
                      Warning in install.packages :
                      ‘lib = “C:/Program Files/R/R-3.6.2/library”‘ is not writable
                      Warning in install.packages :
                      cannot create dir ‘C:\Users\glf19\OneDrive\??’, reason ‘Invalid argument’
                      Error in install.packages : unable to create ‘C:/Users/glf19/OneDrive/??/R/win-library/3.6’
                      > devtools::install_url(‘http://brainsuite.org/wp-content/uploads/2018/05/bssr_0.2.1.tar.gz’)
                      Error in loadNamespace(name) : there is no package called ‘devtools’
                      >

                    • #2057
                      glf19780301
                      Participant

                        I have installed RStudio and Rtools, I don’t know what the problem is, how to correct it

                      • #2058
                        shjoshi
                        Moderator

                          I think this seems to be some permissions issue on your computer for installing packages.
                          This is not related to bssr, but concerns your R installation.

                          I would suggest getting devtools installed first.
                          Since you are using Windows, have you installed Rtools?
                          I would just run install.packages(‘devtools’) first and see if you have any error.

                          S

                        • #2059
                          shjoshi
                          Moderator

                            Can you post the output of
                            install.packages(‘devtools’)?

                          • #2060
                            glf19780301
                            Participant

                              I want to compare the surface and area of each area of the whole brain.How to compare and count “surface area”?

                            • #2061
                              shjoshi
                              Moderator

                                Each individual subject folder contains a file <subjid>.roiwise.stats.txt, where <subjid> should be replaced by your subject Id name.

                                This stats.txt file contains the surface areas for the cortical ROIs mapped by SVREG.

                                To get the index of the ROI and the label, you should look at
                                BrainSuiteXX/svreg/BCI-DNI_brain_atlas/brainsuite_labeldescription.xml
                                or
                                BrainSuiteXX/svreg/BrainSuiteAtlas1/brainsuite_labeldescription.xml
                                depending on the atlas used.

                                To perform a statistical ROI analysis in bssr, look at slide 42 onwards.

                              • #2062
                                glf19780301
                                Participant

                                  > install.packages(‘devtools’)
                                  Error: unexpected input in “install.packages(‘”

                                  Hi Shjoshi,
                                  Thank you very much for your guidance.When I input “> install.packages(‘devtools’) “,it will shows “Error: unexpected input in “install.packages(‘”

                                • #2063
                                  glf19780301
                                  Participant

                                    Could you tell me the specific installation steps of “R”? I can uninstall and start again

                                  • #2064
                                    shjoshi
                                    Moderator

                                      I suspect your R and RStudio installation must be fine. Before you try uninstall/reinstall, can you try running R as an Administrator?

                                      Look at this link,

                                      Then post the output of .libPaths().

                                      If you ran RStudio as an administrator, try installing devtools.
                                      install.packages("devtools")

                                      Note devtools is enclosed in double quotes.
                                      When typing the install.package command above, please type it manually in your R console, since copying pasting from the web may garble some characters.

                                    • #2065
                                      glf19780301
                                      Participant

                                        How do I run R in administrator mode?Can you take a look at my computer remotely?

                                      • #2066
                                        glf19780301
                                        Participant

                                          In my computer, my account is the administrator.

                                        • #2067
                                          shjoshi
                                          Moderator

                                            Ok. Now, when running R, can you find your R installation, right click on it and select “Open as Administrator”?

                                            Look at this too.
                                            https://community.rstudio.com/t/error-after-r-update-lib-c-program-files-r-r-3-5-0-library-is-not-writable/7947/4

                                          • #2068
                                            glf19780301
                                            Participant

                                              * DONE (bssr)
                                              > library(bssr)
                                              Finding BrainSuite atlas file paths…Done. Valid BrainSuite installation found at C:/Program Files/BrainSuite19a
                                              > get_brainsuite_install_path()
                                              [1] “C:/Program Files/BrainSuite19a”
                                              > get_brainsuite_install_path()
                                              [1] “C:/Program Files/BrainSuite19a”
                                              > ?bssr
                                              >

                                              Is that ok? Thank you.

                                            • #2069
                                              glf19780301
                                              Participant

                                                I have two groups of patients, divided into HP group (39 patients) and HC group (40 volunteers), which have all completed the process of svreg.I want to compare the volume, thickness and surface area of each brain region of the two groups.

                                                How to type the command?

                                                1.bss_data <- load_data(type=,…)
                                                2.bss_model <- bss_anova(main_effect=, covariates=, bss_data)
                                                bss_model <- bss_corr(corr_var=, bss_data)
                                                bss_model <- bss_ttest(group_var=, bss_data, paired=)

                                              • #2070
                                                glf19780301
                                                Participant

                                                  Load data

                                                  bss_data <- load_data(type=,…)

                                                  How do I edit the command in this step?

                                                • #2071
                                                  glf19780301
                                                  Participant

                                                    Can you recommend me a few articles published using Brainsuite?

                                                  • #2072
                                                    shjoshi
                                                    Moderator

                                                      Great that you were able to get this installed.
                                                      To run the cortical thickness analysis, look at slide 27. It shows the 3 steps needed to load the data, execute the model, and save the statistical outputs.

                                                      For the first step, you will need 2 things.
                                                      i. Location of the subject directory (subjdir variable in load_bss_data)
                                                      ii. A demographics csv file (csv variable in load_bss_data, which contains columns such as age, sex, and group. Your group will be HC or HP (for now set it as 1 or 0).

                                                      Step 1. Then call the command given in Step 1 by setting the appropriate variables. For e.g. hemi="left" for left hemisphere. If you haven’t smoothed your data, don’t pass the smooth variable.

                                                      Step 2. Call the command in Step 2 as bss_model <- bss_anova(main_effect = "group", covariates = "age", bss_data = bss_data)

                                                      Step 3. Call save_bss_out(bss_data, bss_model, outdir="/path/to/some/output directory")

                                                      To compare volumes, I’d suggest running the TBM analysis. See slide 34 for the steps. The functions are same as above. The options change.

                                                      To compare surface areas, I’d suggest running the ROI analysis. (Slide 43).

                                                    • #2073
                                                      shjoshi
                                                      Moderator

                                                        BTW, just FYI we are getting ready to release the next version of BrainSuite and the statistical tools. This version of bssr will also have functionality that will create html reports. This will be released in 1-2 weeks.

                                                      • #2074
                                                        glf19780301
                                                        Participant

                                                          > load_bss_data
                                                          function (type = “cbm”, subjdir = “”, csv = “”,
                                                          hemi = “left”, smooth = 0, roiids = 0, roimeas = “gmthickness”,
                                                          measure = “”, atlas = “”, eddy = TRUE)
                                                          {
                                                          valid_types <- c(“cbm”, “tbm”, “roi”, “dbm”,
                                                          “nca”)
                                                          if (!type %in% valid_types)
                                                          stop(sprintf(“Valid data types are %s.”, paste(valid_types,
                                                          collapse = “, “)), call. = FALSE)
                                                          switch(type, cbm = {
                                                          bss_data <- load_cbm_data(subjdir = subjdir, csv = csv,
                                                          hemi = hemi, smooth = smooth)
                                                          }, tbm = {
                                                          bss_data <- load_tbm_data(subjdir = subjdir, csv = csv,
                                                          smooth = smooth, atlas = atlas)
                                                          }, dbm = {
                                                          bss_data <- load_dbm_data(subjdir = subjdir, csv = csv,
                                                          measure = measure, smooth = smooth, atlas = atlas,
                                                          eddy = eddy)
                                                          }, roi = {
                                                          bss_data <- load_roi_data(subjdir, csv, roiids, roimeas)
                                                          })
                                                          return(bss_data)

                                                          I did this step, is it right?

                                                        • #2075
                                                          shjoshi
                                                          Moderator

                                                            You should call the function like this:
                                                            bss_data <- load_bss_data(type="cbm", subjdir = "/path/to/the/top level subject directory", csv = "/path/to/your/demographics.csv", hemi="left")

                                                          • #2076
                                                            glf19780301
                                                            Participant

                                                              I am very glad to get your help. I have a preliminary contact with R language. Please give me your detailed guidance.

                                                            • #2077
                                                              shjoshi
                                                              Moderator

                                                                No problems. It will be a good idea to go step by step. Perhaps only focus on one type of analysis first, get familiar with the command structure, and then move to other analysis.
                                                                Loading data first, then running a basic model, and saving the output.

                                                              • #2767
                                                                yemishade
                                                                Participant

                                                                  Hi everyone, I am somewhat new to using Brainsuite and I hope I would find a solution to my problems here. I have a group of MRI scans and I am interested in obtaining the tissue probability ratio (for Grey matter, white matter and CSF) of voxels in every anatomical region in each hemisphere labelled in AAL2. So I wonder if Brainsuite could be used in processing MRI such that I would be able to identify the Tissue probability ration in voxels in every Region.

                                                                  Also, how can I label the anatomical regions of the brain using Brainsuite?

                                                                • #2769
                                                                  Anand Joshi
                                                                  Moderator

                                                                    Sure, Please have a look at
                                                                    http://brainsuite.org/processing/surfaceextraction/pvc/

                                                                    I would go through the brainsuite tutorial which

                                                                    Also for labeling
                                                                    http://brainsuite.org/quickstart/svreg/

                                                                    If you go through getting stated sequence in brainsuite, that will give you good overview of the software as well as answer your questions in more detail. Please feel free to reach out if you still have any questions.
                                                                    For that click on Getting Started on http://brainsuite.org/ and the first 3 subsections would be very useful.

                                                                    • This reply was modified 4 years, 4 months ago by Anand Joshi.
                                                                Viewing 30 reply threads
                                                                • You must be logged in to reply to this topic.