[Tommaso Salvitti]

Below the dataset, I've used

* Example generated by -dataex-. To install: ssc install dataex
clear
input str3 Litter str6 Strain str3 Strain_old byte Pnd double(Peso Temperatura)
"1B"  "B6_2"   "B6"   2   .    .
"2B"  "B6_2"   "B6"   2 1.1 27.6
"3B"  "B6_2"   "B6"   2 1.2 26.8
"4B"  "B6_2"   "B6"   2 2.6 29.3
"5B"  "B6_2"   "B6"   2 1.3 30.3
"6B"  "B6_2"   "B6"   2 1.5 27.7
"7B"  "B6_2"   "B6"   2 1.3 29.5
"8B"  "B6_2"   "B6"   2 1.3 28.1
"1B"  "B6_4"   "B6"   4 1.8 30.5
"2B"  "B6_4"   "B6"   4 1.6 29.6
"3B"  "B6_4"   "B6"   4 1.6 30.6
"4B"  "B6_4"   "B6"   4 1.7 30.4
"5B"  "B6_4"   "B6"   4 1.9 31.2
"6B"  "B6_4"   "B6"   4 2.3 32.1
"7B"  "B6_4"   "B6"   4 2.2 30.8
"8B"  "B6_4"   "B6"   4   2 31.8
"1B"  "B6_6"   "B6"   6 2.6 31.4
"2B"  "B6_6"   "B6"   6 2.2 31.4
"3B"  "B6_6"   "B6"   6 2.3 31.5
"4B"  "B6_6"   "B6"   6 2.6 30.2
"5B"  "B6_6"   "B6"   6 2.6 31.2
"6B"  "B6_6"   "B6"   6 3.3 32.5
"7B"  "B6_6"   "B6"   6 3.2 31.7
"8B"  "B6_6"   "B6"   6 2.9 32.5
"1B"  "B6_8"   "B6"   8 3.5 33.1
"2B"  "B6_8"   "B6"   8 2.9 31.7
"3B"  "B6_8"   "B6"   8 3.2 30.1
"4B"  "B6_8"   "B6"   8 3.3 32.7
"5B"  "B6_8"   "B6"   8 3.3 33.3
"6B"  "B6_8"   "B6"   8 4.5 31.8
"7B"  "B6_8"   "B6"   8 4.2 32.9
"8B"  "B6_8"   "B6"   8   4 32.4
"1B"  "B6_12"  "B6"  12 5.9 33.4
"2B"  "B6_12"  "B6"  12   5 32.9
"3B"  "B6_12"  "B6"  12 5.3 30.2
"4B"  "B6_12"  "B6"  12 4.7 33.5
"5B"  "B6_12"  "B6"  12 4.4 33.8
"6B"  "B6_12"  "B6"  12 6.7 34.4
"7B"  "B6_12"  "B6"  12 6.1 33.8
"8B"  "B6_12"  "B6"  12 5.9 34.3
"1F"  "FVB_2"  "FVB"  2 1.8 31.8
"2F"  "FVB_2"  "FVB"  2 1.6 29.7
"3F"  "FVB_2"  "FVB"  2 1.6 30.4
"4F"  "FVB_2"  "FVB"  2 1.8 30.3
"5F"  "FVB_2"  "FVB"  2 1.8 30.3
"6F"  "FVB_2"  "FVB"  2   2 29.7
"7F"  "FVB_2"  "FVB"  2 1.5 30.9
"8F"  "FVB_2"  "FVB"  2 1.8 30.8
"9F"  "FVB_2"  "FVB"  2 1.8 30.5
"10F" "FVB_2"  "FVB"  2 2.6 31.1
"1F"  "FVB_4"  "FVB"  4 2.6   31
"2F"  "FVB_4"  "FVB"  4 2.4 31.5
"3F"  "FVB_4"  "FVB"  4 2.1 31.3
"4F"  "FVB_4"  "FVB"  4 2.8 30.9
"5F"  "FVB_4"  "FVB"  4 2.6 31.1
"6F"  "FVB_4"  "FVB"  4 2.9 31.7
"7F"  "FVB_4"  "FVB"  4 2.2 30.8
"8F"  "FVB_4"  "FVB"  4 2.7 31.6
"9F"  "FVB_4"  "FVB"  4 2.6 30.8
"10F" "FVB_4"  "FVB"  4 3.6 31.3
"1F"  "FVB_6"  "FVB"  6 3.6 31.7
"2F"  "FVB_6"  "FVB"  6 3.4 32.1
"3F"  "FVB_6"  "FVB"  6 2.7 30.8
"4F"  "FVB_6"  "FVB"  6 3.8 32.5
"5F"  "FVB_6"  "FVB"  6 3.5 31.4
"6F"  "FVB_6"  "FVB"  6   4 32.6
"7F"  "FVB_6"  "FVB"  6 3.1 32.3
"8F"  "FVB_6"  "FVB"  6 3.5 32.8
"9F"  "FVB_6"  "FVB"  6 3.6 32.2
"10F" "FVB_6"  "FVB"  6 4.6 33.6
"1F"  "FVB_8"  "FVB"  8 4.3 33.4
"2F"  "FVB_8"  "FVB"  8 4.5 33.4
"3F"  "FVB_8"  "FVB"  8   3 30.8
"4F"  "FVB_8"  "FVB"  8   5 32.2
"5F"  "FVB_8"  "FVB"  8 4.2 32.1
"6F"  "FVB_8"  "FVB"  8 4.8 33.1
"7F"  "FVB_8"  "FVB"  8   4   32
"8F"  "FVB_8"  "FVB"  8 4.3 32.9
"9F"  "FVB_8"  "FVB"  8 4.5 33.5
"10F" "FVB_8"  "FVB"  8 5.6 33.8
"1F"  "FVB_12" "FVB" 12 5.7 33.9
"2F"  "FVB_12" "FVB" 12 6.1   34
"3F"  "FVB_12" "FVB" 12 3.5 31.5
"4F"  "FVB_12" "FVB" 12 6.2 34.1
"5F"  "FVB_12" "FVB" 12 6.1 33.4
"6F"  "FVB_12" "FVB" 12 6.4 34.6
"7F"  "FVB_12" "FVB" 12 5.2 33.7
"8F"  "FVB_12" "FVB" 12 5.8 34.3
"9F"  "FVB_12" "FVB" 12   6 34.6
"10F" "FVB_12" "FVB" 12 7.1 34.4
"1C"  "CD1_2"  "CD1"  2   2 30.7
"2C"  "CD1_2"  "CD1"  2 1.6 30.2
"3C"  "CD1_2"  "CD1"  2 1.8 30.9
"4C"  "CD1_2"  "CD1"  2   2 29.3
"5C"  "CD1_2"  "CD1"  2 1.8 28.7
"6C"  "CD1_2"  "CD1"  2 2.2 29.9
"7C"  "CD1_2"  "CD1"  2 2.6 30.2
"8C"  "CD1_2"  "CD1"  2 1.7 30.2
"9C"  "CD1_2"  "CD1"  2 2.1 31.1
"10C" "CD1_2"  "CD1"  2 2.7 30.7
end

the code used...

label def order 1 B6 2 FVB 3 CD1 
encode Strain_old, gen(order) label(order) 
sort Pnd order 
gen axis = sum(Pnd != Pnd[_n-1]) + sum(order != order[_n-1]) - 1 

labmask axis, values(Strain_old)
set scheme s1color

bysort Strain: egen med = median(Temperatura)
by Strain: egen lqt = pctile(Temperatura), p(25)
by Strain: egen uqt = pctile(Temperatura), p(75)



local colours gray blue green  

foreach value in B6 FVB CD1 {
    gettoken colour colours : colours 
    local call `call' || rbar lqt med axis if Strain_old == "`value'", barw(0.6) fcolor(none) lcolor(`colour') ///
    || rbar uqt med axis if Strain_old == "`value'", barw(0.6) fcolor(none) lcolor(`colour') ///
    || scatter Peso axis if Strain_old == "`value'", mc(`colour') ms(Oh)
}

twoway `call' xla(1/3 5/7 9/11 13/15 17/19, labsize(small) noticks valuelabel) legend(off) ytitle(Temperature °C) xtitle("") yla(, ang(h)) xli(4(4)16, lw(vthin) lc(gs12)) /// 
xmla(2 "2" 6 "4" 10 "6" 14 "8" 18 "12", tlength(*7) tlc(none) labsize(small)) yla(, ang(h))

I don't know where the mistakes is.

Thank you Nick Cox

[Nick Cox]

 
 scatter Peso axis if Strain_old == "`value'", mc(`colour') ms(Oh)

should be

 
 scatter Temperatura axis if Strain_old == "`value'", mc(`colour') ms(Oh)

[Tommaso Salvitti]

Having been on the pc for 15 hours, I could not see the error .. Thanks a lot Nick Cox

[Tommaso Salvitti]

Good morning Nick Cox and to everybody


below the dataset

* Example generated by -dataex-. To install: ssc install dataex
clear
input byte Subject str7 Strain str1 Pnd str5 Strain_old double Complex
 1 "FF_a"    "a" "FF"      .546331467862883
 2 "FF_a"    "a" "FF"       .39509966728639
 3 "FF_a"    "a" "FF"     .4077292966966801
 4 "FF_a"    "a" "FF"     .5569435640666627
 6 "FF_a"    "a" "FF"     .5694439971635508
 7 "FF_a"    "a" "FF"     .7668754084914023
 9 "FF_a"    "a" "FF"      .529592245288984
10 "FF_a"    "a" "FF"    .38977499661159526
11 "FF_a"    "a" "FF"    .39932506022848785
12 "FF_a"    "a" "FF"     .4405106630046985
14 "FF_a"    "a" "FF"     .3636165520400109
 2 "FF_b"    "b" "FF"     .3574576712639885
 3 "FF_b"    "b" "FF"     .4434051436852918
 4 "FF_b"    "b" "FF"     .4433325360143212
 5 "FF_b"    "b" "FF"     .4182362002972204
 6 "FF_b"    "b" "FF"     .4827678165336439
 7 "FF_b"    "b" "FF"    .38158529505582245
 8 "FF_b"    "b" "FF"     .5090677539682013
 9 "FF_b"    "b" "FF"     .4452743791553812
10 "FF_b"    "b" "FF"     .4665977347107954
11 "FF_b"    "b" "FF"    .23522146450886783
12 "FF_b"    "b" "FF"    .46746811664627574
13 "FF_b"    "b" "FF"     .5486848472063541
 1 "FF_c"    "c" "FF"     .7978107858441659
 2 "FF_c"    "c" "FF"     .7255134367282434
 3 "FF_c"    "c" "FF"     .6843173508124704
 4 "FF_c"    "c" "FF"     .7379227955206507
 5 "FF_c"    "c" "FF"     .5273927071330573
 7 "FF_c"    "c" "FF"     .7284914666386457
 8 "FF_c"    "c" "FF"    1.2165012365433783
10 "FF_c"    "c" "FF"     .5639426413606289
11 "FF_c"    "c" "FF"    1.0451232299410538
12 "FF_c"    "c" "FF"     .8758913890207222
13 "FF_c"    "c" "FF"     .6881810322886441
14 "FF_c"    "c" "FF"     .7944895734394164
 1 "aleF_a"  "a" "aleF"   .4642699836551305
 3 "aleF_a"  "a" "aleF"   .4229081000494186
 4 "aleF_a"  "a" "aleF"   .5321739109333623
 5 "aleF_a"  "a" "aleF"   .4837161463625206
 7 "aleF_a"  "a" "aleF"   .4793442903881181
 9 "aleF_a"  "a" "aleF"   .4985965399981203
10 "aleF_a"  "a" "aleF"   .7853981633974484
11 "aleF_a"  "a" "aleF"   .3939347372399995
12 "aleF_a"  "a" "aleF"    .595165379714826
14 "aleF_a"  "a" "aleF"   .4876750483210819
 1 "aleF_b"  "b" "aleF"   .5373842651993277
 2 "aleF_b"  "b" "aleF"   .7262684537097149
 3 "aleF_b"  "b" "aleF"   .7172864682103272
 4 "aleF_b"  "b" "aleF"   .6171464342684478
 5 "aleF_b"  "b" "aleF"   .7553253905587914
 6 "aleF_b"  "b" "aleF"   .8351880119340882
 7 "aleF_b"  "b" "aleF"   .3639322076101682
 8 "aleF_b"  "b" "aleF"   .8832436898300681
10 "aleF_b"  "b" "aleF"   .9697814053066289
11 "aleF_b"  "b" "aleF"    .614624423286827
12 "aleF_b"  "b" "aleF"   .6135877788059929
13 "aleF_b"  "b" "aleF"   .7345082070708817
 1 "aleF_c"  "c" "aleF"   1.041738426714682
 2 "aleF_c"  "c" "aleF"  1.0529381928213029
 4 "aleF_c"  "c" "aleF"   1.024091424415769
 5 "aleF_c"  "c" "aleF"    .730904380384288
 6 "aleF_c"  "c" "aleF"   .8792995563393209
 7 "aleF_c"  "c" "aleF"   1.036021646008127
 8 "aleF_c"  "c" "aleF"   .8942245151603748
 9 "aleF_c"  "c" "aleF"   .9353315637689095
10 "aleF_c"  "c" "aleF"   .8835108575164353
11 "aleF_c"  "c" "aleF"   .7901004270698991
12 "aleF_c"  "c" "aleF"    .646448444729729
13 "aleF_c"  "c" "aleF"  1.0697033135295395
 1 ".Pups_a" "a" ".Pups"  .5050535370773326
 1 ".Pups_a" "a" ".Pups"   .544451375681954
 1 ".Pups_a" "a" ".Pups"  .1399046578562445
 1 ".Pups_a" "a" ".Pups"   .310797463107646
 1 ".Pups_a" "a" ".Pups"  .2876560949628697
 1 ".Pups_a" "a" ".Pups" .49772785297418687
 1 ".Pups_a" "a" ".Pups" .24990671744613882
 1 ".Pups_a" "a" ".Pups"  .4820012402222069
 1 ".Pups_b" "b" ".Pups"   .376101832364541
 1 ".Pups_b" "b" ".Pups" .43921522994069007
 1 ".Pups_b" "b" ".Pups"  .7312953020520645
 1 ".Pups_b" "b" ".Pups"  .5772024410326464
 1 ".Pups_b" "b" ".Pups"  .6011135433083147
 1 ".Pups_b" "b" ".Pups"  .5871421345546757
 1 ".Pups_b" "b" ".Pups"    .76272758825471
 1 ".Pups_b" "b" ".Pups" .46834850397523203
 1 ".Pups_b" "b" ".Pups"   .716973743737547
 1 ".Pups_b" "b" ".Pups"  .6021712195327014
 1 ".Pups_c" "c" ".Pups"  .5395297632415174
 1 ".Pups_c" "c" ".Pups"  .4143409943426507
 1 ".Pups_c" "c" ".Pups"  .6540422639407324
 1 ".Pups_c" "c" ".Pups"  .5256190648164841
 1 ".Pups_c" "c" ".Pups"   .505367122655699
 1 ".Pups_c" "c" ".Pups"  .6492885003690742
 1 ".Pups_c" "c" ".Pups"  .7726368975579525
 1 ".Pups_c" "c" ".Pups"  .5784077328865905
 1 ".Pups_c" "c" ".Pups"  .3813347725502865
 1 ".Pups_c" "c" ".Pups"  .5904671674641082
end

i use the following code

sort Contesto
by Contesto: egen med = median(Complex)
by Contesto: egen lqt = pctile(Complex), p(25)
by Contesto: egen uqt = pctile(Complex), p(75)


label def Contesto 1 ".Pups_a" 2 ".Pups_b" 3 ".Pups_c" 5 "aleF_a"  6 "aleF_b" 7 "aleF_c" 9 "FF_a" 10 "FF_b" 11 "FF_c"
encode Contesto, gen(axis) label(Contesto)

local colours grey*0.1 blue*0.1 green*0.1 grey*0.3 blue*0.3 green*0.3 grey*0.5 blue*0.5 green*0.5

local Colours grey blue green grey blue green grey blue green

local values 1 2 3 5 6 7 9 10 11

set scheme s1color

forval j = 1/9 {
    local c : word `j' of `colours'
    local C : word `j' of `Colours'
    local v : word `j' of `values'
    
    local call `call' || rbar lqt med axis if axis == `v', barw(0.6) fcolor(`c') lcolor(`C') ///
    || rbar uqt med axis if axis == `v', barw(0.6) fcolor(`c') lcolor(`C')                   ///
    || scatter Complex axis if axis == `v', mc(`C') ms(Oh)
}

twoway `call' xla(1 2 3 5 6 7 9 10 11, noticks valuelabel labsize(small)) legend(off) ytitle(Complex) xtitle("") yla(, ang(h)) xli(4 8, lw(vthin) lc(gs12))

and I obtain the grapgh that i attached name "Complex"

I'm in trouble because I have to reproduce graphics without the colored box (only the outline). How to do? in practice I should get graphs like this which I am attaching called "proof". Can you help me please? Thanks in advanced

[Tommaso Salvitti]

Nick Cox
I did it thanks a lot. However, I was unable to delete the post where I was asking for help (#34)

[Nick Cox]

Pleased you got what you want, but it's time to read the FAQ Advice. https://www.statalist.org/forums/help#closure explains that you cannot delete posts.