1:07 p.m.
Hi there,
I have a function (WIP) to produce input for a boxplot side-by-side comparison of
cross-section statistics of two samples (see below). The list seems to pass the names of
the series, but does not work. If I remove the function body, everything works.
<hansl>
function matrix MatOrderDist (list ls, list lsB[null])
#A function to create data for order stats boxplots (cross-section)
scalar two = 0 #no
scalar entries = nelem(ls)
if exists(lsB)
if nelem(lsB)==entries
two = 1 #yes
endif
endif
strings coln = array(6)
coln[1] = "min"
coln[2] = "lq"
coln[3] = "med"
coln[4] = "uq"
coln[5] = "max"
coln[6] = "avg"
if two
matrix to_plot = zeros(entries,12)
strings colnf = array(12)
loop i=1..12
colnf[i]=coln[i]
colnf[i+1]=coln[i]
i+=1
endloop
else
matrix to_plot = zeros(entries,6)
strings colnf = array(6)
loop i=1..6
colnf[i]=coln[i]
endloop
endif
colnames(to_plot,colnf)
scalar r = 0
scalar offset = 0
if two
offset = 1
endif
loop foreach item ls -q
r+=1
to_plot[r,1] = min($item)
to_plot[r+offset*r,2] = quantile($item,.25)
to_plot[r+offset*r,3] = quantile($item,.5)
to_plot[r+offset*r,4] = quantile($item,.75)
to_plot[r+offset*r,5] = max($item)
to_plot[r+offset*r,6] = mean($item)
endloop
if two
loop foreach item lsB -q
r+=1
to_plot[r+offset*(r-1),1] = min($item)
to_plot[r+offset*(r-1),2] = quantile($item,.25)
to_plot[r+offset*(r-1),3] = quantile($item,.5)
to_plot[r+offset*(r-1),4] = quantile($item,.75)
to_plot[r+offset*(r-1),5] = max($item)
to_plot[r+offset*(r-1),6] = mean($item)
endloop
endif
#boxplot --matrix=to_plot --output=display {set title "Cross-Section Distribution
of Summary Statistics (Original Programme)"; }
return to_plot
end function
</hansl>
1:16 p.m.
Am 05.01.2018 um 14:07 schrieb Schaff, Frederik:
list seems to pass the names of the series, but does not work. If I
remove the function body, everything works.
Yes, see the "Foreach loop" section in the loop chapter 12 of the manual:
"Note that if you use this idiom inside a function (see chapter 13),
looping across a list that has been supplied to the function as an
argument, it is necessary to use the syntax listname.$i
to reference the list-member variables."
(Of course, $i must be replaced with what the loop index is in your case.)
function matrix MatOrderDist (list ls, list lsB[null])
if exists(lsB)
Hey, we had this before, don't use exists() with list function args. It
will always return 1, even for null. Use nelem() directly (with lists
only!).
loop foreach item ls -q
to_plot[r,1] = min($item)
OK, this would need to be ls.$item instead.
loop foreach item lsB -q
to_plot[r+offset*(r-1),1] = min($item)
And this: lsB.$item.
hth,
sven
1:25 p.m.
Many thanks!
I used exists and nelem... because the exists doesn't matter yet but might in future.
Anyway, there are other bugs. I wasn't aware of the forach problem. Thanks!
(And there are still other bugs in my code example...)
Frederik
-----Ursprüngliche Nachricht-----
Von: gretl-users-bounces(a)lists.wfu.edu [mailto:gretl-users-
bounces(a)lists.wfu.edu] Im Auftrag von Sven Schreiber
Gesendet: Freitag, 5. Januar 2018 14:17
An: gretl-users(a)lists.wfu.edu
Betreff: Re: [Gretl-users] passing list to function, loop
Am 05.01.2018 um 14:07 schrieb Schaff, Frederik:
> list seems to pass the names of the series, but does not work. If I
> remove the function body, everything works.
Yes, see the "Foreach loop" section in the loop chapter 12 of the manual:
"Note that if you use this idiom inside a function (see chapter 13), looping
across a list that has been supplied to the function as an argument, it is
necessary to use the syntax listname.$i to reference the list-member
variables."
(Of course, $i must be replaced with what the loop index is in your case.)
> function matrix MatOrderDist (list ls, list lsB[null])
> if exists(lsB)
Hey, we had this before, don't use exists() with list function args. It will
always return 1, even for null. Use nelem() directly (with lists only!).
> loop foreach item ls -q
> to_plot[r,1] = min($item)
OK, this would need to be ls.$item instead.
> loop foreach item lsB -q
> to_plot[r+offset*(r-1),1] = min($item)
And this: lsB.$item.
hth,
sven
_______________________________________________
Gretl-users mailing list
Gretl-users(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-users
2529
days inactive
2529
days old
2 comments
2 participants
participants (2)
-
Schaff, Frederik -
Sven Schreiber