查看: 1748|回复: 1
打印 上一主题 下一主题

[EEG/ERP] erplab插件学习问题汇总 [复制链接]

Rank: 2

水晶
47
心级
38
精华
0
主题
2
帖子
24
跳转到指定楼层
楼主
发表于 2019-4-30 10:10:13 |只看该作者 |倒序浏览
erplab插件非常好用,各位可以去试试,这个是它官方得分析脚本,可以根据自己的需要更改,很方便,不懂得可以在线问我。


附件: 你需要登录才可以下载或查看附件。没有帐号?注册
转发到微博

Rank: 3Rank: 3

水晶
4223
心级
369
精华
0
主题
10
帖子
277
沙发
发表于 2019-5-4 08:49:58 |只看该作者
? ? % CLEAR MEMORY AND THE COMMAND WINDOW
? ? CLEAR;
? ? CLC;


? ? % INITIALIZE THE ALLERP STRUCTURE AND CURRENTERP
? ? ALLERP? ???= BUILDERPSTRUCT([]);
? ? CURRENTERP = 0;


? ? % THIS DEFINES THE SET OF SUBJECTS
? ? SUBJECT_LIST = {'S1', 'S2', 'S3', 'S4', 'S5', 'S6'};
? ? NSUBJ? ?? ???= LENGTH(SUBJECT_LIST); % NUMBER OF SUBJECTS



? ? % PATH TO THE PARENT FOLDER, WHICH CONTAINS THE DATA FOLDERS FOR ALL SUBJECTS
? ? HOME_PATH??= '/USERS/LUCK/DOCUMENTS/SOFTWARE_DEVELOPMENT/ERPLAB_TOOLBOX/TEST_DATA/';


? ? % SET THE SAVE_EVERYTHING VARIABLE TO 1 TO SAVE ALL OF THE INTERMEDIATE FILES TO THE HARD DRIVE
? ? % SET TO 0 TO SAVE ONLY THE INITIAL AND FINAL DATASET AND ERPSET FOR EACH SUBJECT
? ? SAVE_EVERYTHING??= 1;


? ? % SET THE PLOT_PDFS VARIABLE TO 1 TO CREATE PDF FILES WITH THE WAVEFORMS
? ? % FOR EACH SUBJECT (SET TO 0 IF YOU DON'T WANT TO CREATE THE PDF FILES).
? ? PLOT_PDFS = 1;


? ? % LOOP THROUGH ALL SUBJECTS
? ? FOR S=1:NSUBJ

? ?? ???FPRINTF('\N******\NPROCESSING SUBJECT %S\N******\N\N', SUBJECT_LIST{S});

? ?? ???% PATH TO THE FOLDER CONTAINING THE CURRENT SUBJECT'S DATA
? ?? ???DATA_PATH??= [HOME_PATH SUBJECT_LIST{S} '/'];

? ?? ???% CHECK TO MAKE SURE THE DATASET FILE EXISTS
? ?? ???% INITIAL FILENAME = PATH PLUS SUBJECT# PLUS _EEG.SET
? ?? ???SNAME? ?? ?= [DATA_PATH SUBJECT_LIST{S} '_EEG.SET'];

? ?? ???IF EXIST(SNAME, 'FILE')<=0
? ?? ?? ?? ?? ? FPRINTF('\N *** WARNING: %S DOES NOT EXIST *** \N', SNAME);
? ?? ?? ?? ?? ? FPRINTF('\N *** SKIP ALL PROCESSING FOR THIS SUBJECT *** \N\N');
? ?? ???ELSE

? ?? ?? ?? ?%% LOAD DATA
? ?? ?? ?? ?% LOAD ORIGINAL DATASET
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: LOADING DATASET ****\N\N\N', SUBJECT_LIST{S});
? ?? ?? ?? ?EEG = POP_LOADSET('FILENAME', [SUBJECT_LIST{S} '_EEG.SET'], 'FILEPATH', DATA_PATH);


? ?? ?? ?? ?%% CHANNEL LOCATIONS
? ?? ?? ?? ?% ADD THE CHANNEL LOCATIONS
? ?? ?? ?? ?% WE'RE ASSUMING THE FILE 'STANDARD-10-5-CAP385.ELP' IS SOMEWHERE
? ?? ?? ?? ?% IN THE PATH.??THIS CAN BE COPIED FROM
? ?? ?? ?? ?% PLUGINS/DIPFIT2.2/STANDARD_BESA/ INSIDE THE EEGLAB
? ?? ?? ?? ?% FOLDER.
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: ADDING CHANNEL LOCATION INFO ****\N\N\N', SUBJECT_LIST{S});
? ?? ?? ?? ?EEG = POP_CHANEDIT(EEG, 'LOOKUP','STANDARD-10-5-CAP385.ELP');

? ?? ?? ?? ?% SAVE DATASET WITH _CHAN SUFFIX INSTEAD OF _EEG
? ?? ?? ?? ?EEG.SETNAME = [SUBJECT_LIST{S} '_CHAN']; % NAME FOR THE DATASET MENU

? ?? ?? ?? ?IF (SAVE_EVERYTHING)

? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);

? ?? ?? ?? ?END


? ?? ?? ?? ?%% EVENTLIST
? ?? ?? ?? ?% CREATE EVENTLIST AND SAVE (POP_EDITEVENTLIST ADDS _ELIST SUFFIX)
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: CREATING EVENTLIST ****\N\N\N', SUBJECT_LIST{S});
? ?? ?? ?? ?EEG = POP_CREABASICEVENTLIST( EEG , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'ALPHANUMERICCLEANING', 'ON'? ?? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'BOUNDARYNUMERIC'? ???, { -99 }? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'BOUNDARYSTRING'? ?? ?, { 'BOUNDARY' } );
? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_ELIST']; % NAME FOR THE DATASET MENU

? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);
? ?? ?? ?? ?END


? ?? ?? ?? ?%% HIGH-PASS FILTER
? ?? ?? ?? ?% HIGH-PASS FILTER THE EEG
? ?? ?? ?? ?% CHANNELS = 1 TO 16; HIGH-PASS CUTOFF AT 0.1 HZ;
? ?? ?? ?? ?% NO LOWPASS FILTER; ORDER OF THE FILTER = 2.
? ?? ?? ?? ?% TYPE OF FILTER = "BUTTERWORTH"; REMOVE DC OFFSET; FILTER
? ?? ?? ?? ?% "BETWEEN" BOUNDARY EVENTS
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: HIGH-PASS FILTERING EEG AT 0.1 HZ ****\N\N\N', SUBJECT_LIST{S});? ?? ?? ?? ???

? ?? ?? ?? ?EEG??= POP_BASICFILTER( EEG,??1:16 , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'BOUNDARY'??, 'BOUNDARY', ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'CUTOFF'? ? , 0.1? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'DESIGN'? ? , 'BUTTER'??, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'FILTER'? ? , 'HIGHPASS', ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'ORDER'? ???,??2? ?? ???);

? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_HPFILT'];
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);? ?? ?? ?? ???
? ?? ?? ?? ?END



? ?? ?? ?? ?%% RE-REFERENCE
? ?? ?? ?? ?% USE CHANNEL OPERATIONS TO INSERT A BIPOLAR EOG CHANNEL
? ?? ?? ?? ?% AND RE-REFERENCE TO THE AVERAGE OF THE EARLOBES
? ?? ?? ?? ?% EQUATIONS ARE STORED IN 'CHANOPS_REREF_BIVEOG.TXT', WHICH
? ?? ?? ?? ?% MUST BE IN THE HOME DIRECTORY FOR THE EXPERIMENT.
? ?? ?? ?? ?% SAVE OUTPUT WITH _REF SUFFIX.
? ?? ?? ?? ?%

? ?? ?? ?? ?EEG? ?? ?? ?= POP_EEGCHANOPERATOR(EEG, [HOME_PATH 'CHANOPS_REREF_BIVEOG.TXT']);
? ?? ?? ?? ?
? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_REF'];
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);
? ?? ?? ?? ?END



? ?? ?? ?? ?%% BIN
? ?? ?? ?? ?% USE BINLISTER TO SORT THE BINS AND SAVE WITH _BINS SUFFIX
? ?? ?? ?? ?% WE ARE ASSUMING THAT 'BINLISTER_DEMO_1.TXT' IS PRESENT IN THE
? ?? ?? ?? ?% HOME FOLDER.
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: RUNNING BINLISTER ****\N\N\N', SUBJECT_LIST{S});? ?? ?

? ?? ?? ?? ?EEG? ?? ?? ?= POP_BINLISTER( EEG , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'BDF'? ? , [HOME_PATH 'BINLISTER_DEMO_1.TXT'], ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'INDEXEL', 1? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'SENDEL2', 'EEG' , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'VOUTPUT', 'EEG' );

? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_BINS'];
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);
? ?? ?? ?? ?END



? ?? ?? ?? ?%% EPOCH
? ?? ?? ?? ?% EXTRACTS BIN-BASED EPOCHS (200 MS PRE-STIM, 800 MS POST-STIM.
? ?? ?? ?? ?% BASELINE CORRECTION BY PRE-STIM WINDOW)
? ?? ?? ?? ?% THEN SAVE WITH _BE SUFFIX
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: BIN-BASED EPOCHING ****\N\N\N', SUBJECT_LIST{S});

? ?? ?? ?? ?EEG? ?? ?? ?= POP_EPOCHBIN( EEG , [-200.0??800.0],??'PRE');

? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_EPOCHS'];
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? EEG = POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);
? ?? ?? ?? ?END



? ?? ?? ?? ?%% 2 ROUNDS OF ARTIFACT DETECTION
? ?? ?? ?? ?% THEN EXPORT EVENTLIST JUST FOR REVIEW
? ?? ?? ?? ?% SAVE THE PROCESSED EEG TO DISK BECAUSE THE NEXT STEP WILL BE AVERAGING
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: ARTIFACT DETECTION (MOVING WINDOW PEAK-TO-PEAK AND STEP FUNCTION) ****\N\N\N', SUBJECT_LIST{S});? ?? ?? ?? ???

? ?? ?? ?? ?% ARTIFACT DETECTION - RD. 1
? ?? ?? ?? ?%??MOVING WINDOW. TEST WINDOW = [-200
? ?? ?? ?? ?% 798]; THRESHOLD = 100 UV; WINDOW WIDTH = 200 MS;
? ?? ?? ?? ?% WINDOW STEP = 50 MS; CHANNELS = 1 TO 17; FLAGS TO BE ACTIVATED = 1 & 4
? ?? ?? ?? ?%
? ?? ?? ?? ?EEG = POP_ARTMWPPTH( EEG , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'CHANNEL'? ???,??1:17? ?? ?, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'FLAG'? ?? ???, [ 1 4]? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'THRESHOLD'? ?,??100? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'TWINDOW'? ???, [ -200 798], ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'WINDOWSIZE'??,??200? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ???'WINDOWSTEP'??,??50? ?? ???);

? ?? ?? ?? ?%% ARTIFACT DETECTION - RD. 2
? ?? ?? ?? ?% STEP-LIKE ARTIFACTS IN THE BIPOLAR
? ?? ?? ?? ?% VEOG CHANNEL (CHANNEL 14, CREATED EARLIER WITH CHANNEL OPERATIONS)
? ?? ?? ?? ?% THRESHOLD = 30 UV; WINDOW WIDTH = 400 MS;
? ?? ?? ?? ?% WINDOW STEP = 10 MS; FLAGS TO BE ACTIVATED = 1 & 3
? ?? ?? ?? ?EEG = POP_ARTSTEP( EEG , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'CHANNEL'? ?,??14? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'FLAG'? ?? ?, [ 1 3]? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'THRESHOLD' ,??30? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'TWINDOW'? ?, [ -200 798], ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'WINDOWSIZE',??400? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'WINDOWSTEP',??10? ?? ???);

? ?? ?? ?? ?EEG.SETNAME = [EEG.SETNAME '_AR'];
? ?? ?? ?? ?EEG? ?? ?? ?= POP_SAVESET(EEG, 'FILENAME', [EEG.SETNAME '.SET'], 'FILEPATH', DATA_PATH);
? ?? ?? ?? ?EEG? ?? ?? ?= POP_EXPORTEEGEVENTLIST(EEG, 'FILENAME', [DATA_PATH SUBJECT_LIST{S} '_EVENTLIST_AR.TXT']);



? ?? ?? ?? ?% REPORT PERCENTAGE OF REJECTED TRIALS (COLLAPSED ACROSS ALL BINS)
? ?? ?? ?? ?ARTIFACT_PROPORTION = GETARDETECTION(EEG);
? ?? ?? ?? ?FPRINTF('%S: PERCENTAGE OF REJECTED TRIALS WAS %1.2F\N', SUBJECT_LIST{S}, ARTIFACT_PROPORTION);

? ?? ?

? ?? ?? ?? ?%% AVERAGE
? ?? ?? ?? ?% ONLY GOOD TRIALS.??INCLUDE STANDARD DEVIATION.??SAVE TO DISK.
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: AVERAGING ****\N\N\N', SUBJECT_LIST{S});? ?? ?? ?? ???

? ?? ?? ?? ?ERP? ?? ?? ?= POP_AVERAGER( EEG, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'CRITERION'? ?? ?, 'GOOD', ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'DSINDEX'? ?? ???, 6, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'EXCLUDEBOUNDARY', 'ON', ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'SEM'? ?? ?? ?? ?, 'ON');
? ?? ?? ?? ?ERP.ERPNAME = [SUBJECT_LIST{S} '_ERPS'];??% NAME FOR ERPSET MENU
? ?? ?? ?? ?POP_SAVEMYERP(ERP, 'ERPNAME', ERP.ERPNAME, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', DATA_PATH, 'WARNING', 'OFF');



? ?? ?? ?? ?%% FILTER ERP
? ?? ?? ?? ?% CHANNELS = 1 TO 17; NO HIGH-PASS;
? ?? ?? ?? ?% LOWPASS CUTOFF AT 30 HZ; ORDER OF THE FILTER = 2.
? ?? ?? ?? ?% TYPE OF FILTER = "BUTTERWORTH"; DO NOT REMOVE DC OFFSET
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: LOW-PASS FILTERING ERP AT 30 HZ ****\N\N\N', SUBJECT_LIST{S});? ?? ?? ?? ???

? ?? ?? ?? ?ERP = POP_FILTERP( ERP,1:17 , 'CUTOFF',30, 'DESIGN', 'BUTTER', 'FILTER', 'LOWPASS', 'ORDER',2 );

? ?? ?? ?? ?ERP.ERPNAME = [ERP.ERPNAME '_30HZ'];??% NAME FOR ERPSET MENU
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? POP_SAVEMYERP(ERP, 'ERPNAME', ERP.ERPNAME, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', DATA_PATH, 'WARNING', 'OFF');
? ?? ?? ?? ?END



? ?? ?? ?? ?%% BIN OPERATIONS
? ?? ?? ?? ?% CREATE A DIFFERENCE WAVE AND SAVE WITH _DIFF SUFFIX
? ?? ?? ?? ?% DO THIS ON THE UNFILTERED DATA, SO FIRST RELOAD UNFILTERED FILE
? ?? ?? ?? ?% THEN DO A SECOND ROUND OF BIN OPERATIONS AND SAVE WITH _PLUS SUFFIX
? ?? ?? ?? ?%
? ?? ?? ?? ?FPRINTF('\N\N\N**** %S: BIN OPERATIONS (TWO PASSES) ****\N\N\N', SUBJECT_LIST{S});? ?? ?? ?? ???

? ?? ?? ?? ?FNAME = [SUBJECT_LIST{S} '_ERPS.ERP'];??% RE-CREATE FILENAME FOR UNFILTERED ERP
? ?? ?? ?? ?ERP? ?= POP_LOADERP( 'FILENAME', FNAME, 'FILEPATH', DATA_PATH );? ?% LOAD THE FILE??


? ?? ?? ?? ?%% DIFFERENCE WAVE
? ?? ?? ?? ?% NOW MAKE THE DIFFERENCE WAVE, DIRECTLY SPECIFYING THE
? ?? ?? ?? ?% EQUATION THAT MODIFIES THE EXISTING ERPSET

? ?? ?? ?? ?ERP? ?? ?? ?= POP_BINOPERATOR( ERP, {'B3= B2-B1 LABEL RARE MINUS FREQUENT DIFFERENCE WAVE' });

? ?? ?? ?? ?ERP.ERPNAME = [ERP.ERPNAME '_DIFF'];??% NAME FOR ERPSET MENU
? ?? ?? ?? ?IF (SAVE_EVERYTHING)
? ?? ?? ?? ?? ? POP_SAVEMYERP(ERP, 'ERPNAME', ERP.ERPNAME, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', DATA_PATH, 'WARNING', 'OFF');
? ?? ?? ?? ?END

? ?? ?? ?? ?% NOW WE WILL DO BIN OPERATIONS USING A SET OF EQUATIONS
? ?? ?? ?? ?% STORED IN THE FILE 'BIN_EQUATIONS.TXT', WHICH MUST BE IN
? ?? ?? ?? ?% THE HOME FOLDER FOR THE EXPERIMENT

? ?? ?? ?? ?ERP? ?? ?? ?= POP_BINOPERATOR( ERP, [HOME_PATH 'BIN_EQUATIONS.TXT']);

? ?? ?? ?? ?ERP.ERPNAME = [ERP.ERPNAME '_PLUS'];??% NAME FOR ERPSET MENU
? ?? ?? ?? ?POP_SAVEMYERP(ERP, 'ERPNAME', ERP.ERPNAME, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', DATA_PATH, 'WARNING', 'OFF');

? ?? ?

? ?? ?? ?? ?% SAVE THIS FINAL ERP IN THE ALLERP STRUCTURE.??THIS IS NOT
? ?? ?? ?? ?% NECESSARY UNLESS YOU WANT TO SEE THE ERPS IN THE GUI OR IF YOU
? ?? ?? ?? ?% WANT TO ACCESS THEM WITH ANOTHER FUNCTION (E.G., POP_GAVERAGER)

? ?? ?? ?? ?CURRENTERP? ?? ?? ?= CURRENTERP + 1;
? ?? ?? ?? ?ALLERP(CURRENTERP) = ERP;

? ?? ?? ?? ?IF (PLOT_PDFS)
? ?? ?? ?? ?? ? POP_PLOTERPS(ERP, [1 2], 1:17);
? ?? ?? ?? ?? ? POP_EXPORTERPLABFIGURE(ERP, 'FILEPATH', DATA_PATH, 'FORMAT', 'PDF', 'TAG', {'ERP_FIGURE' 'SCALP_FIGURE'});
? ?? ?? ?? ?END



? ?? ???END % END OF THE "IF/ELSE" STATEMENT THAT MAKES SURE THE FILE EXISTS

? ? END % END OF LOOPING THROUGH ALL SUBJECTS


? ? %% GRAND AVERAGE
? ? % MAKE A GRAND AVERAGE. THE FINAL ERP FROM EACH SUBJECT WAS SAVED IN
? ? % ALLERP, AND WE HAVE NSUBJ SUBJECTS, SO THE INDICES OF THE ERPS TO BE AVERAGED
? ? % TOGETHER ARE 1:NSUBJ
? ? % WE'LL ALSO CREATE A FILTERED VERSION AND SAVE IT

? ? ERP = POP_GAVERAGER( ALLERP? ?? ?? ?, ...
? ?? ?? ?? ?? ?? ?? ?? ?'ERPSETS'? ?? ? , [ 1:NSUBJ], ...
? ?? ?? ?? ?? ?? ?? ?? ?'EXCLUDENULLBIN', 'ON', ...
? ?? ?? ?? ?? ?? ?? ?? ?'SEM'? ?? ?? ???, 'ON' ); %'CRITERION', 100 IS LEFT OVER FROM A PREVIOUS VERSION

? ? ERP.ERPNAME? ?? ???= 'GRAND_AVG';??% NAME FOR ERPSET MENU
? ? ERP? ?? ?? ?? ?? ? = POP_SAVEMYERP(ERP, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', HOME_PATH, 'WARNING', 'OFF');

? ? CURRENTERP? ?? ?? ?= CURRENTERP + 1;
? ? ALLERP(CURRENTERP) = ERP;

? ? ERP? ?? ?? ?? ?? ? = POP_FILTERP( ERP, 1:17 , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'CUTOFF' , 30? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'DESIGN' , 'BUTTER' , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'FILTER' , 'LOWPASS', ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ? 'ORDER'??, 2? ?? ???);
? ? ERP.ERPNAME? ?? ???= [ERP.ERPNAME '_30HZ'];??% NAME FOR ERPSET MENU
? ? ERP? ?? ?? ?? ?? ? = POP_SAVEMYERP(ERP, 'FILENAME', [ERP.ERPNAME '.ERP'], 'FILEPATH', HOME_PATH, 'WARNING', 'OFF');
? ? CURRENTERP? ?? ?? ?= CURRENTERP + 1;
? ? ALLERP(CURRENTERP) = ERP;


? ? % MEASURE THE MEAN AMPLITUDE FROM 300-600 MS IN BINS 2 AND 3, CHANNELS 11-13.
? ? % SAVE THE RESULTS IN A VARIABLE NAMED "VALUES" AND IN A FILE NAMED
? ? % "MEASURES.TXT" IN THE HOME FOLDER FOR THE EXPERIMENT.
? ? VALUES = POP_GETERPVALUES( ALLERP, [300 600], [2 3], 11:13? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'BASELINE'??, 'PRE'? ?? ?? ?? ?? ?? ?? ?, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'ERPSETS'? ?, [1:NSUBJ]? ?? ?? ?? ?? ???, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'FILENAME'??, [HOME_PATH 'MEASURES.TXT'], ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'FILENAME'??, 'MEASUREMENT'? ?? ?? ?? ? , ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'MEASURE'? ?, 'MEANBL'? ?? ?? ?? ?? ?? ?, ...
? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?'RESOLUTION', 1? ?? ?? ?? ?? ?? ?? ?? ? );


? ? FPRINTF('\N\N\N**** FINISHED ****\N\N\N');??
已有 1 人评分水晶 收起 理由
tobeornot + 1 赞一个!

总评分:?水晶 + 1? ?查看全部评分

已有之事,后必再有!
您需要登录后才可以回帖 登录 | 注册

bottom

Powered by Discuz! X2

? 2001-2011 Template By Yeei. Comsenz Inc.

回顶部