Carpenter plots of the bayesian estimator

>> t_reaccion=bayestim_medias_nosave(10^5,10,1000,2,0,1,10^-3,.995);
>> hist(t_reaccion,0:1000)

>> t_reaccion2histog_carpenter(t_reaccion,20,’r.-‘);

Advertisement

Mistake creating the system matrix

In previous posts, I made a mistake: I used the transposed of the chemical connectivity matrix. Fortunately, it seems there is not a fundamental difference:

Transposed,

>> M=conectividad2matrizsistema(todas.A_chem’,todas.A_ej,todas.GABA,1,1,30);
>> solucion=matrizsistema2solucion(M,V0);
>> plot(solucion.autovalores,’.’)

Not transposed:

>> M=conectividad2matrizsistema(todas.A_chem,todas.A_ej,todas.GABA,1,1,30);
solucion=matrizsistema2solucion(M,V0);
plot(solucion.autovalores,’.’)

Excitation of neurons connected to “the nose”

sensoriales_nariz=find(todas.S(:,1));
>> V0=zeros(279,1);
>> V0(sensoriales_nariz)=1;
>> solucion=matrizsistema2solucion(M,V0);
>> V=solucion2voltajes(solucion,0:.001:.5);
>> [s,orden]=sort(todas.pos_real);
>> imagesc(V(orden,:))

Note that neurons do not have the usual order. They are ordered by position in the worm.

>> imagesc(log10(abs(V(orden,:))))

Longer times:

>> V=solucion2voltajes(solucion,0:.005:2);

I do not see any propagation…