=====  VBMEG ꡼Ρ

С0.5-0.b.0
CVSversion_0_5_0_b_0

=====  CVS

CVSROOT
  :pserver:username@lobule2:/home/cnb/rhayashi/cvs_repository/vbmeg_genesis


===== С(0.5-0.a.0)νɲùܰ

ѹܰ

()ĥǡեޥåѹѥ᡼ѹؿ󥿡եѹ
      줫ɽ

1.brain
  1-1.Ǿǥåؿεǽɲ
2.basis
  2-1.basisեθߴб
3.bayes
  3-1.job_vb ǤΥѥ᥿åˡ()
  3-2.֤MEGǡΥץ꾮ɤΥåɲ
  3-3.Υǡå»
4.current
  4-1.currentեѿѹ()
  4-2.current_parm.overlap_modeɲ()
5.Х
  5-1.vb_load_resultΥХ
  5-2.proj_runΥХ
  5-3.plot functionΥХ
6.GUI
  6-1.job_plot_currentmapεǽɲ
  6-2.GUIΥХ
  6-3.GUIΥ֥
7.¾ѹ
  7-1.err_cell, err_arrayλѹ()

ɲùܰ

1.current
  1-1.ήɤ߹ߴؿload_currentɲ
2.GUI
  2-1.濴, Ⱦ¤MRIGUIɲ
  2-2.ֹ碌ץ켰(GUI)ɲ
  2-3.job_select_areaɲ
3.¾ɲ
  3-1.vb_versionɲ
  3-2.matlab_version˥Сӵǽɲ

===== ܺ

ѹܾܺ

1.brain
  1-1.Ǿǥåؿεǽɲ
    MRI饤Ťͽ񤭤륰դ䤷
    ե롧
      check_model/check_brain_model.m 
2.basis
  2-1.basisեθߴб
    ΥСver 0.1)Ǻäեɤ߹褦˽
    ե롧
      common/loadfunc/load_basis_info.m, load_basis.m
3.bayes
  3-1.job_vb ǤΥѥ᥿åˡ()
    ؿ󥿥եѹ
      vb_check_variable(B, Cov, Gact, Gall, ix)
                 
      vb_check_variable(B, Cov, Gact, Gall)
    ե롧
      job/job_vb.m,
      estimation/bayes/vb_set_vb_parm.m, vb_check_variable.m
  3-2.֤MEGǡΥץ꾮ɤΥåɲ
    ե롧
      estimation/bayes/get_megdata.m
  3-3.Υǡå»
    ե롧
       Modified Files:
         job/job_vb.m
       Added Files:
         estimation/bayes/vb_check_variable.m
4.current
  4-1.currentեѿѹ()
    ʬѿ
      ver 0.41
        Varact(Nact,Nwindow) : ˷׻
      ver 0.5
        Varact(Nact,Nsample) : ƻǷ׻

    ήλ֥ǥå¸
      Tsample = Time sample index for 'Jact/Jbck'
        Jact( :, t, :) 
       Time sample in MEG data :
         t_meg = Tsample(t) <-> bexp(:, t_meg,:)
       Time (ms)   in MEG data : 
         t_ms  = (t_meg - Pretriger)*(1000/SampleFreq)
  4-2.current_parm.overlap_modeɲ()
    ڥȥ׻Ѥ˽Ťʤ礦ʿѤˡƻǷ׻
      ήñ˷礷ƽϤ⡼ɤɲä
       = 0 : current is averaged over overlapped time window
       = 1 : current is not averaged for overlapped window
             current time series of each time windows 
             are concatenated sequentially for spectral analysis

    ؿ󥿥եѹ
       [weight ,Tindex, Nindex] = overlapped_timewindow_weight
                      
       [weight ,Tindex, Nindex, Nsample] = overlapped_timewindow_weight
    ե롧
        estimation/bayes
          vbmeg_current_reconstruct_tr.m vbmeg_current_reconstruct.m 
          overlapped_timewindow_weight.m
        template/
          set_current_parm_test.m set_current_parm.m
        job/
          job_current_tr.m job_current.m 
5.Х
  5-1.vb_load_resultΥХ
      Ver0.4Ȥθߴ˴ؤХν
    ե롧
      estimate/bayes/vb_load_result.m
  5-2.proj_runΥХ
      job_currentμ¹ԹԤˡߥ󤬤ʤäΤǡͤ
      ̤̤ɽԶ礬ä
    ե롧
      template/proj_run.m
  5-3.plot functionΥХ
     ʲδؿˤĤơԤindex˴ؤХäΤǽ
    ե롧
      plotfunc/plot_meg_2D.m plot_brain_inflate.m plot_brain.m 
6.GUI
  6-1.job_plot_currentmapεǽɲ
      پɽ򥹥饤ǥåȥܥåǲ̤򸫤ʤ
      ѹ뵡ǽɲä
    ե롧
       Modified Files:
         job/job_plot_currentmap/
            job_plot_cmap.fig job_plot_cmap_plotbrain.m 
            job_plot_cmap_print_spatialpattern.m 
            job_plot_cmap_print_temporalpattern.m 
            job_plot_cmap_update_spatialpattern.m job_plot_currentmap.m 
       Added Files:
         job/job_plot_currentmap/
            job_plot_cmap_change_Jratio.m 
  6-2.bayes_parm_editorΥХ
      Execܥ򲡤ȼ¹Υѥ᡼åؿǥ顼ȯ
    ե롧
      gui/bayes_parm_editor/bayes_parm_editor.m
  6-3.GUIΥ֥
      ֲˤ륳ݡͥȤΥեưνä
      Τǽ
    ե롧
       gui/bayes_parm_editor/
          bayes_parm_editor.fig, bayes_parm_editor_advanced.fig
7.¾ѹ
  7-1.err_cell, err_arrayλѹ()
    ʲΤ褦˻ѹ
        [err , Sx, Sy, indx ] = err_cell(x,y)

            if size of cell array x & y is different, 
              err = -1 is returned, and error massage is printed
              Sx, Sy = size of x & y
              indx = []
            if size of each element in the cell is different, 
              err = -2 is returned, and error massage is printed
              indx   = element index with different 
              Sx, Sy = size of elements x{indx} & y{indx}

        [err, Sx, Sy] = err_array(x,y)

            if size of x & y is different, 
             err = -1 is returned, and error massage is printed
             Sx, Sy = size of x & y

    ե롧
       common/utility/err_cell.m, err_array.m

ɲùܾܺ
1.current
  1-1.ήɤ߹ߴؿload_currentɲ
ե롧
        common/loadfunc/load_current.m
    ؿ͡
        % USAGE:
        %  [Jinfo,Jact,Jbck,Varact] = load_current(currfile)
        %
        % --- INPUT:
        %  currfile - current file name
        %
        % --- OUTPUT:
        % Jinfo : structure of current info
        %    .ix_act = Vertex index corresponding to active current
        %    .ix_bck = Vertex index corresponding to background current
        %    .Lact   = # of current direction at one vertex in active region
        %    .Lbck   = # of current direction at one vertex in background
        %    .Tsample    = Time sample index for 'Jact/Jbck'
        %    .SampleFreq = Sample Frequency [Hz]
        %    .Pretriger  = Pretriger period
        % 
        % Jact    : current in focal region
        %   Jact(Nact,Nsample)          for trial_average = ON 
        %   Jact(Nact,Nsample,Ntrials)  for trial_average = OFF
        %     Nact     = Lact * Nvact,  Nvact = length(ix_act)
        %     Nsample  : # of time sample, 
        %     Ntrials  : # of trials in all session]
        %
        %   Jact( n + Nvact*(i-1), t, :) 
        %    current at the vertex 'ix_act(n)', 'i-th' direction, time index 't'
        %
        %   Time sample in MEG data 
        %     t_meg = Tsample(t)
        %   Time (ms)   in MEG data : 
        %     t_ms  = (t_meg - Pretriger)*(1000/SampleFreq)
        %
        %
        % Jbck     : background current
        %   Jbck(Nbck,Nsample)          for trial_average = ON 
        %   Jbck(Nbck,Nsample,Ntrials)  for trial_average = OFF

        % Varact      : posterior variance of act_area
        %   Varact(Nact,Nsample)
        %     When active area is empty,
        %     posterior variance for background area is calculated
2.GUI
  2-1.濴, Ⱦ¤MRIGUIɲ
     MRIaxial, coronal, sagittalγ̤Ǻ㣳
     åƥܥ򲡲ȵ濴Ⱦ¤򥰥ե˳ǧ
     ͤȤƼ뤳ȤǤ롣
     ե롧
       mri_toolbox/sphere_head/*
  2-2.ֹ碌ץ켰(GUI)ɲ
     vivid, fastscan, MEG, MRIκɸϤΰֹ碌򤹤
     ץ켰ɲ
     ե롧
       gui/meg_location/*
  2-3.job_select_areaɲ
     ե롧
       job/job_select_area/*
3.¾ɲ
  3-1.vb_versionɲ
     ե롧
       common/utility/vb_version.m
  3-2.matlab_version˥Сӵǽɲ
    ؿ󥿡եɲá
      [result] = matlab_version(comparative_operator, ver_str)
      [result] = matlab_version(ver_str)
