// This function plot Histogram of Radial distance of
// Branch Point and Terminal Point
// Outputs: Historgram Plots
// Inputs: vector_distance_Apic_BP_b
// vector_distance_Basal_BP_b
// vector_distance_Apic_TP_b
// vector_distance_Basal_TP_b
// written by Jose Gomez, March 2006, jfcgomez@ull.es
objref Hbox_Histb
objref hist1b,hist2b,hist3b, hist4b, hist5b, hist6b, g1b,g2b,g3b,v1,v4
hist1b=vector_distance_Apic_BP_b.c
hist2b=vector_distance_Basal_BP_b.c
hist3b=vector_distance_Apic_TP_b.c
hist4b=vector_distance_Basal_TP_b.c
hist5b=L_TP_apical_b.c
hist6b=L_TP_basal_b.c
proc hist_graphb() {
$o1=$o1.histogram(0,$3,50)
$o1=$o1.reverse
v1=new Vector($o1.size())
v1.indgen(-($3+50),50)
$o4=$o4.histogram(0,$5,50)
v4=new Vector($o4.size())
v4.indgen(0,50)
$o1=$o1.append($o4)
v1=v1.append(v4)
$o1.plot($o2,v1)
}
Hbox_Histb=new HBox()
Hbox_Histb.intercept(1)
g1b=new Graph()
g1b.size(-600,350,0,15)
g2b=new Graph()
g2b.size(-650,400,0,15)
g3b=new Graph()
g3b.size(-1500,750,0,15)
hist_graphb(hist1b,g1b,550,hist2b,250)
g1b.label(0.1,0.95, "BP Radial from soma(um) (Pruning)")
g1b.label(0.1,0.75,"Apical")
g1b.label(0.90,0.75,"Basal")
hist_graphb(hist3b,g2b,600,hist4b,400)
g2b.label(0.1,0.95, "TP Radial from soma(um) (Pruning)")
g2b.label(0.1,0.75,"Apical")
g2b.label(0.90,0.75,"Basal")
hist_graph(hist5b,g3b,1500,hist6b,700)
g3b.label(0.1,0.95, "TP Lenght to soma(um) (Pruning)")
g3b.label(0.1,0.75,"Apical")
g3b.label(0.90,0.75,"Basal")
Hbox_Histb.intercept(0)
Hbox_Histb.map()