I am having trouble subdiving into sub-groups my data. The first thing I was trying to do is selecting group of data every 30 degrees from 0 to 360. It makes 12 groups of data. The secong thing I would like to do is from this group of data do another selection. For each 30 degrees I would like to divide into 5 groups. Thus have 12 Groups of data and from each group have 5 sub groups. I would like to do, for every 30 degrees residual groups (azimuth), the median of the residuals every 50m (epicenter distance). For example for the group of data with residuals that are between 0-30 degrees I would like to do the median of these residuals with epicentral distance from 0-50, 50-100,..., 250-300. This is the code I tried but it says that the indexe exceeds the matrix, I tried other things but have the same problem:
N_pieces3=6;
cheese_angle3=300/N_pieces3;
N_pieces=12;
cheese_angle=360/N_pieces;
for nps=1:N_pieces
resp1=resp(azmp1>=cheese_angle*(nps-1) & azmp1<=cheese_angle*nps);
depip1=depip(azmp1>=cheese_angle*(nps-1) & azmp1<=cheese_angle*nps);
for c=1:N_pieces3
resp3=resp1(depip>=cheese_angle3*(c-1) & depip<=cheese_angle3*c);
depip3=depip1(depip>=cheese_angle3*(c-1) & depip<=cheese_angle3*c );
angle3(c)=cheese_angle3*((c-1)+.5);
Med(nps)=median(resp3)-medAs;
end
angle2=cheese_angle3*((nps-1)+.5);
end
At the end I would like to do a 3D plotting with the residuals median depending on the azimuth and the epicentral distance. If possible a 3D polar scatter. Thank you.
