// This function plot Histogram of Radial distance of // Branch Point and Terminal Point // Outputs: Historgram Plots // Inputs: vector_distance_Apic_BP_a // vector_distance_Basal_BP_a // vector_distance_Apic_TP_a // vector_distance_Basal_TP_a // written by Jose Gomez, March 2006, jfcgomez@ull.es objref Hbox_Hista objref hist1,hist2,hist3, hist4, hist5, hist6, g1,g2,g3,v1,v4 hist1=vector_distance_Apic_BP_a.c hist2=vector_distance_Basal_BP_a.c hist3=vector_distance_Apic_TP_a.c hist4=vector_distance_Basal_TP_a.c hist5=L_TP_apical_a.c hist6=L_TP_basal_a.c proc hist_graph() { \$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_Hista=new HBox() Hbox_Hista.intercept(1) g1=new Graph() g1.size(-600,350,0,15) g2=new Graph() g2.size(-650,400,0,15) g3=new Graph() g3.size(-1500,750,0,15) hist_graph(hist1,g1,550,hist2,250) g1.label(0.1,0.95, "BP Radial from soma(um)") g1.label(0.1,0.75,"Apical") g1.label(0.90,0.75,"Basal") hist_graph(hist3,g2,600,hist4,400) g2.label(0.1,0.95, "TP Radial from soma(um)") g2.label(0.1,0.75,"Apical") g2.label(0.90,0.75,"Basal") hist_graph(hist5,g3,1500,hist6,700) g3.label(0.1,0.95, "TP Lenght to soma(um)") g3.label(0.1,0.75,"Apical") g3.label(0.90,0.75,"Basal") Hbox_Hista.intercept(0) Hbox_Hista.map()