Re: [Gretl-users] aggregate by

Hi, I'd like to point out, that a 'good' 'aggregate by' function would do such a job as described below. Unfortunately I'm stuck with more then bivariate interaction but I'm sure GRETL can do this, this question is how! It's even more complex since the R function is not that good either, cos e.g. for MAVOVA one needs the variance-covariance-structure of the 'interaction series'. That's why in GRETL lists should be processed which makes the stuff not easier. I worked around with the new '^' operator and have some promising results, but a build-in function would be the best. Cheers Leon <hansl> open pizza4.gdt series ages = age<=20 ? 0 : (age>20 && age <=40) ? 1 : 2 set echo off foreign language=R --send-data --quiet library(plyr) ddply(gretldata[,c(1,2,4,6,8)], ~ female + college + ages, numcolwise(mean)) end foreign <hansl> While the bivariate case in rather easy as in <hansl> open pizza4.gdt set echo off series ages = age<=20 ? 0 : (age>20 && age <=40) ? 1 : 2 list G1 = 2 4 list G2 = 2 8 list G3 = 4 8 list G = G1 || G2 scalar times = 2 # should also be possible to be set to nelem(G) #set up of all combination list loop for i=2..times -q string listName = "inter$i" list @listName = null endloop list part1 = null list part2 = G # bivariate interaction loop foreach i G-q part1 = G.$i part2 -= G.$i inter2 += part2^part1 part2 += G.$i endloop <hansl> Am 20.02.2013 17:45, schrieb Riccardo (Jack) Lucchetti: