using NeuroAnalyzer
eeg = load("files/eeg.hdf");[ Info: Loaded: EEG (24 × 308480 × 1; 1204.996 s)NeuroAnalyzer tutorials: EEG study
Load data:
Create the study object:
eeg1 = NeuroAnalyzer.filter(eeg, ch="all", fprototype=:butterworth, ftype=:bs, cutoff=(45, 55), order=4)
eeg2 = NeuroAnalyzer.filter(eeg, ch="all", fprototype=:butterworth, ftype=:bs, cutoff=(45, 55), order=8)
eeg3 = NeuroAnalyzer.filter(eeg, ch="all", fprototype=:butterworth, ftype=:bs, cutoff=(45, 55), order=12)
my_study = create_study([eeg1, eeg2, eeg3], [:g1, :g2, :g3])[ Info: Signal should be tapered prior to filtering to reduce edge artifacts[ Info: Signal should be tapered prior to filtering to reduce edge artifacts
[ Info: Signal should be tapered prior to filtering to reduce edge artifactsNeuroAnalyzer.STUDY(Dict{Symbol, Any}(), NeuroAnalyzer.NEURO[NeuroAnalyzer.NEURO(NeuroAnalyzer.HEADER(Dict{Symbol, Any}(:weight => -1, :id => "", :middle_name => "", :height => -1, :head_circumference => -1, :handedness => "", :last_name => "528004 SIT 52, 20220831-122227-{d589f756-53fc-4f1b-915d-6e3b8c1560ad}", :first_name => ""), Dict{Symbol, Any}(:epoch_id => "", :channel_type => ["eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg" … "eeg", "eeg", "ref", "ref", "eeg", "eeg", "eeg", "eog", "eog", "ecg"], :label => ["Fp1", "Fp2", "F3", "F4", "C3", "C4", "P3", "P4", "O1", "O2" … "T5", "T6", "A1", "A2", "Fz", "Cz", "Pz", "EOG1", "EOG2", "ECG"], :prefiltering => ["HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz" … "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz"], :gain => [0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218 … 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218], :data_type => "eeg", :recording_notes => "", :recording_date => "31.08.22", :sampling_rate => 256, :file_type => "EDF"…), Dict(:name => "", :design => "", :notes => "")), [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [1.114917356466603 75.29824297841112 … -4.432133031516375 1.2064258588154577; -0.8253654797087169 -63.13851838895876 … -2.2396139047108643 0.8194639519631813; … ; 0.7832653446794672 47.76730930462464 … -5.046593171529754 1.5755523935393645; 4.782830129089447 295.536318127934 … -12.297635221061574 3.3449404833300456;;;], Dict{Any, Any}(), 0×5 DataFrame
Row │ id start length description channel
│ String Float64 Float64 String Int64
─────┴────────────────────────────────────────────────, 23×9 DataFrame
Row │ label loc_radius loc_theta loc_x loc_y loc_z loc_radius_sp ⋯
│ String Float64 Float64 Float64 Float64 Float64 Float64 ⋯
─────┼──────────────────────────────────────────────────────────────────────────
1 │ Fp1 1.0 108.0 -0.31 0.95 -0.03 1.0 ⋯
2 │ Fp2 1.0 72.0 0.31 0.95 -0.03 1.0
3 │ F7 1.0 144.0 -0.81 0.59 -0.03 1.0
4 │ F3 0.65 129.0 -0.55 0.67 0.5 1.0
5 │ Fz 0.51 90.0 0.0 0.72 0.7 1.0 ⋯
6 │ F4 0.65 51.0 0.55 0.67 0.5 1.0
7 │ F8 1.0 36.0 0.81 0.59 -0.03 1.0
8 │ T3 1.0 180.0 -1.0 0.0 -0.03 1.0
9 │ C3 0.51 180.0 -0.72 0.0 0.7 1.0 ⋯
10 │ Cz 0.0 0.0 0.0 0.0 1.0 1.0
11 │ C4 0.51 0.0 0.72 0.0 0.7 1.0
⋮ │ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱
14 │ P3 0.65 231.0 -0.55 -0.67 0.5 1.0
15 │ Pz 0.51 270.0 0.0 -0.72 0.7 1.0 ⋯
16 │ P4 0.65 309.0 0.55 -0.67 0.5 1.0
17 │ T6 1.0 324.0 0.81 -0.59 -0.03 1.0
18 │ O1 1.0 252.0 -0.31 -0.95 -0.03 1.0
19 │ O2 1.0 288.0 0.31 -0.95 -0.03 1.0 ⋯
20 │ A1 1.0 192.0 -0.92 -0.23 -0.55 1.1
21 │ A2 1.0 -12.0 0.92 -0.23 -0.55 1.1
22 │ EOG1 1.01 149.0 -0.87 0.51 -0.37 1.0
23 │ EOG2 1.01 31.0 0.87 0.51 -0.37 1.0 ⋯
3 columns and 2 rows omitted, ["reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=iirnotch, ftype=nothing, cutoff=50, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=fir, ftype=hp, cutoff=0.1, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=butterworth, ftype=bs, cutoff=(45, 55), order=4, rp=-1, rs=-1, dir=twopass, w=nothing)"]), NeuroAnalyzer.NEURO(NeuroAnalyzer.HEADER(Dict{Symbol, Any}(:weight => -1, :id => "", :middle_name => "", :height => -1, :head_circumference => -1, :handedness => "", :last_name => "528004 SIT 52, 20220831-122227-{d589f756-53fc-4f1b-915d-6e3b8c1560ad}", :first_name => ""), Dict{Symbol, Any}(:epoch_id => "", :channel_type => ["eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg" … "eeg", "eeg", "ref", "ref", "eeg", "eeg", "eeg", "eog", "eog", "ecg"], :label => ["Fp1", "Fp2", "F3", "F4", "C3", "C4", "P3", "P4", "O1", "O2" … "T5", "T6", "A1", "A2", "Fz", "Cz", "Pz", "EOG1", "EOG2", "ECG"], :prefiltering => ["HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz" … "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz"], :gain => [0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218 … 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218], :data_type => "eeg", :recording_notes => "", :recording_date => "31.08.22", :sampling_rate => 256, :file_type => "EDF"…), Dict(:name => "", :design => "", :notes => "")), [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [0.25761757652603884 75.02557458287775 … -4.402358830484311 1.2507798103615004; -0.00648267199423088 -63.055303120179644 … -1.8981353613387768 0.7411694941290131; … ; 0.021472732615564965 47.001287084617495 … -4.737109149054424 1.6963137085885034; 1.3765932029514865 291.81355675341587 … -12.212119660314992 3.4958531301050333;;;], Dict{Any, Any}(), 0×5 DataFrame
Row │ id start length description channel
│ String Float64 Float64 String Int64
─────┴────────────────────────────────────────────────, 23×9 DataFrame
Row │ label loc_radius loc_theta loc_x loc_y loc_z loc_radius_sp ⋯
│ String Float64 Float64 Float64 Float64 Float64 Float64 ⋯
─────┼──────────────────────────────────────────────────────────────────────────
1 │ Fp1 1.0 108.0 -0.31 0.95 -0.03 1.0 ⋯
2 │ Fp2 1.0 72.0 0.31 0.95 -0.03 1.0
3 │ F7 1.0 144.0 -0.81 0.59 -0.03 1.0
4 │ F3 0.65 129.0 -0.55 0.67 0.5 1.0
5 │ Fz 0.51 90.0 0.0 0.72 0.7 1.0 ⋯
6 │ F4 0.65 51.0 0.55 0.67 0.5 1.0
7 │ F8 1.0 36.0 0.81 0.59 -0.03 1.0
8 │ T3 1.0 180.0 -1.0 0.0 -0.03 1.0
9 │ C3 0.51 180.0 -0.72 0.0 0.7 1.0 ⋯
10 │ Cz 0.0 0.0 0.0 0.0 1.0 1.0
11 │ C4 0.51 0.0 0.72 0.0 0.7 1.0
⋮ │ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱
14 │ P3 0.65 231.0 -0.55 -0.67 0.5 1.0
15 │ Pz 0.51 270.0 0.0 -0.72 0.7 1.0 ⋯
16 │ P4 0.65 309.0 0.55 -0.67 0.5 1.0
17 │ T6 1.0 324.0 0.81 -0.59 -0.03 1.0
18 │ O1 1.0 252.0 -0.31 -0.95 -0.03 1.0
19 │ O2 1.0 288.0 0.31 -0.95 -0.03 1.0 ⋯
20 │ A1 1.0 192.0 -0.92 -0.23 -0.55 1.1
21 │ A2 1.0 -12.0 0.92 -0.23 -0.55 1.1
22 │ EOG1 1.01 149.0 -0.87 0.51 -0.37 1.0
23 │ EOG2 1.01 31.0 0.87 0.51 -0.37 1.0 ⋯
3 columns and 2 rows omitted, ["reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=iirnotch, ftype=nothing, cutoff=50, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=fir, ftype=hp, cutoff=0.1, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=butterworth, ftype=bs, cutoff=(45, 55), order=8, rp=-1, rs=-1, dir=twopass, w=nothing)"]), NeuroAnalyzer.NEURO(NeuroAnalyzer.HEADER(Dict{Symbol, Any}(:weight => -1, :id => "", :middle_name => "", :height => -1, :head_circumference => -1, :handedness => "", :last_name => "528004 SIT 52, 20220831-122227-{d589f756-53fc-4f1b-915d-6e3b8c1560ad}", :first_name => ""), Dict{Symbol, Any}(:epoch_id => "", :channel_type => ["eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg", "eeg" … "eeg", "eeg", "ref", "ref", "eeg", "eeg", "eeg", "eog", "eog", "ecg"], :label => ["Fp1", "Fp2", "F3", "F4", "C3", "C4", "P3", "P4", "O1", "O2" … "T5", "T6", "A1", "A2", "Fz", "Cz", "Pz", "EOG1", "EOG2", "ECG"], :prefiltering => ["HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz" … "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz", "HP:0,18Hz LP:104,0Hz"], :gain => [0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218 … 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218, 0.17935713740749218], :data_type => "eeg", :recording_notes => "", :recording_date => "31.08.22", :sampling_rate => 256, :file_type => "EDF"…), Dict(:name => "", :design => "", :notes => "")), [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [0.0, 0.004, 0.008, 0.012, 0.016, 0.02, 0.023, 0.027, 0.031, 0.035 … 1204.961, 1204.965, 1204.969, 1204.973, 1204.977, 1204.98, 1204.984, 1204.988, 1204.992, 1204.996], [0.06094496337149765 74.93292884405291 … -5.225093520647854 -0.27751334594983834; 0.09735467535359228 -63.030616113869144 … -3.1502355625758622 -0.4391759055258552; … ; -0.07207904585735293 46.85408423958956 … -6.2359703368102375 -0.25575037573053055; -0.17364433008544577 290.09956303513735 … -14.289241169636936 -0.5187725999065581;;;], Dict{Any, Any}(), 0×5 DataFrame
Row │ id start length description channel
│ String Float64 Float64 String Int64
─────┴────────────────────────────────────────────────, 23×9 DataFrame
Row │ label loc_radius loc_theta loc_x loc_y loc_z loc_radius_sp ⋯
│ String Float64 Float64 Float64 Float64 Float64 Float64 ⋯
─────┼──────────────────────────────────────────────────────────────────────────
1 │ Fp1 1.0 108.0 -0.31 0.95 -0.03 1.0 ⋯
2 │ Fp2 1.0 72.0 0.31 0.95 -0.03 1.0
3 │ F7 1.0 144.0 -0.81 0.59 -0.03 1.0
4 │ F3 0.65 129.0 -0.55 0.67 0.5 1.0
5 │ Fz 0.51 90.0 0.0 0.72 0.7 1.0 ⋯
6 │ F4 0.65 51.0 0.55 0.67 0.5 1.0
7 │ F8 1.0 36.0 0.81 0.59 -0.03 1.0
8 │ T3 1.0 180.0 -1.0 0.0 -0.03 1.0
9 │ C3 0.51 180.0 -0.72 0.0 0.7 1.0 ⋯
10 │ Cz 0.0 0.0 0.0 0.0 1.0 1.0
11 │ C4 0.51 0.0 0.72 0.0 0.7 1.0
⋮ │ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱
14 │ P3 0.65 231.0 -0.55 -0.67 0.5 1.0
15 │ Pz 0.51 270.0 0.0 -0.72 0.7 1.0 ⋯
16 │ P4 0.65 309.0 0.55 -0.67 0.5 1.0
17 │ T6 1.0 324.0 0.81 -0.59 -0.03 1.0
18 │ O1 1.0 252.0 -0.31 -0.95 -0.03 1.0
19 │ O2 1.0 288.0 0.31 -0.95 -0.03 1.0 ⋯
20 │ A1 1.0 192.0 -0.92 -0.23 -0.55 1.1
21 │ A2 1.0 -12.0 0.92 -0.23 -0.55 1.1
22 │ EOG1 1.01 149.0 -0.87 0.51 -0.37 1.0
23 │ EOG2 1.01 31.0 0.87 0.51 -0.37 1.0 ⋯
3 columns and 2 rows omitted, ["reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=iirnotch, ftype=nothing, cutoff=50, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=fir, ftype=hp, cutoff=0.1, order=8, rp=-1, rs=-1, dir=twopass, w=nothing)", "reset_components(OBJ)", "filter(OBJ, ch=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], fprototype=butterworth, ftype=bs, cutoff=(45, 55), order=12, rp=-1, rs=-1, dir=twopass, w=nothing)"])], [:g1, :g2, :g3])
(!) Objects in the study objects must have the same size and sampling rate.
Save the study object:
save_study(my_study, file_name="my_study.hdf")Load the study object:
my_study = load_study(file_name="my_study.hdf")(!) Study is saved as HDF5 file. File name extension must be .hdf.