Sorry for the mangled page lay-out. I try again to make the macro readable.
-ftr

Hello,
I did not find a function that produces a Box-Cox transformation in PSPP. On Raynald's SPSS macro pages a SPSS macro with syntax <http://www.spsstools.net/en/syntax/syntax-index/compute/box-cox-transformation/> has been published is published that does the job.

* Box-Cox transformation for all 31 values of lambda between -2 to 1 (increments of .1).
* Raynald Levesque 2003/11/08.

GET FILE='c:\\program files\\spss\\employee data.sav'.
 COMPUTE var1=salary.
VECTOR lam(31) /xl(31).
LOOP idx=1 TO 31.
- COMPUTE lam(idx)=-2.1 + idx * .1.
- DO IF lam(idx)=0.
-     COMPUTE xl(idx)=LN(var1).
- ELSE.
-     COMPUTE xl(idx)=(var1**lam(idx) - 1)/lam(idx).
- END IF.
END LOOP.
EXECUTE.

I used this Box-Cox transformation macro <http://www.spsstools.net/en/syntax/syntax-index/compute/box-cox-transformation/>which produces a series of values for the variable to be transformed called xl for values of lambda above -2.1. The original variable in the test data set, salary, shows a skew of 2.12 and a kurtosis of above 5.

After having run the macro I don't see how to continue. How to decide which is the best lambda value and thus, the best transformation?

I  then did a frequency of all transformed xl variables to find the transformed , i.e. xl, variable with the skew and kurtosis the closest to 0. I also tried to visually find the best transformation.
FREQUENCIES xl1 to xl31
    /FORMAT= NOTABLE
    /STATISTICS=SKEWNESS KURTOSIS.
EXAMINE
VARIABLES= xl1 to xl31
/PLOT= NPPLOT.

I then graphed the skew for each lambda value in a spreadsheet and got skewness values of close to 0 for lambda around -1.

Are these the appropriate ways ?
In my view the results show that xl10 and xl11 show a good combination of skew and kurtosis and are together close to 0.

My questions:
1/ Are there other ways in PSPP to run a box-cox transformation ?
2/ Which is the best way to find the best transformation ?


If not what do you recommend how to decide about the appropriate transformed variable ?


Thanks in advance,

ftr



_______________________________________________
Pspp-users mailing list
Pspp-users@gnu.org
https://lists.gnu.org/mailman/listinfo/pspp-users

Reply via email to