Confocal images and boxplots from my results in strokelitude

on Tuesday, May 15th, 2018 12:26 | by

Confocal image MAX stack of one of the brains at 20x

 

 

and at 40x

 

In this link we have a video of a 3D stainning pattern          zoomed_CC

 

In addition I add here teh boxplots from the final results of the Ping Pong ball setup with these experiments

Print Friendly, PDF & Email

Stainning c105;;c232

on Monday, May 14th, 2018 11:22 | by

 

The first figure shows each of the central complex ring neurons types (Martín Pena et al., 2014). The c105-G4 targets the R1 neurons and the c232-G4 targets the R2 and the R4d neurons

This is the c105-G4 stainning from Martín Pena et al., 2014

232-G4 stainning from Kahsai et al., 2012

 

Axel stainning from c232-G4 alone

Axel stainning from c105-G4 alone

Axel stainning from both drivers together. I would say it really contains both driver lines.

This are both driver lines together as well from Axel. To me it seems that only c105 is present

This are my stainnings at the fluorescence microscope (no confocal). This is to show that in all of the 10-12 brains I have looked at, they all had the c232 pattern present

In addition, they had many more neurons outside from the central complex which I believe belong to the c105-G4 line. This is my only proof to show that c105 is also present, since the R1 neurons seem to be hidden when R2 neurons are stained.

I was also looking to the youtube video you have online, Björn. To me it seems I can only see the R1 ring neuron from the c105

Print Friendly, PDF & Email

EMD with ICA to one sample of torque trace

on Monday, March 19th, 2018 2:35 | by

Trace segment from Maye et al. 2007 in the uniform arena

 

Trace after filtering by selecting the first 8 IMFs (intrinsic mode functions) from EMD (empirical mode decomposition). Since the signal should be quite clean I do not take out the first IMFs. The last IMFs, however, are too slow and change the baseline to much

Some examples of the IMFs obtained from EMD:

Since this is separating behavior adapting to the data intrinsic time scales I am now thinking of analysing with the SMAP algorithm to see if the behavior is more or less nonlinear at certain time scales.

In addition, I have thought of using ICA (Independent component analysis), that is an algorithm famous for the blind source separation problem by extracting the most independent signals from the input signals (in this case the IMFs which are different time scales of the behavior). So the ICs should consist of mixes of different time scales that are correlated together and thus belong (but not necessarily) to the same action/movement module/… Here a few ICs (from 10). My idea is that muscles might coordinate independently between ICs but coordinated within ICs. However to prove that is not that easy I guess

Print Friendly, PDF & Email

quality test before strokelitude experiments

on Monday, January 29th, 2018 12:40 | by

Some traces from this week just so that you have an idea how do they look like. To me they are not the optimal traces I expected. But one can see some signal there. I will start the screen hoping to get enough good traces without too much work.

what do you think is the best quality control for accepting a trace for the analysis or not. I was thinking the 3D mapping gives a good hint but without quantification.

In addition I was writting to Andi Straw to solve the issue of running two cameras in the same computers, but now answer

Print Friendly, PDF & Email

Better analysis of slope dependence on flight length and further insights in the SMAP code

on Monday, December 4th, 2017 1:02 | by

This is another way of showing how the slope of the SMAP analysis varies with length. I have choped the time series in 4 chunks and saw what was the slope for the chunks and for the whole time series. I think this is the best statistical way of doing it, it should not depend on what line was tested. I cannot see any effect.

 

Here below is just to show what I found out in the code. What I thought that theta was controlling for nonlinearity in the model, for me it seems rather a control for under- overfitting in the model. So I might need to read the paper and see what do they say, and if it is the same as what I found out in the code.

 

Print Friendly, PDF & Email

Do not agree with Sathish?

on Monday, November 27th, 2017 2:43 | by

So this is the final graph assuming that 20Hz is the sampling rate. Sathish was not sure what it was and he said he will check.

 

 

 

To have a better overview how the length of the flight compares with the slope obtained from the SMAP. There is no big correlation whatsoever with around 100 flies. Sathish found this in his thesis with data from seventy something WTB, but to me it seems like an anecdotal result.

Here the same as above, just for showing the fit line.

Print Friendly, PDF & Email

All sathish data analysis

on Monday, November 20th, 2017 1:47 | by

This are the results of all of the flies analyzed from Sathish. I just need to know the frequency of the acquisition to exclude the ones that are below 6 minutes.

 

 

 

Print Friendly, PDF & Email

Sathish data

on Monday, November 6th, 2017 1:37 | by

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.

Print Friendly, PDF & Email

Recurrence quantitative analysis

on Monday, October 30th, 2017 11:43 | by

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.

 

Automat

 

One stripe

Openloop

Uniform

Print Friendly, PDF & Email

SMAP results

on Monday, October 9th, 2017 2:34 | by

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.

In addition, I have done some animations of the attractors that I have posted on slack because of size.

Print Friendly, PDF & Email