Skip to content

Updated old naming of EMGHist* to newer MotionHist*. #427

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
bli40 wants to merge 1 commit into
base: master
Choose a base branch
from
Open

Updated old naming of EMGHist* to newer MotionHist*. #427

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
46 changes: 23 additions & 23 deletions detectors/detectStates/SleepScoreMaster/ClusterStates_MakeFigure.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,10 +21,8 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
[zFFTspec,mu,sig] = zscore(log10(swFFTspec)');
if sum(isinf(log10(thFFTspec(:))))==0
[~,mu_th,sig_th] = zscore(log10(thFFTspec)');
disp('im in the if')
else %For Theta over PSS (ThIRASA)
[~,mu_th,sig_th] = zscore(thFFTspec);
disp('im in the else')
[~,mu_th,sig_th] = zscore((thFFTspec)');
end

viewwin =[t_clus(1) t_clus(end)];
Expand All @@ -36,28 +34,30 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
axis xy
set(gca,'YTick',(log2([1 2 4 8 16 32 64 128])))
set(gca,'YTickLabel',{'1','2','4','8','16','32','64','128'})
clim([3.5 6.5])
clim(double([min(mu)-2*max(sig) max(mu)+2*max(sig)]))
% clim([3.5 6.5])
clim(double([min(mu)-1.5*max(sig) max(mu)+1.5*max(sig)]))
xlim(viewwin)
colorbar('east')
ylim([log2(swFFTfreqs(1)) log2(swFFTfreqs(end))+0.2])
set(gca,'XTickLabel',{})
ylabel({'swLFP','f (Hz)'})
title([recordingname,': State Scoring Results']);

subplot(8,1,3)
imagesc(t_clus,log2(thFFTfreqs),log10(thFFTspec))
axis xy
set(gca,'YTick',(log2([1 2 4 8 16 32 64 128])))
set(gca,'YTickLabel',{'1','2','4','8','16','32','64','128'})
caxis(double([min(mu_th)-2*max(sig_th) max(mu_th)+2*max(sig_th)]))
% clim([3.5 6.5])
clim(double([min(mu_th)-1.5*max(sig_th) max(mu_th)+1.5*max(sig_th)]))
xlim(viewwin)
%colorbar('east')
% colorbar('east')
ylim([log2(thFFTfreqs(1)) log2(thFFTfreqs(end))+0.2])
ylabel({'thLFP','f (Hz)'})
set(gca,'XTickLabel',{})

subplot(8,1,4)
%plot(t_clus,-IDX,'LineWidth',2)
% plot(t_clus,-IDX,'LineWidth',2)
hold on
plot(states{1}',-1*ones(size(states{1}))','k','LineWidth',8)
plot(states{2}',-2*ones(size(states{2}))','b','LineWidth',8)
Expand All @@ -71,7 +71,7 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
subplot(6,1,4)
hold on
plot(t_clus,broadbandSlowWave,'k')
%plot(synchtimes',thresh*ones(size(synchtimes))','r')
% plot(synchtimes',thresh*ones(size(synchtimes))','r')
ylabel('SW')
box on
ylim([0 1])
Expand All @@ -81,7 +81,7 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
subplot(6,1,5)
hold on
plot(t_clus,thratio,'k')
%plot(synchtimes',thresh*ones(size(synchtimes))','r')
% plot(synchtimes',thresh*ones(size(synchtimes))','r')
ylabel('Theta')
box on
ylim([0 1])
Expand All @@ -91,7 +91,7 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
subplot(6,1,6)
hold on
plot(t_clus,EMG,'k')
%plot(synchtimes',thresh*ones(size(synchtimes))','r')
% plot(synchtimes',thresh*ones(size(synchtimes))','r')
ylabel('EMG')
box on
ylim([0 1])
Expand Down Expand Up @@ -125,14 +125,14 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
xlabel('Broadband Slow Wave')
title('Step 1: Broadband for NREM')


subplot(3,2,3)
subplot(3,2,3)
hold on
bar(EMGhistbins(EMGhistbins>EMGthresh),EMGhist(EMGhistbins>EMGthresh),'FaceColor','k','barwidth',0.9,'linewidth',1)
bar(EMGhistbins(EMGhistbins<=EMGthresh),EMGhist(EMGhistbins<=EMGthresh),'FaceColor',0.9*[1 1 1],'barwidth',0.9,'linewidth',1)
plot([EMGthresh EMGthresh],[0 max(EMGhist)],'r','LineWidth',1)
bar(MotionHistBins(MotionHistBins>MotionThresh),MotionHist(MotionHistBins>MotionThresh),'FaceColor','k','barwidth',0.9,'linewidth',1)
bar(MotionHistBins(MotionHistBins<=MotionThresh),MotionHist(MotionHistBins<=MotionThresh),'FaceColor',0.9*[1 1 1],'barwidth',0.9,'linewidth',1)
plot([MotionThresh MotionThresh],[0 max(MotionHist)],'r','LineWidth',1)
xlabel('EMG')
title('Step 2: EMG for Muscle Tone')

subplot(3,2,5)
hold on
bar(THhistbins(THhistbins>=THthresh),THhist(THhistbins>=THthresh),'FaceColor','r','barwidth',0.9,'linewidth',1)
Expand All @@ -149,11 +149,11 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
% hold on
plot(broadbandSlowWave(IDX==3,1),EMG(IDX==3),'b.','markersize',0.1)
hold on
plot(broadbandSlowWave(EMG>EMGthresh & IDX==1,1),EMG(EMG>EMGthresh & IDX==1),'k.','markersize',0.1)
plot(broadbandSlowWave(EMG<EMGthresh & IDX==1|IDX==5,1),EMG(EMG<EMGthresh & IDX==1|IDX==5),...
plot(broadbandSlowWave(EMG>MotionThresh & IDX==1,1),EMG(EMG>MotionThresh & IDX==1),'k.','markersize',0.1)
plot(broadbandSlowWave(EMG<MotionThresh & IDX==1|IDX==5,1),EMG(EMG<MotionThresh & IDX==1|IDX==5),...
'.','Color',0.8*[1 1 1],'markersize',0.1)
plot(swthresh*[1 1],get(gca,'ylim'),'r','LineWidth',1)
plot(swthresh*[0 1],EMGthresh*[1 1],'r','LineWidth',1)
plot(swthresh*[0 1],MotionThresh*[1 1],'r','LineWidth',1)
xlabel('Broadband SW');ylabel('EMG')

%[map,C] = hist3([thratio(NREMtimes==0,1),EMG(NREMtimes==0)],[40 40]);
Expand All @@ -167,8 +167,8 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
hold on
plot(thratio(IDX==5,1),EMG(IDX==5,1),'r.','markersize',0.1)
xlabel('Narrowband Theta');ylabel('EMG')
plot(THthresh*[1 1],EMGthresh*[0 1],'r','LineWidth',1)
plot([0 1],EMGthresh*[1 1],'r','LineWidth',1)
plot(THthresh*[1 1],MotionThresh*[0 1],'r','LineWidth',1)
plot([0 1],MotionThresh*[1 1],'r','LineWidth',1)

saveas(gcf,[figloc,recordingname,'_SSCluster2D'],'jpeg')
%saveas(gcf,['/Users/dlevenstein/Code Library/SleepScoreDevelopment/StateScoreFigures/','ThetaEMGExample'],'jpeg')
Expand Down Expand Up @@ -200,8 +200,8 @@ function ClusterStates_MakeFigure(SleepState,basePath,noprompts)
title('Step 1: Broadband for NREM')
subplot(3,3,4)
hold on
bar(EMGhistbins,EMGhist,'FaceColor','none','barwidth',0.9,'linewidth',2)
plot([EMGthresh EMGthresh],[0 max(EMGhist)],'r')
bar(MotionHistBins,MotionHist,'FaceColor','none','barwidth',0.9,'linewidth',2)
plot([MotionThresh MotionThresh],[0 max(MotionHist)],'r')
xlabel('EMG')
title('Step 2: EMG for Muscle Tone')
subplot(3,3,7)
Expand Down
15 changes: 8 additions & 7 deletions detectors/detectStates/SleepScoreMaster/SleepScoreMaster.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -302,15 +302,16 @@
save(bz_sleepstatepath,'SleepState','-v7.3');

%% MAKE THE STATE SCORE OUTPUT FIGURE
%ClusterStates_MakeFigure(stateintervals,stateIDX,figloc,SleepScoreMetrics,StatePlotMaterials);
ClusterStates_MakeFigure(SleepState,basePath,noPrompts);
% ClusterStates_MakeFigure(stateintervals,stateIDX,figloc,SleepScoreMetrics,StatePlotMaterials);
% ClusterStates_MakeFigure_km(SleepState,basePath,noPrompts);

try
ClusterStates_MakeFigure_km(SleepState,basePath,noPrompts);
disp('Figures Saved to StateScoreFigures')
catch
disp('Figure making error')
end
% try
% ClusterStates_MakeFigure(SleepState,basePath,noPrompts);
% disp('Figures Saved to StateScoreFigures')
% catch
% disp('Figure making error')
% end

%% JOIN STATES INTO EPISODES

Expand Down