|
|
@@ -1,34 +0,0 @@
|
|
|
-function [bin]=doBin(NL_x,NL_y,xbin)
|
|
|
-%this subfunction is doing the binning for nonLinearity
|
|
|
-%get different bins
|
|
|
-
|
|
|
-%% 1. get equal amount of observation inside a bin
|
|
|
-%standard score:
|
|
|
-x=prctile(NL_x,linspace(0,100,xbin+1)); %you want equal amount of observations in each bin!And not one bin of only e.g. 10 observations.
|
|
|
-%I want to now where the percentil are
|
|
|
-%prctile gives you now where the edges of the corresponding percentile rank
|
|
|
-%are.
|
|
|
-
|
|
|
-x(end)=x(end)+x(end)/100; %the last bin should be a bit bigger so that the last observation
|
|
|
-%can also go inside it
|
|
|
-[observations,indx]=histc(NL_x,x);
|
|
|
-
|
|
|
-%little trick so that you have 40 equal bin and not one bin at the end
|
|
|
-%that is empty or just has one value
|
|
|
-observations=observations(1:end-1);
|
|
|
-%make average in y axis
|
|
|
-NL_ytemp=cell(1,length(observations));
|
|
|
-for mm=1:length(NL_ytemp)
|
|
|
- NL_ytemp{mm}=NL_y(indx==mm);
|
|
|
-end
|
|
|
-%make average in x axis
|
|
|
-x_AVG=cell(1,length(observations));
|
|
|
-for mm=1:length(x_AVG)
|
|
|
- x_AVG{mm}=NL_x(indx==mm);
|
|
|
-end
|
|
|
-
|
|
|
-%get the average values per bin
|
|
|
-bin.NL_ybin=cellfun(@nanmean,NL_ytemp);
|
|
|
-bin.NL_xbin=cellfun(@nanmean,x_AVG);
|
|
|
-
|
|
|
-end
|
