小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
初学者实在不懂，大家看看到底哪里有问题
&& demopath='/User/yng/Documents/demo';
Undefined function or variable 'EEMs'.
这个EEMs是文件夹，要打开它里面csv格式文件，但直接就出现这样的警告。我用的是matlab R2013a for mac。我正在试运行这个程序。下面是所下载程序的全部，在pc平台上的。红色是正在运行出问题的部分。程序中demopath='D:\Collaborations\PARAFAC tutorial\drEEM parts\demo'; 用右斜杠，但cd()用的是左斜杠，这是不有问题？
% Create an EEM dataset: using the drEEM and N-way toolboxes for MATLAB
% This file shows example code for assembling a multi-way EEM dataset from
% raw (instrument) data files, including applying spectral corrections
% and correction for inner filter effects and/or dilution.
% It accompanies the tutorial in Appendix A of the reference below
% This mfile is part of the drEEM toolbox. Please cite the toolbox
% as follows:
% Murphy K.R., Stedmon C.A., Graeber D. and R. Bro, Fluorescence
%& &&&spectroscopy and multi-way techniques. PARAFAC, Anal. Methods, 2013,
%& &&&DOI:10.160e.
% Copyright (C) 2013 Kathleen R. Murphy
% The University of New South Wales
% Dept Civil and Environmental Engineering
% Water Research Center
% UNSW 2052
% $ Version 0.1.0 $ September 2013 $ First Release
%% Set up toolboxes and demo on MATLAB path (edit paths as required)
%Toolboxes
addpath('D:\PARAFAC\drEEM tutorial\dreem') %edit for correct file destination
addpath('D:\PARAFAC\drEEM tutorial\dreem\nway')& & %edit for correct file destination
%Demo - put a copy of the demo dataset in a writeable folder, e.g.
demopath='D:\Collaborations\PARAFAC tutorial\drEEM parts\demo'; %edit for correct file destination，正在运行
%% Stage 1:& &read in raw data files
cd()正在运行
filetype=1;ext = 'csv';RangeIn='A1..AU105';headers=;display_opt=0;outdat=2;
=readineems(filetype,ext,RangeIn,headers,display_opt,outdat);
Ex=Exmat(1,:); %Since all files have the same excitation wavelengths
Em=Emmat(:,1); %Since all files have the same emission wavelengths
filetype=1;ext = 'csv';RangeIn='A1..AU105';headers=;display_opt=0;outdat=2;
=readineems(filetype,ext,RangeIn,headers,display_opt,outdat);
Exb=Exmat_b(1,:); %Since all files have the same excitation wavelengths
Emb=Emmat_b(:,1); %Since all files have the same emission wavelengths
filetype='R275';ext = 'csv';RangeIn='A2..B167';display_opt=0;outdat=2;
=readinscans(filetype,ext,RangeIn,display_opt,outdat);
RamEx=275; %note landa && 350, so we will need to calculate RamOpt below
filetype='Abs';ext = 'csv';RangeIn='A1..B521';display_opt=0;outdat=0;
=readinscans(filetype,ext,RangeIn,display_opt,outdat);
Excor=csvread('xc06se06n.csv');
Emcor=csvread('mcorrs_4nm.csv');
%Alternative: use xlsread if the file contains text
%Excor=xlsread('xc06se06n.csv','xc06se06n','A1..B75');
%Emcor=xlsread('mcorrs_4nm.csv','mcorrs_4nm','A1..B104');
filetype='QSuv275';ext = 'csv';RangeIn='A2..B75';display_opt=0.5;outdat=0;
=readinscans(filetype,ext,RangeIn,display_opt,outdat);
filetype='QSuv350';ext = 'csv';RangeIn='A2..B75';display_opt=0.5;outdat=0;
=readinscans(filetype,ext,RangeIn,display_opt,outdat);
%% Align data according to records held in a Sample Log
cd(demopath)
% Import sample log
filename='D:\Collaborations\PARAFAC tutorial\drEEM parts\demo\SampleLog_PortSurveyDemo.csv';
strings=; %specify which columns in the log contain text
SampleLog=readlogfile(filename,strings);
% Align data according to the SampleLog. The EEMs are listed in SampleLog.Log_EEMfile
% first align raw data from SampleLog with filelist_eem
sites=& && &alignds(SampleLog,{'EEMfile',filelist_eem},{'Site'}); %Site names are listed in SampleLog.Site
cruises=& & alignds(SampleLog,{'EEMfile',filelist_eem},{'Cruise'}); %Cruise names are listed in SampleLog.Cruise
dates=& && &alignds(SampleLog,{'EEMfile',filelist_eem},{'AnalDate'});
Q=& && && & alignds(SampleLog,{'EEMfile',filelist_eem},{'QS_Slope'});
sampleID=& &alignds(SampleLog,{'EEMfile',filelist_eem},{'SampID'});
replicates= alignds(SampleLog,{'EEMfile',filelist_eem},{'Rep'});
dilfac=& &&&alignds(SampleLog,{'EEMfile',filelist_eem},{'dilutionfactor'});
% next align previously loaded datasets with filelist_eem (ABS, Raman, Blanks, etc).
Sabs=& &&&alignds(SampleLog,{'EEMfile',filelist_eem},{'ABSfile',filelist_abs,S_abs}); %ABS data files are listed in SampleLog.ABSfile
B= alignds(SampleLog,{'EEMfile',filelist_eem},{'BlankFile',filelist_b,X_b}); %blank EEM iles are listed in SampleLog.BlankFile,Blist is used later
Sr=& && & alignds(SampleLog,{'EEMfile',filelist_eem},{'RamanFile',filelist_R,S_R});
Sqsuv275= alignds(SampleLog,{'EEMfile',filelist_eem},{'Qw',filelist_qsuv275,S_qsuv275});
Sqsuv350= alignds(SampleLog,{'EEMfile',filelist_eem},{'Qw',filelist_qsuv350,S_qsuv350});
%% Correct the EEMs
% Attach Wavelength headers to scans
A=;& && && & %add wavelengths to Absorbance scans
Qw275=; %2D Matrix of Raman scans matched to QS dilution series
Qw350=; %2D Matrix of Raman scans matched to QS dilution series
% Eliminate data below or above the wavelength range of Ex and Em correction files
Em_in=Em(Em&=682);
X_in=X(:,Em&=682,:);&&
B_in=B(:,Em&=682,:);
%Turn X_in into a dataset, DS;
DS=assembledataset(X_in,Ex,Em_in,'AU','filelist',filelist_eem,) %#ok&*NOPTS&
% Method 1 - using Raman scans at Ex=275 nm. First we need to locate an optimal
% integration range that will capture the peak at 275 nm
W=;& && && && &&&%2D Matrix of matched (275 nm) Raman scans
RamMat=;& && &&&%2D Matrix of unmatched Raman scans (corresponding to filelist_R)
= ramanintegrationrange(RamMat,filelist_R,RamEx,,0.1);
%From IRmed =&&we choose an integration range of 290-320 nm
%noting that no correction factors exist for Em&290.
W_in=W(:,wave_R&=290);& && && && && &
Qw_in=Qw275(:,wave_qsuv275&=290);& && && && && &
=fdomcorrect(DS,DS.Ex,DS.Em,Emcor,Excor,W_in,RamOpt,A,B_in,,Q,Qw_in);
% Method 2 - using Raman scans at Ex=350 nm extracted from the MilliQ blanks.
%First determine the optimal Raman integration range using the blanks
%in X_b which contain no duplicate EEMs. The best range is 374-426 nm (IRmed)
%however, we have no QS blank measurements below 375 nm, so choose 375 as the
%lower range.
W350 = ;& &%Extract the water Raman scans at Ex=350
= ramanintegrationrange(W350,filelist_b,350,,0.1)
RamOpt=; %calculated Raman integration range: 374-426nm
RamOpt= %default Raman integration range: 381-426nm
W = ;& &%Extract the water Raman scans at Ex=350
=fdomcorrect(DS,DS.Ex,DS.Em,Emcor,Excor,W,RamOpt,A,B_in,,Q,Qw);
%Calculate the QS/RU conversion factors for the three QS dilution series
%and check the ratios are stable
QSonRU=mean(unique(QS_RU2));
fprintf('\n\nQS/RU conversion factors by Method 2:\n')
disp(unique(QS_RU2))
fprintf('...having a mean value of: \n')
disp(QSonRU) % (88.96-90.05 with the default Raman integration range)
% Method 3 - Select the 'Best available' corrected datasets
XcQS=XcQS1;& && &&&%Method 1
XcRU=XcQS1/QSonRU; %Method 1 divided by mean QS/RU ratio from Method 2
% Optional - correct for dilution if necessary
% If samples were diluted before measuring EEMs and Abs scans, divide the
% corrected EEMs by the sorted dilution factors (dilfac NOT df).
XcQS_df1=undilute(XcQS,dilfac);
XcRU_df1=undilute(XcRU,dilfac);
fprintf('\n\n');
fprintf('& && & sampleID& && & diluted& && &undiluted\n');
f350_450=;
disp(f350_450); %compare diluted vs undiluted at 350/450
%% Save the results
save PortSurveyDataset.mat XcQS XcRU Ex Em_in filelist_eem sites replicates cruises dates sampleID
xlswrite('OutData_FDOM.xlsx',cellstr(filelist_eem),'filelist')
cd(demopath)
save('CorrData_FDOM.mat','filelist_eem','IR','IRmed','IRdiff','RamEx','-append')
xlswrite('OutData_FDOM.xlsx',cellstr(filelist_eem),'filelist')
save('CorrData_FDOM.mat','filelist_eem','QSonRU','-append')
xlswrite('OutData_FDOM.xlsx',cellstr(filelist_eem),'filelist')
%%% END OF DATA CORRECTION TUTORIAL
% This mfile is part of the drEEM toolbox. Please cite the toolbox
% as follows:
% Murphy K.R., Stedmon C.A., Graeber D. and R. Bro, 2013, Fluorescence
%& &&&spectroscopy and multi-way techniques. PARAFAC, Anal. Methods, 5,
%& &&&,& &DOI:10.160e.
% drEEM_demo_020: Copyright (C) 2014 Kathleen R. Murphy
% $ Version 0.2.0 $ May 2014 $ Second Release
% Copyright (C) 2014 KR Murphy,
% Water Environment Technology
% Chalmers University of Technology
% Department of Civil and Environmental Engineering
% Water Environment Technology
% Sven Hultins gata 8
% 412 96 G�teborg
% $ Version 0.2.0 $ May 2014 $ Second Release
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研分类排行榜
作者：木木
作者：何杨
作者：何杨
作者：苗元威
站长们都信赖的交易中介平台，安全、高效、便捷。
中小企业网络营销第一平台，解决流量、转化率、品牌。
高收益、移动广告、弹窗CPM、点击CPC、包月广告。
提供公众号出售、求购、代售等交易中介以及增值服务。
增值电信业务经营许可证：苏B2-
编辑***：5
A5创业网 版权所有.
扫一扫关注最新创业资讯