on Monday, November 6th, 2017 1:37 | by Christian Rohrsen
I have analyzed the c105+c232 > tnt data from Sathish. There are a total of 43 flies, although many of them only flew for a few minutes, and therefore should be discarded. Below all of the individual fly scores. I need to analyse now the other groups. This is btw the modified data set, whatever that means for Sathish.
This is a video of the projection from the torque data from Maye et al. 2007. The spikes are not that well sorted in this case as in the Strokelitude. I guess this is because the spikes do not look so smooth.
on Monday, October 30th, 2017 11:43 | by Christian Rohrsen
This is an example of a recurrence plot analysis. In the first graph is shown in single point in time in the optimal embedding dimension and the distance to the other points. For the recurrence plot analysis it is needed to put a threshold to make it binary. This is the second graph. From this second graph one can count many parameters like determinism, laminarity and so on. From what I see, the plots from the Strokelitude as well as Bjoern´s flight simulator in Maye et al 2007 show similar pattern (kind of crosses with vertical and horizontal lines).
This is a measure of the Recurrence Quantitative Analysis of different groups. Recurrence threshold is a tricky and to some extent subjective measure, so this is why I tried two different ones.
DET: recurrence points that form a diagonal line of minimal length, the more diagonal, the more deterministic.
LMAX: Max diagonal line length or divergence. Sometimes considered as an estimator of max. Lyapunov exponent
ENT: Shannon entropy reflects the complexity of the system
TND: info about stationarity (trend)
LAM: Laminarity is related to laminar phases in the system (intermittency). It is tallied as vertical lines over a threshold.
TT: Trapping time, measuring the average length of vertical lines. Related to laminarity.
on Monday, October 16th, 2017 2:26 | by Christian Rohrsen
Video of the attractor projected from a chunk of flight trace. Here one can see the difference in the trajectories of up- and down spikes. So here is another way of spike sorting by the way :).
Coloured traces depending where the locate in the attractor projection. Green means that they lie outside and when they stay with the PCs around zero its in red. The fly1r above does not look so clean as the v8 fly below (the cleanest measure I have)
on Monday, October 9th, 2017 2:34 | by Christian Rohrsen
These are the results of the SMAP for the TNTxWTB. I also have done a few for the c105;;c232xWTB but there is not much to say. I would say that the cleanest lines show a bigger slope, but prone to subjectiveness.
on Monday, October 2nd, 2017 5:35 | by Christian Rohrsen
This is what I showed about one fly from this line showing the attractor.
This graph is what I forgot to show in the lab meeting. There are the 6 best traces from this same line. All of them selected ad hoc subjectively. The three best of them to my eyes are exactly the three above in the graph (v8-this one is the one shown in the picture with the attractors-,v4,v2). What does this means? the ones with better traces (subjectively) have higher offset in the phi, this means that they are more predictable overall (maybe because more resolution?). In addition, they show higher slopes which means more nonlinearity.
on Monday, April 4th, 2016 3:04 | by Christian Rohrsen
In the last meeting Björn proposed to do correlations of cumulative increasing bins. He said to do that taking the zeroth point (last library point where prediction is still not done) and use it for having a potential 1 of correlation coefficient at the beginning. I could not do that because I didnt save the zeroth points, and this will be a bit tedious and confusing considering that many flies were tested and probably the order is not 100% known. Thus, I just did the bins skipping this zeroth point. After all, we should see something similar with this one. First two graphs: c105;;c232>TNT (first and second prediction point), second: WTBxTNT, third: WTBxc105;;c232.
Examples of how each of the flies look like. So they are basically cumulative bins with each single fly (each in different colour). Just to have a hint how does the singularity looks like. Second thing I did is normalize the to have a range from -1 to 1 all of them (I have to double check the range in the script) and also setting them at a starting point of zero. I did this because we do not want to have differences in the correlation coefficient due to a different offset of the values of the wing beat and neither because of the starting point (if the fly was already flying to the right full gas, then it could be that it has an influence in the following prediction).
c105;;c232 –> first at starting at zero without normalizing and then with normalizing. The next is just the RMSE (not so important).
on Monday, March 21st, 2016 1:13 | by Christian Rohrsen
To see if there is an exponential decay in the prediction of the fly traces we did correlations of bins of 40 data points. We have 4 graphs (the last one merged) which consist of predictions at two different points with two different number of neighbours used for the prediction. So we have for each group sucesively: prediction at the first prediction point with the first number of neigbours, then the same with different number of neighbours. The last two are two different numbers of neighbours for the second prediction point. In the order: c105;;c232>TNT, TNTxWTB, c105;;c232xWTB
on Monday, March 14th, 2016 1:48 | by Christian Rohrsen
This is now the results from trying to predict the fly behavior doing ensembles of two predictions for the next 200 data points at two different points of the traces.
From what we see here, there is no “flattening” in the prediction of the fly when the neurons under c105 and c232 are targeted by TNT. This is done with around 14/15 flies for each group with two predictions in each ensemble of the two starting points. That makes a total of 15flies x 3 groups x 2 starting points for prediction x 2 predictions per ensemble = 180 prediction traces. Now I´m trying to calculate it by making correlations of bins in the prediction-observed for the same fly
on Monday, February 15th, 2016 2:25 | by Christian Rohrsen
This is the prediction analysis of flies in the platform under a 20min experiment under dark conditions. The number of experiments change drastically among groups because of technical problems: WTBxTNT is 4, WTBxc105;;c232 is 22, for the experimental line is 6 (c105;;c232>TNT), for the platform without flies is 10. I show the root mean squares and the correlation coefficient for each group.
This is the experimental group: c105;;c232>TNT.
The group without flies on the platform. I expect here to get a very good predictability overall:
on Monday, February 8th, 2016 12:11 | by Christian Rohrsen
These were done the week before last one but I could not upload it last time, so here are they. These are the results for predicting 4-5 flies of each group. Just one prediction from the middle of the time series for the 40 data points ahead in the future. First group is c105;;c232>TNT.
This is c105;;c232 x WTB
As I saw that the graph had so much zig-zag I told Pablo to make a bigger number of tested flies and this is what he is presenting today.
In addition I did analyze other parameters which are all saved under a PDF file below (Strokelitude). This contains some parameters with doubtfull processing which I still don´t trust so I have to find a better way for the calculation of it.